Merge branch 'jolivet/redundant-Cholesky-sbaij' [PR #1968]
* jolivet/redundant-Cholesky-sbaij: PCCHOLESKY in PCREDUNDANT with SBAIJ matrices.
PCREDUNDANT does not currently play nice with SBAIJ mat
Merge branch 'jolivet/redundant-Cholesky-sbaij' [PR #1968]
* jolivet/redundant-Cholesky-sbaij: PCCHOLESKY in PCREDUNDANT with SBAIJ matrices.
PCREDUNDANT does not currently play nice with SBAIJ matrices out-of-the-box and you end up with something like the following.
[0]PETSC ERROR: #1 MatGetFactor() line 4411 in src/mat/interface/matrix.c [0]PETSC ERROR: #2 PCSetUp_LU() line 88 in src/ksp/pc/impls/factor/lu/lu.c [0]PETSC ERROR: #3 PCSetUp() line 894 in src/ksp/pc/interface/precon.c [0]PETSC ERROR: #4 KSPSetUp() line 377 in src/ksp/ksp/interface/itfunc.c [0]PETSC ERROR: #5 PCSetUp_Redundant() line 179 in src/ksp/pc/impls/redundant/redundant.c
This is an attempt to fix this by switching to PCCHOLESKY instead of PCLU for symmetric matrices.
show more ...
|