Bug fixes for ustruct Guccione model and follower pressure load

[aabrown100-git]

Current situation

Slow nonlinear convergence for in ustruct/LV_Guccione_active. See issue: #197

Release Notes

Fix typos in Guccione material model used for ustruct, and follower pressure load implementation for ustruct. These are now consistent with the corresponding lines in svFSI (Fortran)

I made these changes to match the following lines in svFSI (Fortran) code
https://github.com/SimVascular/svFSI/blob/7bdb74cdd517310d6bf32c0bb7d42072aeac8658/Code/Source/svFSI/MATMODELS.f#L809

https://github.com/SimVascular/svFSI/blob/7bdb74cdd517310d6bf32c0bb7d42072aeac8658/Code/Source/svFSI/MATMODELS.f#L812

https://github.com/SimVascular/svFSI/blob/7bdb74cdd517310d6bf32c0bb7d42072aeac8658/Code/Source/svFSI/USTRUCT.f#L1320

Unfortunately, after making these changes, the convergence for ustruct/LV_Guccione_active is still not good. There must be another bug in the tangent somewhere, but I've gone through the Guccione model in get_pk2cc_dev() several times and can't find anything.

Testing

Warning: these changes produced different results for ustruct/LV_Guccione_active than the existing reference solution. I updated the reference solution for this test.

Changes do not significantly improve convergence of ustruct/LV_Guccione_active. There are likely other bugs. However, I think these are clearly bugs, so it makes sense to fix them now.

Code of Conduct & Contributing Guidelines

• I agree to follow the Code of Conduct and Contributing Guidelines.

[codecov]

Codecov Report

All modified and coverable lines are covered by tests ✅

Project coverage is 62.24%. Comparing base (cf3422d) to head (93ba1ef).

Additional details and impacted files

@@           Coverage Diff           @@
##             main     #201   +/-   ##
=======================================
  Coverage   62.24%   62.24%           
=======================================
  Files         103      103           
  Lines       26964    26964           
=======================================
  Hits        16785    16785           
  Misses      10179    10179           

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

[mrp089]

Thank you for digging into this, @aabrown100-git! It makes sense that you had to change the reference solution since the computation of the stress changed.

We're down to 17 Newton iterations from 19. Every bug counts! I'll let @ktbolt review this.

[mrp089]

@ktbolt, can you review this PR? Thank you!

[ktbolt]

Looks good!