我正在尝试从教程中重现一个示例,但我遇到了未定义的网格:
using DiffEqBase
using DiffEqPDEBase
f(t,x,u) = ones(size(x,1)) - .5u
u0_func(x) = zeros(size(x,1))
tspan = (0.0,1.0)
dx = 1//2^(3)
dt = 1//2^(7)
# THE FOLLOWING LINE IS FAILING
mesh = parabolic_squaremesh([0 1 0 1],dx,dt,tspan,:neumann)
u0 = u0_func(mesh.node)
prob = HeatProblem(u0,f,mesh)
sol = solve(prob,FEMDiffEqHeatImplicitEuler())
代码在我尝试使用错误消息创建网格的那一行失败:*UndefVarError: parabolic_squaremesh not defined top-level scope at test.jl:22 所有包都安装没有错误。但是,我无法安装 似乎已贬值的使用 FiniteElementDiffEq 。