Duh, I was calling TaoGetKSP _before_ TaoSetFromOptions so that tao_type was still lmvm...
Blaise > > Hmm, which solver are you using? I fixed this for one solver already. Most > solvers do return the KSP at this stage. For solvers that do not use a KSP I > would prefer that it never be created. > > Barry > > On Jul 30, 2014, at 6:19 PM, Blaise A Bourdin <[email protected]> wrote: > >> Responding to myself to clarify my question: >> >> This is with petsc-maint (3.5) up to date >> >> bourdin@galerkin:tutorials (maint)$ git show >> commit 587bed1b5ce23d03610a5e586bde2d77cb223937 >> Merge: ac2cab2 4856ba2 >> Author: Satish Balay <[email protected]> >> Date: Wed Jul 30 09:53:57 2014 -0500 >> >> Merge branch 'barry/check-mpi-include' into maint >> >> I see in the release notes of 3.5.1 a reference to TAOSetUp()being needed, >> but calling TaoSetUp after TAOSetFromOptions but before TaoGetKSP does not >> make any difference. >> >> Breakpoint 1.1, main (argc=3, argv=0x7fff5fbfe9f8) at >> /opt/HPC/petsc-3.5/src/tao/unconstrained/examples/tutorials/eptorsion1.c:148 >> 148 ierr = TaoSetFromOptions(tao);CHKERRQ(ierr); >> (idb) n >> 150 ierr = TaoSetUp(tao); >> (idb) >> 151 ierr = TaoGetKSP(tao,&ksp); >> (idb) >> 153 ierr = KSPView(ksp,PETSC_VIEWER_STDOUT_WORLD); >> (idb) print tao >> $1 = (struct _p_Tao *) 0x10885c860 >> (idb) print ksp >> $2 = (struct _p_KSP *) 0x0 >> >> >> On Jul 30, 2014, at 6:09 PM, Blaise A Bourdin <[email protected]> wrote: >> >>> Hi, >>> >>> TAOGetKSP returns a null pointer when called before TAOSolve (swap the >>> calls in eptorsion1.c for instance). Is this expected? I want to override >>> the default KSP in my application. >>> >>> Blaise >>> -- >>> Department of Mathematics and Center for Computation & Technology >>> Louisiana State University, Baton Rouge, LA 70803, USA >>> Tel. +1 (225) 578 1612, Fax +1 (225) 578 4276 >>> http://www.math.lsu.edu/~bourdin >>> >>> >>> >>> >>> >>> >>> >> >> -- >> Department of Mathematics and Center for Computation & Technology >> Louisiana State University, Baton Rouge, LA 70803, USA >> Tel. +1 (225) 578 1612, Fax +1 (225) 578 4276 >> http://www.math.lsu.edu/~bourdin >> >> >> >> >> >> >> > -- Department of Mathematics and Center for Computation & Technology Louisiana State University, Baton Rouge, LA 70803, USA Tel. +1 (225) 578 1612, Fax +1 (225) 578 4276 http://www.math.lsu.edu/~bourdin
