Dear,
 
I used to save my Vec results to vts files. Unfortunately this is not working anymore. My Code:

// Creating 3DDM before
ierr = DMDACreateNaturalVector(dynamic_Maxwell, &dynamic_fields_global);CHKERRQ(ierr);
ierr = VecSet(dynamic_fields_global,some_number);CHKERRQ(ierr);

PetscViewer viewer;
ierr = PetscViewerVTKOpen(PETSC_COMM_WORLD, "test.vts", FILE_MODE_WRITE, &viewer);CHKERRQ(ierr);
ierr = VecView(dynamic_fields_global, viewer);CHKERRQ(ierr);
ierr = PetscViewerDestroy(&viewer);CHKERRQ(ierr);

Saving to ASCII still works:
ierr = PetscViewerASCIIOpen(PETSC_COMM_WORLD,"test.vts",&viewer);CHKERRQ(ierr);

I don't get any errors.

 

king regrads,
Grischa

Reply via email to