Убедитесь, что решение удовлетворяет PDE: FiPy

Скажем, у меня есть решение следующей реализации уравнения теплопроводности с использованием FiPy:

from fipy import CellVariable, FaceVariable, TransientTerm, DiffusionTerm, Grid1D
from fipy.tools import numerix

# Set Up Grid and Time Step
nx, dx = 100, 0.01
mesh = Grid1D(dx, nx)

deltaT = 0.9 * (dx**2) / 2
steps = 1000

# Construct Solution Variable
u = CellVariable('u', mesh=mesh, rank=0)

# # Set Dirchelet Boundary Conditions
u.constrain(value=0, where=mesh.facesLeft)
u.constrain(value=0, where=mesh.facesLeft)

# # Set Initial Condition
X = mesh.cellCenters[0]
gaussInit = numerix.exp(-(X-0.5)**2 / (2*0.125**2))
u.setValue(gaussInit)

# Construct PDE
eq = TransientTerm() == DiffusionTerm(coeff=1.0)

# Solve the PDE
T = np.array([0])
Solution = u.value

time = 0
time_iter = 0
while time_iter < steps:
    time += deltaT
    eqn.solve(dt=dt, var=u)

    T = np.append(T, time)
    Solution = np.vstack((Solution, u.value))

Теперь, когда у меня есть решение УЧП Solution на каждом временном шаге T, я хочу посмотреть, удовлетворяют ли эти данные УЧП. Как мне это сделать?

** Редактировать **

Рассмотрим следующий код:

eqn = (TransientTerm(coeff=rho)
       + ConvectionTerm(coeff=v)
       == DiffusionTerm(coeff=D)
       + ImplicitSourceTerm(coeff=S1)
       + S0)

data = np.load('EquationSolution.npz')
t = data['t']
x_cell, x_face = data['x_cell'], data['x_face']
phi_cell, phi_face = data['phi_cell'], data['phi_face']
# _cell means cell-centered
# _face means face-centered
# phi_cell: solution at all times t and all locations x_cell

Я хочу убедиться, что мои загруженные данные решают приведенное выше уравнение eqn. Я не хочу проверять это, решая УЧП, поскольку это будет слишком дорого. Скорее, я хочу как-то использовать eqn.residualVectorAndNorm, чтобы посмотреть, решит ли phi_celleqn. Как мне лучше всего это сделать?