Diff two solutions on two meshes

[tangqi]

Is there a simple way to diff two solutions on two non conforming meshes? For instance, one is a refined solution on a uniform mesh and one is an AMR solution. I found the field-diff app. The example is very simple, but will that work in parallel AMR meshes (quad)? Thanks.

If that's not possible, other suggestions are also welcome. I am trying to do convergence study for problems with no exact solutions.

Qi

[vladotomov]

Hi Qi,

Yes field-diff should work with AMR meshes, but we don't have a parallel version yet. You can use the PrintAsOne and SaveAsOne methods of ParMesh and ParGridFunction and use field-diff on the resulting files.

Vladimir

[tangqi]

Hi, I try to compare a full mesh vs an amr mesh. It is not working for me. See the attached solutions.
archive.tar.gz

Thanks,
Qi

PS, the amr is non-conforming but it is saved as "MFEM mesh v1.0". But I do not think it is the reason.

[tangqi]

ping @vladotomov @tzanio

[vladotomov]

Hi Qi,

Thank you for reporting. I fixed some things and was able to run your comparison with the branch gslib-fix, the PR is #2052.

[tangqi]

Vladimir, thanks for the quick fix. It looks great!