pde_blivet.c running reading blivet_nr.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=1.100000, y=1.000000, z=1.100000 2 x=1.900000, y=1.000000, z=1.100000 3 x=1.900000, y=1.000000, z=1.900000 4 x=1.100000, y=1.000000, z=1.900000 5 x=3.100000, y=1.000000, z=1.100000 6 x=3.900000, y=1.000000, z=1.100000 7 x=3.900000, y=1.000000, z=1.900000 8 x=3.100000, y=1.000000, z=1.900000 9 x=5.100000, y=1.000000, z=1.100000 10 x=5.900000, y=1.000000, z=1.000000 11 x=5.900000, y=1.000000, z=1.900000 12 x=5.100000, y=1.000000, z=1.900000 13 x=2.000000, y=4.000000, z=1.000000 14 x=3.000000, y=4.000000, z=1.000000 15 x=3.000000, y=4.000000, z=2.000000 16 x=2.000000, y=4.000000, z=2.000000 17 x=4.000000, y=4.000000, z=1.000000 18 x=5.000000, y=4.000000, z=1.000000 19 x=5.000000, y=4.000000, z=2.000000 20 x=4.000000, y=4.000000, z=2.000000 21 x=1.000000, y=5.000000, z=1.000000 22 x=6.000000, y=5.000000, z=1.000000 23 x=6.000000, y=5.000000, z=2.000000 24 x=1.000000, y=5.000000, z=2.000000 25 x=1.000000, y=4.000000, z=1.000000 26 x=6.000000, y=4.000000, z=1.000000 27 x=6.000000, y=4.000000, z=2.000000 28 x=1.000000, y=4.000000, z=2.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.597195 2 0.717356 3 0.819192 4 0.717356 5 0.862404 6 0.932039 7 0.977865 8 0.932039 9 0.991458 10 0.999958 11 0.982154 12 0.999574 13 0.985450 14 0.999574 15 0.973848 16 0.999574 17 0.973848 18 0.909297 19 0.808496 20 0.909297 21 0.985450 22 0.675463 23 0.515501 24 0.999574 25 0.932039 26 0.808496 27 0.675463 28 0.985450 ub[0]= 0.597195 ub[1]= 0.717356 ub[2]= 0.819192 ub[3]= 0.717356 ub[4]= 0.862404 ub[5]= 0.932039 ub[6]= 0.977865 ub[7]= 0.932039 ub[8]= 0.991458 ub[9]= 0.999958 ub[10]= 0.982154 ub[11]= 0.999574 ub[12]= 0.985450 ub[13]= 0.999574 ub[14]= 0.973848 ub[15]= 0.999574 ub[16]= 0.973848 ub[17]= 0.909297 ub[18]= 0.808496 ub[19]= 0.909297 ub[20]= 0.985450 ub[21]= 0.675463 ub[22]= 0.515501 ub[23]= 0.999574 ub[24]= 0.932039 ub[25]= 0.808496 ub[26]= 0.675463 ub[27]= 0.985450 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) 8 redundant row (singular) 9 redundant row (singular) 3 C[0]=-0.306826 x^0 y^0 z^0 C[1]=0.291416 x^1 y^0 z^0 C[2]=0.517505 x^0 y^1 z^0 C[3]=0.300482 x^0 y^0 z^1 C[4]=-0.017968 x^2 y^0 z^0 C[5]=-0.039123 x^1 y^1 z^0 C[6]=-0.037248 x^1 y^0 z^1 C[7]=-0.093969 x^0 y^2 z^0 C[8]=-0.045749 x^0 y^1 z^1 C[9]=-0.019717 x^0 y^0 z^2 C[10]=-0.000008 x^3 y^0 z^0 C[11]=-0.000249 x^2 y^1 z^0 C[12]=0.000068 x^2 y^0 z^1 C[13]=0.000555 x^1 y^2 z^0 C[14]=0.000118 x^1 y^1 z^1 C[15]=0.000000 x^1 y^0 z^2 C[16]=0.007239 x^0 y^3 z^0 C[17]=0.001661 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.0038252 DOF fit check maxerr=0.109915 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=1.100000, yg=1.000000, zg=1.100000 j=1, xg=1.900000, yg=1.000000, zg=1.100000 j=2, xg=1.900000, yg=1.000000, zg=1.900000 j=3, xg=1.100000, yg=1.000000, zg=1.900000 j=4, xg=3.100000, yg=1.000000, zg=1.100000 j=5, xg=3.900000, yg=1.000000, zg=1.100000 j=6, xg=3.900000, yg=1.000000, zg=1.900000 j=7, xg=3.100000, yg=1.000000, zg=1.900000 j=8, xg=5.100000, yg=1.000000, zg=1.100000 j=9, xg=5.900000, yg=1.000000, zg=1.000000 j=10, xg=5.900000, yg=1.000000, zg=1.900000 j=11, xg=5.100000, yg=1.000000, zg=1.900000 j=12, xg=2.000000, yg=4.000000, zg=1.000000 j=13, xg=3.000000, yg=4.000000, zg=1.000000 j=14, xg=3.000000, yg=4.000000, zg=2.000000 j=15, xg=2.000000, yg=4.000000, zg=2.000000 j=16, xg=4.000000, yg=4.000000, zg=1.000000 j=17, xg=5.000000, yg=4.000000, zg=1.000000 j=18, xg=5.000000, yg=4.000000, zg=2.000000 j=19, xg=4.000000, yg=4.000000, zg=2.000000 j=20, xg=1.000000, yg=5.000000, zg=1.000000 j=21, xg=6.000000, yg=5.000000, zg=1.000000 j=22, xg=6.000000, yg=5.000000, zg=2.000000 j=23, xg=1.000000, yg=5.000000, zg=2.000000 j=24, xg=1.000000, yg=4.000000, zg=1.000000 j=25, xg=6.000000, yg=4.000000, zg=1.000000 j=26, xg=6.000000, yg=4.000000, zg=2.000000 j=27, xg=1.000000, yg=4.000000, zg=2.000000 j=28, xg=0.500000, yg=0.500000, zg=0.500000 0 x=0.500000, y=0.500000, z=0.500000 , U=0.274605, ue=0.295520, err=0.0209157 1 x=2.500000, y=0.500000, z=0.500000 , U=0.677952, ue=0.644218, err=0.0337344 2 x=4.500000, y=0.500000, z=0.500000 , U=0.962951, ue=0.891207, err=0.0717437 3 x=0.500000, y=1.500000, z=0.500000 , U=0.479777, ue=0.479426, err=0.000351108 4 x=2.500000, y=1.500000, z=0.500000 , U=0.804876, ue=0.783327, err=0.0215492 5 x=4.500000, y=1.500000, z=0.500000 , U=1.004213, ue=0.963558, err=0.0406548 6 x=0.500000, y=2.500000, z=0.500000 , U=0.651767, ue=0.644218, err=0.00754916 7 x=2.500000, y=2.500000, z=0.500000 , U=0.900110, ue=0.891207, err=0.00890255 8 x=4.500000, y=2.500000, z=0.500000 , U=1.015346, ue=0.997495, err=0.0178512 9 x=0.500000, y=3.500000, z=0.500000 , U=0.789961, ue=0.783327, err=0.00663399 10 x=2.500000, y=3.500000, z=0.500000 , U=0.963246, ue=0.963558, err=0.000312355 11 x=4.500000, y=3.500000, z=0.500000 , U=0.996168, ue=0.991665, err=0.00450343 12 x=1.500000, y=4.500000, z=0.500000 , U=0.962523, ue=0.963558, err=0.00103478 13 x=3.500000, y=4.500000, z=0.500000 , U=0.988675, ue=0.991665, err=0.00299007 computed values of DOF maxerr=0.0717437 pde_blivet.c ends