Commits

Lisandro Dalcin  committed 161b52b

Fix bug in initialization of quadrature rules

  • Participants
  • Parent commits 3e14cf9

Comments (0)

Files changed (1)

File src/petiga.c

   }
 
   for (i=0; i<3; i++)
-    if (!iga->collocation && (i >= iga->dim || iga->rule[i]->nqp < 1))
-      {ierr = IGARuleInit(iga->rule[i],iga->axis[i]->p + 1);CHKERRQ(ierr);}
+    if (!iga->collocation) {
+      if (i >= iga->dim || iga->rule[i]->nqp < 1)
+        {ierr = IGARuleInit(iga->rule[i],iga->axis[i]->p + 1);CHKERRQ(ierr);}
+    }
     else
       {ierr = IGARuleReset(iga->rule[i]);CHKERRQ(ierr);}