[petsc-dev] -snes_mf_operator fails with TS

[Barry Smith]

   It is extremely frustrating that something like -snes_mf_operator which has been around since the beginning of time doesn’t work with TS

   Barry

On Nov 11, 2013, at 10:46 AM, Jay J. Billings <billingsjj at ornl.gov> wrote:

> Barry,
> 
> Adding the -snes_mf_operator option fails with the following output
> 
> 0]PETSC ERROR: --------------------- Error Message 
> ------------------------------------
> [0]PETSC ERROR: No support for this operation for this object type!
> [0]PETSC ERROR: Not written for this matrix type!
> [0]PETSC ERROR: 
> ------------------------------------------------------------------------
> [0]PETSC ERROR: Petsc Development GIT revision: 
> a965ca046084fa53248a41da989a0a62cb6266ea  GIT Date: 2013-11-10 14:25:12 
> -0600
> [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: ./xolotl on a arch-linux2-c-opt named 
> antecessor.ornl.gov by bkj Mon Nov 11 11:45:31 2013
> [0]PETSC ERROR: Libraries linked from /opt/petsc-latest_mpich-3.0.1/lib
> [0]PETSC ERROR: Configure run at Mon Nov 11 09:36:40 2013
> [0]PETSC ERROR: Configure options --prefix=/opt/petsc-latest_mpich-3.0.1 
> --with-cc=mpicc --with-cxx=mpicxx --with-fc=mpif77 --with-debugging=no 
> --download-f-blas-lapack=1 --FOPTFLAGS= --with-shared-libraries=1 
> --download-hypre=yes --with-debugging=0
> [0]PETSC ERROR: 
> ------------------------------------------------------------------------
> [0]PETSC ERROR: MatDuplicate() line 4125 in src/mat/interface/matrix.c
> [0]PETSC ERROR: TSGetRHSMats_Private() line 583 in src/ts/interface/ts.c
> [0]PETSC ERROR: TSComputeIJacobian() line 760 in src/ts/interface/ts.c
> [0]PETSC ERROR: SNESTSFormJacobian_ARKIMEX() line 1057 in 
> src/ts/impls/arkimex/arkimex.c
> [0]PETSC ERROR: SNESTSFormJacobian() line 3543 in src/ts/interface/ts.c
> [0]PETSC ERROR: SNESComputeJacobian() line 2241 in src/snes/interface/snes.c
> [0]PETSC ERROR: SNESSolve_NEWTONLS() line 231 in src/snes/impls/ls/ls.c
> [0]PETSC ERROR: SNESSolve() line 3809 in src/snes/interface/snes.c
> [0]PETSC ERROR: TSStep_ARKIMEX() line 777 in src/ts/impls/arkimex/arkimex.c
> [0]PETSC ERROR: TSStep() line 2603 in src/ts/interface/ts.c
> [0]PETSC ERROR: TSSolve() line 2728 in src/ts/interface/ts.c
> [0]PETSC ERROR: checkPetscError() line 61 in 
> /home/bkj/research/xolotl/xolotl_workspace/xolotl-Source at xolotl/xolotlSolver/PetscSolver.cpp
> 
> 
> Jay
> 
> ------------------------------------------------------------------------------
> November Webinars for C, C++, Fortran Developers
> Accelerate application performance with scalable programming models. Explore
> techniques for threading, error checking, porting, and tuning. Get the most 
> from the latest Intel processors and coprocessors. See abstracts and register
> http://pubads.g.doubleclick.net/gampad/clk?id=60136231&iu=/4140/ostg.clktrk
> _______________________________________________
> Xolotl-psi-development mailing list
> Xolotl-psi-development at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/xolotl-psi-development