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	2421	June 1, 2022
How to apply boundary conditions to left-hand matrix Errors dolfinx	6	268	June 25, 2024
Adding rows of Matrix Linear Algebra	2	775	July 21, 2020
Some questions about Matrix and Vectors dolfinx	2	211	October 5, 2023
Invert DG block diagonal matrix General	4	618	September 21, 2023

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

Related topics