Hermite (Délèze) tetrahedron interpolation with 16 parameters Integration of cubic terms on a random tetrahedron. Check integration of the derivatives. x1,y1,z1 = 6.769952D+00 3.599294D+00 3.299880D+00 x2,y2,z2 = 8.474969D-01 8.960828D+00 6.568827D+00 x3,y3,z3 = 9.324702D+00 4.006439D+00 4.156038D+00 x4,y4,z4 = 6.968066D+00 6.569350D+00 8.480313D+00 Nb. of integration points = 5 Int (x) = 4.27797480D+01 Int (0.5d(x**2)/dx) = 4.27797480D+01 Int (x**2)) = 2.69720366D+02 Int (0.33*d(x**3)/dx) = 2.59189956D+02 Int (y**2)) = 2.46094820D+02 Int (0.33*d(y**3)/dy) = 2.27831045D+02 Int (z**2) = 2.32488097D+02 Int (0.33*d(z**3)/dz) = 2.26737824D+02 Nb. of integration points = 11 Int (x) = 4.27797492D+01 Int (0.5d(x**2)/dx) = 4.27797492D+01 Int (x**2)) = 2.69720375D+02 Int (0.33*d(x**3)/dx) = 2.66862001D+02 Int (y**2)) = 2.46094817D+02 Int (0.33*d(y**3)/dy) = 2.38356230D+02 Int (z**2) = 2.32488094D+02 Int (0.33*d(z**3)/dz) = 2.32290620D+02 Nb. of integration points = 14 Int (x) = 4.27797500D+01 Int (0.5d(x**2)/dx) = 4.27797500D+01 Int (x**2)) = 2.69720380D+02 Int (0.33*d(x**3)/dx) = 2.65420527D+02 Int (y**2)) = 2.46094822D+02 Int (0.33*d(y**3)/dy) = 2.44432569D+02 Int (z**2) = 2.32488099D+02 Int (0.33*d(z**3)/dz) = 2.32310094D+02 Nb. of integration points = 15 Int (x) = 4.27797495D+01 Int (0.5d(x**2)/dx) = 4.27797495D+01 Int (x**2)) = 2.69720376D+02 Int (0.33*d(x**3)/dx) = 2.65051171D+02 Int (y**2)) = 2.46094822D+02 Int (0.33*d(y**3)/dy) = 2.43721216D+02 Int (z**2) = 2.32488099D+02 Int (0.33*d(z**3)/dz) = 2.32139806D+02 Nb. of integration points = 24 Int (x) = 4.27797512D+01 Int (0.5d(x**2)/dx) = 4.27797512D+01 Int (x**2)) = 2.69720387D+02 Int (0.33*d(x**3)/dx) = 2.66908222D+02 Int (y**2)) = 2.46094832D+02 Int (0.33*d(y**3)/dy) = 2.44125968D+02 Int (z**2) = 2.32488108D+02 Int (0.33*d(z**3)/dz) = 2.32844251D+02 Nb. of integration points = 31 Int (x) = 4.27797451D+01 Int (0.5d(x**2)/dx) = 4.27797451D+01 Int (x**2)) = 2.69720349D+02 Int (0.33*d(x**3)/dx) = 2.71115397D+02 Int (y**2)) = 2.46094811D+02 Int (0.33*d(y**3)/dy) = 2.45626337D+02 Int (z**2) = 2.32488089D+02 Int (0.33*d(z**3)/dz) = 2.34695612D+02