<div dir="ltr"><div>Thanks Jed, sorry I new using Petsc.<br><br></div>I am runing my program  with default options :<br><div>mpirun -np 1 ./a.out<br>runing the program with a few points (10) the interpolation is aceptable,   <br>
<br>using the monitor options ie.<br>mpirun -np 1 ./a.out -ksp_converged_reason -ksp_monitor_true_residual<br>the output is the next<br> 0 KSP preconditioned resid norm 7.415666989515e+05 true resid norm 1.510533668194e+04 ||r(i)||/||b|| 1.000000000000e+00<br>
  1 KSP preconditioned resid norm 4.935801698308e+02 true resid norm 1.721762948997e+01 ||r(i)||/||b|| 1.139837519183e-03<br>  2 KSP preconditioned resid norm 1.038493428107e-01 true resid norm 9.406615966569e+00 ||r(i)||/||b|| 6.227346112594e-04<br>
Linear solve converged due to CONVERGED_RTOL iterations 2<br><br></div><div>But if encrease the number of points to interpolate, the result is incorrect and the time of convergece encrease, the output in this case is: <br>
</div><div><br>.....<br>138 KSP preconditioned resid norm 4.897817095518e-06 true resid norm 5.455886375264e-08 ||r(i)||/||b|| 4.852601039969e-10<br>139 KSP preconditioned resid norm 3.535973415166e-06 true resid norm 3.717900878947e-08 ||r(i)||/||b|| 3.306793512687e-10<br>
140 KSP preconditioned resid norm 2.450794339359e-06 true resid norm 3.009597005096e-08 ||r(i)||/||b|| 2.676810430480e-10<br>Linear solve converged due to CONVERGED_RTOL iterations 140<br>[0] time (sec) 89.940000<br><br>the  code where the solution is implemented is the next:<br>
KSP ksp;<br>KSPCreate(PETSC_COMM_WORLD,&ksp);<br>KSPSetOperators(ksp,A,A,DIFFERENT_NONZERO_PATTERN);<br>KSPSetTolerances(ksp,1.e-5/(N+nt+1),1.e-50,PETSC_DEFAULT,PETSC_DEFAULT);<br>KSPSetFromOptions(ksp);<br>KSPSolve(ksp,b,x);<br>
<br></div><div>Let me know if you need another information. <br><br></div><div>Regard Alex <br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/1/31 Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>"result is incorrect" is not helpful. Please include diagnostics. Are you using an iterative or direct solver? Start here if you don't know what to check:</div>
<br><a href="http://scicomp.stackexchange.com/questions/513/why-is-my-iterative-linear-solver-not-converging" target="_blank">http://scicomp.stackexchange.com/questions/513/why-is-my-iterative-linear-solver-not-converging</a><div>
<div class="h5"><br>
<div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Jan 31, 2013 at 9:58 AM, Alex Salazar <span dir="ltr"><<a href="mailto:alejandrocof@gmail.com" target="_blank">alejandrocof@gmail.com</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello, everyone <br><br>I want to get the solution of algrebraic system (with a  dense zero-diagonal simetric matrix) of a interpolation problem (using radial basis functions) with Petsc.<br>

<br> For a few points (around 10 elements) using KSP I get the correct solution,  however when  I increase the number of points,  the result is incorrect.<br>
<br>Any Idea ?<br><br>Regards Alejandro. <br><br></div>
</blockquote></div><br></div></div></div></div>
</blockquote></div><br></div>