[petsc-users] SNES test with a nested matrix

Norihiro Watanabe norihiro.w at gmail.com
Thu Apr 17 09:31:32 CDT 2014


Hi, I wonder whether "-snes_type test" option can work also with a nested
matrix. I got the following error in my program and am not sure if there is
a problem in my codes or just PETSc doesn't support it.


Testing hand-coded Jacobian, if the ratio is
O(1.e-8), the hand-coded Jacobian is probably correct.
Run with -snes_test_display to show difference
of hand-coded and finite difference Jacobian.
[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: No support for this operation for this object type!
[0]PETSC ERROR: Mat type nest!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.4.3, Oct, 15, 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: Libraries linked from
/tools/petsc/petsc-3.4.3/linux-gnu-debug/lib
[0]PETSC ERROR: Configure run at Wed Jan 15 15:48:37 2014
[0]PETSC ERROR: Configure options PETSC_ARCH=linux-gnu-debug
--download-f2cblaslapack=1 -with-debugging=1 --download-superlu_dist
--download-hypre=1 --download-ml=1 --download-parmetis --download-metis
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: MatSetValues() line 1085 in
/tools/petsc/petsc-3.4.3/src/mat/interface/matrix.c
[0]PETSC ERROR: SNESComputeJacobianDefault() line 121 in
/tools/petsc/petsc-3.4.3/src/snes/interface/snesj.c
[0]PETSC ERROR: SNESSolve_Test() line 64 in
/tools/petsc/petsc-3.4.3/src/snes/impls/test/snestest.c
[0]PETSC ERROR: SNESSolve() line 3636 in
/tools/petsc/petsc-3.4.3/src/snes/interface/snes.c


Thank you in advance,
Nori
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140417/298123ef/attachment.html>


More information about the petsc-users mailing list