<分区>
这对我来说有点奇怪，我不确定如何正确地为问题命名。我有以下 MWE，它简单地生成一个坐标点列表 (x,t) 并执行一些检查以查看它们是否位于用户规定的边界上。特别是，如果 x[i] == 1.0 和 t[i] != 0.0 那么程序应该打印一条声明。我似乎无法弄清楚为什么从未在此处输入 if 条件。我已经打印出值对 x[i] 和 t[i] 以验证确实存在满足条件的对...
#Load Modules
import numpy as np
import math, random
from pylab import meshgrid
# Create the arrays x and t on an evenly spaced cartesian grid
N = 10
xa = -1.0;
xb = 1.0;
ta = 0.0;
tb = 0.4;
xin = np.arange(xa, xb+0.00001, (xb-xa)/N).reshape((N+1,1))
tin = np.arange(ta, tb+0.00001, (tb-ta)/N).reshape((N+1,1))
X_tmp,T_tmp = meshgrid(xin,tin)
x = np.reshape(X_tmp,((N+1)**2,1))
t = np.reshape(T_tmp,((N+1)**2,1))
# create boundary flags
for i in range(0,(N+1)**2):
if (x[i] == xb and t[i] != ta):
print("We are on the right-side boundary")