fem_check21a_tri.c running Given 3 ux + 2 uy + u = 2 x + 3 y + 13 Analytic solution u(x,y) = 1 + 2 x + 3 y triquad ready for integration over triangles about to read triangles from 22_tss.tri 1 10 12 index=0, last=1 index=2, last=4 index=5, last=7 tri 0 has vertices 1 10 12 1 12 14 index=0, last=1 index=2, last=4 index=5, last=7 tri 1 has vertices 1 12 14 2 15 17 index=0, last=1 index=2, last=4 index=5, last=7 tri 2 has vertices 2 15 17 4 18 19 index=0, last=1 index=2, last=4 index=5, last=7 tri 3 has vertices 4 18 19 6 11 20 index=0, last=1 index=2, last=4 index=5, last=7 tri 4 has vertices 6 11 20 6 20 22 index=0, last=1 index=2, last=4 index=5, last=7 tri 5 has vertices 6 20 22 6 22 24 index=0, last=1 index=2, last=4 index=5, last=7 tri 6 has vertices 6 22 24 8 13 25 index=0, last=1 index=2, last=4 index=5, last=7 tri 7 has vertices 8 13 25 8 25 27 index=0, last=1 index=2, last=4 index=5, last=7 tri 8 has vertices 8 25 27 9 28 29 index=0, last=1 index=2, last=4 index=5, last=7 tri 9 has vertices 9 28 29 5 16 30 index=0, last=1 index=2, last=4 index=5, last=7 tri 10 has vertices 5 16 30 5 30 21 index=0, last=1 index=2, last=4 index=5, last=7 tri 11 has vertices 5 30 21 10 4 11 index=0, last=2 index=3, last=4 index=5, last=7 tri 12 has vertices 10 4 11 12 10 11 index=0, last=2 index=3, last=5 index=6, last=8 tri 13 has vertices 12 10 11 12 11 6 index=0, last=2 index=3, last=5 index=6, last=7 tri 14 has vertices 12 11 6 12 6 13 index=0, last=2 index=3, last=4 index=5, last=7 tri 15 has vertices 12 6 13 14 12 13 index=0, last=2 index=3, last=5 index=6, last=8 tri 16 has vertices 14 12 13 14 13 8 index=0, last=2 index=3, last=5 index=6, last=7 tri 17 has vertices 14 13 8 15 9 16 index=0, last=2 index=3, last=4 index=5, last=7 tri 18 has vertices 15 9 16 17 15 16 index=0, last=2 index=3, last=5 index=6, last=8 tri 19 has vertices 17 15 16 17 16 5 index=0, last=2 index=3, last=5 index=6, last=7 tri 20 has vertices 17 16 5 18 2 17 index=0, last=2 index=3, last=4 index=5, last=7 tri 21 has vertices 18 2 17 19 18 17 index=0, last=2 index=3, last=5 index=6, last=8 tri 22 has vertices 19 18 17 19 17 5 index=0, last=2 index=3, last=5 index=6, last=7 tri 23 has vertices 19 17 5 11 4 19 index=0, last=2 index=3, last=4 index=5, last=7 tri 24 has vertices 11 4 19 20 11 19 index=0, last=2 index=3, last=5 index=6, last=8 tri 25 has vertices 20 11 19 20 19 5 index=0, last=2 index=3, last=5 index=6, last=7 tri 26 has vertices 20 19 5 20 5 21 index=0, last=2 index=3, last=4 index=5, last=7 tri 27 has vertices 20 5 21 22 20 21 index=0, last=2 index=3, last=5 index=6, last=8 tri 28 has vertices 22 20 21 22 21 0 index=0, last=2 index=3, last=5 index=6, last=7 tri 29 has vertices 22 21 0 22 0 23 index=0, last=2 index=3, last=4 index=5, last=7 tri 30 has vertices 22 0 23 24 22 23 index=0, last=2 index=3, last=5 index=6, last=8 tri 31 has vertices 24 22 23 24 23 7 index=0, last=2 index=3, last=5 index=6, last=7 tri 32 has vertices 24 23 7 13 6 24 index=0, last=2 index=3, last=4 index=5, last=7 tri 33 has vertices 13 6 24 25 13 24 index=0, last=2 index=3, last=5 index=6, last=8 tri 34 has vertices 25 13 24 25 24 7 index=0, last=2 index=3, last=5 index=6, last=7 tri 35 has vertices 25 24 7 25 7 26 index=0, last=2 index=3, last=4 index=5, last=7 tri 36 has vertices 25 7 26 27 25 26 index=0, last=2 index=3, last=5 index=6, last=8 tri 37 has vertices 27 25 26 27 26 3 index=0, last=2 index=3, last=5 index=6, last=7 tri 38 has vertices 27 26 3 28 3 26 index=0, last=2 index=3, last=4 index=5, last=7 tri 39 has vertices 28 3 26 29 28 26 index=0, last=2 index=3, last=5 index=6, last=8 tri 40 has vertices 29 28 26 29 26 7 index=0, last=2 index=3, last=5 index=6, last=7 tri 41 has vertices 29 26 7 16 9 29 index=0, last=2 index=3, last=4 index=5, last=7 tri 42 has vertices 16 9 29 30 16 29 index=0, last=2 index=3, last=5 index=6, last=8 tri 43 has vertices 30 16 29 30 29 7 index=0, last=2 index=3, last=5 index=6, last=7 tri 44 has vertices 30 29 7 30 7 23 index=0, last=2 index=3, last=4 index=5, last=7 tri 45 has vertices 30 7 23 21 30 23 index=0, last=2 index=3, last=5 index=6, last=8 tri 46 has vertices 21 30 23 21 23 0 index=0, last=2 index=3, last=5 index=6, last=7 tri 47 has vertices 21 23 0 48 triangles read from 22_tss.tri subtracting minvert=0 from all vertices, using base zero finding unique of nold=144 found unique of n=31 about to read boundary from 22_tss.bound 1 10 index=0, last=1 index=2, last=4 boundary segment nbound=0 has vertices 1, 10 2 15 index=0, last=1 index=2, last=4 boundary segment nbound=1 has vertices 2, 15 3 27 index=0, last=1 index=2, last=4 boundary segment nbound=2 has vertices 3, 27 4 18 index=0, last=1 index=2, last=4 boundary segment nbound=3 has vertices 4, 18 8 14 index=0, last=1 index=2, last=4 boundary segment nbound=4 has vertices 8, 14 9 28 index=0, last=1 index=2, last=4 boundary segment nbound=5 has vertices 9, 28 10 4 index=0, last=2 index=3, last=4 boundary segment nbound=6 has vertices 10, 4 14 1 index=0, last=2 index=3, last=4 boundary segment nbound=7 has vertices 14, 1 15 9 index=0, last=2 index=3, last=4 boundary segment nbound=8 has vertices 15, 9 18 2 index=0, last=2 index=3, last=4 boundary segment nbound=9 has vertices 18, 2 27 8 index=0, last=2 index=3, last=4 boundary segment nbound=10 has vertices 27, 8 28 3 index=0, last=2 index=3, last=4 boundary segment nbound=11 has vertices 28, 3 read Dirichlet boundaries from 22_tss.bound subtracting minvert=0 from all vertices, using base zero finding unique of nold=24 found unique of n=12 freevert na=31, nb=12 freevert nc free = 19 unique boundary 1 unique boundary 2 unique boundary 3 unique boundary 4 unique boundary 8 unique boundary 9 unique boundary 10 unique boundary 14 unique boundary 15 unique boundary 18 unique boundary 27 unique boundary 28 number of free vertices is 19 free vertex 0 free vertex 5 free vertex 6 free vertex 7 free vertex 11 free vertex 12 free vertex 13 free vertex 16 free vertex 17 free vertex 19 free vertex 20 free vertex 21 free vertex 22 free vertex 23 free vertex 24 free vertex 25 free vertex 26 free vertex 29 free vertex 30 about to read coordinates from 22_tss.coord 0.000000 0.000000 index=0, last=8 index=9, last=17 coordinate 0 at 0.0, 0.0 -1.000000 -1.000000 index=0, last=9 index=10, last=19 coordinate 1 at -1.0, -1.0 0.000000 1.000000 index=0, last=8 index=9, last=17 coordinate 2 at 0.0, 1.0 1.000000 -1.000000 index=0, last=8 index=9, last=18 coordinate 3 at 1.0, -1.0 -0.500000 0.000000 index=0, last=9 index=10, last=18 coordinate 4 at -0.5, 0.0 0.000000 0.500000 index=0, last=8 index=9, last=17 coordinate 5 at 0.0, 0.5 -0.500000 -0.500000 index=0, last=9 index=10, last=19 coordinate 6 at -0.5, -0.5 0.500000 -0.500000 index=0, last=8 index=9, last=18 coordinate 7 at 0.5, -0.5 0.000000 -1.000000 index=0, last=8 index=9, last=18 coordinate 8 at 0.0, -1.0 0.500000 0.000000 index=0, last=8 index=9, last=17 coordinate 9 at 0.5, 0.0 -0.750000 -0.500000 index=0, last=9 index=10, last=19 coordinate 10 at -0.75, -0.5 -0.500000 -0.250000 index=0, last=9 index=10, last=19 coordinate 11 at -0.5, -0.25 -0.750000 -0.750000 index=0, last=9 index=10, last=19 coordinate 12 at -0.75, -0.75 -0.250000 -0.750000 index=0, last=9 index=10, last=19 coordinate 13 at -0.25, -0.75 -0.500000 -1.000000 index=0, last=9 index=10, last=19 coordinate 14 at -0.5, -1.0 0.250000 0.500000 index=0, last=8 index=9, last=17 coordinate 15 at 0.25, 0.5 0.250000 0.250000 index=0, last=8 index=9, last=17 coordinate 16 at 0.25, 0.25 0.000000 0.750000 index=0, last=8 index=9, last=17 coordinate 17 at 0.0, 0.75 -0.250000 0.500000 index=0, last=9 index=10, last=18 coordinate 18 at -0.25, 0.5 -0.250000 0.250000 index=0, last=9 index=10, last=18 coordinate 19 at -0.25, 0.25 -0.250000 0.000000 index=0, last=9 index=10, last=18 coordinate 20 at -0.25, 0.0 0.000000 0.250000 index=0, last=8 index=9, last=17 coordinate 21 at 0.0, 0.25 -0.250000 -0.250000 index=0, last=9 index=10, last=19 coordinate 22 at -0.25, -0.25 0.250000 -0.250000 index=0, last=8 index=9, last=18 coordinate 23 at 0.25, -0.25 0.000000 -0.500000 index=0, last=8 index=9, last=18 coordinate 24 at 0.0, -0.5 0.250000 -0.750000 index=0, last=8 index=9, last=18 coordinate 25 at 0.25, -0.75 0.750000 -0.750000 index=0, last=8 index=9, last=18 coordinate 26 at 0.75, -0.75 0.500000 -1.000000 index=0, last=8 index=9, last=18 coordinate 27 at 0.5, -1.0 0.750000 -0.500000 index=0, last=8 index=9, last=18 coordinate 28 at 0.75, -0.5 0.500000 -0.250000 index=0, last=8 index=9, last=18 coordinate 29 at 0.5, -0.25 0.250000 0.000000 index=0, last=8 index=9, last=17 coordinate 30 at 0.25, 0.0 coordinates read from 22_tss.coord xmin=-1.0, xmax=1.0, ymin=-1.0, ymax=1.0 nvert=31, nbound=12, nuniqueb=12, nfree=19, ntri=48 compute global stiffness matrix computing local stiffness matrix for triangle 0 nodes i1=1, i2=10, i3=12 coord (-1.0,-1.0) (-0.75,-0.5) (-0.75,-0.75) cm matrix of am cm = ym solve for cm -3.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm 0.0 -4.0 4.0 cm matrix of am cm = ym solve for cm 4.0 8.0 -4.0 tri_int1a running with np=3, nv=9 integral=0.1718749999999995, at i=12, j=12 intg fg =0.09505208333333307, at i=12 integral=-0.12239583333333298, at i=1, j=12 integral=-0.03906249999999988, at i=10, j=12 finished k=0 computing local stiffness matrix for triangle 1 nodes i1=1, i2=12, i3=14 coord (-1.0,-1.0) (-0.75,-0.75) (-0.5,-1.0) cm matrix of am cm = ym solve for cm -3.0 -2.0 -2.0 cm matrix of am cm = ym solve for cm 4.0 0.0 4.0 cm matrix of am cm = ym solve for cm 0.0 2.0 -2.0 tri_int1a running with np=3, nv=9 integral=0.177083333333333, at i=12, j=12 intg fg =0.18489583333333298, at i=12 integral=-0.2031249999999996, at i=1, j=12 integral=0.04687499999999991, at i=14, j=12 finished k=1 computing local stiffness matrix for triangle 2 nodes i1=2, i2=15, i3=17 coord (0.0,1.0) (0.25,0.5) (0.0,0.75) cm matrix of am cm = ym solve for cm -3.0 4.0 4.0 cm matrix of am cm = ym solve for cm 0.0 4.0 0.0 cm matrix of am cm = ym solve for cm 4.0 -8.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.3281249999999992, at i=17, j=17 intg fg =0.1601562499999996, at i=17 integral=0.21093749999999947, at i=2, j=17 integral=0.12760416666666635, at i=15, j=17 finished k=2 computing local stiffness matrix for triangle 3 nodes i1=4, i2=18, i3=19 coord (-0.5,0.0) (-0.25,0.5) (-0.25,0.25) cm matrix of am cm = ym solve for cm -1.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm -2.0 -4.0 4.0 cm matrix of am cm = ym solve for cm 4.0 8.0 -4.0 tri_int1a running with np=3, nv=9 integral=0.17187499999999958, at i=19, j=19 intg fg =0.1367187499999997, at i=19 integral=-0.12239583333333305, at i=4, j=19 integral=-0.03906249999999991, at i=18, j=19 finished k=3 computing local stiffness matrix for triangle 4 nodes i1=6, i2=11, i3=20 coord (-0.5,-0.5) (-0.5,-0.25) (-0.25,0.0) cm matrix of am cm = ym solve for cm 1.0 4.0 -4.0 cm matrix of am cm = ym solve for cm -2.0 -8.0 4.0 cm matrix of am cm = ym solve for cm 2.0 4.0 0.0 tri_int1a running with np=3, nv=9 integral=0.04687499999999992, at i=6, j=6 intg fg =0.11653645833333309, at i=6 integral=-0.1640624999999997, at i=11, j=6 integral=0.1276041666666664, at i=20, j=6 integral=-0.16145833333333304, at i=11, j=11 intg fg =0.11848958333333309, at i=11 integral=0.04427083333333326, at i=6, j=11 integral=0.1276041666666664, at i=20, j=11 integral=0.13020833333333304, at i=20, j=20 intg fg =0.12174479166666638, at i=20 integral=0.04427083333333324, at i=6, j=20 integral=-0.1640624999999996, at i=11, j=20 finished k=4 computing local stiffness matrix for triangle 5 nodes i1=6, i2=20, i3=22 coord (-0.5,-0.5) (-0.25,0.0) (-0.25,-0.25) cm matrix of am cm = ym solve for cm -1.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm 0.0 -4.0 4.0 cm matrix of am cm = ym solve for cm 2.0 8.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.11979166666666645, at i=6, j=6 intg fg =0.11783854166666645, at i=6 integral=-0.03906249999999992, at i=20, j=6 integral=0.169270833333333, at i=22, j=6 integral=-0.03645833333333327, at i=20, j=20 intg fg =0.12304687499999979, at i=20 integral=-0.12239583333333315, at i=6, j=20 integral=0.16927083333333307, at i=22, j=20 integral=0.17187499999999956, at i=22, j=22 intg fg =0.1210937499999997, at i=22 integral=-0.12239583333333304, at i=6, j=22 integral=-0.0390624999999999, at i=20, j=22 finished k=5 computing local stiffness matrix for triangle 6 nodes i1=6, i2=22, i3=24 coord (-0.5,-0.5) (-0.25,-0.25) (0.0,-0.5) cm matrix of am cm = ym solve for cm -1.0 -2.0 -2.0 cm matrix of am cm = ym solve for cm 2.0 0.0 4.0 cm matrix of am cm = ym solve for cm 0.0 2.0 -2.0 tri_int1a running with np=3, nv=9 integral=-0.19791666666666632, at i=6, j=6 intg fg =0.23046874999999956, at i=6 integral=0.1718749999999997, at i=22, j=6 integral=0.04687499999999991, at i=24, j=6 integral=0.17708333333333298, at i=22, j=22 intg fg =0.2369791666666662, at i=22 integral=-0.20312499999999967, at i=6, j=22 integral=0.04687499999999992, at i=24, j=22 integral=0.05208333333333322, at i=24, j=24 intg fg =0.2356770833333328, at i=24 integral=-0.20312499999999956, at i=6, j=24 integral=0.1718749999999996, at i=22, j=24 finished k=6 computing local stiffness matrix for triangle 7 nodes i1=8, i2=13, i3=25 coord (0.0,-1.0) (-0.25,-0.75) (0.25,-0.75) cm matrix of am cm = ym solve for cm -3.0 0.0 -4.0 cm matrix of am cm = ym solve for cm 2.0 -2.0 2.0 cm matrix of am cm = ym solve for cm 2.0 2.0 2.0 tri_int1a running with np=3, nv=9 integral=-0.031249999999999938, at i=13, j=13 intg fg =0.2174479166666662, at i=13 integral=-0.16145833333333298, at i=8, j=13 integral=0.2135416666666662, at i=25, j=13 integral=0.21874999999999944, at i=25, j=25 intg fg =0.22265624999999944, at i=25 integral=-0.16145833333333293, at i=8, j=25 integral=-0.036458333333333245, at i=13, j=25 finished k=7 computing local stiffness matrix for triangle 8 nodes i1=8, i2=25, i3=27 coord (0.0,-1.0) (0.25,-0.75) (0.5,-1.0) cm matrix of am cm = ym solve for cm -1.0 -2.0 -2.0 cm matrix of am cm = ym solve for cm 4.0 0.0 4.0 cm matrix of am cm = ym solve for cm -2.0 2.0 -2.0 tri_int1a running with np=3, nv=9 integral=0.177083333333333, at i=25, j=25 intg fg =0.2265624999999995, at i=25 integral=-0.2031249999999996, at i=8, j=25 integral=0.04687499999999991, at i=27, j=25 finished k=8 computing local stiffness matrix for triangle 9 nodes i1=9, i2=28, i3=29 coord (0.5,0.0) (0.75,-0.5) (0.5,-0.25) cm matrix of am cm = ym solve for cm -1.0 4.0 4.0 cm matrix of am cm = ym solve for cm -2.0 4.0 0.0 cm matrix of am cm = ym solve for cm 4.0 -8.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.32812499999999933, at i=29, j=29 intg fg =0.13932291666666638, at i=29 integral=0.21093749999999956, at i=9, j=29 integral=0.1276041666666664, at i=28, j=29 finished k=9 computing local stiffness matrix for triangle 10 nodes i1=5, i2=16, i3=30 coord (0.0,0.5) (0.25,0.25) (0.25,0.0) cm matrix of am cm = ym solve for cm 1.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm -2.0 8.0 4.0 cm matrix of am cm = ym solve for cm 2.0 -4.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.11979166666666644, at i=5, j=5 intg fg =0.14778645833333307, at i=5 integral=0.33593749999999933, at i=16, j=5 integral=-0.20572916666666627, at i=30, j=5 integral=0.338541666666666, at i=16, j=16 intg fg =0.1471354166666664, at i=16 integral=-0.1223958333333331, at i=5, j=16 integral=-0.20572916666666627, at i=30, j=16 integral=-0.20312499999999953, at i=30, j=30 intg fg =0.14518229166666632, at i=30 integral=-0.12239583333333305, at i=5, j=30 integral=0.3359374999999992, at i=16, j=30 finished k=10 computing local stiffness matrix for triangle 11 nodes i1=5, i2=30, i3=21 coord (0.0,0.5) (0.25,0.0) (0.0,0.25) cm matrix of am cm = ym solve for cm -1.0 4.0 4.0 cm matrix of am cm = ym solve for cm 0.0 4.0 0.0 cm matrix of am cm = ym solve for cm 2.0 -8.0 -4.0 tri_int1a running with np=3, nv=9 integral=0.2135416666666662, at i=5, j=5 intg fg =0.14648437499999972, at i=5 integral=0.1276041666666664, at i=30, j=5 integral=-0.33072916666666596, at i=21, j=5 integral=0.1302083333333331, at i=30, j=30 intg fg =0.14388020833333307, at i=30 integral=0.2109374999999996, at i=5, j=30 integral=-0.3307291666666661, at i=21, j=30 integral=-0.32812499999999933, at i=21, j=21 intg fg =0.1445312499999997, at i=21 integral=0.21093749999999956, at i=5, j=21 integral=0.1276041666666664, at i=30, j=21 finished k=11 computing local stiffness matrix for triangle 12 nodes i1=10, i2=4, i3=11 coord (-0.75,-0.5) (-0.5,0.0) (-0.5,-0.25) cm matrix of am cm = ym solve for cm -2.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm -1.0 -4.0 4.0 cm matrix of am cm = ym solve for cm 4.0 8.0 -4.0 tri_int1a running with np=3, nv=9 integral=0.17187499999999942, at i=11, j=11 intg fg =0.1158854166666663, at i=11 integral=-0.12239583333333293, at i=10, j=11 integral=-0.03906249999999987, at i=4, j=11 finished k=12 computing local stiffness matrix for triangle 13 nodes i1=12, i2=10, i3=11 coord (-0.75,-0.75) (-0.75,-0.5) (-0.5,-0.25) cm matrix of am cm = ym solve for cm 1.0 4.0 -4.0 cm matrix of am cm = ym solve for cm -3.0 -8.0 4.0 cm matrix of am cm = ym solve for cm 3.0 4.0 0.0 tri_int1a running with np=3, nv=9 integral=0.04687499999999992, at i=12, j=12 intg fg =0.1035156249999998, at i=12 integral=-0.1640624999999997, at i=10, j=12 integral=0.12760416666666644, at i=11, j=12 integral=0.13020833333333304, at i=11, j=11 intg fg =0.10872395833333309, at i=11 integral=0.04427083333333324, at i=12, j=11 integral=-0.16406249999999964, at i=10, j=11 finished k=13 computing local stiffness matrix for triangle 14 nodes i1=12, i2=11, i3=6 coord (-0.75,-0.75) (-0.5,-0.25) (-0.5,-0.5) cm matrix of am cm = ym solve for cm -2.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm 0.0 -4.0 4.0 cm matrix of am cm = ym solve for cm 3.0 8.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.11979166666666649, at i=12, j=12 intg fg =0.10481770833333318, at i=12 integral=-0.03906249999999994, at i=11, j=12 integral=0.16927083333333307, at i=6, j=12 integral=-0.03645833333333328, at i=11, j=11 intg fg =0.1100260416666665, at i=11 integral=-0.12239583333333316, at i=12, j=11 integral=0.1692708333333331, at i=6, j=11 integral=0.17187499999999944, at i=6, j=6 intg fg =0.10807291666666632, at i=6 integral=-0.12239583333333295, at i=12, j=6 integral=-0.039062499999999875, at i=11, j=6 finished k=14 computing local stiffness matrix for triangle 15 nodes i1=12, i2=6, i3=13 coord (-0.75,-0.75) (-0.5,-0.5) (-0.25,-0.75) cm matrix of am cm = ym solve for cm -2.0 -2.0 -2.0 cm matrix of am cm = ym solve for cm 3.0 0.0 4.0 cm matrix of am cm = ym solve for cm 0.0 2.0 -2.0 tri_int1a running with np=3, nv=9 integral=-0.1979166666666663, at i=12, j=12 intg fg =0.20442708333333293, at i=12 integral=0.17187499999999972, at i=6, j=12 integral=0.04687499999999991, at i=13, j=12 integral=0.177083333333333, at i=6, j=6 intg fg =0.2109374999999996, at i=6 integral=-0.2031249999999996, at i=12, j=6 integral=0.04687499999999991, at i=13, j=6 integral=0.05208333333333322, at i=13, j=13 intg fg =0.2096354166666662, at i=13 integral=-0.20312499999999953, at i=12, j=13 integral=0.17187499999999964, at i=6, j=13 finished k=15 computing local stiffness matrix for triangle 16 nodes i1=14, i2=12, i3=13 coord (-0.5,-1.0) (-0.75,-0.75) (-0.25,-0.75) cm matrix of am cm = ym solve for cm -3.0 0.0 -4.0 cm matrix of am cm = ym solve for cm 1.0 -2.0 2.0 cm matrix of am cm = ym solve for cm 3.0 2.0 2.0 tri_int1a running with np=3, nv=9 integral=-0.03124999999999994, at i=12, j=12 intg fg =0.19661458333333295, at i=12 integral=-0.16145833333333304, at i=14, j=12 integral=0.21354166666666624, at i=13, j=12 integral=0.21874999999999944, at i=13, j=13 intg fg =0.20182291666666616, at i=13 integral=-0.16145833333333293, at i=14, j=13 integral=-0.03645833333333324, at i=12, j=13 finished k=16 computing local stiffness matrix for triangle 17 nodes i1=14, i2=13, i3=8 coord (-0.5,-1.0) (-0.25,-0.75) (0.0,-1.0) cm matrix of am cm = ym solve for cm -2.0 -2.0 -2.0 cm matrix of am cm = ym solve for cm 4.0 0.0 4.0 cm matrix of am cm = ym solve for cm -1.0 2.0 -2.0 tri_int1a running with np=3, nv=9 integral=0.177083333333333, at i=13, j=13 intg fg =0.20572916666666627, at i=13 integral=-0.2031249999999996, at i=14, j=13 integral=0.04687499999999991, at i=8, j=13 finished k=17 computing local stiffness matrix for triangle 18 nodes i1=15, i2=9, i3=16 coord (0.25,0.5) (0.5,0.0) (0.25,0.25) cm matrix of am cm = ym solve for cm -2.0 4.0 4.0 cm matrix of am cm = ym solve for cm -1.0 4.0 0.0 cm matrix of am cm = ym solve for cm 4.0 -8.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.32812499999999933, at i=16, j=16 intg fg =0.149739583333333, at i=16 integral=0.21093749999999958, at i=15, j=16 integral=0.1276041666666664, at i=9, j=16 finished k=18 computing local stiffness matrix for triangle 19 nodes i1=17, i2=15, i3=16 coord (0.0,0.75) (0.25,0.5) (0.25,0.25) cm matrix of am cm = ym solve for cm 1.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm -3.0 8.0 4.0 cm matrix of am cm = ym solve for cm 3.0 -4.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.11979166666666644, at i=17, j=17 intg fg =0.15559895833333304, at i=17 integral=0.33593749999999933, at i=15, j=17 integral=-0.20572916666666627, at i=16, j=17 integral=-0.20312499999999958, at i=16, j=16 intg fg =0.15299479166666635, at i=16 integral=-0.12239583333333309, at i=17, j=16 integral=0.3359374999999993, at i=15, j=16 finished k=19 computing local stiffness matrix for triangle 20 nodes i1=17, i2=16, i3=5 coord (0.0,0.75) (0.25,0.25) (0.0,0.5) cm matrix of am cm = ym solve for cm -2.0 4.0 4.0 cm matrix of am cm = ym solve for cm 0.0 4.0 0.0 cm matrix of am cm = ym solve for cm 3.0 -8.0 -4.0 tri_int1a running with np=3, nv=9 integral=0.2135416666666662, at i=17, j=17 intg fg =0.1542968749999997, at i=17 integral=0.1276041666666664, at i=16, j=17 integral=-0.330729166666666, at i=5, j=17 integral=0.1302083333333331, at i=16, j=16 intg fg =0.15169270833333307, at i=16 integral=0.2109374999999996, at i=17, j=16 integral=-0.3307291666666661, at i=5, j=16 integral=-0.32812499999999933, at i=5, j=5 intg fg =0.15234374999999967, at i=5 integral=0.21093749999999953, at i=17, j=5 integral=0.1276041666666664, at i=16, j=5 finished k=20 computing local stiffness matrix for triangle 21 nodes i1=18, i2=2, i3=17 coord (-0.25,0.5) (0.0,1.0) (0.0,0.75) cm matrix of am cm = ym solve for cm 0.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm -3.0 -4.0 4.0 cm matrix of am cm = ym solve for cm 4.0 8.0 -4.0 tri_int1a running with np=3, nv=9 integral=0.1718749999999997, at i=17, j=17 intg fg =0.15755208333333307, at i=17 integral=-0.12239583333333312, at i=18, j=17 integral=-0.03906249999999993, at i=2, j=17 finished k=21 computing local stiffness matrix for triangle 22 nodes i1=19, i2=18, i3=17 coord (-0.25,0.25) (-0.25,0.5) (0.0,0.75) cm matrix of am cm = ym solve for cm 3.0 4.0 -4.0 cm matrix of am cm = ym solve for cm -3.0 -8.0 4.0 cm matrix of am cm = ym solve for cm 1.0 4.0 0.0 tri_int1a running with np=3, nv=9 integral=0.04687499999999993, at i=19, j=19 intg fg =0.1451822916666664, at i=19 integral=-0.16406249999999972, at i=18, j=19 integral=0.12760416666666646, at i=17, j=19 integral=0.13020833333333304, at i=17, j=17 intg fg =0.15039062499999967, at i=17 integral=0.04427083333333324, at i=19, j=17 integral=-0.1640624999999996, at i=18, j=17 finished k=22 computing local stiffness matrix for triangle 23 nodes i1=19, i2=17, i3=5 coord (-0.25,0.25) (0.0,0.75) (0.0,0.5) cm matrix of am cm = ym solve for cm 0.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm -2.0 -4.0 4.0 cm matrix of am cm = ym solve for cm 3.0 8.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.11979166666666644, at i=19, j=19 intg fg =0.1464843749999997, at i=19 integral=-0.03906249999999992, at i=17, j=19 integral=0.169270833333333, at i=5, j=19 integral=-0.036458333333333266, at i=17, j=17 intg fg =0.151692708333333, at i=17 integral=-0.12239583333333309, at i=19, j=17 integral=0.169270833333333, at i=5, j=17 integral=0.17187499999999967, at i=5, j=5 intg fg =0.14973958333333304, at i=5 integral=-0.12239583333333309, at i=19, j=5 integral=-0.039062499999999924, at i=17, j=5 finished k=23 computing local stiffness matrix for triangle 24 nodes i1=11, i2=4, i3=19 coord (-0.5,-0.25) (-0.5,0.0) (-0.25,0.25) cm matrix of am cm = ym solve for cm 2.0 4.0 -4.0 cm matrix of am cm = ym solve for cm -3.0 -8.0 4.0 cm matrix of am cm = ym solve for cm 2.0 4.0 0.0 tri_int1a running with np=3, nv=9 integral=0.04687499999999992, at i=11, j=11 intg fg =0.12434895833333308, at i=11 integral=-0.16406249999999967, at i=4, j=11 integral=0.1276041666666664, at i=19, j=11 integral=0.13020833333333304, at i=19, j=19 intg fg =0.12955729166666635, at i=19 integral=0.04427083333333324, at i=11, j=19 integral=-0.1640624999999996, at i=4, j=19 finished k=24 computing local stiffness matrix for triangle 25 nodes i1=20, i2=11, i3=19 coord (-0.25,0.0) (-0.5,-0.25) (-0.25,0.25) cm matrix of am cm = ym solve for cm 3.0 8.0 -4.0 cm matrix of am cm = ym solve for cm -1.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm -1.0 -4.0 4.0 tri_int1a running with np=3, nv=9 integral=0.1718749999999997, at i=20, j=20 intg fg =0.12890624999999975, at i=20 integral=-0.12239583333333312, at i=11, j=20 integral=-0.039062499999999924, at i=19, j=20 integral=-0.11979166666666646, at i=11, j=11 intg fg =0.1256510416666664, at i=11 integral=0.16927083333333304, at i=20, j=11 integral=-0.039062499999999924, at i=19, j=11 integral=-0.036458333333333245, at i=19, j=19 intg fg =0.1308593749999997, at i=19 integral=0.16927083333333295, at i=20, j=19 integral=-0.12239583333333307, at i=11, j=19 finished k=25 computing local stiffness matrix for triangle 26 nodes i1=20, i2=19, i3=5 coord (-0.25,0.0) (-0.25,0.25) (0.0,0.5) cm matrix of am cm = ym solve for cm 2.0 4.0 -4.0 cm matrix of am cm = ym solve for cm -2.0 -8.0 4.0 cm matrix of am cm = ym solve for cm 1.0 4.0 0.0 tri_int1a running with np=3, nv=9 integral=0.046874999999999924, at i=20, j=20 intg fg =0.13736979166666638, at i=20 integral=-0.1640624999999997, at i=19, j=20 integral=0.12760416666666644, at i=5, j=20 integral=-0.16145833333333304, at i=19, j=19 intg fg =0.1393229166666664, at i=19 integral=0.04427083333333326, at i=20, j=19 integral=0.12760416666666644, at i=5, j=19 integral=0.13020833333333304, at i=5, j=5 intg fg =0.1425781249999997, at i=5 integral=0.04427083333333324, at i=20, j=5 integral=-0.1640624999999996, at i=19, j=5 finished k=26 computing local stiffness matrix for triangle 27 nodes i1=20, i2=5, i3=21 coord (-0.25,0.0) (0.0,0.5) (0.0,0.25) cm matrix of am cm = ym solve for cm 0.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm -1.0 -4.0 4.0 cm matrix of am cm = ym solve for cm 2.0 8.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.11979166666666644, at i=20, j=20 intg fg =0.13867187499999975, at i=20 integral=-0.03906249999999992, at i=5, j=20 integral=0.169270833333333, at i=21, j=20 integral=-0.03645833333333326, at i=5, j=5 intg fg =0.14388020833333307, at i=5 integral=-0.1223958333333331, at i=20, j=5 integral=0.169270833333333, at i=21, j=5 integral=0.17187499999999964, at i=21, j=21 intg fg =0.14192708333333304, at i=21 integral=-0.12239583333333309, at i=20, j=21 integral=-0.03906249999999992, at i=5, j=21 finished k=27 computing local stiffness matrix for triangle 28 nodes i1=22, i2=20, i3=21 coord (-0.25,-0.25) (-0.25,0.0) (0.0,0.25) cm matrix of am cm = ym solve for cm 1.0 4.0 -4.0 cm matrix of am cm = ym solve for cm -1.0 -8.0 4.0 cm matrix of am cm = ym solve for cm 1.0 4.0 0.0 tri_int1a running with np=3, nv=9 integral=0.04687499999999992, at i=22, j=22 intg fg =0.1295572916666664, at i=22 integral=-0.16406249999999967, at i=20, j=22 integral=0.1276041666666664, at i=21, j=22 integral=-0.16145833333333307, at i=20, j=20 intg fg =0.13151041666666644, at i=20 integral=0.04427083333333326, at i=22, j=20 integral=0.12760416666666644, at i=21, j=20 integral=0.13020833333333304, at i=21, j=21 intg fg =0.1347656249999997, at i=21 integral=0.04427083333333324, at i=22, j=21 integral=-0.1640624999999996, at i=20, j=21 finished k=28 computing local stiffness matrix for triangle 29 nodes i1=22, i2=21, i3=0 coord (-0.25,-0.25) (0.0,0.25) (0.0,0.0) cm matrix of am cm = ym solve for cm 0.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm 0.0 -4.0 4.0 cm matrix of am cm = ym solve for cm 1.0 8.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.11979166666666644, at i=22, j=22 intg fg =0.13085937499999975, at i=22 integral=-0.03906249999999992, at i=21, j=22 integral=0.169270833333333, at i=0, j=22 integral=-0.036458333333333266, at i=21, j=21 intg fg =0.13606770833333307, at i=21 integral=-0.1223958333333331, at i=22, j=21 integral=0.169270833333333, at i=0, j=21 integral=0.17187499999999964, at i=0, j=0 intg fg =0.13411458333333304, at i=0 integral=-0.12239583333333307, at i=22, j=0 integral=-0.03906249999999992, at i=21, j=0 finished k=29 computing local stiffness matrix for triangle 30 nodes i1=22, i2=0, i3=23 coord (-0.25,-0.25) (0.0,0.0) (0.25,-0.25) cm matrix of am cm = ym solve for cm 0.0 -2.0 -2.0 cm matrix of am cm = ym solve for cm 1.0 0.0 4.0 cm matrix of am cm = ym solve for cm 0.0 2.0 -2.0 tri_int1a running with np=3, nv=9 integral=-0.19791666666666632, at i=22, j=22 intg fg =0.25651041666666613, at i=22 integral=0.17187499999999964, at i=0, j=22 integral=0.0468749999999999, at i=23, j=22 integral=0.17708333333333298, at i=0, j=0 intg fg =0.2630208333333328, at i=0 integral=-0.20312499999999967, at i=22, j=0 integral=0.04687499999999992, at i=23, j=0 integral=0.05208333333333322, at i=23, j=23 intg fg =0.26171874999999944, at i=23 integral=-0.20312499999999956, at i=22, j=23 integral=0.1718749999999996, at i=0, j=23 finished k=30 computing local stiffness matrix for triangle 31 nodes i1=24, i2=22, i3=23 coord (0.0,-0.5) (-0.25,-0.25) (0.25,-0.25) cm matrix of am cm = ym solve for cm -1.0 0.0 -4.0 cm matrix of am cm = ym solve for cm 1.0 -2.0 2.0 cm matrix of am cm = ym solve for cm 1.0 2.0 2.0 tri_int1a running with np=3, nv=9 integral=-0.15624999999999967, at i=24, j=24 intg fg =0.24739583333333284, at i=24 integral=-0.03645833333333325, at i=22, j=24 integral=0.21354166666666624, at i=23, j=24 integral=-0.031249999999999944, at i=22, j=22 intg fg =0.24869791666666619, at i=22 integral=-0.16145833333333304, at i=24, j=22 integral=0.21354166666666627, at i=23, j=22 integral=0.2187499999999995, at i=23, j=23 intg fg =0.25390624999999944, at i=23 integral=-0.16145833333333295, at i=24, j=23 integral=-0.03645833333333325, at i=22, j=23 finished k=31 computing local stiffness matrix for triangle 32 nodes i1=24, i2=23, i3=7 coord (0.0,-0.5) (0.25,-0.25) (0.5,-0.5) cm matrix of am cm = ym solve for cm 0.0 -2.0 -2.0 cm matrix of am cm = ym solve for cm 2.0 0.0 4.0 cm matrix of am cm = ym solve for cm -1.0 2.0 -2.0 tri_int1a running with np=3, nv=9 integral=-0.19791666666666632, at i=24, j=24 intg fg =0.25130208333333287, at i=24 integral=0.17187499999999967, at i=23, j=24 integral=0.0468749999999999, at i=7, j=24 integral=0.17708333333333298, at i=23, j=23 intg fg =0.2578124999999995, at i=23 integral=-0.20312499999999967, at i=24, j=23 integral=0.04687499999999992, at i=7, j=23 integral=0.05208333333333322, at i=7, j=7 intg fg =0.25651041666666613, at i=7 integral=-0.20312499999999956, at i=24, j=7 integral=0.1718749999999996, at i=23, j=7 finished k=32 computing local stiffness matrix for triangle 33 nodes i1=13, i2=6, i3=24 coord (-0.25,-0.75) (-0.5,-0.5) (0.0,-0.5) cm matrix of am cm = ym solve for cm -2.0 0.0 -4.0 cm matrix of am cm = ym solve for cm 1.0 -2.0 2.0 cm matrix of am cm = ym solve for cm 2.0 2.0 2.0 tri_int1a running with np=3, nv=9 integral=-0.15624999999999972, at i=13, j=13 intg fg =0.2213541666666663, at i=13 integral=-0.036458333333333266, at i=6, j=13 integral=0.21354166666666635, at i=24, j=13 integral=-0.031249999999999938, at i=6, j=6 intg fg =0.22265624999999956, at i=6 integral=-0.16145833333333298, at i=13, j=6 integral=0.2135416666666662, at i=24, j=6 integral=0.21874999999999944, at i=24, j=24 intg fg =0.22786458333333276, at i=24 integral=-0.16145833333333293, at i=13, j=24 integral=-0.036458333333333245, at i=6, j=24 finished k=33 computing local stiffness matrix for triangle 34 nodes i1=25, i2=13, i3=24 coord (0.25,-0.75) (-0.25,-0.75) (0.0,-0.5) cm matrix of am cm = ym solve for cm -1.0 2.0 -2.0 cm matrix of am cm = ym solve for cm -1.0 -2.0 -2.0 cm matrix of am cm = ym solve for cm 3.0 0.0 4.0 tri_int1a running with np=3, nv=9 integral=0.052083333333333225, at i=25, j=25 intg fg =0.23046874999999958, at i=25 integral=-0.20312499999999964, at i=13, j=25 integral=0.1718749999999997, at i=24, j=25 integral=-0.1979166666666663, at i=13, j=13 intg fg =0.2252604166666663, at i=13 integral=0.04687499999999993, at i=25, j=13 integral=0.1718749999999997, at i=24, j=13 integral=0.17708333333333293, at i=24, j=24 intg fg =0.23177083333333284, at i=24 integral=0.046874999999999896, at i=25, j=24 integral=-0.20312499999999953, at i=13, j=24 finished k=34 computing local stiffness matrix for triangle 35 nodes i1=25, i2=24, i3=7 coord (0.25,-0.75) (0.0,-0.5) (0.5,-0.5) cm matrix of am cm = ym solve for cm -2.0 0.0 -4.0 cm matrix of am cm = ym solve for cm 2.0 -2.0 2.0 cm matrix of am cm = ym solve for cm 1.0 2.0 2.0 tri_int1a running with np=3, nv=9 integral=-0.15624999999999972, at i=25, j=25 intg fg =0.2421874999999996, at i=25 integral=-0.036458333333333266, at i=24, j=25 integral=0.21354166666666635, at i=7, j=25 integral=-0.03124999999999994, at i=24, j=24 intg fg =0.24348958333333284, at i=24 integral=-0.16145833333333298, at i=25, j=24 integral=0.2135416666666662, at i=7, j=24 integral=0.21874999999999944, at i=7, j=7 intg fg =0.24869791666666607, at i=7 integral=-0.16145833333333293, at i=25, j=7 integral=-0.036458333333333245, at i=24, j=7 finished k=35 computing local stiffness matrix for triangle 36 nodes i1=25, i2=7, i3=26 coord (0.25,-0.75) (0.5,-0.5) (0.75,-0.75) cm matrix of am cm = ym solve for cm 0.0 -2.0 -2.0 cm matrix of am cm = ym solve for cm 3.0 0.0 4.0 cm matrix of am cm = ym solve for cm -2.0 2.0 -2.0 tri_int1a running with np=3, nv=9 integral=-0.19791666666666632, at i=25, j=25 intg fg =0.24609374999999956, at i=25 integral=0.17187499999999972, at i=7, j=25 integral=0.04687499999999991, at i=26, j=25 integral=0.177083333333333, at i=7, j=7 intg fg =0.2526041666666662, at i=7 integral=-0.2031249999999996, at i=25, j=7 integral=0.04687499999999991, at i=26, j=7 integral=0.05208333333333321, at i=26, j=26 intg fg =0.2513020833333327, at i=26 integral=-0.20312499999999953, at i=25, j=26 integral=0.17187499999999958, at i=7, j=26 finished k=36 computing local stiffness matrix for triangle 37 nodes i1=27, i2=25, i3=26 coord (0.5,-1.0) (0.25,-0.75) (0.75,-0.75) cm matrix of am cm = ym solve for cm -3.0 0.0 -4.0 cm matrix of am cm = ym solve for cm 3.0 -2.0 2.0 cm matrix of am cm = ym solve for cm 1.0 2.0 2.0 tri_int1a running with np=3, nv=9 integral=-0.031249999999999938, at i=25, j=25 intg fg =0.23828124999999953, at i=25 integral=-0.16145833333333298, at i=27, j=25 integral=0.2135416666666662, at i=26, j=25 integral=0.21874999999999944, at i=26, j=26 intg fg =0.24348958333333273, at i=26 integral=-0.16145833333333293, at i=27, j=26 integral=-0.036458333333333245, at i=25, j=26 finished k=37 computing local stiffness matrix for triangle 38 nodes i1=27, i2=26, i3=3 coord (0.5,-1.0) (0.75,-0.75) (1.0,-1.0) cm matrix of am cm = ym solve for cm 0.0 -2.0 -2.0 cm matrix of am cm = ym solve for cm 4.0 0.0 4.0 cm matrix of am cm = ym solve for cm -3.0 2.0 -2.0 tri_int1a running with np=3, nv=9 integral=0.177083333333333, at i=26, j=26 intg fg =0.24739583333333284, at i=26 integral=-0.2031249999999996, at i=27, j=26 integral=0.04687499999999991, at i=3, j=26 finished k=38 computing local stiffness matrix for triangle 39 nodes i1=28, i2=3, i3=26 coord (0.75,-0.5) (1.0,-1.0) (0.75,-0.75) cm matrix of am cm = ym solve for cm 0.0 4.0 4.0 cm matrix of am cm = ym solve for cm -3.0 4.0 0.0 cm matrix of am cm = ym solve for cm 4.0 -8.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.3281249999999992, at i=26, j=26 intg fg =0.12890624999999967, at i=26 integral=0.2109374999999995, at i=28, j=26 integral=0.12760416666666635, at i=3, j=26 finished k=39 computing local stiffness matrix for triangle 40 nodes i1=29, i2=28, i3=26 coord (0.5,-0.25) (0.75,-0.5) (0.75,-0.75) cm matrix of am cm = ym solve for cm 3.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm -3.0 8.0 4.0 cm matrix of am cm = ym solve for cm 1.0 -4.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.11979166666666649, at i=29, j=29 intg fg =0.13476562499999978, at i=29 integral=0.33593749999999944, at i=28, j=29 integral=-0.20572916666666635, at i=26, j=29 integral=-0.20312499999999956, at i=26, j=26 intg fg =0.13216145833333307, at i=26 integral=-0.12239583333333308, at i=29, j=26 integral=0.33593749999999933, at i=28, j=26 finished k=40 computing local stiffness matrix for triangle 41 nodes i1=29, i2=26, i3=7 coord (0.5,-0.25) (0.75,-0.75) (0.5,-0.5) cm matrix of am cm = ym solve for cm 0.0 4.0 4.0 cm matrix of am cm = ym solve for cm -2.0 4.0 0.0 cm matrix of am cm = ym solve for cm 3.0 -8.0 -4.0 tri_int1a running with np=3, nv=9 integral=0.21354166666666624, at i=29, j=29 intg fg =0.13346354166666638, at i=29 integral=0.12760416666666644, at i=26, j=29 integral=-0.330729166666666, at i=7, j=29 integral=0.1302083333333331, at i=26, j=26 intg fg =0.13085937499999975, at i=26 integral=0.21093749999999958, at i=29, j=26 integral=-0.3307291666666661, at i=7, j=26 integral=-0.32812499999999933, at i=7, j=7 intg fg =0.13151041666666638, at i=7 integral=0.21093749999999953, at i=29, j=7 integral=0.1276041666666664, at i=26, j=7 finished k=41 computing local stiffness matrix for triangle 42 nodes i1=16, i2=9, i3=29 coord (0.25,0.25) (0.5,0.0) (0.5,-0.25) cm matrix of am cm = ym solve for cm 2.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm -3.0 8.0 4.0 cm matrix of am cm = ym solve for cm 2.0 -4.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.11979166666666645, at i=16, j=16 intg fg =0.14518229166666638, at i=16 integral=0.33593749999999933, at i=9, j=16 integral=-0.20572916666666624, at i=29, j=16 integral=-0.20312499999999956, at i=29, j=29 intg fg =0.1425781249999997, at i=29 integral=-0.12239583333333308, at i=16, j=29 integral=0.3359374999999993, at i=9, j=29 finished k=42 computing local stiffness matrix for triangle 43 nodes i1=30, i2=16, i3=29 coord (0.25,0.0) (0.25,0.25) (0.5,-0.25) cm matrix of am cm = ym solve for cm 3.0 -8.0 -4.0 cm matrix of am cm = ym solve for cm -1.0 4.0 4.0 cm matrix of am cm = ym solve for cm -1.0 4.0 0.0 tri_int1a running with np=3, nv=9 integral=-0.3281249999999994, at i=30, j=30 intg fg =0.14192708333333307, at i=30 integral=0.2109374999999996, at i=16, j=30 integral=0.12760416666666644, at i=29, j=30 integral=0.21354166666666627, at i=16, j=16 intg fg =0.14388020833333307, at i=16 integral=-0.3307291666666661, at i=30, j=16 integral=0.1276041666666664, at i=29, j=16 integral=0.13020833333333304, at i=29, j=29 intg fg =0.14127604166666635, at i=29 integral=-0.33072916666666596, at i=30, j=29 integral=0.21093749999999953, at i=16, j=29 finished k=43 computing local stiffness matrix for triangle 44 nodes i1=30, i2=29, i3=7 coord (0.25,0.0) (0.5,-0.25) (0.5,-0.5) cm matrix of am cm = ym solve for cm 2.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm -2.0 8.0 4.0 cm matrix of am cm = ym solve for cm 1.0 -4.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.11979166666666645, at i=30, j=30 intg fg =0.13736979166666644, at i=30 integral=0.33593749999999933, at i=29, j=30 integral=-0.20572916666666624, at i=7, j=30 integral=0.33854166666666596, at i=29, j=29 intg fg =0.13671874999999972, at i=29 integral=-0.1223958333333331, at i=30, j=29 integral=-0.20572916666666627, at i=7, j=29 integral=-0.20312499999999956, at i=7, j=7 intg fg =0.1347656249999997, at i=7 integral=-0.12239583333333307, at i=30, j=7 integral=0.3359374999999992, at i=29, j=7 finished k=44 computing local stiffness matrix for triangle 45 nodes i1=30, i2=7, i3=23 coord (0.25,0.0) (0.5,-0.5) (0.25,-0.25) cm matrix of am cm = ym solve for cm 0.0 4.0 4.0 cm matrix of am cm = ym solve for cm -1.0 4.0 0.0 cm matrix of am cm = ym solve for cm 2.0 -8.0 -4.0 tri_int1a running with np=3, nv=9 integral=0.2135416666666662, at i=30, j=30 intg fg =0.13606770833333307, at i=30 integral=0.12760416666666644, at i=7, j=30 integral=-0.33072916666666596, at i=23, j=30 integral=0.13020833333333307, at i=7, j=7 intg fg =0.1334635416666664, at i=7 integral=0.2109374999999996, at i=30, j=7 integral=-0.330729166666666, at i=23, j=7 integral=-0.32812499999999933, at i=23, j=23 intg fg =0.13411458333333304, at i=23 integral=0.21093749999999956, at i=30, j=23 integral=0.1276041666666664, at i=7, j=23 finished k=45 computing local stiffness matrix for triangle 46 nodes i1=21, i2=30, i3=23 coord (0.0,0.25) (0.25,0.0) (0.25,-0.25) cm matrix of am cm = ym solve for cm 1.0 -4.0 -0.0 cm matrix of am cm = ym solve for cm -1.0 8.0 4.0 cm matrix of am cm = ym solve for cm 1.0 -4.0 -4.0 tri_int1a running with np=3, nv=9 integral=-0.11979166666666644, at i=21, j=21 intg fg =0.13997395833333307, at i=21 integral=0.33593749999999933, at i=30, j=21 integral=-0.20572916666666627, at i=23, j=21 integral=0.33854166666666596, at i=30, j=30 intg fg =0.1393229166666664, at i=30 integral=-0.1223958333333331, at i=21, j=30 integral=-0.20572916666666624, at i=23, j=30 integral=-0.20312499999999956, at i=23, j=23 intg fg =0.13736979166666635, at i=23 integral=-0.12239583333333307, at i=21, j=23 integral=0.3359374999999993, at i=30, j=23 finished k=46 computing local stiffness matrix for triangle 47 nodes i1=21, i2=23, i3=0 coord (0.0,0.25) (0.25,-0.25) (0.0,0.0) cm matrix of am cm = ym solve for cm 0.0 4.0 4.0 cm matrix of am cm = ym solve for cm 0.0 4.0 0.0 cm matrix of am cm = ym solve for cm 1.0 -8.0 -4.0 tri_int1a running with np=3, nv=9 integral=0.21354166666666627, at i=21, j=21 intg fg =0.13867187499999972, at i=21 integral=0.12760416666666644, at i=23, j=21 integral=-0.330729166666666, at i=0, j=21 integral=0.1302083333333331, at i=23, j=23 intg fg =0.13606770833333307, at i=23 integral=0.2109374999999996, at i=21, j=23 integral=-0.3307291666666661, at i=0, j=23 integral=-0.3281249999999992, at i=0, j=0 intg fg =0.1367187499999997, at i=0 integral=0.21093749999999953, at i=21, j=0 integral=0.12760416666666638, at i=23, j=0 finished k=47 tri 0 has vertices 1, 10, 12 tri 1 has vertices 1, 12, 14 tri 2 has vertices 2, 15, 17 tri 3 has vertices 4, 18, 19 tri 4 has vertices 6, 11, 20 tri 5 has vertices 6, 20, 22 tri 6 has vertices 6, 22, 24 tri 7 has vertices 8, 13, 25 tri 8 has vertices 8, 25, 27 tri 9 has vertices 9, 28, 29 tri 10 has vertices 5, 16, 30 tri 11 has vertices 5, 30, 21 tri 12 has vertices 10, 4, 11 tri 13 has vertices 12, 10, 11 tri 14 has vertices 12, 11, 6 tri 15 has vertices 12, 6, 13 tri 16 has vertices 14, 12, 13 tri 17 has vertices 14, 13, 8 tri 18 has vertices 15, 9, 16 tri 19 has vertices 17, 15, 16 tri 20 has vertices 17, 16, 5 tri 21 has vertices 18, 2, 17 tri 22 has vertices 19, 18, 17 tri 23 has vertices 19, 17, 5 tri 24 has vertices 11, 4, 19 tri 25 has vertices 20, 11, 19 tri 26 has vertices 20, 19, 5 tri 27 has vertices 20, 5, 21 tri 28 has vertices 22, 20, 21 tri 29 has vertices 22, 21, 0 tri 30 has vertices 22, 0, 23 tri 31 has vertices 24, 22, 23 tri 32 has vertices 24, 23, 7 tri 33 has vertices 13, 6, 24 tri 34 has vertices 25, 13, 24 tri 35 has vertices 25, 24, 7 tri 36 has vertices 25, 7, 26 tri 37 has vertices 27, 25, 26 tri 38 has vertices 27, 26, 3 tri 39 has vertices 28, 3, 26 tri 40 has vertices 29, 28, 26 tri 41 has vertices 29, 26, 7 tri 42 has vertices 16, 9, 29 tri 43 has vertices 30, 16, 29 tri 44 has vertices 30, 29, 7 tri 45 has vertices 30, 7, 23 tri 46 has vertices 21, 30, 23 tri 47 has vertices 21, 23, 0 vertices, coordinates and analytic values coordinate 0 at 0.0, 0.0, uana=1.0 coordinate 1 at -1.0, -1.0, uana=-4.0 coordinate 2 at 0.0, 1.0, uana=4.0 coordinate 3 at 1.0, -1.0, uana=0.0 coordinate 4 at -0.5, 0.0, uana=0.0 coordinate 5 at 0.0, 0.5, uana=2.5 coordinate 6 at -0.5, -0.5, uana=-1.5 coordinate 7 at 0.5, -0.5, uana=0.5 coordinate 8 at 0.0, -1.0, uana=-2.0 coordinate 9 at 0.5, 0.0, uana=2.0 coordinate 10 at -0.75, -0.5, uana=-2.0 coordinate 11 at -0.5, -0.25, uana=-0.75 coordinate 12 at -0.75, -0.75, uana=-2.75 coordinate 13 at -0.25, -0.75, uana=-1.75 coordinate 14 at -0.5, -1.0, uana=-3.0 coordinate 15 at 0.25, 0.5, uana=3.0 coordinate 16 at 0.25, 0.25, uana=2.25 coordinate 17 at 0.0, 0.75, uana=3.25 coordinate 18 at -0.25, 0.5, uana=2.0 coordinate 19 at -0.25, 0.25, uana=1.25 coordinate 20 at -0.25, 0.0, uana=0.5 coordinate 21 at 0.0, 0.25, uana=1.75 coordinate 22 at -0.25, -0.25, uana=-0.25 coordinate 23 at 0.25, -0.25, uana=0.75 coordinate 24 at 0.0, -0.5, uana=-0.5 coordinate 25 at 0.25, -0.75, uana=-0.75 coordinate 26 at 0.75, -0.75, uana=0.25 coordinate 27 at 0.5, -1.0, uana=-1.0 coordinate 28 at 0.75, -0.5, uana=1.0 coordinate 29 at 0.5, -0.25, uana=1.25 coordinate 30 at 0.25, 0.0, uana=1.5 k computed stiffness matrix, if debug K(0,0)=0.02083333333333337 K(0,1)=0.0 K(0,2)=0.0 K(0,3)=0.0 K(0,4)=0.0 K(0,5)=0.0 K(0,6)=0.0 K(0,7)=0.0 K(0,8)=0.0 K(0,9)=0.0 K(0,10)=0.0 K(0,11)=0.0 K(0,12)=0.0 K(0,13)=0.0 K(0,14)=0.0 K(0,15)=0.0 K(0,16)=0.0 K(0,17)=0.0 K(0,18)=0.0 K(0,19)=0.0 K(0,20)=0.0 K(0,21)=0.1718749999999996 K(0,22)=-0.3255208333333327 K(0,23)=0.1744791666666663 K(0,24)=0.0 K(0,25)=0.0 K(0,26)=0.0 K(0,27)=0.0 K(0,28)=0.0 K(0,29)=0.0 K(0,30)=0.0 K(1,0)=0.0 K(1,1)=1.0 K(1,2)=0.0 K(1,3)=0.0 K(1,4)=0.0 K(1,5)=0.0 K(1,6)=0.0 K(1,7)=0.0 K(1,8)=0.0 K(1,9)=0.0 K(1,10)=0.0 K(1,11)=0.0 K(1,12)=0.0 K(1,13)=0.0 K(1,14)=0.0 K(1,15)=0.0 K(1,16)=0.0 K(1,17)=0.0 K(1,18)=0.0 K(1,19)=0.0 K(1,20)=0.0 K(1,21)=0.0 K(1,22)=0.0 K(1,23)=0.0 K(1,24)=0.0 K(1,25)=0.0 K(1,26)=0.0 K(1,27)=0.0 K(1,28)=0.0 K(1,29)=0.0 K(1,30)=0.0 K(2,0)=0.0 K(2,1)=0.0 K(2,2)=1.0 K(2,3)=0.0 K(2,4)=0.0 K(2,5)=0.0 K(2,6)=0.0 K(2,7)=0.0 K(2,8)=0.0 K(2,9)=0.0 K(2,10)=0.0 K(2,11)=0.0 K(2,12)=0.0 K(2,13)=0.0 K(2,14)=0.0 K(2,15)=0.0 K(2,16)=0.0 K(2,17)=0.0 K(2,18)=0.0 K(2,19)=0.0 K(2,20)=0.0 K(2,21)=0.0 K(2,22)=0.0 K(2,23)=0.0 K(2,24)=0.0 K(2,25)=0.0 K(2,26)=0.0 K(2,27)=0.0 K(2,28)=0.0 K(2,29)=0.0 K(2,30)=0.0 K(3,0)=0.0 K(3,1)=0.0 K(3,2)=0.0 K(3,3)=1.0 K(3,4)=0.0 K(3,5)=0.0 K(3,6)=0.0 K(3,7)=0.0 K(3,8)=0.0 K(3,9)=0.0 K(3,10)=0.0 K(3,11)=0.0 K(3,12)=0.0 K(3,13)=0.0 K(3,14)=0.0 K(3,15)=0.0 K(3,16)=0.0 K(3,17)=0.0 K(3,18)=0.0 K(3,19)=0.0 K(3,20)=0.0 K(3,21)=0.0 K(3,22)=0.0 K(3,23)=0.0 K(3,24)=0.0 K(3,25)=0.0 K(3,26)=0.0 K(3,27)=0.0 K(3,28)=0.0 K(3,29)=0.0 K(3,30)=0.0 K(4,0)=0.0 K(4,1)=0.0 K(4,2)=0.0 K(4,3)=0.0 K(4,4)=1.0 K(4,5)=0.0 K(4,6)=0.0 K(4,7)=0.0 K(4,8)=0.0 K(4,9)=0.0 K(4,10)=0.0 K(4,11)=0.0 K(4,12)=0.0 K(4,13)=0.0 K(4,14)=0.0 K(4,15)=0.0 K(4,16)=0.0 K(4,17)=0.0 K(4,18)=0.0 K(4,19)=0.0 K(4,20)=0.0 K(4,21)=0.0 K(4,22)=0.0 K(4,23)=0.0 K(4,24)=0.0 K(4,25)=0.0 K(4,26)=0.0 K(4,27)=0.0 K(4,28)=0.0 K(4,29)=0.0 K(4,30)=0.0 K(5,0)=0.0 K(5,1)=0.0 K(5,2)=0.0 K(5,3)=0.0 K(5,4)=0.0 K(5,5)=0.03124999999999989 K(5,6)=0.0 K(5,7)=0.0 K(5,8)=0.0 K(5,9)=0.0 K(5,10)=0.0 K(5,11)=0.0 K(5,12)=0.0 K(5,13)=0.0 K(5,14)=0.0 K(5,15)=0.0 K(5,16)=0.46354166666666574 K(5,17)=0.1718749999999996 K(5,18)=0.0 K(5,19)=-0.2864583333333327 K(5,20)=-0.07812499999999986 K(5,21)=-0.16145833333333295 K(5,22)=0.0 K(5,23)=0.0 K(5,24)=0.0 K(5,25)=0.0 K(5,26)=0.0 K(5,27)=0.0 K(5,28)=0.0 K(5,29)=0.0 K(5,30)=-0.07812499999999986 K(6,0)=0.0 K(6,1)=0.0 K(6,2)=0.0 K(6,3)=0.0 K(6,4)=0.0 K(6,5)=0.0 K(6,6)=0.046874999999999646 K(6,7)=0.0 K(6,8)=0.0 K(6,9)=0.0 K(6,10)=0.0 K(6,11)=-0.20312499999999956 K(6,12)=-0.3255208333333326 K(6,13)=-0.11458333333333307 K(6,14)=0.0 K(6,15)=0.0 K(6,16)=0.0 K(6,17)=0.0 K(6,18)=0.0 K(6,19)=0.0 K(6,20)=0.08854166666666649 K(6,21)=0.0 K(6,22)=0.3411458333333327 K(6,23)=0.0 K(6,24)=0.26041666666666613 K(6,25)=0.0 K(6,26)=0.0 K(6,27)=0.0 K(6,28)=0.0 K(6,29)=0.0 K(6,30)=0.0 K(7,0)=0.0 K(7,1)=0.0 K(7,2)=0.0 K(7,3)=0.0 K(7,4)=0.0 K(7,5)=0.0 K(7,6)=0.0 K(7,7)=0.046874999999999806 K(7,8)=0.0 K(7,9)=0.0 K(7,10)=0.0 K(7,11)=0.0 K(7,12)=0.0 K(7,13)=0.0 K(7,14)=0.0 K(7,15)=0.0 K(7,16)=0.0 K(7,17)=0.0 K(7,18)=0.0 K(7,19)=0.0 K(7,20)=0.0 K(7,21)=0.0 K(7,22)=0.0 K(7,23)=-0.1588541666666664 K(7,24)=-0.23958333333333282 K(7,25)=-0.36458333333333254 K(7,26)=0.17447916666666632 K(7,27)=0.0 K(7,28)=0.0 K(7,29)=0.5468749999999988 K(7,30)=0.08854166666666655 K(8,0)=0.0 K(8,1)=0.0 K(8,2)=0.0 K(8,3)=0.0 K(8,4)=0.0 K(8,5)=0.0 K(8,6)=0.0 K(8,7)=0.0 K(8,8)=1.0 K(8,9)=0.0 K(8,10)=0.0 K(8,11)=0.0 K(8,12)=0.0 K(8,13)=0.0 K(8,14)=0.0 K(8,15)=0.0 K(8,16)=0.0 K(8,17)=0.0 K(8,18)=0.0 K(8,19)=0.0 K(8,20)=0.0 K(8,21)=0.0 K(8,22)=0.0 K(8,23)=0.0 K(8,24)=0.0 K(8,25)=0.0 K(8,26)=0.0 K(8,27)=0.0 K(8,28)=0.0 K(8,29)=0.0 K(8,30)=0.0 K(9,0)=0.0 K(9,1)=0.0 K(9,2)=0.0 K(9,3)=0.0 K(9,4)=0.0 K(9,5)=0.0 K(9,6)=0.0 K(9,7)=0.0 K(9,8)=0.0 K(9,9)=1.0 K(9,10)=0.0 K(9,11)=0.0 K(9,12)=0.0 K(9,13)=0.0 K(9,14)=0.0 K(9,15)=0.0 K(9,16)=0.0 K(9,17)=0.0 K(9,18)=0.0 K(9,19)=0.0 K(9,20)=0.0 K(9,21)=0.0 K(9,22)=0.0 K(9,23)=0.0 K(9,24)=0.0 K(9,25)=0.0 K(9,26)=0.0 K(9,27)=0.0 K(9,28)=0.0 K(9,29)=0.0 K(9,30)=0.0 K(10,0)=0.0 K(10,1)=0.0 K(10,2)=0.0 K(10,3)=0.0 K(10,4)=0.0 K(10,5)=0.0 K(10,6)=0.0 K(10,7)=0.0 K(10,8)=0.0 K(10,9)=0.0 K(10,10)=1.0 K(10,11)=0.0 K(10,12)=0.0 K(10,13)=0.0 K(10,14)=0.0 K(10,15)=0.0 K(10,16)=0.0 K(10,17)=0.0 K(10,18)=0.0 K(10,19)=0.0 K(10,20)=0.0 K(10,21)=0.0 K(10,22)=0.0 K(10,23)=0.0 K(10,24)=0.0 K(10,25)=0.0 K(10,26)=0.0 K(10,27)=0.0 K(10,28)=0.0 K(10,29)=0.0 K(10,30)=0.0 K(11,0)=0.0 K(11,1)=0.0 K(11,2)=0.0 K(11,3)=0.0 K(11,4)=-0.20312499999999953 K(11,5)=0.0 K(11,6)=0.21354166666666635 K(11,7)=0.0 K(11,8)=0.0 K(11,9)=0.0 K(11,10)=-0.2864583333333326 K(11,11)=0.031249999999999584 K(11,12)=-0.07812499999999992 K(11,13)=0.0 K(11,14)=0.0 K(11,15)=0.0 K(11,16)=0.0 K(11,17)=0.0 K(11,18)=0.0 K(11,19)=0.08854166666666649 K(11,20)=0.29687499999999944 K(11,21)=0.0 K(11,22)=0.0 K(11,23)=0.0 K(11,24)=0.0 K(11,25)=0.0 K(11,26)=0.0 K(11,27)=0.0 K(11,28)=0.0 K(11,29)=0.0 K(11,30)=0.0 K(12,0)=0.0 K(12,1)=-0.3255208333333326 K(12,2)=0.0 K(12,3)=0.0 K(12,4)=0.0 K(12,5)=0.0 K(12,6)=0.3411458333333328 K(12,7)=0.0 K(12,8)=0.0 K(12,9)=0.0 K(12,10)=-0.20312499999999958 K(12,11)=0.08854166666666649 K(12,12)=0.04687499999999961 K(12,13)=0.26041666666666613 K(12,14)=-0.11458333333333312 K(12,15)=0.0 K(12,16)=0.0 K(12,17)=0.0 K(12,18)=0.0 K(12,19)=0.0 K(12,20)=0.0 K(12,21)=0.0 K(12,22)=0.0 K(12,23)=0.0 K(12,24)=0.0 K(12,25)=0.0 K(12,26)=0.0 K(12,27)=0.0 K(12,28)=0.0 K(12,29)=0.0 K(12,30)=0.0 K(13,0)=0.0 K(13,1)=0.0 K(13,2)=0.0 K(13,3)=0.0 K(13,4)=0.0 K(13,5)=0.0 K(13,6)=0.13541666666666638 K(13,7)=0.0 K(13,8)=-0.11458333333333307 K(13,9)=0.0 K(13,10)=0.0 K(13,11)=0.0 K(13,12)=-0.23958333333333276 K(13,13)=0.06249999999999972 K(13,14)=-0.36458333333333254 K(13,15)=0.0 K(13,16)=0.0 K(13,17)=0.0 K(13,18)=0.0 K(13,19)=0.0 K(13,20)=0.0 K(13,21)=0.0 K(13,22)=0.0 K(13,23)=0.0 K(13,24)=0.3854166666666661 K(13,25)=0.26041666666666613 K(13,26)=0.0 K(13,27)=0.0 K(13,28)=0.0 K(13,29)=0.0 K(13,30)=0.0 K(14,0)=0.0 K(14,1)=0.0 K(14,2)=0.0 K(14,3)=0.0 K(14,4)=0.0 K(14,5)=0.0 K(14,6)=0.0 K(14,7)=0.0 K(14,8)=0.0 K(14,9)=0.0 K(14,10)=0.0 K(14,11)=0.0 K(14,12)=0.0 K(14,13)=0.0 K(14,14)=1.0 K(14,15)=0.0 K(14,16)=0.0 K(14,17)=0.0 K(14,18)=0.0 K(14,19)=0.0 K(14,20)=0.0 K(14,21)=0.0 K(14,22)=0.0 K(14,23)=0.0 K(14,24)=0.0 K(14,25)=0.0 K(14,26)=0.0 K(14,27)=0.0 K(14,28)=0.0 K(14,29)=0.0 K(14,30)=0.0 K(15,0)=0.0 K(15,1)=0.0 K(15,2)=0.0 K(15,3)=0.0 K(15,4)=0.0 K(15,5)=0.0 K(15,6)=0.0 K(15,7)=0.0 K(15,8)=0.0 K(15,9)=0.0 K(15,10)=0.0 K(15,11)=0.0 K(15,12)=0.0 K(15,13)=0.0 K(15,14)=0.0 K(15,15)=1.0 K(15,16)=0.0 K(15,17)=0.0 K(15,18)=0.0 K(15,19)=0.0 K(15,20)=0.0 K(15,21)=0.0 K(15,22)=0.0 K(15,23)=0.0 K(15,24)=0.0 K(15,25)=0.0 K(15,26)=0.0 K(15,27)=0.0 K(15,28)=0.0 K(15,29)=0.0 K(15,30)=0.0 K(16,0)=0.0 K(16,1)=0.0 K(16,2)=0.0 K(16,3)=0.0 K(16,4)=0.0 K(16,5)=-0.45312499999999917 K(16,6)=0.0 K(16,7)=0.0 K(16,8)=0.0 K(16,9)=0.46354166666666574 K(16,10)=0.0 K(16,11)=0.0 K(16,12)=0.0 K(16,13)=0.0 K(16,14)=0.0 K(16,15)=0.5468749999999989 K(16,16)=0.03125000000000003 K(16,17)=0.08854166666666652 K(16,18)=0.0 K(16,19)=0.0 K(16,20)=0.0 K(16,21)=0.0 K(16,22)=0.0 K(16,23)=0.0 K(16,24)=0.0 K(16,25)=0.0 K(16,26)=0.0 K(16,27)=0.0 K(16,28)=0.0 K(16,29)=-0.07812499999999983 K(16,30)=-0.5364583333333324 K(17,0)=0.0 K(17,1)=0.0 K(17,2)=0.17187499999999956 K(17,3)=0.0 K(17,4)=0.0 K(17,5)=-0.161458333333333 K(17,6)=0.0 K(17,7)=0.0 K(17,8)=0.0 K(17,9)=0.0 K(17,10)=0.0 K(17,11)=0.0 K(17,12)=0.0 K(17,13)=0.0 K(17,14)=0.0 K(17,15)=0.4635416666666657 K(17,16)=-0.07812499999999986 K(17,17)=0.03125000000000005 K(17,18)=-0.2864583333333327 K(17,19)=-0.07812499999999986 K(17,20)=0.0 K(17,21)=0.0 K(17,22)=0.0 K(17,23)=0.0 K(17,24)=0.0 K(17,25)=0.0 K(17,26)=0.0 K(17,27)=0.0 K(17,28)=0.0 K(17,29)=0.0 K(17,30)=0.0 K(18,0)=0.0 K(18,1)=0.0 K(18,2)=0.0 K(18,3)=0.0 K(18,4)=0.0 K(18,5)=0.0 K(18,6)=0.0 K(18,7)=0.0 K(18,8)=0.0 K(18,9)=0.0 K(18,10)=0.0 K(18,11)=0.0 K(18,12)=0.0 K(18,13)=0.0 K(18,14)=0.0 K(18,15)=0.0 K(18,16)=0.0 K(18,17)=0.0 K(18,18)=1.0 K(18,19)=0.0 K(18,20)=0.0 K(18,21)=0.0 K(18,22)=0.0 K(18,23)=0.0 K(18,24)=0.0 K(18,25)=0.0 K(18,26)=0.0 K(18,27)=0.0 K(18,28)=0.0 K(18,29)=0.0 K(18,30)=0.0 K(19,0)=0.0 K(19,1)=0.0 K(19,2)=0.0 K(19,3)=0.0 K(19,4)=-0.28645833333333265 K(19,5)=0.29687499999999944 K(19,6)=0.0 K(19,7)=0.0 K(19,8)=0.0 K(19,9)=0.0 K(19,10)=0.0 K(19,11)=-0.07812499999999983 K(19,12)=0.0 K(19,13)=0.0 K(19,14)=0.0 K(19,15)=0.0 K(19,16)=0.0 K(19,17)=0.08854166666666655 K(19,18)=-0.20312499999999964 K(19,19)=0.031249999999999833 K(19,20)=0.2135416666666662 K(19,21)=0.0 K(19,22)=0.0 K(19,23)=0.0 K(19,24)=0.0 K(19,25)=0.0 K(19,26)=0.0 K(19,27)=0.0 K(19,28)=0.0 K(19,29)=0.0 K(19,30)=0.0 K(20,0)=0.0 K(20,1)=0.0 K(20,2)=0.0 K(20,3)=0.0 K(20,4)=0.0 K(20,5)=0.08854166666666652 K(20,6)=-0.07812499999999992 K(20,7)=0.0 K(20,8)=0.0 K(20,9)=0.0 K(20,10)=0.0 K(20,11)=-0.2864583333333327 K(20,12)=0.0 K(20,13)=0.0 K(20,14)=0.0 K(20,15)=0.0 K(20,16)=0.0 K(20,17)=0.0 K(20,18)=0.0 K(20,19)=-0.2031249999999996 K(20,20)=0.03124999999999989 K(20,21)=0.29687499999999944 K(20,22)=0.21354166666666632 K(20,23)=0.0 K(20,24)=0.0 K(20,25)=0.0 K(20,26)=0.0 K(20,27)=0.0 K(20,28)=0.0 K(20,29)=0.0 K(20,30)=0.0 K(21,0)=-0.161458333333333 K(21,1)=0.0 K(21,2)=0.0 K(21,3)=0.0 K(21,4)=0.0 K(21,5)=0.17187499999999964 K(21,6)=0.0 K(21,7)=0.0 K(21,8)=0.0 K(21,9)=0.0 K(21,10)=0.0 K(21,11)=0.0 K(21,12)=0.0 K(21,13)=0.0 K(21,14)=0.0 K(21,15)=0.0 K(21,16)=0.0 K(21,17)=0.0 K(21,18)=0.0 K(21,19)=0.0 K(21,20)=-0.2864583333333327 K(21,21)=0.031249999999999917 K(21,22)=-0.07812499999999986 K(21,23)=-0.07812499999999983 K(21,24)=0.0 K(21,25)=0.0 K(21,26)=0.0 K(21,27)=0.0 K(21,28)=0.0 K(21,29)=0.0 K(21,30)=0.46354166666666574 K(22,0)=0.34114583333333265 K(22,1)=0.0 K(22,2)=0.0 K(22,3)=0.0 K(22,4)=0.0 K(22,5)=0.0 K(22,6)=-0.3255208333333327 K(22,7)=0.0 K(22,8)=0.0 K(22,9)=0.0 K(22,10)=0.0 K(22,11)=0.0 K(22,12)=0.0 K(22,13)=0.0 K(22,14)=0.0 K(22,15)=0.0 K(22,16)=0.0 K(22,17)=0.0 K(22,18)=0.0 K(22,19)=0.0 K(22,20)=-0.20312499999999956 K(22,21)=0.08854166666666649 K(22,22)=0.046874999999999806 K(22,23)=0.2604166666666662 K(22,24)=-0.11458333333333312 K(22,25)=0.0 K(22,26)=0.0 K(22,27)=0.0 K(22,28)=0.0 K(22,29)=0.0 K(22,30)=0.0 K(23,0)=-0.15885416666666646 K(23,1)=0.0 K(23,2)=0.0 K(23,3)=0.0 K(23,4)=0.0 K(23,5)=0.0 K(23,6)=0.0 K(23,7)=0.17447916666666632 K(23,8)=0.0 K(23,9)=0.0 K(23,10)=0.0 K(23,11)=0.0 K(23,12)=0.0 K(23,13)=0.0 K(23,14)=0.0 K(23,15)=0.0 K(23,16)=0.0 K(23,17)=0.0 K(23,18)=0.0 K(23,19)=0.0 K(23,20)=0.0 K(23,21)=0.08854166666666655 K(23,22)=-0.23958333333333282 K(23,23)=0.04687499999999989 K(23,24)=-0.3645833333333326 K(23,25)=0.0 K(23,26)=0.0 K(23,27)=0.0 K(23,28)=0.0 K(23,29)=0.0 K(23,30)=0.5468749999999989 K(24,0)=0.0 K(24,1)=0.0 K(24,2)=0.0 K(24,3)=0.0 K(24,4)=0.0 K(24,5)=0.0 K(24,6)=-0.23958333333333282 K(24,7)=0.26041666666666613 K(24,8)=0.0 K(24,9)=0.0 K(24,10)=0.0 K(24,11)=0.0 K(24,12)=0.0 K(24,13)=-0.3645833333333325 K(24,14)=0.0 K(24,15)=0.0 K(24,16)=0.0 K(24,17)=0.0 K(24,18)=0.0 K(24,19)=0.0 K(24,20)=0.0 K(24,21)=0.0 K(24,22)=0.13541666666666635 K(24,23)=0.3854166666666659 K(24,24)=0.06249999999999967 K(24,25)=-0.11458333333333309 K(24,26)=0.0 K(24,27)=0.0 K(24,28)=0.0 K(24,29)=0.0 K(24,30)=0.0 K(25,0)=0.0 K(25,1)=0.0 K(25,2)=0.0 K(25,3)=0.0 K(25,4)=0.0 K(25,5)=0.0 K(25,6)=0.0 K(25,7)=0.3854166666666661 K(25,8)=-0.36458333333333254 K(25,9)=0.0 K(25,10)=0.0 K(25,11)=0.0 K(25,12)=0.0 K(25,13)=-0.23958333333333287 K(25,14)=0.0 K(25,15)=0.0 K(25,16)=0.0 K(25,17)=0.0 K(25,18)=0.0 K(25,19)=0.0 K(25,20)=0.0 K(25,21)=0.0 K(25,22)=0.0 K(25,23)=0.0 K(25,24)=0.13541666666666644 K(25,25)=0.0624999999999997 K(25,26)=0.26041666666666613 K(25,27)=-0.11458333333333307 K(25,28)=0.0 K(25,29)=0.0 K(25,30)=0.0 K(26,0)=0.0 K(26,1)=0.0 K(26,2)=0.0 K(26,3)=0.17447916666666627 K(26,4)=0.0 K(26,5)=0.0 K(26,6)=0.0 K(26,7)=-0.1588541666666665 K(26,8)=0.0 K(26,9)=0.0 K(26,10)=0.0 K(26,11)=0.0 K(26,12)=0.0 K(26,13)=0.0 K(26,14)=0.0 K(26,15)=0.0 K(26,16)=0.0 K(26,17)=0.0 K(26,18)=0.0 K(26,19)=0.0 K(26,20)=0.0 K(26,21)=0.0 K(26,22)=0.0 K(26,23)=0.0 K(26,24)=0.0 K(26,25)=-0.23958333333333276 K(26,26)=0.046874999999999944 K(26,27)=-0.36458333333333254 K(26,28)=0.5468749999999989 K(26,29)=0.0885416666666665 K(26,30)=0.0 K(27,0)=0.0 K(27,1)=0.0 K(27,2)=0.0 K(27,3)=0.0 K(27,4)=0.0 K(27,5)=0.0 K(27,6)=0.0 K(27,7)=0.0 K(27,8)=0.0 K(27,9)=0.0 K(27,10)=0.0 K(27,11)=0.0 K(27,12)=0.0 K(27,13)=0.0 K(27,14)=0.0 K(27,15)=0.0 K(27,16)=0.0 K(27,17)=0.0 K(27,18)=0.0 K(27,19)=0.0 K(27,20)=0.0 K(27,21)=0.0 K(27,22)=0.0 K(27,23)=0.0 K(27,24)=0.0 K(27,25)=0.0 K(27,26)=0.0 K(27,27)=1.0 K(27,28)=0.0 K(27,29)=0.0 K(27,30)=0.0 K(28,0)=0.0 K(28,1)=0.0 K(28,2)=0.0 K(28,3)=0.0 K(28,4)=0.0 K(28,5)=0.0 K(28,6)=0.0 K(28,7)=0.0 K(28,8)=0.0 K(28,9)=0.0 K(28,10)=0.0 K(28,11)=0.0 K(28,12)=0.0 K(28,13)=0.0 K(28,14)=0.0 K(28,15)=0.0 K(28,16)=0.0 K(28,17)=0.0 K(28,18)=0.0 K(28,19)=0.0 K(28,20)=0.0 K(28,21)=0.0 K(28,22)=0.0 K(28,23)=0.0 K(28,24)=0.0 K(28,25)=0.0 K(28,26)=0.0 K(28,27)=0.0 K(28,28)=1.0 K(28,29)=0.0 K(28,30)=0.0 K(29,0)=0.0 K(29,1)=0.0 K(29,2)=0.0 K(29,3)=0.0 K(29,4)=0.0 K(29,5)=0.0 K(29,6)=0.0 K(29,7)=-0.5364583333333323 K(29,8)=0.0 K(29,9)=0.5468749999999989 K(29,10)=0.0 K(29,11)=0.0 K(29,12)=0.0 K(29,13)=0.0 K(29,14)=0.0 K(29,15)=0.0 K(29,16)=0.08854166666666645 K(29,17)=0.0 K(29,18)=0.0 K(29,19)=0.0 K(29,20)=0.0 K(29,21)=0.0 K(29,22)=0.0 K(29,23)=0.0 K(29,24)=0.0 K(29,25)=0.0 K(29,26)=-0.07812499999999992 K(29,27)=0.0 K(29,28)=0.46354166666666585 K(29,29)=0.031249999999999833 K(29,30)=-0.45312499999999906 K(30,0)=0.0 K(30,1)=0.0 K(30,2)=0.0 K(30,3)=0.0 K(30,4)=0.0 K(30,5)=0.08854166666666656 K(30,6)=0.0 K(30,7)=-0.0781249999999998 K(30,8)=0.0 K(30,9)=0.0 K(30,10)=0.0 K(30,11)=0.0 K(30,12)=0.0 K(30,13)=0.0 K(30,14)=0.0 K(30,15)=0.0 K(30,16)=0.5468749999999989 K(30,17)=0.0 K(30,18)=0.0 K(30,19)=0.0 K(30,20)=0.0 K(30,21)=-0.45312499999999917 K(30,22)=0.0 K(30,23)=-0.5364583333333321 K(30,24)=0.0 K(30,25)=0.0 K(30,26)=0.0 K(30,27)=0.0 K(30,28)=0.0 K(30,29)=0.46354166666666574 K(30,30)=0.03124999999999989 f computed forcing function and boundary, if debug F(0)=0.5338541666666655, ug[0]=0.0 F(1)=-4.0, ug[1]=0.0 F(2)=4.0, ug[2]=0.0 F(3)=0.0, ug[3]=0.0 F(4)=0.0, ug[4]=0.0 F(5)=0.8828124999999982, ug[5]=0.0 F(6)=1.0065104166666645, ug[6]=0.0 F(7)=1.1575520833333308, ug[7]=0.0 F(8)=-2.0, ug[8]=0.0 F(9)=2.0, ug[9]=0.0 F(10)=-2.0, ug[10]=0.0 F(11)=0.7031249999999984, ug[11]=0.0 F(12)=0.8893229166666649, ug[12]=0.0 F(13)=1.2812499999999976, ug[13]=0.0 F(14)=-3.0, ug[14]=0.0 F(15)=3.0, ug[15]=0.0 F(16)=0.8906249999999983, ug[16]=0.0 F(17)=0.9296874999999981, ug[17]=0.0 F(18)=2.0, ug[18]=0.0 F(19)=0.8281249999999982, ug[19]=0.0 F(20)=0.7812499999999986, ug[20]=0.0 F(21)=0.8359374999999982, ug[21]=0.0 F(22)=1.1236979166666643, ug[22]=0.0 F(23)=1.1809895833333308, ug[23]=0.0 F(24)=1.437499999999997, ug[24]=0.0 F(25)=1.4062499999999973, ug[25]=0.0 F(26)=1.1341145833333308, ug[26]=0.0 F(27)=-1.0, ug[27]=0.0 F(28)=1.0, ug[28]=0.0 F(29)=0.8281249999999984, ug[29]=0.0 F(30)=0.8437499999999983, ug[30]=0.0 ug computed Galerkin, Ua analytic, error ug[0]=0.9999999999999967, Ua=1.0, err=-3.3306690738754696E-15 ug[1]=-4.0, Ua=-4.0, err=0.0 ug[2]=4.0, Ua=4.0, err=0.0 ug[3]=0.0, Ua=0.0, err=0.0 ug[4]=0.0, Ua=0.0, err=0.0 ug[5]=2.499999999999998, Ua=2.5, err=-2.220446049250313E-15 ug[6]=-1.500000000000003, Ua=-1.5, err=-3.1086244689504383E-15 ug[7]=0.49999999999999867, Ua=0.5, err=-1.3322676295501878E-15 ug[8]=-2.0, Ua=-2.0, err=0.0 ug[9]=2.0, Ua=2.0, err=0.0 ug[10]=-2.0, Ua=-2.0, err=0.0 ug[11]=-0.75, Ua=-0.75, err=0.0 ug[12]=-2.7499999999999987, Ua=-2.75, err=1.3322676295501878E-15 ug[13]=-1.7499999999999996, Ua=-1.75, err=4.440892098500626E-16 ug[14]=-3.0, Ua=-3.0, err=0.0 ug[15]=3.0, Ua=3.0, err=0.0 ug[16]=2.2500000000000004, Ua=2.25, err=4.440892098500626E-16 ug[17]=3.25, Ua=3.25, err=0.0 ug[18]=2.0, Ua=2.0, err=0.0 ug[19]=1.2499999999999993, Ua=1.25, err=-6.661338147750939E-16 ug[20]=0.5000000000000011, Ua=0.5, err=1.1102230246251565E-15 ug[21]=1.75, Ua=1.75, err=0.0 ug[22]=-0.24999999999999956, Ua=-0.25, err=4.440892098500626E-16 ug[23]=0.75, Ua=0.75, err=0.0 ug[24]=-0.49999999999999933, Ua=-0.5, err=6.661338147750939E-16 ug[25]=-0.7499999999999986, Ua=-0.75, err=1.4432899320127035E-15 ug[26]=0.2500000000000009, Ua=0.25, err=8.881784197001252E-16 ug[27]=-1.0, Ua=-1.0, err=0.0 ug[28]=1.0, Ua=1.0, err=0.0 ug[29]=1.2500000000000004, Ua=1.25, err=4.440892098500626E-16 ug[30]=1.5000000000000004, Ua=1.5, err=4.440892098500626E-16 maxerr=3.3306690738754696E-15, avgerr=5.9092515826822845E-16