<div dir="ltr">On Wed, Mar 27, 2013 at 3:13 PM, Anil . <span dir="ltr"><<a href="mailto:dasans@gmail.com" target="_blank">dasans@gmail.com</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Matt,<div><br></div><div>Petsc Options I use to run are....</div><div>mpirun -np 4 ./reader -pc_type asm -sub_pc_type lu -sub_mat_type dense -ksp_monitor -ksp_rtol 1e-13 -ksp_max_it 100 -vecscatter_alltoall -log_summary<br>
</div></div></blockquote><div><br></div><div style>1) Always run with -ksp_view.</div><div style><br></div><div style>2) The relative tolerance is probably too tight, but that is secondary</div><div style><br></div><div style>
3) Something is really wrong here. I am guessing something in the input not what you want. If the interaction is</div><div style> truly short range, you would see significant drop in the residual on the first iteration. First, take a look at the</div>
<div style> matrix using -mat_view draw:: -draw_pause -1. It should be banded.</div><div style><br></div><div style> Matt</div><div style><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div>
</div><div>Attached is the output..It also contains my petsc configuration</div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Mar 27, 2013 at 7:59 AM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@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"><div class="im"><div>On Wed, Mar 27, 2013 at 11:58 AM, Anil . <span dir="ltr"><<a href="mailto:dasans@gmail.com" target="_blank">dasans@gmail.com</a>></span> wrote:<br>
</div></div><div class="gmail_extra"><div class="gmail_quote"><div class="im"><div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Matt,<div><br></div><div>I am having around 3481 particles that are placed in an unstructured manner.</div><div>Attached is the image showing the distribution.</div></div></blockquote><div><br></div></div>
</div><div>
Show me your PETSc options, and try playing with the number of blocks. If you look</div><div>at the PetRBF paper, we give guidance for choosing the sizes.</div><div><br></div><div> Matt</div><div><div class="h5"><div><div>
<div> </div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Mar 26, 2013 at 11:45 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@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"><div>On Mon, Mar 25, 2013 at 10:38 PM, Anil . <span dir="ltr"><<a href="mailto:dasans@gmail.com" target="_blank">dasans@gmail.com</a>></span> wrote:<br></div><div class="gmail_extra"><div class="gmail_quote">
<div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">1) Could not find the petrbf mailing list<div>2) Petrbf runs perfectly</div><div>3) Attached is the output with -ksp_view -ksp_monitor</div><div><br></div><div>Just point me in the right direction. Issues might be very basic as I am starting to use Petsc</div>
</div></blockquote><div><br></div></div><div>This output is a little strange. Some partitions have 0 entries. I am guessing this problem is very</div><div>small. For PeRBF, it does turn out to be optimal to use small blocks, but the block size depends</div>
<div>on your interaction scale. Right now you have 75 blocks, which might be too many for your small</div><div>problem.</div><div><br></div><div> Matt</div><div><div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="gmail_extra"><div class="gmail_quote">On Sat, Mar 23, 2013 at 2:18 AM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@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"><div><div>On Fri, Mar 22, 2013 at 10:58 PM, Anil . <span dir="ltr"><<a href="mailto:dasans@gmail.com" target="_blank">dasans@gmail.com</a>></span> wrote:<br>
</div></div><div class="gmail_extra"><div class="gmail_quote"><div><div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Hi,<div><br></div><div>I have a text file containing N rows.</div><div>Each row with x,y,omega values.</div><div>I am trying to interpolate this data onto a regular grid using petrbf</div><div>But the KSP does not converge and am not able to find the reason.</div>
<div><br></div><div>The code is available with the text files at</div><div><a href="https://www.dropbox.com/s/cypuwugbxo07kx0/rbf-interpolation.tar.gz" target="_blank">https://www.dropbox.com/s/cypuwugbxo07kx0/rbf-interpolation.tar.gz</a><br>
</div><div><br></div><div>I am very new to petsc and any direction how o proceed would be helpful.</div></div></blockquote><div><br></div></div></div><div>1) Did you mail the petrbf list?</div><div><br></div><div>2) Could you run the petrbf examples?</div>
<div><br></div><div>3) We cannot tell anything about convergence without the output of -ksp_view -ksp_monitor.</div><div><br></div><div> Matt</div><div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><span><font color="#888888"><div>-- <br><div>Sincerely</div>Anil Das P V<br>
</div></font></span></div>
</blockquote></div></div><span><font color="#888888"><br><br clear="all"><span><font color="#888888"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</font></span></font></span></div></div><span><font color="#888888">
</font></span></blockquote></div><span><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div>Sincerely</div>Anil Das P V<br>
</font></span></div>
</blockquote></div></div></div><div><div><br><br clear="all"><span><font color="#888888"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</font></span></div></div></div></div><span><font color="#888888">
</font></span></blockquote></div><span><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div>Sincerely</div>Anil Das P V<br>
</font></span></div>
</blockquote></div></div></div></div></div><div><div class="h5"><div><div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div></div></div></div></div>
</blockquote></div><div><div class="h5"><br><br clear="all"><div><br></div>-- <br><div>Sincerely</div>Anil Das P V<br>
</div></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>