fem_check33_la.c running Given uxxx(x,y,z) + 2 uyyy(x,y,z) + 3 uzzz(x,y,z) + 4 uxxy(x,y,z) + 5 uxxz(x,y,z) + 6 uyyx(x,y,z) + 7 uyyz(x,y,z) + 8 uzzx(x,y,z) + 9 uzzy(x,y,z) + 10 uxx(x,y,z) + 11 uyy(x,y,z) + 12 uzz(x,y,z) + 13 uxy(x,y,z) + 14 uxz(x,y,z) + 15 uyz(x,y,z) + 16 ux(x,y,z) + 17 uy(x,y,z) + 18 uz(x,y,z) + 19 u(x,y,z) = F(x,y,z) F(x,y,z) = 406.0 + 208.0*x*y*z*z + 240.0*x*x*y*z + 88.0*x*x*z*z + 136.0*x*x*y*z*z + 128.0*x*y*y*z*z + 217.0*x + 302.0*y + 80.0*y*y*z*z + 144.0*x*x*y*y*z + 96.0*x*x*y*y + 224.0*x*y*y*z + 76.0*x*x*y*y*z*z + 144.0*x*x*y + 80.0*y*y*z + 64.0*y*z*z + 133.0*y*z + 128.0*x*y*y + 112.0*x*x*z + 96.0*x*z*z + 114.0*x*z + 95.0*x*y + 431.0*z + 19.0*x*x*x*x + 38.0*y*y*y*y + 57.0*z*z*z*z + 64.0*x*x*x + 136.0*y*y*y + 216.0*z*z*z + 120.0*x*x + 264.0*y*y + 432.0*z*z xmin<=x<=xmax ymin<=y<=ymax zmin<=z<=zmax Boundaries Analytic solution u(x,y,z) = x^4 + 2 y^4 + 3 z^4 + 4 x^2 y^2 z^2 + 5 x y + 6 x z + 7 y z + 8 xmin=0, xmax=1, ymin=0, ymax=1, zmin=0, zmax=1 nx=5, ny=5, nz=5 time start = 0.000000 seconds x grid and analytic solution at ymin,zmin i=0, Ua( 0.000, 0.000, 0.000)= 8.0000 i=1, Ua( 0.250, 0.000, 0.000)= 8.0039 i=2, Ua( 0.500, 0.000, 0.000)= 8.0625 i=3, Ua( 0.750, 0.000, 0.000)= 8.3164 i=4, Ua( 1.000, 0.000, 0.000)= 9.0000 y grid and analytic solution at xmin, zmin ii=0, Ua( 0.000, 0.000, 0.000)= 8.0000 ii=1, Ua( 0.000, 0.250, 0.000)= 8.0078 ii=2, Ua( 0.000, 0.500, 0.000)= 8.1250 ii=3, Ua( 0.000, 0.750, 0.000)= 8.6328 ii=4, Ua( 0.000, 1.000, 0.000)=10.0000 z grid and analytic solution at xmin,ymin iii=0, Ua( 0.000, 0.000, 0.000)= 8.0000 iii=1, Ua( 0.000, 0.000, 0.250)= 8.0117 iii=2, Ua( 0.000, 0.000, 0.500)= 8.1875 iii=3, Ua( 0.000, 0.000, 0.750)= 8.9492 iii=4, Ua( 0.000, 0.000, 1.000)=11.0000 time now = 0.000000 seconds, time for previous section = 0.000000 seconds time now = 0.000000 seconds, time for previous section = 0.000000 seconds boundary i=0,x=0, ii=0,y=0, iii=0,z=0 is 8 boundary i=0,x=0, ii=1,y=0.25, iii=0,z=0 is 8.00781 boundary i=0,x=0, ii=2,y=0.5, iii=0,z=0 is 8.125 boundary i=0,x=0, ii=3,y=0.75, iii=0,z=0 is 8.63281 boundary i=0,x=0, ii=4,y=1, iii=0,z=0 is 10 boundary i=0,x=0, ii=0,y=0, iii=4,z=1 is 11 boundary i=0,x=0, ii=1,y=0.25, iii=4,z=1 is 12.7578 boundary i=0,x=0, ii=2,y=0.5, iii=4,z=1 is 14.625 boundary i=0,x=0, ii=3,y=0.75, iii=4,z=1 is 16.8828 boundary i=0,x=0, ii=4,y=1, iii=4,z=1 is 20 boundary i=0,x=0, ii=0,y=0, iii=0,z=0 is 8 boundary i=0,x=0, ii=0,y=0, iii=1,z=0.25 is 8.01172 boundary i=0,x=0, ii=0,y=0, iii=2,z=0.5 is 8.1875 boundary i=0,x=0, ii=0,y=0, iii=3,z=0.75 is 8.94922 boundary i=0,x=0, ii=0,y=0, iii=4,z=1 is 11 boundary i=0,x=0, ii=4,y=1, iii=0,z=0 is 10 boundary i=0,x=0, ii=4,y=1, iii=1,z=0.25 is 11.7617 boundary i=0,x=0, ii=4,y=1, iii=2,z=0.5 is 13.6875 boundary i=0,x=0, ii=4,y=1, iii=3,z=0.75 is 16.1992 boundary i=0,x=0, ii=4,y=1, iii=4,z=1 is 20 boundary i=1,x=0.25, ii=0,y=0, iii=0,z=0 is 8.00391 boundary i=1,x=0.25, ii=1,y=0.25, iii=0,z=0 is 8.32422 boundary i=1,x=0.25, ii=2,y=0.5, iii=0,z=0 is 8.75391 boundary i=1,x=0.25, ii=3,y=0.75, iii=0,z=0 is 9.57422 boundary i=1,x=0.25, ii=4,y=1, iii=0,z=0 is 11.2539 boundary i=1,x=0.25, ii=0,y=0, iii=4,z=1 is 12.5039 boundary i=1,x=0.25, ii=1,y=0.25, iii=4,z=1 is 14.5898 boundary i=1,x=0.25, ii=2,y=0.5, iii=4,z=1 is 16.8164 boundary i=1,x=0.25, ii=3,y=0.75, iii=4,z=1 is 19.4648 boundary i=1,x=0.25, ii=4,y=1, iii=4,z=1 is 23.0039 boundary i=1,x=0.25, ii=0,y=0, iii=0,z=0 is 8.00391 boundary i=1,x=0.25, ii=0,y=0, iii=1,z=0.25 is 8.39062 boundary i=1,x=0.25, ii=0,y=0, iii=2,z=0.5 is 8.94141 boundary i=1,x=0.25, ii=0,y=0, iii=3,z=0.75 is 10.0781 boundary i=1,x=0.25, ii=0,y=0, iii=4,z=1 is 12.5039 boundary i=1,x=0.25, ii=4,y=1, iii=0,z=0 is 11.2539 boundary i=1,x=0.25, ii=4,y=1, iii=1,z=0.25 is 13.4062 boundary i=1,x=0.25, ii=4,y=1, iii=2,z=0.5 is 15.7539 boundary i=1,x=0.25, ii=4,y=1, iii=3,z=0.75 is 18.7188 boundary i=1,x=0.25, ii=4,y=1, iii=4,z=1 is 23.0039 boundary i=2,x=0.5, ii=0,y=0, iii=0,z=0 is 8.0625 boundary i=2,x=0.5, ii=1,y=0.25, iii=0,z=0 is 8.69531 boundary i=2,x=0.5, ii=2,y=0.5, iii=0,z=0 is 9.4375 boundary i=2,x=0.5, ii=3,y=0.75, iii=0,z=0 is 10.5703 boundary i=2,x=0.5, ii=4,y=1, iii=0,z=0 is 12.5625 boundary i=2,x=0.5, ii=0,y=0, iii=4,z=1 is 14.0625 boundary i=2,x=0.5, ii=1,y=0.25, iii=4,z=1 is 16.5078 boundary i=2,x=0.5, ii=2,y=0.5, iii=4,z=1 is 19.1875 boundary i=2,x=0.5, ii=3,y=0.75, iii=4,z=1 is 22.3828 boundary i=2,x=0.5, ii=4,y=1, iii=4,z=1 is 26.5625 boundary i=2,x=0.5, ii=0,y=0, iii=0,z=0 is 8.0625 boundary i=2,x=0.5, ii=0,y=0, iii=1,z=0.25 is 8.82422 boundary i=2,x=0.5, ii=0,y=0, iii=2,z=0.5 is 9.75 boundary i=2,x=0.5, ii=0,y=0, iii=3,z=0.75 is 11.2617 boundary i=2,x=0.5, ii=0,y=0, iii=4,z=1 is 14.0625 boundary i=2,x=0.5, ii=4,y=1, iii=0,z=0 is 12.5625 boundary i=2,x=0.5, ii=4,y=1, iii=1,z=0.25 is 15.1367 boundary i=2,x=0.5, ii=4,y=1, iii=2,z=0.5 is 18 boundary i=2,x=0.5, ii=4,y=1, iii=3,z=0.75 is 21.5742 boundary i=2,x=0.5, ii=4,y=1, iii=4,z=1 is 26.5625 boundary i=3,x=0.75, ii=0,y=0, iii=0,z=0 is 8.31641 boundary i=3,x=0.75, ii=1,y=0.25, iii=0,z=0 is 9.26172 boundary i=3,x=0.75, ii=2,y=0.5, iii=0,z=0 is 10.3164 boundary i=3,x=0.75, ii=3,y=0.75, iii=0,z=0 is 11.7617 boundary i=3,x=0.75, ii=4,y=1, iii=0,z=0 is 14.0664 boundary i=3,x=0.75, ii=0,y=0, iii=4,z=1 is 15.8164 boundary i=3,x=0.75, ii=1,y=0.25, iii=4,z=1 is 18.6523 boundary i=3,x=0.75, ii=2,y=0.5, iii=4,z=1 is 21.8789 boundary i=3,x=0.75, ii=3,y=0.75, iii=4,z=1 is 25.7773 boundary i=3,x=0.75, ii=4,y=1, iii=4,z=1 is 30.8164 boundary i=3,x=0.75, ii=0,y=0, iii=0,z=0 is 8.31641 boundary i=3,x=0.75, ii=0,y=0, iii=1,z=0.25 is 9.45312 boundary i=3,x=0.75, ii=0,y=0, iii=2,z=0.5 is 10.7539 boundary i=3,x=0.75, ii=0,y=0, iii=3,z=0.75 is 12.6406 boundary i=3,x=0.75, ii=0,y=0, iii=4,z=1 is 15.8164 boundary i=3,x=0.75, ii=4,y=1, iii=0,z=0 is 14.0664 boundary i=3,x=0.75, ii=4,y=1, iii=1,z=0.25 is 17.0938 boundary i=3,x=0.75, ii=4,y=1, iii=2,z=0.5 is 20.5664 boundary i=3,x=0.75, ii=4,y=1, iii=3,z=0.75 is 24.9062 boundary i=3,x=0.75, ii=4,y=1, iii=4,z=1 is 30.8164 boundary i=4,x=1, ii=0,y=0, iii=0,z=0 is 9 boundary i=4,x=1, ii=1,y=0.25, iii=0,z=0 is 10.2578 boundary i=4,x=1, ii=2,y=0.5, iii=0,z=0 is 11.625 boundary i=4,x=1, ii=3,y=0.75, iii=0,z=0 is 13.3828 boundary i=4,x=1, ii=4,y=1, iii=0,z=0 is 16 boundary i=4,x=1, ii=0,y=0, iii=4,z=1 is 18 boundary i=4,x=1, ii=1,y=0.25, iii=4,z=1 is 21.2578 boundary i=4,x=1, ii=2,y=0.5, iii=4,z=1 is 25.125 boundary i=4,x=1, ii=3,y=0.75, iii=4,z=1 is 29.8828 boundary i=4,x=1, ii=4,y=1, iii=4,z=1 is 36 boundary i=4,x=1, ii=0,y=0, iii=0,z=0 is 9 boundary i=4,x=1, ii=0,y=0, iii=1,z=0.25 is 10.5117 boundary i=4,x=1, ii=0,y=0, iii=2,z=0.5 is 12.1875 boundary i=4,x=1, ii=0,y=0, iii=3,z=0.75 is 14.4492 boundary i=4,x=1, ii=0,y=0, iii=4,z=1 is 18 boundary i=4,x=1, ii=4,y=1, iii=0,z=0 is 16 boundary i=4,x=1, ii=4,y=1, iii=1,z=0.25 is 19.5117 boundary i=4,x=1, ii=4,y=1, iii=2,z=0.5 is 23.6875 boundary i=4,x=1, ii=4,y=1, iii=3,z=0.75 is 28.9492 boundary i=4,x=1, ii=4,y=1, iii=4,z=1 is 36 boundary i=0,x=0, ii=0,y=0, iii=0,z=0 is 8 boundary i=1,x=0.25, ii=0,y=0, iii=0,z=0 is 8.00391 boundary i=2,x=0.5, ii=0,y=0, iii=0,z=0 is 8.0625 boundary i=3,x=0.75, ii=0,y=0, iii=0,z=0 is 8.31641 boundary i=4,x=1, ii=0,y=0, iii=0,z=0 is 9 boundary i=0,x=0, ii=0,y=0, iii=4,z=1 is 11 boundary i=1,x=0.25, ii=0,y=0, iii=4,z=1 is 12.5039 boundary i=2,x=0.5, ii=0,y=0, iii=4,z=1 is 14.0625 boundary i=3,x=0.75, ii=0,y=0, iii=4,z=1 is 15.8164 boundary i=4,x=1, ii=0,y=0, iii=4,z=1 is 18 boundary i=0,x=0, ii=0,y=0, iii=0,z=0 is 8 boundary i=0,x=0, ii=0,y=0, iii=1,z=0.25 is 8.01172 boundary i=0,x=0, ii=0,y=0, iii=2,z=0.5 is 8.1875 boundary i=0,x=0, ii=0,y=0, iii=3,z=0.75 is 8.94922 boundary i=0,x=0, ii=0,y=0, iii=4,z=1 is 11 boundary i=4,x=1, ii=0,y=0, iii=0,z=0 is 9 boundary i=4,x=1, ii=0,y=0, iii=1,z=0.25 is 10.5117 boundary i=4,x=1, ii=0,y=0, iii=2,z=0.5 is 12.1875 boundary i=4,x=1, ii=0,y=0, iii=3,z=0.75 is 14.4492 boundary i=4,x=1, ii=0,y=0, iii=4,z=1 is 18 boundary i=0,x=0, ii=1,y=0.25, iii=0,z=0 is 8.00781 boundary i=1,x=0.25, ii=1,y=0.25, iii=0,z=0 is 8.32422 boundary i=2,x=0.5, ii=1,y=0.25, iii=0,z=0 is 8.69531 boundary i=3,x=0.75, ii=1,y=0.25, iii=0,z=0 is 9.26172 boundary i=4,x=1, ii=1,y=0.25, iii=0,z=0 is 10.2578 boundary i=0,x=0, ii=1,y=0.25, iii=4,z=1 is 12.7578 boundary i=1,x=0.25, ii=1,y=0.25, iii=4,z=1 is 14.5898 boundary i=2,x=0.5, ii=1,y=0.25, iii=4,z=1 is 16.5078 boundary i=3,x=0.75, ii=1,y=0.25, iii=4,z=1 is 18.6523 boundary i=4,x=1, ii=1,y=0.25, iii=4,z=1 is 21.2578 boundary i=0,x=0, ii=1,y=0.25, iii=0,z=0 is 8.00781 boundary i=0,x=0, ii=1,y=0.25, iii=1,z=0.25 is 8.45703 boundary i=0,x=0, ii=1,y=0.25, iii=2,z=0.5 is 9.07031 boundary i=0,x=0, ii=1,y=0.25, iii=3,z=0.75 is 10.2695 boundary i=0,x=0, ii=1,y=0.25, iii=4,z=1 is 12.7578 boundary i=4,x=1, ii=1,y=0.25, iii=0,z=0 is 10.2578 boundary i=4,x=1, ii=1,y=0.25, iii=1,z=0.25 is 12.2227 boundary i=4,x=1, ii=1,y=0.25, iii=2,z=0.5 is 14.3828 boundary i=4,x=1, ii=1,y=0.25, iii=3,z=0.75 is 17.1602 boundary i=4,x=1, ii=1,y=0.25, iii=4,z=1 is 21.2578 boundary i=0,x=0, ii=2,y=0.5, iii=0,z=0 is 8.125 boundary i=1,x=0.25, ii=2,y=0.5, iii=0,z=0 is 8.75391 boundary i=2,x=0.5, ii=2,y=0.5, iii=0,z=0 is 9.4375 boundary i=3,x=0.75, ii=2,y=0.5, iii=0,z=0 is 10.3164 boundary i=4,x=1, ii=2,y=0.5, iii=0,z=0 is 11.625 boundary i=0,x=0, ii=2,y=0.5, iii=4,z=1 is 14.625 boundary i=1,x=0.25, ii=2,y=0.5, iii=4,z=1 is 16.8164 boundary i=2,x=0.5, ii=2,y=0.5, iii=4,z=1 is 19.1875 boundary i=3,x=0.75, ii=2,y=0.5, iii=4,z=1 is 21.8789 boundary i=4,x=1, ii=2,y=0.5, iii=4,z=1 is 25.125 boundary i=0,x=0, ii=2,y=0.5, iii=0,z=0 is 8.125 boundary i=0,x=0, ii=2,y=0.5, iii=1,z=0.25 is 9.01172 boundary i=0,x=0, ii=2,y=0.5, iii=2,z=0.5 is 10.0625 boundary i=0,x=0, ii=2,y=0.5, iii=3,z=0.75 is 11.6992 boundary i=0,x=0, ii=2,y=0.5, iii=4,z=1 is 14.625 boundary i=4,x=1, ii=2,y=0.5, iii=0,z=0 is 11.625 boundary i=4,x=1, ii=2,y=0.5, iii=1,z=0.25 is 14.0742 boundary i=4,x=1, ii=2,y=0.5, iii=2,z=0.5 is 16.8125 boundary i=4,x=1, ii=2,y=0.5, iii=3,z=0.75 is 20.2617 boundary i=4,x=1, ii=2,y=0.5, iii=4,z=1 is 25.125 boundary i=0,x=0, ii=3,y=0.75, iii=0,z=0 is 8.63281 boundary i=1,x=0.25, ii=3,y=0.75, iii=0,z=0 is 9.57422 boundary i=2,x=0.5, ii=3,y=0.75, iii=0,z=0 is 10.5703 boundary i=3,x=0.75, ii=3,y=0.75, iii=0,z=0 is 11.7617 boundary i=4,x=1, ii=3,y=0.75, iii=0,z=0 is 13.3828 boundary i=0,x=0, ii=3,y=0.75, iii=4,z=1 is 16.8828 boundary i=1,x=0.25, ii=3,y=0.75, iii=4,z=1 is 19.4648 boundary i=2,x=0.5, ii=3,y=0.75, iii=4,z=1 is 22.3828 boundary i=3,x=0.75, ii=3,y=0.75, iii=4,z=1 is 25.7773 boundary i=4,x=1, ii=3,y=0.75, iii=4,z=1 is 29.8828 boundary i=0,x=0, ii=3,y=0.75, iii=0,z=0 is 8.63281 boundary i=0,x=0, ii=3,y=0.75, iii=1,z=0.25 is 9.95703 boundary i=0,x=0, ii=3,y=0.75, iii=2,z=0.5 is 11.4453 boundary i=0,x=0, ii=3,y=0.75, iii=3,z=0.75 is 13.5195 boundary i=0,x=0, ii=3,y=0.75, iii=4,z=1 is 16.8828 boundary i=4,x=1, ii=3,y=0.75, iii=0,z=0 is 13.3828 boundary i=4,x=1, ii=3,y=0.75, iii=1,z=0.25 is 16.3477 boundary i=4,x=1, ii=3,y=0.75, iii=2,z=0.5 is 19.7578 boundary i=4,x=1, ii=3,y=0.75, iii=3,z=0.75 is 24.0352 boundary i=4,x=1, ii=3,y=0.75, iii=4,z=1 is 29.8828 boundary i=0,x=0, ii=4,y=1, iii=0,z=0 is 10 boundary i=1,x=0.25, ii=4,y=1, iii=0,z=0 is 11.2539 boundary i=2,x=0.5, ii=4,y=1, iii=0,z=0 is 12.5625 boundary i=3,x=0.75, ii=4,y=1, iii=0,z=0 is 14.0664 boundary i=4,x=1, ii=4,y=1, iii=0,z=0 is 16 boundary i=0,x=0, ii=4,y=1, iii=4,z=1 is 20 boundary i=1,x=0.25, ii=4,y=1, iii=4,z=1 is 23.0039 boundary i=2,x=0.5, ii=4,y=1, iii=4,z=1 is 26.5625 boundary i=3,x=0.75, ii=4,y=1, iii=4,z=1 is 30.8164 boundary i=4,x=1, ii=4,y=1, iii=4,z=1 is 36 boundary i=0,x=0, ii=4,y=1, iii=0,z=0 is 10 boundary i=0,x=0, ii=4,y=1, iii=1,z=0.25 is 11.7617 boundary i=0,x=0, ii=4,y=1, iii=2,z=0.5 is 13.6875 boundary i=0,x=0, ii=4,y=1, iii=3,z=0.75 is 16.1992 boundary i=0,x=0, ii=4,y=1, iii=4,z=1 is 20 boundary i=4,x=1, ii=4,y=1, iii=0,z=0 is 16 boundary i=4,x=1, ii=4,y=1, iii=1,z=0.25 is 19.5117 boundary i=4,x=1, ii=4,y=1, iii=2,z=0.5 is 23.6875 boundary i=4,x=1, ii=4,y=1, iii=3,z=0.75 is 28.9492 boundary i=4,x=1, ii=4,y=1, iii=4,z=1 is 36 boundary i=0,x=0, ii=0,y=0, iii=0,z=0 is 8 boundary i=0,x=0, ii=1,y=0.25, iii=0,z=0 is 8.00781 boundary i=0,x=0, ii=2,y=0.5, iii=0,z=0 is 8.125 boundary i=0,x=0, ii=3,y=0.75, iii=0,z=0 is 8.63281 boundary i=0,x=0, ii=4,y=1, iii=0,z=0 is 10 boundary i=4,x=1, ii=0,y=0, iii=0,z=0 is 9 boundary i=4,x=1, ii=1,y=0.25, iii=0,z=0 is 10.2578 boundary i=4,x=1, ii=2,y=0.5, iii=0,z=0 is 11.625 boundary i=4,x=1, ii=3,y=0.75, iii=0,z=0 is 13.3828 boundary i=4,x=1, ii=4,y=1, iii=0,z=0 is 16 boundary i=0,x=0, ii=0,y=0, iii=0,z=0 is 8 boundary i=1,x=0.25, ii=0,y=0, iii=0,z=0 is 8.00391 boundary i=2,x=0.5, ii=0,y=0, iii=0,z=0 is 8.0625 boundary i=3,x=0.75, ii=0,y=0, iii=0,z=0 is 8.31641 boundary i=4,x=1, ii=0,y=0, iii=0,z=0 is 9 boundary i=0,x=0, ii=4,y=1, iii=0,z=0 is 10 boundary i=1,x=0.25, ii=4,y=1, iii=0,z=0 is 11.2539 boundary i=2,x=0.5, ii=4,y=1, iii=0,z=0 is 12.5625 boundary i=3,x=0.75, ii=4,y=1, iii=0,z=0 is 14.0664 boundary i=4,x=1, ii=4,y=1, iii=0,z=0 is 16 boundary i=0,x=0, ii=0,y=0, iii=1,z=0.25 is 8.01172 boundary i=0,x=0, ii=1,y=0.25, iii=1,z=0.25 is 8.45703 boundary i=0,x=0, ii=2,y=0.5, iii=1,z=0.25 is 9.01172 boundary i=0,x=0, ii=3,y=0.75, iii=1,z=0.25 is 9.95703 boundary i=0,x=0, ii=4,y=1, iii=1,z=0.25 is 11.7617 boundary i=4,x=1, ii=0,y=0, iii=1,z=0.25 is 10.5117 boundary i=4,x=1, ii=1,y=0.25, iii=1,z=0.25 is 12.2227 boundary i=4,x=1, ii=2,y=0.5, iii=1,z=0.25 is 14.0742 boundary i=4,x=1, ii=3,y=0.75, iii=1,z=0.25 is 16.3477 boundary i=4,x=1, ii=4,y=1, iii=1,z=0.25 is 19.5117 boundary i=0,x=0, ii=0,y=0, iii=1,z=0.25 is 8.01172 boundary i=1,x=0.25, ii=0,y=0, iii=1,z=0.25 is 8.39062 boundary i=2,x=0.5, ii=0,y=0, iii=1,z=0.25 is 8.82422 boundary i=3,x=0.75, ii=0,y=0, iii=1,z=0.25 is 9.45312 boundary i=4,x=1, ii=0,y=0, iii=1,z=0.25 is 10.5117 boundary i=0,x=0, ii=4,y=1, iii=1,z=0.25 is 11.7617 boundary i=1,x=0.25, ii=4,y=1, iii=1,z=0.25 is 13.4062 boundary i=2,x=0.5, ii=4,y=1, iii=1,z=0.25 is 15.1367 boundary i=3,x=0.75, ii=4,y=1, iii=1,z=0.25 is 17.0938 boundary i=4,x=1, ii=4,y=1, iii=1,z=0.25 is 19.5117 boundary i=0,x=0, ii=0,y=0, iii=2,z=0.5 is 8.1875 boundary i=0,x=0, ii=1,y=0.25, iii=2,z=0.5 is 9.07031 boundary i=0,x=0, ii=2,y=0.5, iii=2,z=0.5 is 10.0625 boundary i=0,x=0, ii=3,y=0.75, iii=2,z=0.5 is 11.4453 boundary i=0,x=0, ii=4,y=1, iii=2,z=0.5 is 13.6875 boundary i=4,x=1, ii=0,y=0, iii=2,z=0.5 is 12.1875 boundary i=4,x=1, ii=1,y=0.25, iii=2,z=0.5 is 14.3828 boundary i=4,x=1, ii=2,y=0.5, iii=2,z=0.5 is 16.8125 boundary i=4,x=1, ii=3,y=0.75, iii=2,z=0.5 is 19.7578 boundary i=4,x=1, ii=4,y=1, iii=2,z=0.5 is 23.6875 boundary i=0,x=0, ii=0,y=0, iii=2,z=0.5 is 8.1875 boundary i=1,x=0.25, ii=0,y=0, iii=2,z=0.5 is 8.94141 boundary i=2,x=0.5, ii=0,y=0, iii=2,z=0.5 is 9.75 boundary i=3,x=0.75, ii=0,y=0, iii=2,z=0.5 is 10.7539 boundary i=4,x=1, ii=0,y=0, iii=2,z=0.5 is 12.1875 boundary i=0,x=0, ii=4,y=1, iii=2,z=0.5 is 13.6875 boundary i=1,x=0.25, ii=4,y=1, iii=2,z=0.5 is 15.7539 boundary i=2,x=0.5, ii=4,y=1, iii=2,z=0.5 is 18 boundary i=3,x=0.75, ii=4,y=1, iii=2,z=0.5 is 20.5664 boundary i=4,x=1, ii=4,y=1, iii=2,z=0.5 is 23.6875 boundary i=0,x=0, ii=0,y=0, iii=3,z=0.75 is 8.94922 boundary i=0,x=0, ii=1,y=0.25, iii=3,z=0.75 is 10.2695 boundary i=0,x=0, ii=2,y=0.5, iii=3,z=0.75 is 11.6992 boundary i=0,x=0, ii=3,y=0.75, iii=3,z=0.75 is 13.5195 boundary i=0,x=0, ii=4,y=1, iii=3,z=0.75 is 16.1992 boundary i=4,x=1, ii=0,y=0, iii=3,z=0.75 is 14.4492 boundary i=4,x=1, ii=1,y=0.25, iii=3,z=0.75 is 17.1602 boundary i=4,x=1, ii=2,y=0.5, iii=3,z=0.75 is 20.2617 boundary i=4,x=1, ii=3,y=0.75, iii=3,z=0.75 is 24.0352 boundary i=4,x=1, ii=4,y=1, iii=3,z=0.75 is 28.9492 boundary i=0,x=0, ii=0,y=0, iii=3,z=0.75 is 8.94922 boundary i=1,x=0.25, ii=0,y=0, iii=3,z=0.75 is 10.0781 boundary i=2,x=0.5, ii=0,y=0, iii=3,z=0.75 is 11.2617 boundary i=3,x=0.75, ii=0,y=0, iii=3,z=0.75 is 12.6406 boundary i=4,x=1, ii=0,y=0, iii=3,z=0.75 is 14.4492 boundary i=0,x=0, ii=4,y=1, iii=3,z=0.75 is 16.1992 boundary i=1,x=0.25, ii=4,y=1, iii=3,z=0.75 is 18.7188 boundary i=2,x=0.5, ii=4,y=1, iii=3,z=0.75 is 21.5742 boundary i=3,x=0.75, ii=4,y=1, iii=3,z=0.75 is 24.9062 boundary i=4,x=1, ii=4,y=1, iii=3,z=0.75 is 28.9492 boundary i=0,x=0, ii=0,y=0, iii=4,z=1 is 11 boundary i=0,x=0, ii=1,y=0.25, iii=4,z=1 is 12.7578 boundary i=0,x=0, ii=2,y=0.5, iii=4,z=1 is 14.625 boundary i=0,x=0, ii=3,y=0.75, iii=4,z=1 is 16.8828 boundary i=0,x=0, ii=4,y=1, iii=4,z=1 is 20 boundary i=4,x=1, ii=0,y=0, iii=4,z=1 is 18 boundary i=4,x=1, ii=1,y=0.25, iii=4,z=1 is 21.2578 boundary i=4,x=1, ii=2,y=0.5, iii=4,z=1 is 25.125 boundary i=4,x=1, ii=3,y=0.75, iii=4,z=1 is 29.8828 boundary i=4,x=1, ii=4,y=1, iii=4,z=1 is 36 boundary i=0,x=0, ii=0,y=0, iii=4,z=1 is 11 boundary i=1,x=0.25, ii=0,y=0, iii=4,z=1 is 12.5039 boundary i=2,x=0.5, ii=0,y=0, iii=4,z=1 is 14.0625 boundary i=3,x=0.75, ii=0,y=0, iii=4,z=1 is 15.8164 boundary i=4,x=1, ii=0,y=0, iii=4,z=1 is 18 boundary i=0,x=0, ii=4,y=1, iii=4,z=1 is 20 boundary i=1,x=0.25, ii=4,y=1, iii=4,z=1 is 23.0039 boundary i=2,x=0.5, ii=4,y=1, iii=4,z=1 is 26.5625 boundary i=3,x=0.75, ii=4,y=1, iii=4,z=1 is 30.8164 boundary i=4,x=1, ii=4,y=1, iii=4,z=1 is 36 calling gauleg xmin=0, xmax=1, npx=12 calling gauleg ymin=0, ymax=1, npy=12 calling gauleg zmin=0, zmax=1, npz=12 xx[1]=0.00921968, xx[2]=0.0479414, wx[1]=0.0235877, wx[2]=0.0534697 yy[1]=0.00921968, yy[2]=0.0479414, wy[1]=0.0235877, wy[2]=0.0534697 zz[1]=0.00921968, zz[2]=0.0479414, wz[1]=0.0235877, wz[2]=0.0534697 galk(xx[2],yy[2],zz[2])=-5773.5, at i=1,j=1,ii=1,jj=1,iii=1,jjj=1 galf=107.266, at i=1,ii=1,iii=1 time now = 0.000000 seconds, time for previous section = 0.000000 seconds compute stiffness matrix Legendre=3.92786, at i=1,j=0,ii=4,jj=5,iii=4,jjj=5 Legendre=30.8343, at i=1,j=1,ii=4,jj=5,iii=4,jjj=5 Legendre=4.70871, at i=1,j=2,ii=4,jj=5,iii=4,jjj=5 Legendre=-0.402012, at i=1,j=3,ii=4,jj=5,iii=4,jjj=5 Legendre=3.05964, at i=1,j=4,ii=4,jj=5,iii=4,jjj=5 Legendre=-2.17733, at i=2,j=0,ii=4,jj=5,iii=4,jjj=5 Legendre=-9.28356, at i=2,j=1,ii=4,jj=5,iii=4,jjj=5 Legendre=29.8327, at i=2,j=2,ii=4,jj=5,iii=4,jjj=5 Legendre=5.05755, at i=2,j=3,ii=4,jj=5,iii=4,jjj=5 Legendre=-7.63118, at i=2,j=4,ii=4,jj=5,iii=4,jjj=5 Legendre=0.505727, at i=3,j=0,ii=4,jj=5,iii=4,jjj=5 Legendre=6.20171, at i=3,j=1,ii=4,jj=5,iii=4,jjj=5 Legendre=-8.93472, at i=3,j=2,ii=4,jj=5,iii=4,jjj=5 Legendre=30.2141, at i=3,j=3,ii=4,jj=5,iii=4,jjj=5 Legendre=14.1416, at i=3,j=4,ii=4,jj=5,iii=4,jjj=5 time now = 108.120000 seconds, time for previous section = 108.120000 seconds Legendre integration=91.4455, f at i=1, ii=4, iii=4 Legendre integration=36.7115, f at i=2, ii=4, iii=4 Legendre integration=125.973, f at i=3, ii=4, iii=4 k computed stiffness matrix, see above f computed forcing function, see above time now = 108.140000 seconds, time for previous section = 0.020000 seconds solve 125 equations in 125 unknowns equations solved time now = 108.150000 seconds, time for previous section = 0.010000 seconds ug computed Galerkin, Ua analytic, error ug[0,0,0]= 8.000, Ua= 8.000, err=7.63833e-14 ug[0,0,1]= 8.012, Ua= 8.012, err=7.10543e-14 ug[0,0,2]= 8.187, Ua= 8.188, err=6.57252e-14 ug[0,0,3]= 8.949, Ua= 8.949, err=5.32907e-14 ug[0,0,4]= 11.000, Ua= 11.000, err=0 ug[0,1,0]= 8.008, Ua= 8.008, err=1.59872e-14 ug[0,1,1]= 8.457, Ua= 8.457, err=3.55271e-14 ug[0,1,2]= 9.070, Ua= 9.070, err=4.61853e-14 ug[0,1,3]= 10.270, Ua= 10.270, err=7.10543e-15 ug[0,1,4]= 12.758, Ua= 12.758, err=4.44089e-14 ug[0,2,0]= 8.125, Ua= 8.125, err=9.9476e-14 ug[0,2,1]= 9.012, Ua= 9.012, err=2.13163e-14 ug[0,2,2]= 10.062, Ua= 10.062, err=2.13163e-14 ug[0,2,3]= 11.699, Ua= 11.699, err=7.99361e-14 ug[0,2,4]= 14.625, Ua= 14.625, err=1.66978e-13 ug[0,3,0]= 8.633, Ua= 8.633, err=5.15143e-14 ug[0,3,1]= 9.957, Ua= 9.957, err=5.32907e-14 ug[0,3,2]= 11.445, Ua= 11.445, err=1.45661e-13 ug[0,3,3]= 13.520, Ua= 13.520, err=1.24345e-14 ug[0,3,4]= 16.883, Ua= 16.883, err=0 ug[0,4,0]= 10.000, Ua= 10.000, err=0 ug[0,4,1]= 11.762, Ua= 11.762, err=0 ug[0,4,2]= 13.688, Ua= 13.688, err=0 ug[0,4,3]= 16.199, Ua= 16.199, err=0 ug[0,4,4]= 20.000, Ua= 20.000, err=0 ug[1,0,0]= 8.004, Ua= 8.004, err=3.19744e-14 ug[1,0,1]= 8.391, Ua= 8.391, err=4.44089e-14 ug[1,0,2]= 8.941, Ua= 8.941, err=5.86198e-14 ug[1,0,3]= 10.078, Ua= 10.078, err=1.24345e-13 ug[1,0,4]= 12.504, Ua= 12.504, err=0 ug[1,1,0]= 8.324, Ua= 8.324, err=1.59872e-13 ug[1,1,1]= 9.149, Ua= 9.149, err=9.52127e-13 ug[1,1,2]= 10.141, Ua= 10.141, err=7.10543e-13 ug[1,1,3]= 11.720, Ua= 11.720, err=4.83169e-13 ug[1,1,4]= 14.590, Ua= 14.590, err=1.42109e-14 ug[1,2,0]= 8.754, Ua= 8.754, err=6.21725e-14 ug[1,2,1]= 10.020, Ua= 10.020, err=2.41585e-13 ug[1,2,2]= 11.457, Ua= 11.457, err=4.01457e-13 ug[1,2,3]= 13.488, Ua= 13.488, err=3.8014e-13 ug[1,2,4]= 16.816, Ua= 16.816, err=0 ug[1,3,0]= 9.574, Ua= 9.574, err=1.77636e-14 ug[1,3,1]= 11.282, Ua= 11.282, err=2.89546e-13 ug[1,3,2]= 13.172, Ua= 13.172, err=7.19425e-13 ug[1,3,3]= 15.665, Ua= 15.665, err=6.66134e-13 ug[1,3,4]= 19.465, Ua= 19.465, err=6.75016e-14 ug[1,4,0]= 11.254, Ua= 11.254, err=0 ug[1,4,1]= 13.406, Ua= 13.406, err=3.55271e-14 ug[1,4,2]= 15.754, Ua= 15.754, err=2.64677e-13 ug[1,4,3]= 18.719, Ua= 18.719, err=1.77636e-14 ug[1,4,4]= 23.004, Ua= 23.004, err=0 ug[2,0,0]= 8.063, Ua= 8.062, err=7.10543e-15 ug[2,0,1]= 8.824, Ua= 8.824, err=1.77636e-14 ug[2,0,2]= 9.750, Ua= 9.750, err=3.55271e-14 ug[2,0,3]= 11.262, Ua= 11.262, err=3.73035e-14 ug[2,0,4]= 14.062, Ua= 14.062, err=4.44089e-14 ug[2,1,0]= 8.695, Ua= 8.695, err=2.30926e-14 ug[2,1,1]= 9.898, Ua= 9.898, err=1.68932e-12 ug[2,1,2]= 11.273, Ua= 11.273, err=1.71951e-12 ug[2,1,3]= 13.242, Ua= 13.242, err=1.1493e-12 ug[2,1,4]= 16.508, Ua= 16.508, err=2.4869e-14 ug[2,2,0]= 9.438, Ua= 9.438, err=2.30926e-14 ug[2,2,1]= 11.090, Ua= 11.090, err=7.28306e-13 ug[2,2,2]= 12.937, Ua= 12.938, err=6.32383e-13 ug[2,2,3]= 15.402, Ua= 15.402, err=3.28626e-13 ug[2,2,4]= 19.188, Ua= 19.188, err=3.55271e-15 ug[2,3,0]= 10.570, Ua= 10.570, err=1.06581e-13 ug[2,3,1]= 12.680, Ua= 12.680, err=4.15668e-13 ug[2,3,2]= 15.023, Ua= 15.023, err=5.15143e-14 ug[2,3,3]= 18.023, Ua= 18.023, err=4.61853e-13 ug[2,3,4]= 22.383, Ua= 22.383, err=3.55271e-15 ug[2,4,0]= 12.562, Ua= 12.562, err=0 ug[2,4,1]= 15.137, Ua= 15.137, err=5.32907e-15 ug[2,4,2]= 18.000, Ua= 18.000, err=0 ug[2,4,3]= 21.574, Ua= 21.574, err=0 ug[2,4,4]= 26.563, Ua= 26.562, err=2.4869e-14 ug[3,0,0]= 8.316, Ua= 8.316, err=0 ug[3,0,1]= 9.453, Ua= 9.453, err=1.77636e-15 ug[3,0,2]= 10.754, Ua= 10.754, err=5.32907e-15 ug[3,0,3]= 12.641, Ua= 12.641, err=0 ug[3,0,4]= 15.816, Ua= 15.816, err=0 ug[3,1,0]= 9.262, Ua= 9.262, err=0 ug[3,1,1]= 10.845, Ua= 10.845, err=1.05871e-12 ug[3,1,2]= 12.609, Ua= 12.609, err=1.57563e-12 ug[3,1,3]= 14.978, Ua= 14.978, err=1.5401e-12 ug[3,1,4]= 18.652, Ua= 18.652, err=0 ug[3,2,0]= 10.316, Ua= 10.316, err=7.10543e-15 ug[3,2,1]= 12.363, Ua= 12.363, err=8.61533e-13 ug[3,2,2]= 14.645, Ua= 14.645, err=9.50351e-13 ug[3,2,3]= 17.582, Ua= 17.582, err=9.91207e-13 ug[3,2,4]= 21.879, Ua= 21.879, err=0 ug[3,3,0]= 11.762, Ua= 11.762, err=1.77636e-15 ug[3,3,1]= 14.290, Ua= 14.290, err=5.8975e-13 ug[3,3,2]= 17.141, Ua= 17.141, err=4.54747e-13 ug[3,3,3]= 20.735, Ua= 20.735, err=2.02505e-13 ug[3,3,4]= 25.777, Ua= 25.777, err=0 ug[3,4,0]= 14.066, Ua= 14.066, err=0 ug[3,4,1]= 17.094, Ua= 17.094, err=0 ug[3,4,2]= 20.566, Ua= 20.566, err=0 ug[3,4,3]= 24.906, Ua= 24.906, err=0 ug[3,4,4]= 30.816, Ua= 30.816, err=0 ug[4,0,0]= 9.000, Ua= 9.000, err=0 ug[4,0,1]= 10.512, Ua= 10.512, err=0 ug[4,0,2]= 12.188, Ua= 12.188, err=0 ug[4,0,3]= 14.449, Ua= 14.449, err=0 ug[4,0,4]= 18.000, Ua= 18.000, err=0 ug[4,1,0]= 10.258, Ua= 10.258, err=0 ug[4,1,1]= 12.223, Ua= 12.223, err=0 ug[4,1,2]= 14.383, Ua= 14.383, err=0 ug[4,1,3]= 17.160, Ua= 17.160, err=0 ug[4,1,4]= 21.258, Ua= 21.258, err=0 ug[4,2,0]= 11.625, Ua= 11.625, err=0 ug[4,2,1]= 14.074, Ua= 14.074, err=0 ug[4,2,2]= 16.812, Ua= 16.812, err=0 ug[4,2,3]= 20.262, Ua= 20.262, err=0 ug[4,2,4]= 25.125, Ua= 25.125, err=0 ug[4,3,0]= 13.383, Ua= 13.383, err=0 ug[4,3,1]= 16.348, Ua= 16.348, err=0 ug[4,3,2]= 19.758, Ua= 19.758, err=0 ug[4,3,3]= 24.035, Ua= 24.035, err=0 ug[4,3,4]= 29.883, Ua= 29.883, err=0 ug[4,4,0]= 16.000, Ua= 16.000, err=0 ug[4,4,1]= 19.512, Ua= 19.512, err=0 ug[4,4,2]= 23.688, Ua= 23.688, err=0 ug[4,4,3]= 28.949, Ua= 28.949, err=0 ug[4,4,4]= 36.000, Ua= 36.000, err=0 total CPU time = 108.150000 seconds maxerr=1.71951e-12, avgerr=1.80677e-13 end fem_check33_la.c