No need to hard-code velocity assembly for Q1

jedbrown · Jun 6, 2011 · b9fc817 · b9fc817
1 parent d1d2ab2
commit b9fc817
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/fs/tests/stokes.c b/src/fs/tests/stokes.c
@@ -722,7 +722,7 @@ static dErr StokesJacobianAssemble_Velocity(Stokes stk,Mat Ap,Vec Mdiag,Vec gx)
           }
         }
       }
-      err = dFSMatSetValuesBlockedExpanded(stk->fsu,Ap,8,rowcol,8,rowcol,&K[0][0][0][0],ADD_VALUES);dCHK(err);
+      err = dFSMatSetValuesBlockedExpanded(stk->fsu,Ap,P,rowcol,P,rowcol,&K[0][0][0][0],ADD_VALUES);dCHK(err);
       for (dInt i=0; i<P; i++) {
         dScalar Mentry = 0;
         for (dInt q=0; q<Q; q++) Mentry += interp[q][i] * jw[q] * interp[q][i]; /* Integrate the diagonal entry over this element */

diff --git a/src/fs/tests/vht.c b/src/fs/tests/vht.c
@@ -1802,7 +1802,7 @@ static dErr VHTJacobianAssemble_Velocity(VHT vht,Mat Buu,Vec Mdiag,Vec X)
           }
         }
       }
-      err = dFSMatSetValuesBlockedExpanded(vht->fsu,Buu,8,rowcol,8,rowcol,&K[0][0][0][0],ADD_VALUES);dCHK(err);
+      err = dFSMatSetValuesBlockedExpanded(vht->fsu,Buu,P,rowcol,P,rowcol,&K[0][0][0][0],ADD_VALUES);dCHK(err);
       for (dInt i=0; i<P; i++) {
         dScalar Mentry = 0;
         for (dInt q=0; q<Q; q++) Mentry += interp[q][i] * jw[q] * interp[q][i]; /* Integrate the diagonal entry over this element */