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 <[email protected]> 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@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 > [email protected] > https://lists.sourceforge.net/lists/listinfo/xolotl-psi-development
