<div dir="ltr">Thanks for your fix Hong.<br></div><div class="gmail_extra"><br><div class="gmail_quote">2016-05-03 17:48 GMT+02:00 Hong <span dir="ltr"><<a href="mailto:hzhang@mcs.anl.gov" target="_blank">hzhang@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 class="gmail_extra"><div class="gmail_quote">Alejandro :</div><div class="gmail_quote">Fixed</div><div class="gmail_quote"><a href="https://bitbucket.org/petsc/petsc/commits/0043d7fe17a5a8a1e4768ab560023392b3862352" target="_blank">https://bitbucket.org/petsc/petsc/commits/0043d7fe17a5a8a1e4768ab560023392b3862352</a></div><div class="gmail_quote">It will merged to petsc-release after our nightly tests.</div><div class="gmail_quote"><br></div><div class="gmail_quote">Your detailed report is very helpful, thanks!</div><span class="HOEnZb"><font color="#888888"><div class="gmail_quote"><br></div><div class="gmail_quote">Hong</div></font></span><div><div class="h5"><div class="gmail_quote"><br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div>Hi,<br><br></div>After updating petsc from version 3.6.2 to 3.7, I had the following problem:<br><br></div>The result that I get from doing something like MatMult(schurComplement,v1,v2); differs between the petsc versions that I mentioned. I traced back the problem to the function MatGetSubMatrices_SeqSBAIJ(). What I did in order to temporarily solve my problem was to use in petsc 3.7 the code of MatGetSubMatrices_SeqSBAIJ() corresponding to petsc 3.6.2. <div><div><br></div><div></div><div>How I solved the problem is not a real solution, but it allows me to keep working. My question is: is it a bug or do you think that there could be something wrong with my code from which I use petsc?<br><br></div><div>I think it should be a bug. If I run the ksp example 59 which uses FETI-DP (whose implementation also uses something like MatMult(schurComplement,v1,v2)) with petsc 3.7, I get the following unexpected answer<br></div><div><br>mpiexec -n 4 ./ex59 -npx 2 -npy 2 -nex 60 -ney 60 -ksp_max_it 3<br>------------------FETI-DP stats-------------------------------<br>Number of degrees of freedom               :      355 <br>Number of iterations                       :        3 <br>Eigenvalues preconditioned operator        : 1.49e+00 8.86e+00<br>Error betweeen exact and computed solution : 1.91e+00<br>--------------------------------------------------------------<br><br>If I run the same problem with petsc 3.7 but with MatGetSubMatrices_SeqSBAIJ() corresponding to petsc 3.6.2, I get the expected answer<br><br>mpiexec -n 4 ./ex59 -npx 2 -npy 2 -nex 60 -ney 60 -ksp_max_it 3<br>------------------FETI-DP stats-------------------------------<br>Number of degrees of freedom               :      355 <br>Number of iterations                       :        3 <br>Eigenvalues preconditioned operator        : 1.08e+00 2.08e+00<br>Error betweeen exact and computed solution : 3.60e-02<br>--------------------------------------------------------------<br><br></div><div><br></div><div>Thanks,<br></div><div>Alejandro <br></div></div></div>
</blockquote></div><br></div></div></div></div>
</blockquote></div><br></div>