Skip to content

Implement PETSc proxy matrix assembly for 3D freeslip cases

Eugenio D'Ascoli requested to merge edascoli/PETScFreeSlipTestApply3D into master

Adding the function saveProjectNormalOperator3D to EdgeDoFProjectNormalOperator.cpp and P1ProjectNormalOperator.cpp to utilize the PETSc coarse grid solver for Freeslip boundary condition cases in 3D. The function saveProjectNormalOperator3D is added to correctly assemble the PETSc proxy matrix for the 3D normal projection operator for cases with Freeslip boundary conditions. Prior to this commit 3D cases, like a spherical shell, were not supported.

Merge request reports