<br><br><div class="gmail_quote">On Mon, Mar 28, 2011 at 10:44 PM, Sturla Molden <span dir="ltr"><<a href="mailto:sturla@molden.no">sturla@molden.no</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Den 28.03.2011 19:12, skrev Pearu Peterson:<br>
<div class="im">><br>
> FYI, f2py in numpy 1.6.x supports also assumed shape arrays.<br>
<br>
</div>How did you do that? Chasm-interop, C bindings from F03, or marshalling<br>
through explicit-shape?<br></blockquote><div><br>The latter.<br> Basically, if you have<br><br>subroutine foo(a)<br>real a(:)<br>end<br><br>then f2py automatically creates a wrapper subroutine<br><br>subroutine wrapfoo(a, n)<br>
real a(n)<br>integer n<br>!f2py intent(in) :: a<br>!f2py intent(hide) :: n = shape(a,0)<br>interface<br>subroutine foo(a)<br>real a(:)<br>end<br>end interface<br>call foo(a)<br>end<br><br>that can be wrapped with f2py in ordinary way.<br>
 <br></div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Can f2py pass strided memory from NumPy to Fortran?<br>
<div><div class="h5"><br></div></div></blockquote><div><br>No. I haven't thought about it.<br><br>Pearu<br><br></div></div>