pde44e_eq.c running Given uxxxx(x,y,z,t) + uyyyy(x,y,z,t) + uzzzz(x,y,z,t) + utttt(x,y,z,t) + uxyzt(x,y,z,t) = F(x,y,z,t) F(x,y,z,t) = (1 + 7xyzt + 6xxyyzztt +x^3Y^3z^3t^3)xp(x*y*z*t) xmin<=x<=xmax ymin<=y<=ymax zmin<=z<=zmax Boundaries Analytic solution u(x,y,z,t) = exp(x*y*z*t) xmin=0, xmax=1, ymin=0, ymax=1 zmin=0, zmax=1, tmin=0, tmax=1 nx=5, ny=5, nz=5, nt=5 time start = 0.000000 seconds x grid and analytic solution at ymin,zmin,tmin i=0, Ua( 0.000, 0.000, 0.000, 0.000)= 1.0000 i=1, Ua( 0.250, 0.000, 0.000, 0.000)= 1.0000 i=2, Ua( 0.500, 0.000, 0.000, 0.000)= 1.0000 i=3, Ua( 0.750, 0.000, 0.000, 0.000)= 1.0000 i=4, Ua( 1.000, 0.000, 0.000, 0.000)= 1.0000 y grid and analytic solution at xmin,zmin,tmin ii=0, Ua( 0.000, 0.000, 0.000, 0.000)= 1.0000 ii=1, Ua( 0.000, 0.250, 0.000, 0.000)= 1.0000 ii=2, Ua( 0.000, 0.500, 0.000, 0.000)= 1.0000 ii=3, Ua( 0.000, 0.750, 0.000, 0.000)= 1.0000 ii=4, Ua( 0.000, 1.000, 0.000, 0.000)= 1.0000 z grid and analytic solution at xmin,ymin,tmin iii=0, Ua( 0.000, 0.000, 0.000, 0.000)= 1.0000 iii=1, Ua( 0.000, 0.000, 0.250, 0.000)= 1.0000 iii=2, Ua( 0.000, 0.000, 0.500, 0.000)= 1.0000 iii=3, Ua( 0.000, 0.000, 0.750, 0.000)= 1.0000 iii=4, Ua( 0.000, 0.000, 1.000, 0.000)= 1.0000 t grid and analytic solution at xmin,ymin,zmin iiii=0, Ua( 0.000, 0.000, 0.000, 0.000)= 1.0000 iiii=1, Ua( 0.000, 0.000, 0.000, 0.250)= 1.0000 iiii=2, Ua( 0.000, 0.000, 0.000, 0.500)= 1.0000 iiii=3, Ua( 0.000, 0.000, 0.000, 0.750)= 1.0000 iiii=4, Ua( 0.000, 0.000, 0.000, 1.000)= 1.0000 check rderiv(4,nx,0,hx,cxxxx) cxxxx[0]=256.000000 cxxxx[1]=-1024.000000 cxxxx[2]=1536.000000 cxxxx[3]=-1024.000000 cxxxx[4]=256.000000 check rderiv(4,nx,2,hx,cxxxx) cxxxx[0]=256.000000 cxxxx[1]=-1024.000000 cxxxx[2]=1536.000000 cxxxx[3]=-1024.000000 cxxxx[4]=256.000000 put solution in Ua vector time now = 0.000000 seconds, time for previous section = 0.000000 seconds initialize A time now = 0.000000 seconds, time for previous section = 0.000000 seconds compute A matrix A computed stiffness matrix time now = 0.010000 seconds, time for previous section = 0.010000 seconds solve 81 equations in 81 unknowns equations solved time now = 0.020000 seconds, time for previous section = 0.010000 seconds U computed, Ua analytic, error ug[1,1,1,1]= 1.0038, Ua= 1.0039, err=0.000101622 ug[1,1,1,2]= 1.0078, Ua= 1.0078, err=7.89915e-05 ug[1,1,1,3]= 1.0118, Ua= 1.0118, err=1.70143e-05 ug[1,1,2,1]= 1.0078, Ua= 1.0078, err=7.89915e-05 ug[1,1,2,2]= 1.0157, Ua= 1.0157, err=6.06493e-05 ug[1,1,2,3]= 1.0237, Ua= 1.0237, err=1.22358e-05 ug[1,1,3,1]= 1.0118, Ua= 1.0118, err=1.70143e-05 ug[1,1,3,2]= 1.0237, Ua= 1.0237, err=1.22358e-05 ug[1,1,3,3]= 1.0358, Ua= 1.0358, err=1.23721e-06 ug[1,2,1,1]= 1.0078, Ua= 1.0078, err=7.89915e-05 ug[1,2,1,2]= 1.0157, Ua= 1.0157, err=6.06493e-05 ug[1,2,1,3]= 1.0237, Ua= 1.0237, err=1.22358e-05 ug[1,2,2,1]= 1.0157, Ua= 1.0157, err=6.06493e-05 ug[1,2,2,2]= 1.0317, Ua= 1.0317, err=4.68289e-05 ug[1,2,2,3]= 1.0480, Ua= 1.0480, err=9.85253e-06 ug[1,2,3,1]= 1.0237, Ua= 1.0237, err=1.22358e-05 ug[1,2,3,2]= 1.0480, Ua= 1.0480, err=9.85253e-06 ug[1,2,3,3]= 1.0728, Ua= 1.0728, err=2.19015e-06 ug[1,3,1,1]= 1.0118, Ua= 1.0118, err=1.70143e-05 ug[1,3,1,2]= 1.0237, Ua= 1.0237, err=1.22358e-05 ug[1,3,1,3]= 1.0358, Ua= 1.0358, err=1.23721e-06 ug[1,3,2,1]= 1.0237, Ua= 1.0237, err=1.22358e-05 ug[1,3,2,2]= 1.0480, Ua= 1.0480, err=9.85253e-06 ug[1,3,2,3]= 1.0728, Ua= 1.0728, err=2.19015e-06 ug[1,3,3,1]= 1.0358, Ua= 1.0358, err=1.23721e-06 ug[1,3,3,2]= 1.0728, Ua= 1.0728, err=2.19015e-06 ug[1,3,3,3]= 1.1112, Ua= 1.1112, err=1.68402e-06 ug[2,1,1,1]= 1.0078, Ua= 1.0078, err=7.89915e-05 ug[2,1,1,2]= 1.0157, Ua= 1.0157, err=6.06493e-05 ug[2,1,1,3]= 1.0237, Ua= 1.0237, err=1.22358e-05 ug[2,1,2,1]= 1.0157, Ua= 1.0157, err=6.06493e-05 ug[2,1,2,2]= 1.0317, Ua= 1.0317, err=4.68289e-05 ug[2,1,2,3]= 1.0480, Ua= 1.0480, err=9.85253e-06 ug[2,1,3,1]= 1.0237, Ua= 1.0237, err=1.22358e-05 ug[2,1,3,2]= 1.0480, Ua= 1.0480, err=9.85253e-06 ug[2,1,3,3]= 1.0728, Ua= 1.0728, err=2.19015e-06 ug[2,2,1,1]= 1.0157, Ua= 1.0157, err=6.06493e-05 ug[2,2,1,2]= 1.0317, Ua= 1.0317, err=4.68289e-05 ug[2,2,1,3]= 1.0480, Ua= 1.0480, err=9.85253e-06 ug[2,2,2,1]= 1.0317, Ua= 1.0317, err=4.68289e-05 ug[2,2,2,2]= 1.0644, Ua= 1.0645, err=4.4654e-05 ug[2,2,2,3]= 1.0983, Ua= 1.0983, err=2.01776e-05 ug[2,2,3,1]= 1.0480, Ua= 1.0480, err=9.85253e-06 ug[2,2,3,2]= 1.0983, Ua= 1.0983, err=2.01776e-05 ug[2,2,3,3]= 1.1510, Ua= 1.1510, err=1.88294e-05 ug[2,3,1,1]= 1.0237, Ua= 1.0237, err=1.22358e-05 ug[2,3,1,2]= 1.0480, Ua= 1.0480, err=9.85253e-06 ug[2,3,1,3]= 1.0728, Ua= 1.0728, err=2.19015e-06 ug[2,3,2,1]= 1.0480, Ua= 1.0480, err=9.85253e-06 ug[2,3,2,2]= 1.0983, Ua= 1.0983, err=2.01776e-05 ug[2,3,2,3]= 1.1510, Ua= 1.1510, err=1.88294e-05 ug[2,3,3,1]= 1.0728, Ua= 1.0728, err=2.19015e-06 ug[2,3,3,2]= 1.1510, Ua= 1.1510, err=1.88294e-05 ug[2,3,3,3]= 1.2348, Ua= 1.2348, err=1.94772e-05 ug[3,1,1,1]= 1.0118, Ua= 1.0118, err=1.70143e-05 ug[3,1,1,2]= 1.0237, Ua= 1.0237, err=1.22358e-05 ug[3,1,1,3]= 1.0358, Ua= 1.0358, err=1.23721e-06 ug[3,1,2,1]= 1.0237, Ua= 1.0237, err=1.22358e-05 ug[3,1,2,2]= 1.0480, Ua= 1.0480, err=9.85253e-06 ug[3,1,2,3]= 1.0728, Ua= 1.0728, err=2.19015e-06 ug[3,1,3,1]= 1.0358, Ua= 1.0358, err=1.23721e-06 ug[3,1,3,2]= 1.0728, Ua= 1.0728, err=2.19015e-06 ug[3,1,3,3]= 1.1112, Ua= 1.1112, err=1.68402e-06 ug[3,2,1,1]= 1.0237, Ua= 1.0237, err=1.22358e-05 ug[3,2,1,2]= 1.0480, Ua= 1.0480, err=9.85253e-06 ug[3,2,1,3]= 1.0728, Ua= 1.0728, err=2.19015e-06 ug[3,2,2,1]= 1.0480, Ua= 1.0480, err=9.85253e-06 ug[3,2,2,2]= 1.0983, Ua= 1.0983, err=2.01776e-05 ug[3,2,2,3]= 1.1510, Ua= 1.1510, err=1.88294e-05 ug[3,2,3,1]= 1.0728, Ua= 1.0728, err=2.19015e-06 ug[3,2,3,2]= 1.1510, Ua= 1.1510, err=1.88294e-05 ug[3,2,3,3]= 1.2348, Ua= 1.2348, err=1.94772e-05 ug[3,3,1,1]= 1.0358, Ua= 1.0358, err=1.23721e-06 ug[3,3,1,2]= 1.0728, Ua= 1.0728, err=2.19015e-06 ug[3,3,1,3]= 1.1112, Ua= 1.1112, err=1.68402e-06 ug[3,3,2,1]= 1.0728, Ua= 1.0728, err=2.19015e-06 ug[3,3,2,2]= 1.1510, Ua= 1.1510, err=1.88294e-05 ug[3,3,2,3]= 1.2348, Ua= 1.2348, err=1.94772e-05 ug[3,3,3,1]= 1.1112, Ua= 1.1112, err=1.68402e-06 ug[3,3,3,2]= 1.2348, Ua= 1.2348, err=1.94772e-05 ug[3,3,3,3]= 1.3722, Ua= 1.3722, err=1.19436e-05 total CPU time = 0.020000 seconds maxerr=0.000101622, avgerr=2.67288e-06 end pde44e_eq.c