当前位置: 代码迷 >> 综合 >> f2py::Fortran 90 模块中的可分配数组
  详细解决方案

f2py::Fortran 90 模块中的可分配数组

热度:82   发布时间:2024-01-04 01:58:44.0

F2PY has basic support for Fortran 90 module allocatable arrays.

请看示例Fortran 90 文件:

module modreal, allocatable, dimension(:,:) :: b 
containssubroutine foointeger kif (allocated(b)) thenprint*, "b=["do k = 1,size(b,1)print*, b(k,1:size(b,2))enddoprint*, "]"elseprint*, "b is not allocated"endifend subroutine foo
end module mod


在dos窗口编译一下:

f2py -c -m allocarr allocarr.f90


下面演示在python中如何使用:

>>> allocarr.mod.foo()b=[1.000000       2.000000       3.000000    4.000000       5.000000       6.000000    ]
>>> allocarr.mod.b                             # b is Fortran-contiguous
array([[ 1.,  2.,  3.],[ 4.,  5.,  6.]],'f')
>>> allocarr.mod.b = [[1,2,3],[4,5,6],[7,8,9]] # reallocate/initialize b
>>> allocarr.mod.foo()b=[1.000000       2.000000       3.000000    4.000000       5.000000       6.000000    7.000000       8.000000       9.000000    ]
>>> allocarr.mod.b = None                      # deallocate array
>>> allocarr.mod.foo()b is not allocated


更详细的介绍见官方网站:

http://cens.ioc.ee/projects/f2py2e/