Commits

Nathan Collier  committed 393cc01

update convergence test, slowdown observed

  • Participants
  • Parent commits 2322d00

Comments (0)

Files changed (2)

File test/ConvTest.c

 {
   PetscInt nen,dim;
   IGAPointGetSizes(p,0,&nen,0);
-  IGAPointGetDims(p,&dim,0);
+  IGAPointGetDims(p,&dim,0,0);
 
   const PetscReal *N0;
   IGAPointGetBasisFuns(p,0,&N0);
 {
   PetscInt nen,dim,i;
   IGAPointGetSizes(p,0,&nen,0);
-  IGAPointGetDims(p,&dim,0);
+  IGAPointGetDims(p,&dim,0,0);
   
   PetscInt Nb[3] = {0,0,0};
-  Nb[0] = p->parent->BD[0]->nnp;
-  Nb[1] = p->parent->BD[1]->nnp;
-  Nb[2] = p->parent->BD[2]->nnp;
+  Nb[0] = p->parent->parent->axis[0]->nnp;
+  Nb[1] = p->parent->parent->axis[1]->nnp;
+  Nb[2] = p->parent->parent->axis[2]->nnp;
   
   const PetscReal *N0,(*N2)[dim][dim];
   IGAPointGetBasisFuns(p,0,(const PetscReal**)&N0);
   IGAPointFormValue(p,U,&ua);
 
   PetscInt i,dim;
-  IGAPointGetDims(p,&dim,0);
+  IGAPointGetDims(p,&dim,0,0);
   for (i=0; i<dim; i++) ue *= sin(2*CONST*PETSC_PI*p->point[i]);
 
   PetscReal error = PetscAbsScalar(ua-ue);

File test/ConvTest.py

         r = int(round(r))
         assert p+1==r
 dt = time.time()-t0
+print "time %.3f (%.3f)" % (dt,143.299206018)
 assert abs(dt-143.299206018)/143.299206018 < 0.1
-print "time %.3f (%.3f)" % (dt,143.299206018)
+