<html><head></head><body><div style="font-family: Verdana;font-size: 12.0px;"><div>MatGetSubmatrix is in the current PETSc released called MatCreateSubmatrix. Is this correct ?
<div>
<div name="quote" style="margin:10px 5px 5px 10px; padding: 10px 0 10px 10px; border-left:2px solid #C3D9E5; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div style="margin:0 0 10px 0;"> </div>
<div name="quoted-content">
<div>
<div class="gmail_quote">
<div>On Mon, Jul 9, 2018 at 8:38 PM Marius Buerkle <<a href="mailto:mbuerkle@web.de" onclick="parent.window.location.href='mbuerkle@web.de'; return false;" target="_blank">mbuerkle@web.de</a>> wrote:</div>
<blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex;border-left: 1.0px rgb(204,204,204) solid;padding-left: 1.0ex;">
<div>
<div style="font-family: Verdana;font-size: 12.0px;">
<div> I see. What I want to do is to calculate the matrix product C=A*B' between two sparse matrices A=(A11 0 , A21 0) and B=(B11 0 , B21 0) where C will be dense in the end but I just want to calculate some selected entries C_ij of C. At the moment I extract submatricies for the corresponding rows and columns,</div>
</div>
</div>
</blockquote>
<div> </div>
<div>I think that is the right way, but you should only need MatGetSubmatrix for that.</div>
<div> </div>
<div> Thanks,</div>
<div> </div>
<div> Matt</div>
<div> </div>
<blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex;border-left: 1.0px rgb(204,204,204) solid;padding-left: 1.0ex;">
<div>
<div style="font-family: Verdana;font-size: 12.0px;">
<div>so I was wondering if there is a simpler or performancer-wise faster way. I assume there is not such thing as a restricted MatMatMul which just calculated the lets say predefined nonzero entries of C.
<div style="margin: 10.0px 5.0px 5.0px 10.0px;padding: 10.0px 0 10.0px 10.0px;border-left: 2.0px solid rgb(195,217,229);">
<div style="margin: 0 0 10.0px 0;"> </div>
<div><br/>
These are "internal" routines that we rarely expect end users to use since they are specific for particular matrix implementations. As such they are also kind of strange to use from Fortran since none of the AIJ data structures can be made visible to Fortran.<br/>
<br/>
Could you explain why you want them from Fortran and maybe we'll have alternative suggestions on how you can achieve the same effect.<br/>
<br/>
Barry<br/>
<br/>
<br/>
> On Jul 5, 2018, at 3:05 AM, Marius Buerkle <<a href="mailto:mbuerkle@web.de" onclick="parent.window.location.href='mbuerkle@web.de'; return false;" target="_blank">mbuerkle@web.de</a>> wrote:<br/>
><br/>
> or MatMPIAIJGetLocalMatCondensed for that matter.<br/>
><br/>
><br/>
><br/>
> Hi !<br/>
><br/>
> Is MatMPIAIJGetSeqAIJ implemented for fortran?<br/>
><br/>
> best,<br/>
> Marius<br/>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<div> </div>
--
<div class="gmail_signature">
<div>
<div>
<div>
<div>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><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div></div></body></html>