Manufactured solution test

mimesis-inria / caribou

Multi-physics computation library

GNU Lesser General Public License v3.0

29 stars 17 forks source link

Manufactured solution test #116

Closed Sidaty1 closed 2 years ago

Sidaty1 commented 2 years ago

We would like to validate the HyperelasticForceField_FeniCS on a manufactured solution

Ziemnono commented 2 years ago

Works for linear tetrahedron using SaintVenantKirchhoff, a bit of work has still to be done for quadratic tetrahedron (difference due to index changes)

Ziemnono commented 2 years ago

Last results for the quadratic tetrahedron

SOFA relative L2 error at 100% load: 9.582959168931972e-06
FEniCS relative L2 error at 100% load: 9.582959166061108e-06
Difference in relative L2 error at 100% load: 2.8708640147220507e-15

Seems good enough for me :)