Commits

Lisandro Dalcin committed 18a5bb0

Fix IGASave() not writting spatial dimension

Comments (0)

Files changed (1)

     }
   }
   if (iga->geometry) {
+    PetscInt nsd;
+    ierr = IGAGetSpatialDim(iga,&nsd);CHKERRQ(ierr);
+    ierr = PetscViewerBinaryWrite(viewer,&nsd,1,PETSC_INT,PETSC_TRUE);CHKERRQ(ierr);
     ierr = IGASaveGeometry(iga,viewer);CHKERRQ(ierr);
   }
   ierr = PetscViewerFlush(viewer);CHKERRQ(ierr);