Fortran嵌入Python

利用forpy库，可以在Fortran程序内嵌入Python.

program test_forpyuse forpy_modimplicit noneinteger::ierrierr=forpy_initialize()!!---------------!!list!!---------------blocktype(list)::my_listierr=list_create(my_list)ierr=my_list%append(19)ierr=my_list%append("Hello world")ierr=my_list%append(3.14d0)ierr=print_py(my_list)call my_list%destroyend block!!----------------!!dict!!----------------blocktype(dict)::direal::a_valueierr=dict_create(di)ierr=di%setitem("temperature",273.0)ierr=di%setitem("pressure",1013.0)ierr=di%getitem(a_value,"pressure")write(*,*) 'pressure= ',a_valueierr=di%getitem(a_value,'does not exist')if(ierr/=0) thenif(exception_matches(KeyError)) thenwrite(*,*) 'Key not found....'call err_clearelsewrite(*,*) 'Unknown error.....'stopendifendif!default valueierr=di%get(a_value,"volumne",1.0)write(*,*) 'volume= ',a_valuecall di%destroy()endblock!!----------------!!tuple!!----------------blocktype(tuple)::tutype(object)::iteminteger::tu_leninteger::iinteger::int_valuecharacter(len=:),allocatable::str_valueierr=tuple_create(tu,4)ierr=tu%setitem(0,17)ierr=tu%setitem(1,"hello")ierr=tu%setitem(2,23)ierr=tu%setitem(3,"world")ierr=tu%len(tu_len)do i=0,tu_len-1 !!python 0-basedierr=tu%getitem(item,i)if(is_int(item))thenierr=cast(int_value,item)write(*,*) int_valueelse if(is_str(item)) thenierr=cast(str_value,item)write(*,*) str_valueendifend docall tu%destroy()endblock!!----------------------!!import python module !!----------------------blocktype(module_py)::datetimetype(object)::date,today,today_strcharacter(len=:),allocatable::today_fortranierr=import_py(datetime,"datetime")ierr=datetime%getattribute(date,'date')ierr=call_py(today,date,"today")ierr=call_py(today_str,today,'isoformat')ierr=cast(today_fortran,today_str)write(*,*) 'Today is ',today_fortran call datetime%destroy()call date%destroycall today%destroycall today_str%destroyendblock!!---------------------------------!!import my python module !!---------------------------------blocktype(list)::pathstype(tuple)::argstype(dict)::kwargs,xxxtype(module_py)::mymoduletype(object)::return_valuetype(ndarray)::cellsreal,dimension(2,3)::matrixcharacter(len=:),allocatable::return_stringierr=get_sys_path(paths)ierr=paths%append(".")ierr=import_py(mymodule,'mymodule')ierr=tuple_create(args,3)ierr=args%setitem(0,12)ierr=args%setitem(1,'Hi')ierr=args%setitem(2,.true.)ierr=dict_create(kwargs)ierr=kwargs%setitem("message","Hello world")matrix=1111ierr=ndarray_create_nocopy(cells,matrix) !!ndarray(shallow copy)!ierr=kwargs%setitem("cells",cells) !!成员是ndarrayierr=dict_create(xxx)ierr=xxx%setitem("cells",cells)ierr=kwargs%setitem("xxx",xxx) !!成员是dictierr=call_py(return_value,mymodule,"print_args",args,kwargs)ierr=cast(return_string,return_value)write(*,*) return_string matrix(1,1)=99999 !!python端也被改变ierr=call_py(return_value,mymodule,"print_args",args,kwargs)ierr=cast(return_string,return_value)write(*,*) return_string call args%destroycall kwargs%destroycall mymodule%destroycall return_value%destroycall cells%destroycall xxx%destroycall paths%destroyendblock!!-----------------!!array!!-----------------!!deep-copy : creating a numpy array from a Fortran array blockreal,dimension(:,:),allocatable::matrixtype(ndarray)::arrinteger::i,jallocate(matrix(2,3))do i=1,size(matrix,1)do j=1,size(matrix,2)matrix(i,j)=real(i)*jend doenddoierr=ndarray_create(arr,matrix)ierr=print_py(arr)call arr%destroyend block!!creating a numpy wrapper for a Fortran arrayblock!Fortran: asynchronousinteger::i,jreal,asynchronous::matrix(2,3)type(ndarray)::arrforall(i=1:2,j=1:3) matrix(i,j)=real(i)*jierr=ndarray_create_nocopy(arr,matrix)ierr=print_py(arr)matrix(1,1)=1234.0ierr=print_py(arr)call arr%destroyend block!how to return a ndarray from a subroutineblockuse iso_fortran_env,only:real64type(ndarray)::arrinteger::i,jreal(kind=real64),dimension(:,:),pointer::matrixierr=ndarray_create_empty(arr,[2,3],dtype='float64')ierr=arr%get_data(matrix)forall(i=1:2,j=1:3) matrix(i,j)=real(i,kind=real64)*jierr=print_py(arr)call arr%destroyend block!!----------------------------------call forpy_finalizeend program test_forpy!!====================================
!!编译:
!!   gfortran -c forpy_mod.F90
!!   gfortran test_forpy.F90 forpy_mod.o `python3-config --ldflags --embed`
!!
!!

mymodule.py

def print_args(*args,**kwargs):print("Arguments: ",args)print("keyword arguments: ",kwargs)return "Return from mymodule.print_args"

---

资料

文档,API