sparse_abc.c running with user files abc.h abc.c SparseWrite_All nrow=25 xmin=-1, xmax=1, nx=7, hx=0.333333 ymin=-1, ymax=1, ny=7, hy=0.333333 internal cells zeroed build_coef for point ix=2, iy=2 deriv pdxx returns bb=12, aa=-1 16 -30 16 -1 deriv pdyy returns bb=12, aa=-1 16 -30 16 -1 deriv pdx returns bb=12, aa=1 -8 0 8 -1 deriv pdy returns bb=12, aa=1 -8 0 8 -1 pdxx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 -0.75 0.00 0.00 i=1, 0.00 0.00 12.00 0.00 0.00 i=2, 0.00 0.00 -22.50 0.00 0.00 i=3, 0.00 0.00 12.00 0.00 0.00 i=4, 0.00 0.00 -0.75 0.00 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, 0.06 -0.50 0.00 0.50 -0.06 i=1, -0.50 4.00 -0.00 -4.00 0.50 i=2, 0.00 -0.00 0.00 0.00 -0.00 i=3, 0.50 -4.00 0.00 4.00 -0.50 i=4, -0.06 0.50 -0.00 -0.50 0.06 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, -0.75 12.00 -22.50 12.00 -0.75 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.25 0.00 0.00 i=1, 0.00 0.00 -2.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 2.00 0.00 0.00 i=4, 0.00 0.00 -0.25 0.00 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 0.25 -2.00 0.00 2.00 -0.25 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 working row 6, i=2, j=2, bi=0, bj=0 a1=0.303265, b1=1.36627, c1=4.87397, d1=-0.222222, e1=-0.037037 f1=-1.66667, c=-41.222 ij=0, u[ii][ij]=5.46506 ij=5, u[ii][ij]=4.08363 ij=10, u[ii][ij]=-5.46506 ij=15, u[ii][ij]=0.683133 ij=1, u[ii][ij]=58.5618 ij=6, u[ii][ij]=-118.155 ij=11, u[ii][ij]=58.4136 ij=16, u[ii][ij]=-3.64622 ij=2, u[ii][ij]=-5.46506 ij=7, u[ii][ij]=3.19474 ij=12, u[ii][ij]=5.46506 ij=17, u[ii][ij]=-0.683133 ij=3, u[ii][ij]=0.683133 ij=8, u[ii][ij]=-0.171893 ij=13, u[ii][ij]=-0.683133 ij=18, u[ii][ij]=0.0853916 ii=6, ut[ii][ii]=-118.155, ut[ii][cs]=-58.8925 working row 11, i=2, j=3, bi=0, bj=0 a1=0.423241, b1=1.4, c1=4.56694, d1=0.111111, e1=-1.02716e-33 f1=-1, c=-14.0842 ij=0, u[ii][ij]=-0.7 ij=5, u[ii][ij]=5.6 ij=10, u[ii][ij]=4.85667 ij=15, u[ii][ij]=-5.6 ij=20, u[ii][ij]=0.7 ij=1, u[ii][ij]=-3.4252 ij=6, u[ii][ij]=54.8032 ij=11, u[ii][ij]=-113.279 ij=16, u[ii][ij]=54.8032 ij=21, u[ii][ij]=-3.4252 ij=2, u[ii][ij]=0.7 ij=7, u[ii][ij]=-5.6 ij=12, u[ii][ij]=5.30111 ij=17, u[ii][ij]=5.6 ij=22, u[ii][ij]=-0.7 ij=3, u[ii][ij]=-0.0875 ij=8, u[ii][ij]=0.7 ij=13, u[ii][ij]=-0.345208 ij=18, u[ii][ij]=-0.7 ij=23, u[ii][ij]=0.0875 ii=11, ut[ii][ii]=-113.279, ut[ii][cs]=-15.5446 working row 16, i=2, j=4, bi=0, bj=0 a1=0.59068, b1=1.36627, c1=4.20422, d1=0.444444, e1=-0.037037 f1=-0.333333, c=14.0103 ii=16, ut[ii][ii]=-108.219, ut[ii][cs]=56.6188 working row 7, i=3, j=2, bi=0, bj=0 a1=0.358266, b1=1.4, c1=4.30704, d1=-0.333333, e1=-6.16791e-18 f1=-0.666667, c=-19.999 ij=0, u[ii][ij]=-0.7 ij=5, u[ii][ij]=-0.352033 ij=10, u[ii][ij]=0.7 ij=15, u[ii][ij]=-0.0875 ij=1, u[ii][ij]=5.6 ij=6, u[ii][ij]=4.96585 ij=11, u[ii][ij]=-5.6 ij=16, u[ii][ij]=0.7 ij=2, u[ii][ij]=51.6844 ij=7, u[ii][ij]=-105.636 ij=12, u[ii][ij]=51.6844 ij=17, u[ii][ij]=-3.23028 ij=3, u[ii][ij]=-5.6 ij=8, u[ii][ij]=3.63252 ij=13, u[ii][ij]=5.6 ij=18, u[ii][ij]=-0.7 ij=4, u[ii][ij]=0.7 ij=9, u[ii][ij]=-0.185366 ij=14, u[ii][ij]=-0.7 ij=19, u[ii][ij]=0.0875 ii=7, ut[ii][ii]=-105.636, ut[ii][cs]=-24.9793 working row 12, i=3, j=3, bi=0, bj=0 a1=0.5, b1=1.4, c1=4, d1=-5.55112e-17, e1=-1.71057e-49 f1=-2.77556e-16, c=7 ij=0, u[ii][ij]=0.0875 ij=5, u[ii][ij]=-0.7 ij=10, u[ii][ij]=-0.375 ij=15, u[ii][ij]=0.7 ij=20, u[ii][ij]=-0.0875 ij=1, u[ii][ij]=-0.7 ij=6, u[ii][ij]=5.6 ij=11, u[ii][ij]=6 ij=16, u[ii][ij]=-5.6 ij=21, u[ii][ij]=0.7 ij=2, u[ii][ij]=-3 ij=7, u[ii][ij]=48 ij=12, u[ii][ij]=-101.25 ij=17, u[ii][ij]=48 ij=22, u[ii][ij]=-3 ij=3, u[ii][ij]=0.7 ij=8, u[ii][ij]=-5.6 ij=13, u[ii][ij]=6 ij=18, u[ii][ij]=5.6 ij=23, u[ii][ij]=-0.7 ij=4, u[ii][ij]=-0.0875 ij=9, u[ii][ij]=0.7 ij=14, u[ii][ij]=-0.375 ij=19, u[ii][ij]=-0.7 ij=24, u[ii][ij]=0.0875 ii=12, ut[ii][ii]=-101.25, ut[ii][cs]=7 working row 17, i=3, j=4, bi=0, bj=0 a1=0.697806, b1=1.4, c1=3.63728, d1=0.333333, e1=-6.16791e-18 f1=0.666667, c=36.32 ii=17, ut[ii][ii]=-96.8728, ut[ii][cs]=87.9453 working row 8, i=4, j=2, bi=0, bj=0 a1=0.423241, b1=1.36627, c1=3.51581, d1=-0.222222, e1=0.037037 f1=0.333333, c=2.94055 ii=8, ut[ii][ii]=-88.2953, ut[ii][cs]=9.7605 working row 13, i=4, j=3, bi=0, bj=0 a1=0.59068, b1=1.4, c1=3.20878, d1=0.111111, e1=1.02716e-33 f1=1, c=30.2788 ii=13, ut[ii][ii]=-84.4877, ut[ii][cs]=42.5907 working row 18, i=4, j=4, bi=0, bj=0 a1=0.824361, b1=1.36627, c1=2.84605, d1=0.444444, e1=0.037037 f1=1.66667, c=62.0175 ij=6, u[ii][ij]=0.0853916 ij=11, u[ii][ij]=-0.683133 ij=16, u[ii][ij]=-0.507159 ij=21, u[ii][ij]=0.683133 ij=7, u[ii][ij]=-0.683133 ij=12, u[ii][ij]=5.46506 ij=17, u[ii][ij]=9.00344 ij=22, u[ii][ij]=-5.46506 ij=8, u[ii][ij]=-2.12528 ij=13, u[ii][ij]=34.0786 ij=18, u[ii][ij]=-80.9177 ij=23, u[ii][ij]=34.2267 ij=9, u[ii][ij]=0.683133 ij=14, u[ii][ij]=-5.46506 ij=19, u[ii][ij]=10.7812 ij=24, u[ii][ij]=5.46506 ii=18, ut[ii][ii]=-80.9177, ut[ii][cs]=140.198 standard central cells generated build_coef for point ix=1, iy=2 deriv pdxx returns bb=12, aa=11 -20 6 4 -1 deriv pdyy returns bb=12, aa=-1 16 -30 16 -1 deriv pdx returns bb=12, aa=-3 -10 18 -6 1 deriv pdy returns bb=12, aa=1 -8 0 8 -1 pdxx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 8.25 0.00 0.00 i=1, 0.00 0.00 -15.00 0.00 0.00 i=2, 0.00 0.00 4.50 0.00 0.00 i=3, 0.00 0.00 3.00 0.00 0.00 i=4, 0.00 0.00 -0.75 0.00 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, -0.19 1.50 -0.00 -1.50 0.19 i=1, -0.62 5.00 -0.00 -5.00 0.62 i=2, 1.12 -9.00 0.00 9.00 -1.12 i=3, -0.38 3.00 -0.00 -3.00 0.38 i=4, 0.06 -0.50 0.00 0.50 -0.06 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, -0.75 12.00 -22.50 12.00 -0.75 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 -0.75 0.00 0.00 i=1, 0.00 0.00 -2.50 0.00 0.00 i=2, 0.00 0.00 4.50 0.00 0.00 i=3, 0.00 0.00 -1.50 0.00 0.00 i=4, 0.00 0.00 0.25 0.00 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.25 -2.00 0.00 2.00 -0.25 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 working row edge i=1, j=2 working row 5, i=1, j=2, bi=-1, bj=0 a1=0.256709, b1=1.27416, c1=5.2802, d1=0.111111, e1=-0.0740741 f1=-2.66667, c=-57.1145 ij=0, u[ii][ij]=69.8814 ij=5, u[ii][ij]=-125.6 ij=10, u[ii][ij]=56.8435 ij=15, u[ii][ij]=-3.14528 ij=1, u[ii][ij]=-11.4674 ij=6, u[ii][ij]=1.65519 ij=11, u[ii][ij]=11.4674 ij=16, u[ii][ij]=-1.43343 ij=2, u[ii][ij]=3.82247 ij=7, u[ii][ij]=0.603459 ij=12, u[ii][ij]=-3.82247 ij=17, u[ii][ij]=0.477809 ij=3, u[ii][ij]=-0.637079 ij=8, u[ii][ij]=-0.164754 ij=13, u[ii][ij]=0.637079 ij=18, u[ii][ij]=-0.0796348 ii=5, ut[ii][ii]=-125.6, ut[ii][cs]=-71.3669 working row edge i=1, j=3 working row 10, i=1, j=3, bi=-1, bj=0 a1=0.358266, b1=1.4, c1=4.97317, d1=0.444444, e1=-2.05433e-33 f1=-2, c=-30.7048 ij=0, u[ii][ij]=-4.60487 ij=5, u[ii][ij]=66.678 ij=10, u[ii][ij]=-120.381 ij=15, u[ii][ij]=52.678 ij=20, u[ii][ij]=-2.85487 ij=1, u[ii][ij]=1.575 ij=6, u[ii][ij]=-12.6 ij=11, u[ii][ij]=3.6122 ij=16, u[ii][ij]=12.6 ij=21, u[ii][ij]=-1.575 ij=2, u[ii][ij]=-0.525 ij=7, u[ii][ij]=4.2 ij=12, u[ii][ij]=0.40813 ij=17, u[ii][ij]=-4.2 ij=22, u[ii][ij]=0.525 ij=3, u[ii][ij]=0.0875 ij=8, u[ii][ij]=-0.7 ij=13, u[ii][ij]=-0.157588 ij=18, u[ii][ij]=0.7 ij=23, u[ii][ij]=-0.0875 ii=10, ut[ii][ii]=-120.381, ut[ii][cs]=-28.0771 working row edge i=1, j=4 working row 15, i=1, j=4, bi=-1, bj=0 a1=0.5, b1=1.27416, c1=4.61044, d1=0.777778, e1=-0.0740741 f1=-1.33333, c=-3.31897 ii=15, ut[ii][ii]=-114.513, ut[ii][cs]=23.9413 build_coef for point ix=3, iy=2 deriv pdxx returns bb=12, aa=-1 4 6 -20 11 deriv pdyy returns bb=12, aa=-1 16 -30 16 -1 deriv pdx returns bb=12, aa=-1 6 -18 10 3 deriv pdy returns bb=12, aa=1 -8 0 8 -1 pdxx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 -0.75 0.00 0.00 i=1, 0.00 0.00 3.00 0.00 0.00 i=2, 0.00 0.00 4.50 0.00 0.00 i=3, 0.00 0.00 -15.00 0.00 0.00 i=4, 0.00 0.00 8.25 0.00 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, -0.06 0.50 -0.00 -0.50 0.06 i=1, 0.38 -3.00 0.00 3.00 -0.38 i=2, -1.12 9.00 -0.00 -9.00 1.12 i=3, 0.62 -5.00 0.00 5.00 -0.62 i=4, 0.19 -1.50 0.00 1.50 -0.19 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, -0.75 12.00 -22.50 12.00 -0.75 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 -0.25 0.00 0.00 i=1, 0.00 0.00 1.50 0.00 0.00 i=2, 0.00 0.00 -4.50 0.00 0.00 i=3, 0.00 0.00 2.50 0.00 0.00 i=4, 0.00 0.00 0.75 0.00 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.25 -2.00 0.00 2.00 -0.25 i=4, 0.00 0.00 0.00 0.00 0.00 working row edge i=5, j=2 working row 9, i=5, j=2, bi=1, bj=0 a1=0.5, b1=1.27416, c1=2.41157, d1=0.111111, e1=0.0740741 f1=1.33333, c=25.104 ii=9, ut[ii][ii]=-60.1492, ut[ii][cs]=-28.9059 working row edge i=5, j=3 working row 14, i=5, j=3, bi=1, bj=0 a1=0.697806, b1=1.4, c1=2.10453, d1=0.444444, e1=2.05433e-33 f1=2, c=54.4672 ii=14, ut[ii][ii]=-54.708, ut[ii][cs]=-52.6363 working row edge i=5, j=4 working row 19, i=5, j=4, bi=1, bj=0 a1=0.973867, b1=1.27416, c1=1.74181, d1=0.777778, e1=0.0740741 f1=2.66667, c=90.5414 ij=6, u[ii][ij]=-0.0796348 ij=11, u[ii][ij]=0.637079 ij=16, u[ii][ij]=-0.924845 ij=21, u[ii][ij]=-0.637079 ij=7, u[ii][ij]=0.477809 ij=12, u[ii][ij]=-3.82247 ij=17, u[ii][ij]=4.08827 ij=22, u[ii][ij]=3.82247 ij=8, u[ii][ij]=-1.43343 ij=13, u[ii][ij]=11.4674 ij=18, u[ii][ij]=0.882402 ij=23, u[ii][ij]=-11.4674 ij=9, u[ii][ij]=-0.491491 ij=14, u[ii][ij]=14.3828 ij=19, u[ii][ij]=-49.1876 ij=24, u[ii][ij]=27.4207 ii=19, ut[ii][ii]=-49.1876, ut[ii][cs]=-67.9943 build_coef for point ix=2, iy=1 deriv pdxx returns bb=12, aa=-1 16 -30 16 -1 deriv pdyy returns bb=12, aa=11 -20 6 4 -1 deriv pdx returns bb=12, aa=1 -8 0 8 -1 deriv pdy returns bb=12, aa=-3 -10 18 -6 1 pdxx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 -0.75 0.00 0.00 0.00 i=1, 0.00 12.00 0.00 0.00 0.00 i=2, 0.00 -22.50 0.00 0.00 0.00 i=3, 0.00 12.00 0.00 0.00 0.00 i=4, 0.00 -0.75 0.00 0.00 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, -0.19 -0.62 1.12 -0.38 0.06 i=1, 1.50 5.00 -9.00 3.00 -0.50 i=2, -0.00 -0.00 0.00 -0.00 0.00 i=3, -1.50 -5.00 9.00 -3.00 0.50 i=4, 0.19 0.62 -1.12 0.38 -0.06 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 8.25 -15.00 4.50 3.00 -0.75 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.25 0.00 0.00 0.00 i=1, 0.00 -2.00 0.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 2.00 0.00 0.00 0.00 i=4, 0.00 -0.25 0.00 0.00 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, -0.75 -2.50 4.50 -1.50 0.25 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 working row edge i=2, j=1 working row 1, i=2, j=1, bi=0, bj=-1 a1=0.217299, b1=1.27416, c1=5.13387, d1=-0.555556, e1=-0.148148 f1=-2.33333, c=-66.2645 ij=0, u[ii][ij]=10.0895 ij=5, u[ii][ij]=-11.4674 ij=10, u[ii][ij]=3.82247 ij=15, u[ii][ij]=-0.637079 ij=1, u[ii][ij]=-83.8603 ij=6, u[ii][ij]=22.4358 ij=11, u[ii][ij]=15.6238 ij=16, u[ii][ij]=-3.88744 ij=2, u[ii][ij]=-4.87431 ij=7, u[ii][ij]=11.4674 ij=12, u[ii][ij]=-3.82247 ij=17, u[ii][ij]=0.637079 ij=3, u[ii][ij]=0.772263 ij=8, u[ii][ij]=-1.43343 ij=13, u[ii][ij]=0.477809 ij=18, u[ii][ij]=-0.0796348 ii=1, ut[ii][ii]=-83.8603, ut[ii][cs]=66.5411 working row edge i=3, j=1 working row 2, i=3, j=1, bi=0, bj=-1 a1=0.256709, b1=1.4, c1=4.56694, d1=-0.666667, e1=-2.46716e-17 f1=-1.33333, c=-44.6463 ij=0, u[ii][ij]=-1.2342 ij=5, u[ii][ij]=1.575 ij=10, u[ii][ij]=-0.525 ij=15, u[ii][ij]=0.0875 ij=1, u[ii][ij]=11.4138 ij=6, u[ii][ij]=-12.6 ij=11, u[ii][ij]=4.2 ij=16, u[ii][ij]=-0.7 ij=2, u[ii][ij]=-75.6133 ij=7, u[ii][ij]=20.5512 ij=12, u[ii][ij]=13.7008 ij=17, u[ii][ij]=-3.4252 ij=3, u[ii][ij]=-5.25283 ij=8, u[ii][ij]=12.6 ij=13, u[ii][ij]=-4.2 ij=18, u[ii][ij]=0.7 ij=4, u[ii][ij]=0.849135 ij=9, u[ii][ij]=-1.575 ij=14, u[ii][ij]=0.525 ij=19, u[ii][ij]=-0.0875 ii=2, ut[ii][ii]=-75.6133, ut[ii][cs]=-1.71905 working row edge i=4, j=1 working row 3, i=4, j=1, bi=0, bj=-1 a1=0.303265, b1=1.27416, c1=3.77571, d1=-0.555556, e1=0.148148 f1=-0.333333, c=-20.4143 ii=3, ut[ii][ii]=-64.1629, ut[ii][cs]=-23.7704 build_coef for point ix=2, iy=3 deriv pdxx returns bb=12, aa=-1 16 -30 16 -1 deriv pdyy returns bb=12, aa=-1 4 6 -20 11 deriv pdx returns bb=12, aa=1 -8 0 8 -1 deriv pdy returns bb=12, aa=-1 6 -18 10 3 pdxx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 -0.75 0.00 i=1, 0.00 0.00 0.00 12.00 0.00 i=2, 0.00 0.00 0.00 -22.50 0.00 i=3, 0.00 0.00 0.00 12.00 0.00 i=4, 0.00 0.00 0.00 -0.75 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, -0.06 0.38 -1.12 0.62 0.19 i=1, 0.50 -3.00 9.00 -5.00 -1.50 i=2, -0.00 0.00 -0.00 0.00 0.00 i=3, -0.50 3.00 -9.00 5.00 1.50 i=4, 0.06 -0.38 1.12 -0.62 -0.19 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, -0.75 3.00 4.50 -15.00 8.25 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.25 0.00 i=1, 0.00 0.00 0.00 -2.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 0.00 2.00 0.00 i=4, 0.00 0.00 0.00 -0.25 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, -0.25 1.50 -4.50 2.50 0.75 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 working row edge i=2, j=5 working row 21, i=2, j=5, bi=0, bj=1 a1=0.824361, b1=1.27416, c1=3.77571, d1=0.777778, e1=-0.148148 f1=0.333333, c=42.5437 ii=21, ut[ii][ii]=-75.2208, ut[ii][cs]=-350.31 working row edge i=3, j=5 working row 22, i=3, j=5, bi=0, bj=1 a1=0.973867, b1=1.4, c1=3.20878, d1=0.666667, e1=-2.46716e-17 f1=1.33333, c=69.1188 ii=22, ut[ii][ii]=-68.7103, ut[ii][cs]=-396.662 working row edge i=4, j=5 working row 23, i=4, j=5, bi=0, bj=1 a1=1.15049, b1=1.27416, c1=2.41755, d1=0.777778, e1=0.148148 f1=2.33333, c=100.52 ij=6, u[ii][ij]=-0.0796348 ij=11, u[ii][ij]=0.477809 ij=16, u[ii][ij]=-1.43343 ij=21, u[ii][ij]=0.127927 ij=7, u[ii][ij]=0.637079 ij=12, u[ii][ij]=-3.82247 ij=17, u[ii][ij]=11.4674 ij=22, u[ii][ij]=5.87951 ij=8, u[ii][ij]=-1.8502 ij=13, u[ii][ij]=7.47487 ij=18, u[ii][ij]=10.2123 ij=23, u[ii][ij]=-59.4455 ij=9, u[ii][ij]=-0.637079 ij=14, u[ii][ij]=3.82247 ij=19, u[ii][ij]=-11.4674 ij=24, u[ii][ij]=21.7322 ii=23, ut[ii][ii]=-59.4455, ut[ii][cs]=-337.037 four sides, one element in, special case, generated build_coef for point ix=1, iy=1 deriv pdxx returns bb=12, aa=11 -20 6 4 -1 deriv pdyy returns bb=12, aa=11 -20 6 4 -1 deriv pdx returns bb=12, aa=-3 -10 18 -6 1 deriv pdy returns bb=12, aa=-3 -10 18 -6 1 pdxx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 8.25 0.00 0.00 0.00 i=1, 0.00 -15.00 0.00 0.00 0.00 i=2, 0.00 4.50 0.00 0.00 0.00 i=3, 0.00 3.00 0.00 0.00 0.00 i=4, 0.00 -0.75 0.00 0.00 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, 0.56 1.88 -3.38 1.12 -0.19 i=1, 1.88 6.25 -11.25 3.75 -0.62 i=2, -3.38 -11.25 20.25 -6.75 1.12 i=3, 1.12 3.75 -6.75 2.25 -0.38 i=4, -0.19 -0.62 1.12 -0.38 0.06 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 8.25 -15.00 4.50 3.00 -0.75 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 -0.75 0.00 0.00 0.00 i=1, 0.00 -2.50 0.00 0.00 0.00 i=2, 0.00 4.50 0.00 0.00 0.00 i=3, 0.00 -1.50 0.00 0.00 0.00 i=4, 0.00 0.25 0.00 0.00 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, -0.75 -2.50 4.50 -1.50 0.25 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 working row 0, i=1, j=1, bi=-1, bj=-1 a1=0.18394, b1=1.00354, c1=5.5401, d1=-0.222222, e1=-0.296296 f1=-3.33333, c=-80.1985 ij=0, u[ii][ij]=-81.6256 ij=5, u[ii][ij]=12.3073 ij=10, u[ii][ij]=20.828 ij=15, u[ii][ij]=-4.85636 ij=1, u[ii][ij]=-11.4621 ij=6, u[ii][ij]=20.3217 ij=11, u[ii][ij]=-6.77389 ij=16, u[ii][ij]=1.12898 ij=2, u[ii][ij]=4.64843 ij=7, u[ii][ij]=-6.77389 ij=12, u[ii][ij]=2.25796 ij=17, u[ii][ij]=-0.376327 ij=3, u[ii][ij]=-0.820723 ij=8, u[ii][ij]=1.12898 ij=13, u[ii][ij]=-0.376327 ij=18, u[ii][ij]=0.0627212 ii=0, ut[ii][ii]=-81.6256, ut[ii][cs]=213.508 build_coef for point ix=1, iy=3 deriv pdxx returns bb=12, aa=11 -20 6 4 -1 deriv pdyy returns bb=12, aa=-1 4 6 -20 11 deriv pdx returns bb=12, aa=-3 -10 18 -6 1 deriv pdy returns bb=12, aa=-1 6 -18 10 3 pdxx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 8.25 0.00 i=1, 0.00 0.00 0.00 -15.00 0.00 i=2, 0.00 0.00 0.00 4.50 0.00 i=3, 0.00 0.00 0.00 3.00 0.00 i=4, 0.00 0.00 0.00 -0.75 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, 0.19 -1.12 3.38 -1.88 -0.56 i=1, 0.62 -3.75 11.25 -6.25 -1.88 i=2, -1.12 6.75 -20.25 11.25 3.38 i=3, 0.38 -2.25 6.75 -3.75 -1.12 i=4, -0.06 0.38 -1.12 0.62 0.19 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, -0.75 3.00 4.50 -15.00 8.25 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 -0.75 0.00 i=1, 0.00 0.00 0.00 -2.50 0.00 i=2, 0.00 0.00 0.00 4.50 0.00 i=3, 0.00 0.00 0.00 -1.50 0.00 i=4, 0.00 0.00 0.00 0.25 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, -0.25 1.50 -4.50 2.50 0.75 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 0.00 0.00 0.00 0.00 0.00 i=4, 0.00 0.00 0.00 0.00 0.00 working row 20, i=1, j=5, bi=-1, bj=1 a1=0.697806, b1=1.00354, c1=4.18194, d1=1.11111, e1=-0.296296 f1=-0.666667, c=22.8574 ii=20, ut[ii][ii]=-83.6535, ut[ii][cs]=-277.927 build_coef for point ix=3, iy=1 deriv pdxx returns bb=12, aa=-1 4 6 -20 11 deriv pdyy returns bb=12, aa=11 -20 6 4 -1 deriv pdx returns bb=12, aa=-1 6 -18 10 3 deriv pdy returns bb=12, aa=-3 -10 18 -6 1 pdxx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 -0.75 0.00 0.00 0.00 i=1, 0.00 3.00 0.00 0.00 0.00 i=2, 0.00 4.50 0.00 0.00 0.00 i=3, 0.00 -15.00 0.00 0.00 0.00 i=4, 0.00 8.25 0.00 0.00 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, 0.19 0.62 -1.12 0.38 -0.06 i=1, -1.12 -3.75 6.75 -2.25 0.38 i=2, 3.38 11.25 -20.25 6.75 -1.12 i=3, -1.88 -6.25 11.25 -3.75 0.62 i=4, -0.56 -1.88 3.38 -1.12 0.19 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, 8.25 -15.00 4.50 3.00 -0.75 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 -0.25 0.00 0.00 0.00 i=1, 0.00 1.50 0.00 0.00 0.00 i=2, 0.00 -4.50 0.00 0.00 0.00 i=3, 0.00 2.50 0.00 0.00 0.00 i=4, 0.00 0.75 0.00 0.00 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, -0.75 -2.50 4.50 -1.50 0.25 i=4, 0.00 0.00 0.00 0.00 0.00 working row 4, i=5, j=1, bi=1, bj=-1 a1=0.358266, b1=1.00354, c1=2.67147, d1=-0.222222, e1=0.296296 f1=0.666667, c=2.46519 ii=4, ut[ii][ii]=-52.3478, ut[ii][cs]=-51.7228 build_coef for point ix=3, iy=3 deriv pdxx returns bb=12, aa=-1 4 6 -20 11 deriv pdyy returns bb=12, aa=-1 4 6 -20 11 deriv pdx returns bb=12, aa=-1 6 -18 10 3 deriv pdy returns bb=12, aa=-1 6 -18 10 3 pdxx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 -0.75 0.00 i=1, 0.00 0.00 0.00 3.00 0.00 i=2, 0.00 0.00 0.00 4.50 0.00 i=3, 0.00 0.00 0.00 -15.00 0.00 i=4, 0.00 0.00 0.00 8.25 0.00 x pdxy j=0 j=1 j=2 j=3 j=4 y i=0, 0.06 -0.38 1.12 -0.62 -0.19 i=1, -0.38 2.25 -6.75 3.75 1.12 i=2, 1.12 -6.75 20.25 -11.25 -3.38 i=3, -0.62 3.75 -11.25 6.25 1.88 i=4, -0.19 1.12 -3.38 1.88 0.56 pdyy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, -0.75 3.00 4.50 -15.00 8.25 i=4, 0.00 0.00 0.00 0.00 0.00 pdx j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 -0.25 0.00 i=1, 0.00 0.00 0.00 1.50 0.00 i=2, 0.00 0.00 0.00 -4.50 0.00 i=3, 0.00 0.00 0.00 2.50 0.00 i=4, 0.00 0.00 0.00 0.75 0.00 pdy j=0 j=1 j=2 j=3 j=4 y i=0, 0.00 0.00 0.00 0.00 0.00 i=1, 0.00 0.00 0.00 0.00 0.00 i=2, 0.00 0.00 0.00 0.00 0.00 i=3, -0.25 1.50 -4.50 2.50 0.75 i=4, 0.00 0.00 0.00 0.00 0.00 working row 24, i=5, j=5, bi=1, bj=1 a1=1.35914, b1=1.00354, c1=1.31331, d1=1.11111, e1=0.296296 f1=3.33333, c=137.698 ij=6, u[ii][ij]=0.0627212 ij=11, u[ii][ij]=-0.376327 ij=16, u[ii][ij]=1.12898 ij=21, u[ii][ij]=-1.92435 ij=7, u[ii][ij]=-0.376327 ij=12, u[ii][ij]=2.25796 ij=17, u[ii][ij]=-6.77389 ij=22, u[ii][ij]=9.50736 ij=8, u[ii][ij]=1.12898 ij=13, u[ii][ij]=-6.77389 ij=18, u[ii][ij]=20.3217 ij=23, u[ii][ij]=-10.1737 ij=9, u[ii][ij]=-1.68627 ij=14, u[ii][ij]=8.14764 ij=19, u[ii][ij]=-6.71327 ij=24, u[ii][ij]=-26.9627 ii=24, ut[ii][ii]=-26.9627, ut[ii][cs]=-520.864 matrix initialized initial matrix i=0, j=0 .. 25 -81.626 -11.462 4.648 -0.821 0.000 12.307 20.322 -6.774 1.129 0.000 20.828 -6.774 2.258 -0.376 0.000 -4.856 1.129 -0.376 0.063 0.000 0.000 0.000 0.000 0.000 0.000 213.508 i=1, j=0 .. 25 10.089 -83.860 -4.874 0.772 0.000 -11.467 22.436 11.467 -1.433 0.000 3.822 15.624 -3.822 0.478 0.000 -0.637 -3.887 0.637 -0.080 0.000 0.000 0.000 0.000 0.000 0.000 66.541 i=2, j=0 .. 25 -1.234 11.414 -75.613 -5.253 0.849 1.575 -12.600 20.551 12.600 -1.575 -0.525 4.200 13.701 -4.200 0.525 0.087 -0.700 -3.425 0.700 -0.087 0.000 0.000 0.000 0.000 0.000 -1.719 i=3, j=0 .. 25 0.000 -1.163 11.121 -64.163 -3.843 0.000 1.433 -11.467 17.657 11.467 0.000 -0.478 3.822 11.105 -3.822 0.000 0.080 -0.637 -2.795 0.637 0.000 0.000 0.000 0.000 0.000 -23.770 i=4, j=0 .. 25 0.000 0.414 -3.022 13.902 -52.348 0.000 -1.129 6.774 -20.322 24.645 0.000 0.376 -2.258 6.774 3.807 0.000 -0.063 0.376 -1.129 -1.302 0.000 0.000 0.000 0.000 0.000 -51.723 i=5, j=0 .. 25 69.881 -11.467 3.822 -0.637 0.000 -125.600 1.655 0.603 -0.165 0.000 56.843 11.467 -3.822 0.637 0.000 -3.145 -1.433 0.478 -0.080 0.000 0.000 0.000 0.000 0.000 0.000 -71.367 i=6, j=0 .. 25 5.465 58.562 -5.465 0.683 0.000 4.084 -118.155 3.195 -0.172 0.000 -5.465 58.414 5.465 -0.683 0.000 0.683 -3.646 -0.683 0.085 0.000 0.000 0.000 0.000 0.000 0.000 -58.892 i=7, j=0 .. 25 -0.700 5.600 51.684 -5.600 0.700 -0.352 4.966 -105.636 3.633 -0.185 0.700 -5.600 51.684 5.600 -0.700 -0.087 0.700 -3.230 -0.700 0.087 0.000 0.000 0.000 0.000 0.000 -24.979 i=8, j=0 .. 25 0.000 -0.683 5.465 42.116 -5.465 0.000 -0.373 5.523 -88.295 4.634 0.000 0.683 -5.465 42.264 5.465 0.000 -0.085 0.683 -2.646 -0.683 0.000 0.000 0.000 0.000 0.000 9.761 i=9, j=0 .. 25 0.000 0.637 -3.822 11.467 22.420 0.000 -0.403 1.667 1.750 -60.149 0.000 -0.637 3.822 -11.467 35.458 0.000 0.080 -0.478 1.433 -2.624 0.000 0.000 0.000 0.000 0.000 -28.906 i=10, j=0 .. 25 -4.605 1.575 -0.525 0.087 0.000 66.678 -12.600 4.200 -0.700 0.000 -120.381 3.612 0.408 -0.158 0.000 52.678 12.600 -4.200 0.700 0.000 -2.855 -1.575 0.525 -0.087 0.000 -28.077 i=11, j=0 .. 25 -0.700 -3.425 0.700 -0.087 0.000 5.600 54.803 -5.600 0.700 0.000 4.857 -113.279 5.301 -0.345 0.000 -5.600 54.803 5.600 -0.700 0.000 0.700 -3.425 -0.700 0.087 0.000 -15.545 i=12, j=0 .. 25 0.087 -0.700 -3.000 0.700 -0.087 -0.700 5.600 48.000 -5.600 0.700 -0.375 6.000 -101.250 6.000 -0.375 0.700 -5.600 48.000 5.600 -0.700 -0.087 0.700 -3.000 -0.700 0.087 7.000 i=13, j=0 .. 25 0.000 0.087 -0.700 -2.407 0.700 0.000 -0.700 5.600 38.505 -5.600 0.000 -0.415 6.866 -84.488 7.310 0.000 0.700 -5.600 38.505 5.600 0.000 -0.087 0.700 -2.407 -0.700 42.591 i=14, j=0 .. 25 0.000 -0.087 0.525 -1.575 -0.703 0.000 0.700 -4.200 12.600 18.254 0.000 -0.634 2.760 1.140 -54.708 0.000 -0.700 4.200 -12.600 32.254 0.000 0.087 -0.525 1.575 -2.453 -52.636 i=15, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 -4.273 1.433 -0.478 0.080 0.000 61.844 -11.467 3.822 -0.637 0.000 -114.513 5.750 0.333 -0.181 0.000 48.806 11.467 -3.822 0.637 0.000 23.941 i=16, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 -0.683 -3.162 0.683 -0.085 0.000 5.465 50.525 -5.465 0.683 0.000 6.199 -108.219 7.977 -0.554 0.000 -5.465 50.377 5.465 -0.683 0.000 56.619 i=17, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 0.087 -0.700 -2.728 0.700 -0.087 -0.700 5.600 43.647 -5.600 0.700 -0.440 7.707 -96.873 9.040 -0.607 0.700 -5.600 43.647 5.600 -0.700 87.945 i=18, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 0.000 0.085 -0.683 -2.125 0.683 0.000 -0.683 5.465 34.079 -5.465 0.000 -0.507 9.003 -80.918 10.781 0.000 0.683 -5.465 34.227 5.465 140.198 i=19, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 0.000 -0.080 0.478 -1.433 -0.491 0.000 0.637 -3.822 11.467 14.383 0.000 -0.925 4.088 0.882 -49.188 0.000 -0.637 3.822 -11.467 27.421 -67.994 i=20, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 -2.435 -1.129 0.376 -0.063 0.000 8.338 6.774 -2.258 0.376 0.000 31.442 -20.322 6.774 -1.129 0.000 -83.654 19.430 -3.337 0.382 0.000 -277.927 i=21, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 0.637 -2.795 -0.637 0.080 0.000 -3.822 11.105 3.822 -0.478 0.000 11.467 17.657 -11.467 1.433 0.000 1.966 -75.221 17.819 -1.609 0.000 -350.310 i=22, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 -0.087 0.700 -2.407 -0.700 0.087 0.525 -4.200 9.626 4.200 -0.525 -1.575 12.600 14.439 -12.600 1.575 0.311 3.353 -68.710 20.020 -1.772 -396.662 i=23, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 0.000 -0.080 0.637 -1.850 -0.637 0.000 0.478 -3.822 7.475 3.822 0.000 -1.433 11.467 10.212 -11.467 0.000 0.128 5.880 -59.446 21.732 -337.037 i=24, j=0 .. 25 0.000 0.000 0.000 0.000 0.000 0.000 0.063 -0.376 1.129 -1.686 0.000 -0.376 2.258 -6.774 8.148 0.000 1.129 -6.774 20.322 -6.713 0.000 -1.924 9.507 -10.174 -26.963 -520.864 computed solution us[x,y] row i=1, at x=-0.666667, first y=-0.666667 -1.41 1.67 3.70 5.15 6.44 row i=2, at x=-0.333333, first y=-0.666667 1.00 3.85 5.96 7.78 9.74 row i=3, at x=-5.55112e-17, first y=-0.666667 2.74 5.59 8.00 10.41 13.26 row i=4, at x=0.333333, first y=-0.666667 4.04 7.11 10.04 13.26 17.22 row i=5, at x=0.666667, first y=-0.666667 5.11 8.63 12.30 16.56 21.85 exact solution u(x,y) -1.41 1.67 3.70 5.15 6.44 1.00 3.85 5.96 7.78 9.74 2.74 5.59 8.00 10.41 13.26 4.04 7.11 10.04 13.26 17.22 5.11 8.63 12.30 16.56 21.85 difference exact-computed solution u(x,y)-us[x,y] -0.00 -0.00 0.00 0.00 0.00 -0.00 0.00 0.00 0.00 0.00 -0.00 -0.00 0.00 -0.00 0.00 -0.00 -0.00 -0.00 -0.00 0.00 -0.00 -0.00 -0.00 -0.00 -0.00 total error=5.30687e-14, average error=2.12275e-15, max error=7.10543e-15