so assume you want to solve:
int_0^1 grad(u)grad(v) dx = int_0^1 v dx
on [0,1] mesh with let’s say N=5 cells (i.e. 6 nodes). For the basis, we want to use tent functions. I.e. something like this Tent functions
The cells are of width h=1/N, so they are equidistant.
I somehow fail to find out how I can use such a simple 1D mesh.
Please have a look at LaplacianCurve.edp.
Very little adjustements are needed, see below.
func f = 1 ;
macro Grad3(uh) dx(uh) // EOM
// with problem
problem Lap3dL(uLPb,vLPb) = int1d(ThL3)(Grad3(uLPb)'*Grad3(vLPb))
- int1d ( ThL3 ) ( f * vLPb )