Fortran 90 到 Python 数组答案

【问题标题】：Fortran 90 to Python arrayFortran 90 到 Python 数组
【发布时间】：2016-05-05 18:45:41
【问题描述】：

如何将以下 Fortran 声明转换为 Python（numpy 数组）？

INTEGER, ALLOCATABLE::          I(:)
DOUBLE PRECISION, ALLOCATABLE:: A(:)
DOUBLE PRECISION, ALLOCATABLE:: B(:,:)
DIMENSION Z(*)
DIMENSION X(N)

【问题讨论】：

【解决方案1】：

要声明一个整数，只需像这样赋值：

i=1

要声明一个浮点数，像这样赋值：

a=1.0

对于numpy 数组最好提前知道大小：

b=np.zeros((100,100)) # a 100x100 array, initialized with zeros
x=np.ones(1000)   # a 100 array , initialized with ones

如果您不提前知道数组的大小，请分配一个空列表，如下所示：

xlist=[]

然后像这样填充列表（比如在循环内）：

xlist.append(5.34)

为什么这是首选方法，请参阅this SO question。

您可以使用第二个循环制作列表列表，如以下乘法表示例：

aa=[]
for i in range(10):
     a=[]
     for j in range(10):
         a.append(i*j)
     aa.append(a)

当您完成创建列表（任意维度）并想要快速时，将它们转换为像这样的 numpy 数组

my2darray=np.array(aa)

【讨论】：