<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Hi,<br><br>I'm trying to build a differential equation integrator and later a stochastic differential equation integrator.<br><br>I'm having trouble getting f2py to work where the callback itself receives an array from the Fortran routine does some work on it and then passes an array back.  <br><br>For the stoachastic integrator I'll need 2 callbacks both dealing with arrays.<br><br>The idea is the code that never changes (ie the integrator) will be in Fortran and the code that changes (ie the callbacks defining differential equations) will be different for each problem.<br><br>To test the idea I've written basic code which should pass an array back and forth between Python and Fortran if it works right.<br><br>Here is some code which doesn't work properly:-<br><br>SUBROUTINE CallbackTest(dv,v0,Vout,N)<br>    !IMPLICIT
 NONE<br>        <br>cF2PY     intent( hide ):: N<br>    INTEGER:: N, ic<br>                <br>    EXTERNAL:: dv        <br><br>    DOUBLE PRECISION, DIMENSION( N ), INTENT(IN):: v0        <br>    DOUBLE PRECISION, DIMENSION( N ), INTENT(OUT):: Vout<br>            <br>    DOUBLE PRECISION, DIMENSION( N ):: Vnow<br>    DOUBLE PRECISION, DIMENSION( N )::  temp<br>        <br>    Vnow = v0<br>        <br><br>    temp = dv(Vnow, N)<br><br>    DO ic = 1, N<br>       
 Vout( ic ) = temp(ic)<br>    END DO    <br>        <br>END SUBROUTINE CallbackTest<br><br><br><br>When I test it with this python code I find the code just replicates the first term of the array!<br><br><br><br><br>from numpy import *<br>import callback as c<br><br>def dV(v):<br>    print 'in Python dV: V is: ',v<br>    return v.copy()    <br><br>arr = array([2.0, 4.0, 6.0, 8.0])<br><br>print 'Arr is: ', arr<br><br>output = c.CallbackTest(dV, arr)<br><br>print 'Out is: ', output<br><br><br><br><br>Arr is:  [ 2.  4.  6.  8.]<br><br>in Python dV: V is:  [ 2.  4.  6.  8.]<br><br>Out is:  [ 2.  2.  2.  2.]<br><br><br><br>Any ideas how I should do this, and also how do I get the code to work with implicit none not commented out?<br><br>Thanks<br><br>Jon</td></tr></table><br>