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	2084	June 1, 2022
Some questions about Matrix and Vectors dolfinx	2	201	October 5, 2023
Evaluate PETSc matrix for form that has gradient of test function dolfinx	1	595	December 10, 2021
Function parameters	2	494	August 23, 2019
Manipulating vector data of fem.Function dolfinx	2	671	April 30, 2023

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

Related Topics