[petsc-users] Suggestions for code with dof >> 1 ?

[Jed Brown]

Christophe Ortiz <christophe.ortiz at ciemat.es> 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).
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20131023/66a0ca3b/attachment.pgp>