Hello, I’m trying to do a fairly simple thing: Interpolate an expression into a function named “test”, and then take the vector norm via PETSc (in dolfinx from today). The underlying vector has the right entries, however the norm routine always returns the initial value after test.vector.set(0.0), …

Interpolate into function, and fail in PETSc vector norm (dolfinx)

dokken October 13, 2020, 8:28am 2

Add the line:

test.vector.ghostUpdate(addv=PETSc.InsertMode.ADD, mode=PETSc.ScatterMode.REVERSE)

test.vector.assemble()

after interpolation.

Different results in serial and parallel run dolfinx

Topic		Replies	Views
Some questions about solving the nonlinear vector-valued problem dolfinx	7	184	January 10, 2024
Complex valued Neumann b.c dolfinx	3	491	August 5, 2021
Inconsistent results with different dolfinx/petsc versions	0	397	July 28, 2022
Taking function converting it to petsc4py vector, and then putting it back into function?	5	1431	April 25, 2021
Evaluate PETSc matrix for form that has gradient of test function dolfinx	1	563	December 10, 2021