On Thu, May 10, 2012 at 3:45 PM, Chetan Jhurani <span dir="ltr"><<a href="mailto:chetan.jhurani@gmail.com" target="_blank">chetan.jhurani@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi wizards of "->ops->" indirections,<br>
<br>
I'm trying to use bcgs without preconditioning, for now, and<br>
the iterations using "-vec_type cusp -mat_type mpiaijcusp" don't<br>
match serial or non-GPU options.  I've attached the test program<br>
and the 4 outputs (serial/parallel + CPU/GPU).  All this is<br>
with petsc-dev downloaded now and real scalars.<br>
<br>
Only the parallel GPU results are different starting from<br>
third residual norm seen in results.txt.  The other three match<br>
one another.  Am I doing something wrong?<br>
<br>
fbcgs (bcgs with -ksp_bcgs_flexible) works fine with all the<br>
serial/parallel or CPU/GPU options I've tried.<br>
<br>
Let me know if you need the matrix, rhs, and initial guess<br>
binary files that are read in by the test program.<br></blockquote><div><br></div><div>That would be great. This looks like a bug that should be tracked down.</div><div><br></div><div>    MAtt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Thanks,<br>
<br>
Chetan<br>
<br>
</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<br>