test_nuderiv2dgc.c on corner solve Uxx(x,y) + 2*Uyy(x,y) = 14 boundary and solution U(x,y) = x^2+2xy+3y^2+4x+5y+6 B 2 B = boundary P 5 P = DOF B P 1 4 P P P 3 7 8 B B B 0 6 9 i= 0, b= 0, xp=0.000000, yp=0.000000, U=6.000000 i= 1, b= 0, xp=0.000000, yp=1.000000, U=14.000000 i= 2, b= 0, xp=0.000000, yp=2.000000, U=28.000000 i= 3, b= 1, xp=0.500000, yp=0.500000, U=12.000000 i= 4, b= 1, xp=0.500000, yp=1.000000, U=17.250000 i= 5, b= 1, xp=0.500000, yp=1.500000, U=24.000000 i= 6, b= 0, xp=1.000000, yp=0.000000, U=11.000000 i= 7, b= 1, xp=1.000000, yp=0.500000, U=15.250000 i= 8, b= 1, xp=1.500000, yp=0.500000, U=19.000000 i= 9, b= 0, xp=2.000000, yp=0.000000, U=18.000000 build derivative for Point=3 nu2dxx at 3, nupx=7 j= 0, ipx[j]= 3, cxx[j]=-2.000000 j= 1, ipx[j]= 1, cxx[j]=3.000000 j= 2, ipx[j]= 2, cxx[j]=-0.500000 j= 3, ipx[j]= 6, cxx[j]=3.000000 j= 4, ipx[j]= 7, cxx[j]=-8.000000 j= 5, ipx[j]= 8, cxx[j]=6.000000 j= 6, ipx[j]= 9, cxx[j]=-1.500000 nu2dyy at 3, nupy=7 j= 0, ipy[j]= 3, cyy[j]=4.000000 j= 1, ipy[j]= 1, cyy[j]=-2.000000 j= 2, ipy[j]= 2, cyy[j]=1.000000 j= 3, ipy[j]= 6, cyy[j]=2.000000 j= 4, ipy[j]= 7, cyy[j]=-8.000000 j= 5, ipy[j]= 8, cyy[j]=4.000000 j= 6, ipy[j]= 9, cyy[j]=-1.000000 at 3 sum=14.000000, f[i]=14.000000, err=3.552714e-14 build derivative for Point=4 nu2dxx at 4, nupx=7 j= 0, ipx[j]= 4, cxx[j]=-2.666667 j= 1, ipx[j]= 1, cxx[j]=2.666667 j= 2, ipx[j]= 2, cxx[j]=0.000000 j= 3, ipx[j]= 6, cxx[j]=1.333333 j= 4, ipx[j]= 7, cxx[j]=-5.333333 j= 5, ipx[j]= 8, cxx[j]=5.333333 j= 6, ipx[j]= 9, cxx[j]=-1.333333 nu2dyy at 4, nupy=7 j= 0, ipy[j]= 4, cyy[j]=5.333333 j= 1, ipy[j]= 1, cyy[j]=-1.333333 j= 2, ipy[j]= 2, cyy[j]=0.000000 j= 3, ipy[j]= 6, cyy[j]=5.333333 j= 4, ipy[j]= 7, cyy[j]=-13.333333 j= 5, ipy[j]= 8, cyy[j]=5.333333 j= 6, ipy[j]= 9, cyy[j]=-1.333333 at 4 sum=14.000000, f[i]=14.000000, err=6.572520e-14 build derivative for Point=5 nu2dxx at 5, nupx=7 j= 0, ipx[j]= 5, cxx[j]=-1.333333 j= 1, ipx[j]= 1, cxx[j]=2.000000 j= 2, ipx[j]= 2, cxx[j]=0.333333 j= 3, ipx[j]= 6, cxx[j]=2.000000 j= 4, ipx[j]= 7, cxx[j]=-8.000000 j= 5, ipx[j]= 8, cxx[j]=6.666667 j= 6, ipx[j]= 9, cxx[j]=-1.666667 nu2dyy at 5, nupy=7 j= 0, ipy[j]= 5, cyy[j]=2.666667 j= 1, ipy[j]= 1, cyy[j]=-0.000000 j= 2, ipy[j]= 2, cyy[j]=-0.666667 j= 3, ipy[j]= 6, cyy[j]=4.000000 j= 4, ipy[j]= 7, cyy[j]=-8.000000 j= 5, ipy[j]= 8, cyy[j]=2.666667 j= 6, ipy[j]= 9, cyy[j]=-0.666667 at 5 sum=14.000000, f[i]=14.000000, err=6.039613e-14 build derivative for Point=7 nu2dxx at 7, nupx=8 j= 0, ipx[j]= 7, cxx[j]=-8.000000 j= 1, ipx[j]= 1, cxx[j]=0.000000 j= 2, ipx[j]= 2, cxx[j]=-0.000000 j= 3, ipx[j]= 3, cxx[j]=4.000000 j= 4, ipx[j]= 4, cxx[j]=-0.000000 j= 5, ipx[j]= 5, cxx[j]=0.000000 j= 6, ipx[j]= 8, cxx[j]=4.000000 j= 7, ipx[j]= 9, cxx[j]=0.000000 nu2dyy at 7, nupy=8 j= 0, ipy[j]= 7, cyy[j]=0.000000 j= 1, ipy[j]= 1, cyy[j]=0.000000 j= 2, ipy[j]= 2, cyy[j]=-0.000000 j= 3, ipy[j]= 3, cyy[j]=4.000000 j= 4, ipy[j]= 4, cyy[j]=-8.000000 j= 5, ipy[j]= 5, cyy[j]=4.000000 j= 6, ipy[j]= 8, cyy[j]=-0.000000 j= 7, ipy[j]= 9, cyy[j]=0.000000 at 7 sum=14.000000, f[i]=14.000000, err=1.243450e-14 build derivative for Point=8 nu2dxx at 8, nupx=10 j= 0, ipx[j]= 8, cxx[j]=9.333333 j= 1, ipx[j]= 1, cxx[j]=-4.000000 j= 2, ipx[j]= 2, cxx[j]=0.666667 j= 3, ipx[j]= 3, cxx[j]=12.000000 j= 4, ipx[j]= 4, cxx[j]=8.000000 j= 5, ipx[j]= 5, cxx[j]=-2.666667 j= 6, ipx[j]= 6, cxx[j]=2.000000 j= 7, ipx[j]= 7, cxx[j]=-24.000000 j= 8, ipx[j]= 0, cxx[j]=-1.000000 j= 9, ipx[j]= 9, cxx[j]=-0.333333 nu2dyy at 8, nupy=10 j= 0, ipy[j]= 8, cyy[j]=2.666667 j= 1, ipy[j]= 1, cyy[j]=2.000000 j= 2, ipy[j]= 2, cyy[j]=-1.666667 j= 3, ipy[j]= 3, cyy[j]=4.000000 j= 4, ipy[j]= 4, cyy[j]=-8.000000 j= 5, ipy[j]= 5, cyy[j]=6.666667 j= 6, ipy[j]= 6, cyy[j]=4.000000 j= 7, ipy[j]= 7, cyy[j]=-8.000000 j= 8, ipy[j]= 0, cyy[j]=-1.000000 j= 9, ipy[j]= 9, cyy[j]=-0.666667 at 8 sum=14.000000, f[i]=14.000000, err=-1.030287e-13 end test_nuderiv2dgc.c ndof=5