Simple example of f2py

Author: jussienko

demos/interface/fortran/README.md

@@ -0,0 +1,8 @@
+# Building the extension
+
+Fortran source can be built into a C-extension as
+```
+f2py -c add.f90 -m add
+```
+(For Python 3 use f2py3)
+

demos/interface/fortran/add.f90

@@ -0,0 +1,12 @@
+subroutine add(a, b, n)
+
+  implicit none
+
+  ! f2py understands only limited number of kind parameters
+  real(kind=8), intent(inout) :: a(n)
+  real(kind=8), intent(inout) :: b(n)
+  integer :: n
+
+  a = a + b
+
+end subroutine

demos/interface/fortran/example.py

@@ -0,0 +1,8 @@
+import add
+import numpy as np
+
+a = np.random.random(10)
+b = np.ones_like(a)
+print(a)
+add.add(a, b)
+print(a)