Hi all, So in my DMPlex FEM (serial) simulations, I want to output the solution U to a binary file. I have the following:
ierr = PetscViewerBinaryOpen(PETSC_COMM_WORLD,"sol.dat",FILE_MODE_WRITE,&viewer);CHKERRQ(ierr); ierr = VecView(U,viewer);CHKERRQ(ierr); ierr = PetscViewerDestroy(&viewer);CHKERRQ(ierr); The mesh I am solving has 13368 nodes, and 404 of them are constrained using the DMPlexAddBoundary(...) function. However, when I output the vector I only get 12964 values, which I believe is leaving out the 404 constrained ones. How do I get all dofs in the vector U (both free and constrained) printed? Thanks, -- Justin Chang PhD Candidate, Civil Engineering - Computational Sciences University of Houston, Department of Civil and Environmental Engineering Houston, TX 77004 (512) 963-3262
