Feature/povolny1/polar decomposition

GEOS-DEV / LvArray

Portable HPC Containers (C++)

BSD 3-Clause "New" or "Revised" License

47 stars 10 forks source link

Feature/povolny1/polar decomposition #277

Open povolny1 opened 1 year ago

povolny1 commented 1 year ago

This PR implements the polar decomposition for real square 2x2 and 3x3 matrices using Higham's iterative algorithm. This is needed for some constitutive models and particularly a new hypoelastic update that I'm implementing.

povolny1 commented 1 year ago

Unit test is implemented and passing, I will make an accompanying PR for GEOSX soon