Christophe Ortiz <[email protected]> writes: >> That does not sound right. Can you send a reproducible test case (best >> as a patch to PETSc)? >> > > What do you mean by patch to PETSc ? How do I do that ?
See "contributing a small patch" (or any of the other workflows). https://bitbucket.org/petsc/petsc/wiki/Home > Sorry about that. Here is the entire error message: > > [0]PETSC ERROR: No support for this operation for this object type! > [0]PETSC ERROR: Mat type mffd! > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: Petsc Release Version 3.4.1, Jun, 10, 2013 > [0]PETSC ERROR: See docs/changes/index.html for recent updates. > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > [0]PETSC ERROR: See docs/index.html for manual pages. > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: ./diffusion_c on a ifort-icc-nompi-double-debug named > mazinger.ciemat.es by u5751 Wed Oct 23 12:46:04 2013 > [0]PETSC ERROR: Libraries linked from > /home/u5751/petsc-3.4.1/ifort-icc-nompi-double-debug/lib > [0]PETSC ERROR: Configure run at Mon Oct 14 05:43:50 2013 > [0]PETSC ERROR: Configure options --with-mpi=0 --with-fc=ifort > --with-cc=icc --with-debugging=1 --with-scalar-type=real > --with-precision=double --with-blas-lapack-dir=/opt/intel/mkl > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: MatZeroEntries() line 5189 in src/mat/interface/matrix.c > [0]PETSC ERROR: FormIJacobian() line 689 in > src/ts/examples/tutorials/diffusion.c This is your code calling MatZeroEntries on the Amat (first Mat argument). The convention is to assemble into Pmat (the preconditioning matrix; second Mat argument).
pgpixFexq6OBb.pgp
Description: PGP signature
