Dear Hong,
Yes, I am using this SAME_NONZERO_PATTERN to optimize the code,hehe!
Thanks very much!
Best,
Rebecca
On Sep 1, 2011, at 7:35 AM, Hong Zhang wrote:
Rebecca:
Thanks for your kind reply. So if I set up the following:
ierr = DAGetMatrix(DMMGGetDA(dmmg), MATAIJ,
Rebecca:
Thanks for your kind reply. So if I set up the following:
? ? ? ?ierr = DAGetMatrix(DMMGGetDA(dmmg), MATAIJ, jacobian);CHKERRQ(ierr);
? ? ? ?KSP ? ? ksp;
? ? ? ?ierr = SNESGetKSP(DMMGGetSNES(dmmg),ksp);CHKERRQ(ierr);
? ? ? ?ierr =
Hello all,
Is the default nonzero pattern for SuperLU_dist is SAME_NONZERO_PATTERN?
Thanks!
Best,
Rebecca
On Aug 31, 2011, at 4:40 PM, Xuefei (Rebecca) Yuan wrote:
Hello all,
Is the default nonzero pattern for SuperLU_dist is SAME_NONZERO_PATTERN?
There is no default nonzero pattern. You always must set it when you call
KSPSetOperators() or in your SNES FormFunction. The PETSc SUPERLU_dist
Dear Barry,
Thanks for your kind reply. So if I set up the following:
ierr = DAGetMatrix(DMMGGetDA(dmmg), MATAIJ, jacobian);CHKERRQ(ierr);
KSP ksp;
ierr = SNESGetKSP(DMMGGetSNES(dmmg),ksp);CHKERRQ(ierr);
ierr =