KSPSolveTranspose and PCFieldSplit

[Randall Mackie]

When I try to use KSPSolveTranspose with PCFieldSplit, I get the following errors:

[0]PETSC ERROR: PCApplyTranspose() line 483 in src/ksp/pc/interface/precon.c
[0]PETSC ERROR: No support for this operation for this object type!
[0]PETSC ERROR:  !
[0]PETSC ERROR: KSPInitialResidual() line 56 in src/ksp/ksp/interface/itres.c
[0]PETSC ERROR: KSPSolve_BCGS() line 41 in src/ksp/ksp/impls/bcgs/bcgs.c
[0]PETSC ERROR: KSPSolveTranspose() line 562 in src/ksp/ksp/interface/itfunc.c


So then I thought in my code I would just take the transpose of A and call
KSPSolve, but then I get the following errors:

[0]PETSC ERROR: MatGetSubMatrix_SeqAIJ() line 1686 in src/mat/impls/aij/seq/aij.c
[0]PETSC ERROR: Nonconforming object sizes!
[0]PETSC ERROR: Cannot reuse matrix. wrong no of nonzeros!
[0]PETSC ERROR: MatGetSubMatrices_SeqAIJ() line 1799 in src/mat/impls/aij/seq/aij.c
[0]PETSC ERROR: MatGetSubMatrices() line 4828 in src/mat/interface/matrix.c
[0]PETSC ERROR: MatGetSubMatrix() line 5438 in src/mat/interface/matrix.c
[0]PETSC ERROR: PCSetUp_FieldSplit() line 162 in src/ksp/pc/impls/fieldsplit/fieldsplit.c
[0]PETSC ERROR: PCSetUp() line 798 in src/ksp/pc/interface/precon.c
[0]PETSC ERROR: KSPSetUp() line 234 in src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: KSPSolve() line 334 in src/ksp/ksp/interface/itfunc.c


Is there some reason why the PCFIeildSplit is not compatible with KSPSolveTranspose?

I'm solving with the following options:


          -em_ksp_truemonitor \
          -em_ksp_type bcgs \
          -em_pc_type fieldsplit \
          -em_pc_fieldsplit_default \
          -em_pc_fieldsplit_type multiplicative \
          -em_fieldsplit_pc_type bjacobi \
          -em_fieldsplit_sub_pc_type ilu \
          -em_fieldsplit_sub_pc_factor_levels 4 \
          -em_fieldsplit_sub_pc_factor_fill 10 \

Randy

-- 
Randall Mackie
GSY-USA, Inc.
PMB# 643
2261 Market St.,
San Francisco, CA 94114-1600
Tel (415) 469-8649
Fax (415) 469-5044

California Registered Geophysicist
License No. GP 1034