Integral conservation of adapted function?
12 months ago by
it appears that adaptation of a function to a refined mesh does not preserve the integral of the function.
The following code produces two different results for me:
print(str(assemble(u0*dx))) refined_mesh = refine(mesh, cell_markers) adapt(u0,refined_mesh) dx = Measure("dx",domain =refined_mesh) print(str(assemble(u0*dx)))
Is this a typical behavior or am I doing something wrong?
Community: FEniCS Project
9 months ago by
As a follow up question: Is it possible to prevent this, meaning using the integral conservation as constraint for adapting the function?
Please login to add an answer/comment or follow this question.