<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Oct 1, 2013 at 2:01 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> writes:<br>
> I put the example with VecGetArray() into petsc-dev and if you send<br>
> me one with DMDAVecGetArray() I'll add that as well. Since the<br>
> derived type has to appear explicitly in the code I can't put<br>
> something in the library only in examples that people can copy and<br>
> change.<br>
<br>
</div>This is the sad state of Fortran. Alternatively, it should be possible<br>
to call the C DMDAVecGetArray via the ISO C bindings in Fortran 2003.<br>
<br></blockquote><div><br></div><div>I will check that. I have the book fortran 95/2003 by Metcalf and there is a chapter dedicated to iso c bindings.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Yet another example that the best way to call Fortran libraries from<br>
Fortran is to use the ISO C bindings.<br>
</blockquote></div><br></div></div>