The Stokes problem is no longer symmetric when written in conservativ…

…e form with ice velocity != total velocity
jedbrown · May 29, 2011 · 0d00934 · 0d00934
1 parent 8a9c0ed
commit 0d00934
Showing 1 changed file with 0 additions and 1 deletion.
diff --git a/src/fs/tests/vht.c b/src/fs/tests/vht.c
@@ -835,7 +835,6 @@ static dErr VHTGetMatrices(VHT vht,dBool use_jblock,Mat *J,Mat *P)
   err = MatSetOptionsPrefix(Buu,"Buu_");dCHK(err);
   err = MatSetOptionsPrefix(Bpp,"Bpp_");dCHK(err);
   err = MatSetOptionsPrefix(Bee,"Bee_");dCHK(err);
-  err = MatSetOption(Buu,MAT_SYMMETRIC,PETSC_TRUE);dCHK(err);
   err = MatSetOption(Bpp,MAT_SYMMETRIC,PETSC_TRUE);dCHK(err);
   err = MatSetFromOptions(Buu);dCHK(err);
   err = MatSetFromOptions(Bpp);dCHK(err);