Solution u is known but still I am getting an error for H1 norm due to dx(u) term. Why?

exampleofdisk.edp (1.4 KB)

I didn’t bother to diff it but the derivatives apparently are defined by the element
and mesh. The integral just specifies the mesh. I only replaced the derivatives but you can
replace the value term too. Also just a coding note, you use “pow(float,int)” with small int
and only want
an int. I’m not use of the guarantee about getting an int back but floating point.
AFAICT FF does not support “<<” for shift but repeated int multiplies may be worth considering.

exampleofdisk.edp (1.5 KB)

1 Like