Tworzenie python i fortran przyjaciele

Question

Załóżmy, że musimy wywołać funkcję Fortran, która zwraca niektóre wartości, w programie Pythona. I okazało się, że przepisanie kodu FORTRAN w taki sposób:

subroutine pow2(in_x, out_x) 
     implicit none 
     real, intent(in)  :: in_x 
!f2py real, intent(in, out) :: out_x 
     real, intent(out)  :: out_x 
     out_x = in_x ** 2 
     return 
end 

i nazywając go w Pythonie w taki sposób:

import modulename 
a = 2.0 
b = 0.0 
b = modulename.pow2(a, b) 

daje nam wynik roboczą. Czy mogę inaczej nazwać funkcję fortrana, bo myślę, że pierwszy sposób jest trochę niezgrabny?

Answer 1

myślę, że wystarczy zmienić swój podpis funkcji f2py lekko (tak, że out_x jest tylko intent(out) i in_x jest tylko intent(in)):

subroutine pow2(in_x, out_x) 
    implicit none 
    real, intent(in) :: in_x 
    !f2py real, intent(in) :: in_x 
    real, intent(out)  :: out_x 
    !f2py real, intent(out) :: out_x 
    out_x = in_x ** 2 
    return 
end subroutine pow2 

Teraz skompilować:

f2py -m test -c test.f90 

Teraz uruchom:

>>> import test 
>>> test.pow2(3) #only need to pass intent(in) parameters :-) 
9.0 
>>> 

---

Należy zauważyć, że w tym przypadku, f2py jest w stanie poprawnie zeskanować podpis funkcji bez specjalne !f2py komentarze:

!test2.f90 
subroutine pow2(in_x, out_x) 
    implicit none 
    real, intent(in) :: in_x 
    real, intent(out)  :: out_x 
    out_x = in_x ** 2 
    return 
end subroutine pow2 

kompilacji:

f2py -m test2 -c test2.f90 

run:

>>> import test2 
>>> test2.pow2(3) #only need to pass intent(in) parameters :-) 
9.0 
>>> 

Answer 2

Kolejna korzyść dla uniknięcia ding intent(inout) argumenty są (z kilkoma innymi ograniczeniami) wynikową funkcją can be considered PURE. Jest to związane z brakiem efektów ubocznych podjętych przez języki funkcjonalne i opłaca się z lepszą optymalizacją i wykrywaniem błędów przez kompilator Fortranu i jest szczególnie ważne dla autokalaryzacji.

Answer 3

Jeśli używasz ipython, spróbuj magiczne komendy, co sprawia, że ​​w stanie napisać

%%fortran 
subroutine pow2(in_x, out_x) 
    implicit none 
    real, intent(in) :: in_x 
    real, intent(out)  :: out_x 
    out_x = in_x ** 2 
    return 
end subroutine pow2 

a następnie po prostu użyć tej funkcji w kodzie (bez dodatkowego importu).