error3d.fee

READ_MESH $1.msh DIMENSIONS 3 READ_FIELD Temperature as T

INTEGRATE 1 RESULT volume
h = (volume/cells)^(1/3)

T_mms(x,y,z) := $2

INTEGRATE (T(x,y,z)-T_mms(x,y,z))^2 RESULT e_2
error_2 = sqrt(e_2)

FIND_EXTREMA abs(T(x,y,z)-T_mms(x,y,z)) MAX error_inf

PRINT %e h error_inf error_2 log(h) log(error_inf) log(error_2) %g $3 cells nodes volume

WRITE_MESH $1.vtk T T_mms T(x,y,z)-T_mms(x,y,z)