<div dir="ltr">Hi Corrado,<div><br></div><div>The TRON solver should be using the STCG as a default ksp type in order to get trust region information from the KSP. However, you are correct that it should at least be able to run using a different KSP without bailing out. I will get this fixed.</div>
<div><br></div><div>Jason</div><div><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Mar 13, 2014 at 11:17 AM, Corrado Maurini <span dir="ltr"><<a href="mailto:corrado.maurini@upmc.fr" target="_blank">corrado.maurini@upmc.fr</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">



<div style="word-wrap:break-word">
<font face="Menlo"><span style="font-size:11px">Hi,</span></font>
<div><font face="Menlo"><span style="font-size:11px"><br>
</span></font></div>
<div><font face="Menlo"><span style="font-size:11px">Running in petsc/src/tao/bound/examples/tutorials <br>
<br>
./jbearing2 -tao_view -tao_type tron -ksp_type gmres<br>
<br>
I get the error below. With tao 2.2 and petsc 3.4 and the same options, I did not get any error. </span></font>
<div><font face="Menlo"><span style="font-size:11px"><br>
</span></font></div>
<div><font face="Menlo"><span style="font-size:11px">Now the TAO TRON method seems to work only with STCG as KSP. Is there any reason? </span><br>
<br>
<span style="font-size:11px">Corrado</span><br>
<br>
———</font></div>
<div><font face="Menlo"><br>
<span style="font-size:11px"> ./jbearing2 -tao_view -tao_type tron -ksp_type gmres<br>
<br>
---- Journal Bearing Problem SHB-----<br>
mx: 50,  my: 50,  ecc: 0.1<br>
<br>
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>
[0]PETSC ERROR: No support for this operation for this object type<br>
[0]PETSC ERROR: Cannot locate function KSPSTCGSetRadius_C in object<br>
[0]PETSC ERROR: See <a href="http://http://www" target="_blank">http://http://www</a>.<a href="http://mcs.anl.gov/petsc/documentation/faq.html" target="_blank">mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.<br>

[0]PETSC ERROR: Petsc Development GIT revision: v3.4.3-3575-gf3e7060  GIT Date: 2014-03-11 19:24:46 -0500<br>
[0]PETSC ERROR: ./jbearing2 on a arch-darwin-cxx-opt named muspratt by Maurini Thu Mar 13 17:16:30 2014<br>
[0]PETSC ERROR: Configure options --download-blacs=1 --download-hypre=1 --download-metis=1 --download-ml=1 --download-mumps=1 --download-parmetis=1 --download-superlu_dist=1 --download-umfpack=1 --with-c-support=1 --with-clanguage=cxx --with-debugging=0 --with-scalapack-dir=/usr/local/
 --with-scotch-dir=/usr/local/ --with-shared-libraries=1 --with-x11=1 COPTFLAGS=-O2 PETSC_ARCH=arch-darwin-cxx-opt<br>
[0]PETSC ERROR: #1 KSPSTCGSetRadius() line 37 in /opt/HPC/petsc/src/ksp/ksp/impls/cg/stcg/stcg.c<br>
[0]PETSC ERROR: #2 TaoSolve_TRON() line 170 in /opt/HPC/petsc/src/tao/bound/impls/tron/tron.c<br>
[0]PETSC ERROR: #3 TaoSolve() line 192 in /opt/HPC/petsc/src/tao/interface/taosolver.c<br>
[0]PETSC ERROR: #4 main() line 169 in /opt/HPC/petsc/src/tao/bound/examples/tutorials/jbearing2.c<br>
[0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov----------<br>
--------------------------------------------------------------------------<br>
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD<br>
with errorcode 56.<br>
<br>
NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.<br>
You may or may not see output from other processes, depending on<br>
exactly when Open MPI kills them.<br>
--------------------------------------------------------------------------</span></font></div>
</div>
</div>

</blockquote></div><br></div></div>