I have added PCApplyTranspose() for PCCOMPOSITE in petsc-dev. You will need
to switch to petsc-dev
http://www.mcs.anl.gov/petsc/petsc-as/developers/index.html to access it.
Sorry for the inconvenience.
Barry
On Jun 23, 2011, at 4:13 AM, Andrew T Barker wrote:
> I've made a PCCOMPOSIT
Thanks for helping me check the code.
I agree that moving the KSPSetOperators(ksp,M,M,SAME_NONZERO_PATTERN)
outside the iteration loop lifted the crash of superlu_dist. However,
I have the following concerns:
1) If we only call KSPSetOperators once outside the loop, the ksp
solver would be the
Use the KSP routines, there is no reason to use the low-level routines, they
will not be any less efficient.
Barry
On Jun 24, 2011, at 10:20 AM, Xiangdong Liang wrote:
> Thanks for helping me check the code.
>
> I agree that moving the KSPSetOperators(ksp,M,M,SAME_NONZERO_PATTERN)
> out
So, any ideas?
when I call MatMult(A,b,x,info) the ghosts get updated. When I call
KSPSolve(ksp,b,x,info) the ghosts dont get updated (the multiplication
routine is called... I checked that)... It seems that KSP calles
MatMult with different vectors... not the original b and x...
ov/pipermail/petsc-users/attachments/20110624/7e2cabf6/attachment.htm>
