pde_blivet.c running reading blivet.inp pde_read_ucd returned: 1 based indexing n_vertices (boundary)= 28 n_cells (boundary faces)= 20 n_ndata (Dirichlet at vertices)= 1 boundary vertices: 1 x=0.000000, y=0.000000, z=0.000000 2 x=1.000000, y=0.000000, z=0.000000 3 x=1.000000, y=0.000000, z=1.000000 4 x=0.000000, y=0.000000, z=1.000000 5 x=2.000000, y=0.000000, z=0.000000 6 x=3.000000, y=0.000000, z=0.000000 7 x=3.000000, y=0.000000, z=1.000000 8 x=2.000000, y=0.000000, z=1.000000 9 x=4.000000, y=0.000000, z=0.000000 10 x=5.000000, y=0.000000, z=0.000000 11 x=5.000000, y=0.000000, z=1.000000 12 x=4.000000, y=0.000000, z=1.000000 13 x=1.000000, y=3.000000, z=0.000000 14 x=2.000000, y=3.000000, z=0.000000 15 x=2.000000, y=3.000000, z=1.000000 16 x=1.000000, y=3.000000, z=1.000000 17 x=3.000000, y=3.000000, z=0.000000 18 x=4.000000, y=3.000000, z=0.000000 19 x=4.000000, y=3.000000, z=1.000000 20 x=3.000000, y=3.000000, z=1.000000 21 x=0.000000, y=4.000000, z=0.000000 22 x=5.000000, y=4.000000, z=0.000000 23 x=5.000000, y=4.000000, z=1.000000 24 x=0.000000, y=4.000000, z=1.000000 25 x=0.000000, y=3.000000, z=0.000000 26 x=5.000000, y=3.000000, z=0.000000 27 x=5.000000, y=3.000000, z=1.000000 28 x=0.000000, y=3.000000, z=1.000000 boundary cells: 1 matl=0, type=4 1 2 3 4 2 matl=0, type=4 5 6 7 8 3 matl=0, type=4 9 10 11 12 4 matl=0, type=4 13 14 15 16 5 matl=0, type=4 17 18 19 20 6 matl=0, type=4 21 22 23 24 7 matl=0, type=4 1 21 24 4 8 matl=0, type=4 2 13 16 3 9 matl=0, type=4 5 14 15 8 10 matl=0, type=4 6 17 20 7 11 matl=0, type=4 9 18 19 12 12 matl=0, type=4 10 22 23 11 13 matl=0, type=4 1 2 13 25 14 matl=0, type=4 5 6 17 14 15 matl=0, type=4 9 10 26 18 16 matl=0, type=4 25 26 22 21 17 matl=0, type=4 4 3 16 28 18 matl=0, type=4 8 7 20 15 19 matl=0, type=4 12 11 27 19 20 matl=0, type=4 28 27 23 24 Dirichlet boundary: Dirichlet Boundary 1 0.000000 2 0.198669 3 0.389418 4 0.198669 5 0.389418 6 0.564642 7 0.717356 8 0.564642 9 0.717356 10 0.841471 11 0.932039 12 0.841471 13 0.717356 14 0.841471 15 0.932039 16 0.841471 17 0.932039 18 0.985449 19 0.999573 20 0.985449 21 0.717356 22 0.973847 23 0.909297 24 0.841471 25 0.564642 26 0.999573 27 0.973847 28 0.717356 check input by fitting terms used to find fit 0 a^0 * b^0 * c^0 1 a^1 * b^0 * c^0 2 a^0 * b^1 * c^0 3 a^0 * b^0 * c^1 4 a^2 * b^0 * c^0 5 a^1 * b^1 * c^0 6 a^1 * b^0 * c^1 7 a^0 * b^2 * c^0 8 a^0 * b^1 * c^1 9 a^0 * b^0 * c^2 10 a^3 * b^0 * c^0 11 a^2 * b^1 * c^0 12 a^2 * b^0 * c^1 13 a^1 * b^2 * c^0 14 a^1 * b^1 * c^1 15 a^1 * b^0 * c^2 16 a^0 * b^3 * c^0 17 a^0 * b^2 * c^1 18 a^0 * b^1 * c^2 19 a^0 * b^0 * c^3 nnn= 20 polynomial terms nnn= 20 space allocated redundant row (singular) 9 redundant row (singular) 8 redundant row (singular) 2 redundant row (singular) 6 redundant row (singular) 19 C[0]=-0.003065 x^0 y^0 z^0 C[1]=0.208332 x^1 y^0 z^0 C[2]=0.222420 x^0 y^1 z^0 C[3]=0.202860 x^0 y^0 z^1 C[4]=-0.004011 x^2 y^0 z^0 C[5]=-0.010670 x^1 y^1 z^0 C[6]=-0.010508 x^1 y^0 z^1 C[7]=-0.010678 x^0 y^2 z^0 C[8]=-0.017250 x^0 y^1 z^1 C[9]=0.000000 x^0 y^0 z^2 C[10]=-0.000775 x^3 y^0 z^0 C[11]=-0.002382 x^2 y^1 z^0 C[12]=-0.002318 x^2 y^0 z^1 C[13]=-0.001665 x^1 y^2 z^0 C[14]=-0.004225 x^1 y^1 z^1 C[15]=0.000000 x^1 y^0 z^2 C[16]=0.000000 x^0 y^3 z^0 C[17]=-0.000411 x^0 y^2 z^1 C[18]=0.000000 x^0 y^1 z^2 C[19]=0.000000 x^0 y^0 z^3 fit3 maxerr=0.00345447 DOF fit check maxerr=0.0119482 build PDE A, Y and solve nuderiv3dg.c self test completed with 0 errors input to nu3dxx for DOF 0 m=nu3dxx(order=2,npoints=29,point=28, j=0, xg=0.000000, yg=0.000000, zg=0.000000 j=1, xg=1.000000, yg=0.000000, zg=0.000000 j=2, xg=1.000000, yg=0.000000, zg=1.000000 j=3, xg=0.000000, yg=0.000000, zg=1.000000 j=4, xg=2.000000, yg=0.000000, zg=0.000000 j=5, xg=3.000000, yg=0.000000, zg=0.000000 j=6, xg=3.000000, yg=0.000000, zg=1.000000 j=7, xg=2.000000, yg=0.000000, zg=1.000000 j=8, xg=4.000000, yg=0.000000, zg=0.000000 j=9, xg=5.000000, yg=0.000000, zg=0.000000 j=10, xg=5.000000, yg=0.000000, zg=1.000000 j=11, xg=4.000000, yg=0.000000, zg=1.000000 j=12, xg=1.000000, yg=3.000000, zg=0.000000 j=13, xg=2.000000, yg=3.000000, zg=0.000000 j=14, xg=2.000000, yg=3.000000, zg=1.000000 j=15, xg=1.000000, yg=3.000000, zg=1.000000 j=16, xg=3.000000, yg=3.000000, zg=0.000000 j=17, xg=4.000000, yg=3.000000, zg=0.000000 j=18, xg=4.000000, yg=3.000000, zg=1.000000 j=19, xg=3.000000, yg=3.000000, zg=1.000000 j=20, xg=0.000000, yg=4.000000, zg=0.000000 j=21, xg=5.000000, yg=4.000000, zg=0.000000 j=22, xg=5.000000, yg=4.000000, zg=1.000000 j=23, xg=0.000000, yg=4.000000, zg=1.000000 j=24, xg=0.000000, yg=3.000000, zg=0.000000 j=25, xg=5.000000, yg=3.000000, zg=0.000000 j=26, xg=5.000000, yg=3.000000, zg=1.000000 j=27, xg=0.000000, yg=3.000000, zg=1.000000 j=28, xg=0.500000, yg=0.500000, zg=0.500000 0 x=0.500000, y=0.500000, z=0.500000 , U=0.287740, ue=0.295520, err=0.00778031 1 x=2.500000, y=0.500000, z=0.500000 , U=0.595087, ue=0.644218, err=0.0491305 2 x=4.500000, y=0.500000, z=0.500000 , U=0.790418, ue=0.891207, err=0.100789 3 x=0.500000, y=1.500000, z=0.500000 , U=0.474503, ue=0.479426, err=0.0049221 4 x=2.500000, y=1.500000, z=0.500000 , U=0.750697, ue=0.783327, err=0.0326302 5 x=4.500000, y=1.500000, z=0.500000 , U=0.891403, ue=0.963558, err=0.0721551 6 x=0.500000, y=2.500000, z=0.500000 , U=0.640728, ue=0.644218, err=0.00348973 7 x=2.500000, y=2.500000, z=0.500000 , U=0.877953, ue=0.891207, err=0.0132544 8 x=4.500000, y=2.500000, z=0.500000 , U=0.956221, ue=0.997495, err=0.0412744 9 x=0.500000, y=3.500000, z=0.500000 , U=0.786413, ue=0.783327, err=0.00308654 10 x=2.500000, y=3.500000, z=0.500000 , U=0.976856, ue=0.963558, err=0.0132978 11 x=4.500000, y=3.500000, z=0.500000 , U=0.984871, ue=0.991665, err=0.00679428 12 x=1.500000, y=4.500000, z=0.500000 , U=1.003055, ue=0.963558, err=0.0394971 13 x=3.500000, y=4.500000, z=0.500000 , U=1.040282, ue=0.991665, err=0.0486168 computed values of DOF maxerr=0.100789 pde_blivet.c ends