Set symmetry options on Stokes preconditioning matrices

jedbrown · Apr 24, 2011 · fbebec2 · fbebec2
1 parent 00d06c4
commit fbebec2
Showing 1 changed file with 9 additions and 0 deletions.
diff --git a/src/fs/tests/stokes.c b/src/fs/tests/stokes.c
@@ -410,6 +410,15 @@ static dErr StokesGetMatrices(Stokes stk,dBool use_jblock,Mat *J,Mat *Jp)
   err = dFSGetMatrix(stk->fsp,stk->mattype_D,&D);dCHK(err);
   err = MatSetOptionsPrefix(A,"Ap_");dCHK(err);
   err = MatSetOptionsPrefix(D,"Dp_");dCHK(err);
+  err = MatSetOption(A,MAT_SYMMETRIC,PETSC_TRUE);dCHK(err);
+  err = MatSetOption(D,MAT_SYMMETRIC,PETSC_TRUE);dCHK(err);
+  {
+    dBool seqsbaij;
+    err = PetscTypeCompare((PetscObject)A,MATSEQSBAIJ,&seqsbaij);dCHK(err);
+    if (seqsbaij) {err = MatSetOption(A,MAT_IGNORE_LOWER_TRIANGULAR,PETSC_TRUE);dCHK(err);}
+    err = PetscTypeCompare((PetscObject)D,MATSEQSBAIJ,&seqsbaij);dCHK(err);
+    if (seqsbaij) {err = MatSetOption(D,MAT_IGNORE_LOWER_TRIANGULAR,PETSC_TRUE);dCHK(err);}
+  }
   err = MatCreateNest(stk->comm,2,splitis,2,splitis,((Mat[]){A,NULL,NULL,D}),Jp);dCHK(err);
   err = MatSetOptionsPrefix(*Jp,"Jp_");dCHK(err);
   err = MatSetFromOptions(*Jp);dCHK(err);