<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sun, Sep 24, 2017 at 8:57 PM, Adrian Croucher <span dir="ltr"><<a href="mailto:a.croucher@auckland.ac.nz" target="_blank">a.croucher@auckland.ac.nz</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class=""><br>
<br>
On 25/09/17 13:37, Adrian Croucher wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
Actually it looks like DMGetDefaultSF() does already work in Fortran though, and could be used to do the same thing.<br>
<br>
Do you think that would be a reasonable way to do it?<br>
<br>
</blockquote>
<br></span>
Hmm, I just tried it and it looks like the default SF produced by this function doesn't have any overlap. So that won't work for me.<br>
<br>
Guess I might need a Fortran interface for PetscSFSetGraph/PetscSFGetGrap<wbr>h after all. If you could put one it that would be great.<br>
<br>
Another simpler possibility that would work for my case is if there was a DMCopySF() function for copying the point SF from one DM to another. Up to you if you think that would be of enough general usefulness to warrant inclusion.</blockquote><div><br></div><div>If you truly need the exact same SF for your grid, you should be able to use DMGet/SetPointSF() since it will just reference count it for you. Then</div><div>the default SF is created automatically from the point SF. Does this work?</div><div><br></div><div>  Thanks,</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"><div class="HOEnZb"><div class="h5">
- Adrian<br>
<br>
-- <br>
Dr Adrian Croucher<br>
Senior Research Fellow<br>
Department of Engineering Science<br>
University of Auckland, New Zealand<br>
email: <a href="mailto:a.croucher@auckland.ac.nz" target="_blank">a.croucher@auckland.ac.nz</a><br>
tel: <a href="tel:%2B64%20%280%299%20923%204611" value="+6499234611" target="_blank">+64 (0)9 923 4611</a><br>
<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><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><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>