title '1D Accuracy Test - Gaussian'
select
ngrid = 5
stages = 4
errlim = staged(1e-2, 1e-3, 1e-4, 1e-5)
coordinates
cartesian1
Variables
u
definitions
k = 1
w=0.1
u0 = exp(-x^2/w^2)
s = -dxx(u0)
left=point(-1)
right=point(1)
equations
div(K*grad(u)) +s = 0
boundaries
Region 1
start left point value(u)=u0
line to right point value(u)=u0
monitors
elevation(u) from left to right
plots
elevation(u,u0) from left to right
elevation(u-u0) from left to right
elevation(-div(grad(u)),s) from left to right
histories
history(sqrt(integral((u-u0)^2)/integral(1)), errlim)
end