[petsc-users] Outputting cell data in stead of point data while writing .vtr file

[Smit Thijs]

Hi All,

I am outputting several vectors to a .vtr file successfully for viewing in Paraview. At this moment the information is written to point data. How can I change this and make sure the data is written to cell data?

The code I am currently using for outputting:

PetscViewer viewer;

ierr = PetscViewerVTKOpen(PETSC_COMM_WORLD, "test.vtr", FILE_MODE_WRITE, &viewer);
CHKERRQ(ierr);

ierr = DMView(nd, viewer);
CHKERRQ(ierr);

PetscObjectSetName((PetscObject)xPhys,"xPhys");
ierr = VecView(xPhys, viewer);
CHKERRQ(ierr);

PetscObjectSetName((PetscObject)S,"SvonMises");
ierr = VecView(S, viewer);
CHKERRQ(ierr);

ierr = PetscViewerDestroy(&viewer);
CHKERRQ(ierr);

Best regards,

Thijs Smit

PhD Candidate
ETH Zurich
Institute for Biomechanics

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20210311/a0add8e9/attachment.html>