Hi all, I have an assembled matrix with zero rows, where I’d like to insert ones on the diagonal. In older versions there was the A.ident_zeros() function which did exactly that, and fast. Such a function seems not to be implemented in dolfinx, so I tried to find a workaround with PETSc, using: A…

I would suggest the following: import time import dolfinx import numpy as np import ufl from dolfinx import cpp, fem from mpi4py import MPI from petsc4py import PETSc N = 75 degree = 1 mesh = dolfinx.mesh.create_unit_cube(MPI.COMM_WORLD, N, N, N) V = dolfinx.fem.FunctionSpace(mesh, ("CG", degree)…

Modify matrix diagonal -- dolfinx version for A.ident_zeros()

General

CastriMik December 7, 2022, 1:46pm 2

You can find a version of ident_zeros() for DOLFINx in this post.

Topic		Replies	Views
What's the best way to assemble a form? dolfinx	14	2283	June 1, 2022
Some questions about Matrix and Vectors dolfinx	2	209	October 5, 2023
Access matrix in dolfinx I/O assembly , dolfinx	2	63	January 14, 2025
Remove many zero entries from assembled sparse matrix or avoid them Linear Algebra	0	533	June 3, 2021
How to apply boundary conditions to left-hand matrix Errors dolfinx	6	229	June 25, 2024

Modify matrix diagonal -- dolfinx version for A.ident_zeros()

Related topics