Question

我创建了一个2D numpy数组，并在迭代过程中的随机位置引入了“ 1”和“ -1”。现在，我想要包含那些索引和减去索引的索引列表。这是我的代码：

import math
import numpy as np
import itertools

def function(
        250,
        10,
        50):


    crystal = np.zeros((250,250))

    defect_concentration = np.zeros(50)
    for k in range(50):


        i= np.random.randint(0,250-1,size =2)
        while crystal[i[0]][i[1]] != 0:
            i = np.random.randint(0,250-1,size =2)
        crystal[i[0]][i[1]] = -1
        j = np.random.randint(0,250-1,size =2)
        while crystal[j[0]][j[1]] !=0:
            j = np.random.randint(0,250-1,size =2)
        crystal[j[0]][j[1]] = 1 

        b = []
        c = []
        dist = math.sqrt(((i[0]-j[0])**2)+((i[1]-j[1])**2))

        if dist <= 10:
            crystal[i[0]][i[1]] = 0
            crystal[j[0]][j[1]] = 0

        else:    

            for x in range(len(crystal)):
                for y in range(len(crystal)):

                    if crystal[x][y] < 0:
                        b.append((x,y))
                     elif crystal[x][y] != 0:
                        c.append((x,y))

但是，我没有得到预期的输出。

b（1）

-------------------------------------------------- ---------------------------- NameError追溯（最近的呼叫   最后）在（）   ----> 1 b（1）

NameError：名称'b'未定义

Answer 1

这不是您定义函数的方式：

def function(
        250,
        10,
        50):

您可以清理crystal索引：

In [512]: crystal[i[0]][i[1]]=10
In [513]: crystal[i[0]][i[1]]
Out[513]: 10.0

更多惯用的numpy：

In [514]: crystal[i[0], i[1]]
Out[514]: 10.0

甚至

In [515]: crystal[tuple(i)]
Out[515]: 10.0

由于i只有2个数字，因此它们是等效的，但是不久以后[...][...]索引风格就会咬住你。

通过2D numpy数组创建具有'1'和'-1'元素的索引列表

1 个答案: