I am trying to project values on particular points onto a function space so that I can compare it to error values that live on that function space. Fenics provides a function called ‘project’ that accomplishes this, but now I am working with fenicsx and I can’t find the function anywhere (I looked at dolfinx and several of its submodules). Does it still exist and if so, where can I find it?
The problem with the project syntax is the uninformed user will assume it will handle many esoteric cases. For example, incomplete data in f, singular components in V, obscure functions f, reusing the LU factorisation for multiple projections onto the same space, and so on… Each of these user dependent cases should be handled by whatever the user’s numerical problem requires, and not a catch-all project function.
Thanks for the explanation @nate. I am working with @bouncer on doing a rather easy calculation for which previously we used the project function from FEniCS. Let me rephrase the question:
We have a function with a reference solution and a function with the calculated solution:
over the course of the simulation, u_n gets resolved. At some point we want to calculate the relative error between the reference and the calculated solution in the l2 norm. Previously we did this by:
for which we needed the project functionality. Is there a simpler way to calculate the total error in dolfinx? I am guessing what we are trying to do is not the most efficient way.