[petsc-users] Complex vector / matrix to VTK

[Juris Vencels]

Hello,

I want to print vector / matrix to VTK file.

I am running simple tutorial example ex15.c, but since I have compiled 
PETSc with complex scalars (I need them for my program) I get an error 
below. Can I force PETSc to print only real, imaginary (or both) part(s) 
of the solution?

Another question - ksp/ex42.c explicitly writes header information to 
VTK file while snes/ex15.c has only one line - to open VTK file. How 
these two examples are different in terms of writing output to VTK?
src/ksp/ksp/examples/tutorials/ex42.c 
<http://www.mcs.anl.gov/petsc/petsc-current/src/ksp/ksp/examples/tutorials/ex42.c.html>
src/snes/examples/tutorials/ex15.c 
<http://www.mcs.anl.gov/petsc/petsc-current/src/snes/examples/tutorials/ex15.c.html>

Thanks!
Juris

> vencels at bigone:/scratch/petsc-git/src/snes/examples/tutorials$ *./ex15 
> -snes_monitor_short -da_grid_x 20 -da_grid_y 20 -p 1.3 -lambda 1 
> -jtype NEWTON -o ex15.vts*
>   0 SNES Function norm 0.690597
>   1 SNES Function norm 0.437399
>   2 SNES Function norm 0.309415
>   3 SNES Function norm 0.132358
>   4 SNES Function norm 0.0919276
>   5 SNES Function norm 0.0752639
>   6 SNES Function norm 0.0669195
>   7 SNES Function norm 0.0592599
>   8 SNES Function norm 0.0516658
>   9 SNES Function norm 0.0441565
>  10 SNES Function norm 0.0367695
>  11 SNES Function norm 0.0288595
>  12 SNES Function norm 0.0277417
>  13 SNES Function norm 0.0085536
>  14 SNES Function norm 0.00140906
>  15 SNES Function norm 7.35341e-05
>  16 SNES Function norm 3.63992e-07
>  17 SNES Function norm 1.034e-11
> CONVERGED_FNORM_RELATIVE Number of nonlinear iterations = 17
> [0]PETSC ERROR: --------------------- Error Message 
> --------------------------------------------------------------
> *[0]PETSC ERROR: No support for this operation for this object type**
> **[0]PETSC ERROR: Complex values not supported*
> [0]PETSC ERROR: See 
> http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
> [0]PETSC ERROR: Petsc Development GIT revision: v3.5.3-2213-g589572e  
> GIT Date: 2015-03-07 10:56:38 -0600
> [0]PETSC ERROR: ./ex15 on a intelmkl-debug named bigone by vencels Mon 
> Mar 16 10:45:52 2015
> [0]PETSC ERROR: Configure options --with-cc=icc --with-cxx=icpc 
> --with-fc=ifort --with-blas-lapack-dir=/opt/intel/mkl/ 
> --with-fortran-kernels=generic --with-scalar-type=complex --with-mpi=0 
> --download-superlu
> [0]PETSC ERROR: #1 DMDAVTKWriteAll_VTS() line 30 in 
> /scratch/petsc-git/src/dm/impls/da/grvtk.c
> [0]PETSC ERROR: #2 DMDAVTKWriteAll() line 482 in 
> /scratch/petsc-git/src/dm/impls/da/grvtk.c
> [0]PETSC ERROR: #3 PetscViewerFlush_VTK() line 78 in 
> /scratch/petsc-git/src/sys/classes/viewer/impls/vtk/vtkv.c
> [0]PETSC ERROR: #4 PetscViewerFlush() line 30 in 
> /scratch/petsc-git/src/sys/classes/viewer/interface/flush.c
> [0]PETSC ERROR: #5 PetscViewerDestroy() line 101 in 
> /scratch/petsc-git/src/sys/classes/viewer/interface/view.c
> [0]PETSC ERROR: #6 main() line 257 in 
> /scratch/petsc-git/src/snes/examples/tutorials/ex15.c
> [0]PETSC ERROR: PETSc Option Table entries:
> [0]PETSC ERROR: -da_grid_x 20
> [0]PETSC ERROR: -da_grid_y 20
> [0]PETSC ERROR: -jtype NEWTON
> [0]PETSC ERROR: -lambda 1
> [0]PETSC ERROR: -o ex15.vts
> [0]PETSC ERROR: -p 1.3
> [0]PETSC ERROR: -snes_monitor_short
> [0]PETSC ERROR: ----------------End of Error Message -------send 
> entire error message to petsc-maint at mcs.anl.gov----------
> Aborted (core dumped)


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