Forward GaussFamily when using dQUADRATURE_METHOD_SELF

jedbrown · Mar 18, 2011 · 4a2033e · 4a2033e
1 parent ea93ac7
commit 4a2033e
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/src/fs/interface/fsruleset.c b/src/fs/interface/fsruleset.c
@@ -43,6 +43,7 @@ dErr dFSGetPreferredQuadratureRuleSet(dFS fs,dMeshESH set,dEntType etype,dEntTop
   rset->topo = etopo;
   rset->n = ents_s;
   err = dJacobiGetQuadrature(fs->jacobi,method,&quad);dCHK(err);
+  if (method == dQUADRATURE_METHOD_SELF) {err = dQuadratureTensorSetGaussFamily(quad,dGAUSS_LOBATTO);dCHK(err);}
   err = dQuadratureGetRules(quad,ents_s,topo,order,&rset->rules);dCHK(err);
   err = dFree3(ents,topo,order);dCHK(err);
   *ruleset = rset;