<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">Alejandro:</div><div class="gmail_quote">I can reproduce this error. I'm investigating it.</div><div class="gmail_quote"><br></div><div class="gmail_quote">I made minor update of MatGetSubMatrices_SeqSBAIJ() in v3.7 which only removed requirement of isrow=iscol...do not know what leads to this error yet.</div><div class="gmail_quote"><br></div><div class="gmail_quote">Hong</div><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>