fem_check22p_tri.c running uxx(x,y)+uyy(x,y)=f(x,y) u(x,y)=x*x+2*x*y+3*y*y+4*x+5*y+6 f(x,y)=8 triquad ready for integration over triangles about to read triangles from 9p.tri 1 2 3 index=1, last=2 index=3, last=4 index=5, last=6 tri 0 has vertices 1 2 3 1 3 4 index=1, last=2 index=3, last=4 index=5, last=6 tri 1 has vertices 1 3 4 1 4 5 index=1, last=2 index=3, last=4 index=5, last=6 tri 2 has vertices 1 4 5 1 5 6 index=1, last=2 index=3, last=4 index=5, last=6 tri 3 has vertices 1 5 6 1 6 7 index=1, last=2 index=3, last=4 index=5, last=6 tri 4 has vertices 1 6 7 1 7 8 index=1, last=2 index=3, last=4 index=5, last=6 tri 5 has vertices 1 7 8 1 8 9 index=1, last=2 index=3, last=4 index=5, last=6 tri 6 has vertices 1 8 9 1 9 2 index=1, last=2 index=3, last=4 index=5, last=6 tri 7 has vertices 1 9 2 8 triangles read from 9p.tri subtracting minvert=1 from all vertices, using base zero finding unique of nold=24 found unique of n=9 about to read boundary from 9p.bound 2 3 index=1, last=2 index=3, last=4 boundary segment nbound=0 has vertices 2, 3 3 4 index=1, last=2 index=3, last=4 boundary segment nbound=1 has vertices 3, 4 4 5 index=1, last=2 index=3, last=4 boundary segment nbound=2 has vertices 4, 5 5 6 index=1, last=2 index=3, last=4 boundary segment nbound=3 has vertices 5, 6 6 7 index=1, last=2 index=3, last=4 boundary segment nbound=4 has vertices 6, 7 7 8 index=1, last=2 index=3, last=4 boundary segment nbound=5 has vertices 7, 8 8 9 index=1, last=2 index=3, last=4 boundary segment nbound=6 has vertices 8, 9 9 2 index=1, last=2 index=3, last=4 boundary segment nbound=7 has vertices 9, 2 read Dirichlet boundaries from 9p.bound subtracting minvert=1 from all vertices, using base zero finding unique of nold=16 found unique of n=8 freevert na=9, nb=8 freevert nc free = 1 unique boundary 1 unique boundary 2 unique boundary 3 unique boundary 4 unique boundary 5 unique boundary 6 unique boundary 7 unique boundary 8 number of free vertices is 1 free vertex 0 about to read coordinates from 9p.coord 0.5 0.5 index=1, last=4 index=5, last=8 coordinate 0 at 0.5, 0.5 0.0 0.0 index=1, last=4 index=5, last=8 coordinate 1 at 0.0, 0.0 0.0 0.5 index=1, last=4 index=5, last=8 coordinate 2 at 0.0, 0.5 0.0 1.0 index=1, last=4 index=5, last=8 coordinate 3 at 0.0, 1.0 0.5 1.0 index=1, last=4 index=5, last=8 coordinate 4 at 0.5, 1.0 1.0 1.0 index=1, last=4 index=5, last=8 coordinate 5 at 1.0, 1.0 1.0 0.5 index=1, last=4 index=5, last=8 coordinate 6 at 1.0, 0.5 1.0 0.0 index=1, last=4 index=5, last=8 coordinate 7 at 1.0, 0.0 0.5 0.0 index=1, last=4 index=5, last=8 coordinate 8 at 0.5, 0.0 coordinates read from 9p.coord xmin=0.0, xmax=1.0, ymin=0.0, ymax=1.0 nvert=9, nbound=8, nuniqueb=8, nfree=1, ntri=8 compute global stiffness matrix computing local stiffness matrix for triangle 0 nodes i1=0, i2=1, i3=2 coord (0.5,0.5) (0.0,0.0) (0.0,0.5) cm matrix of am cm = ym solve for cm 0.0 -2.0 0.0 8.0 0.0 -0.0 cm matrix of am cm = ym solve for cm 1.0 0.0 -6.0 0.0 0.0 8.0 cm matrix of am cm = ym solve for cm 0.0 2.0 -2.0 8.0 -16.0 8.0 tri_int1a running with np=3, nv=9 i=0, xs=0.3938297308804235, ys=0.405795297023964, fs=8.0 ph1=0.4531553936419275, ph2=-0.11741319744964662, ph3=-0.022785734101997246 gs1=16.0, gs2=16.0, gs3=32.0 i=1, xs=0.39382973088042394, ys=0.446914865440212, fs=8.0 ph1=0.45315539364192925, ph2=-0.0836260170297296, ph3=-0.08362601702972983 gs1=16.0, gs2=16.0, gs3=32.0 i=2, xs=0.3938297308804235, ys=0.4880344338564595, fs=8.0 ph1=0.4531553936419275, ph2=-0.022785734101997468, ph3=-0.11741319744964662 gs1=16.0, gs2=16.0, gs3=32.0 i=3, xs=0.20473343222036747, ys=0.23801046613994947, fs=8.0 ph1=-0.07414063829088052, ph2=0.02512905909755142, ph3=-0.05769517994728385 gs1=16.0, gs2=16.0, gs3=32.0 i=4, xs=0.204733432220367, ys=0.3523667161101835, fs=8.0 ph1=-0.07414063829088113, ph2=-0.12090187568290367, ph3=-0.12090187568290367 gs1=16.0, gs2=16.0, gs3=32.0 i=5, xs=0.204733432220367, ys=0.4667229660804175, fs=8.0 ph1=-0.07414063829088113, ph2=-0.05769517994728468, ph3=0.025129059097551476 gs1=16.0, gs2=16.0, gs3=32.0 i=6, xs=0.04429397975635198, ys=0.09565280716116498, fs=8.0 ph1=-0.07289230637145501, ph2=0.49927883317549826, ph3=-0.08161582159044704 gs1=16.0, gs2=16.0, gs3=32.0 i=7, xs=0.04429397975635199, ys=0.272146989878176, fs=8.0 ph1=-0.07289230637145502, ph2=-0.04037006647103969, ph3=-0.0403700664710398 gs1=16.0, gs2=16.0, gs3=32.0 i=8, xs=0.044293979756351964, ys=0.448641172595187, fs=8.0 ph1=-0.072892306371455, ph2=-0.08161582159044722, ph3=0.4992788331754985 gs1=16.0, gs2=16.0, gs3=32.0 i=0, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=16.0, B=-0.072892306371455, ws=0.013953605120761 integral=-2.0816681711721685E-17, at i=0, j=0 i=0, A=8.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=8.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=8.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=8.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=8.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=8.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=8.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=8.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=8.0, B=-0.072892306371455, ws=0.013953605120761 intg fg =-1.0408340855860843E-17, at i=0 i=0, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=16.0, B=-0.072892306371455, ws=0.013953605120761 integral=-2.0816681711721685E-17, at i=1, j=0 i=0, A=32.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=32.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=32.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=32.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=32.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=32.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=32.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=32.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=32.0, B=-0.072892306371455, ws=0.013953605120761 integral=-4.163336342344337E-17, at i=2, j=0 finished k=0 computing local stiffness matrix for triangle 1 nodes i1=0, i2=2, i3=3 coord (0.5,0.5) (0.0,0.5) (0.0,1.0) cm matrix of am cm = ym solve for cm 0.0 -2.0 0.0 8.0 0.0 -0.0 cm matrix of am cm = ym solve for cm 6.0 -14.0 -14.0 8.0 16.0 8.0 cm matrix of am cm = ym solve for cm 3.0 0.0 -10.0 0.0 0.0 8.0 tri_int1a running with np=3, nv=9 i=0, xs=0.3938297308804235, ys=0.5119655661435405, fs=8.0 ph1=0.4531553936419275, ph2=-0.11741319744964596, ph3=-0.022785734101997246 gs1=16.0, gs2=32.0, gs3=16.0 i=1, xs=0.39382973088042394, ys=0.553085134559788, fs=8.0 ph1=0.45315539364192925, ph2=-0.08362601702972983, ph3=-0.08362601702972983 gs1=16.0, gs2=32.0, gs3=16.0 i=2, xs=0.3938297308804235, ys=0.594204702976036, fs=8.0 ph1=0.4531553936419275, ph2=-0.022785734101997246, ph3=-0.11741319744964684 gs1=16.0, gs2=32.0, gs3=16.0 i=3, xs=0.20473343222036747, ys=0.533277033919582, fs=8.0 ph1=-0.07414063829088052, ph2=0.025129059097550588, ph3=-0.05769517994728357 gs1=16.0, gs2=32.0, gs3=16.0 i=4, xs=0.204733432220367, ys=0.6476332838898164, fs=8.0 ph1=-0.07414063829088113, ph2=-0.12090187568290345, ph3=-0.12090187568290389 gs1=16.0, gs2=32.0, gs3=16.0 i=5, xs=0.204733432220367, ys=0.7619895338600505, fs=8.0 ph1=-0.07414063829088113, ph2=-0.05769517994728446, ph3=0.025129059097551476 gs1=16.0, gs2=32.0, gs3=16.0 i=6, xs=0.04429397975635198, ys=0.551358827404813, fs=8.0 ph1=-0.07289230637145501, ph2=0.49927883317549826, ph3=-0.08161582159044745 gs1=16.0, gs2=32.0, gs3=16.0 i=7, xs=0.04429397975635199, ys=0.727853010121824, fs=8.0 ph1=-0.07289230637145502, ph2=-0.04037006647104047, ph3=-0.04037006647103958 gs1=16.0, gs2=32.0, gs3=16.0 i=8, xs=0.044293979756351964, ys=0.904347192838835, fs=8.0 ph1=-0.072892306371455, ph2=-0.08161582159044656, ph3=0.4992788331754987 gs1=16.0, gs2=32.0, gs3=16.0 i=0, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=16.0, B=-0.072892306371455, ws=0.013953605120761 integral=-2.0816681711721685E-17, at i=0, j=0 i=0, A=8.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=8.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=8.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=8.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=8.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=8.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=8.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=8.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=8.0, B=-0.072892306371455, ws=0.013953605120761 intg fg =-1.0408340855860843E-17, at i=0 i=0, A=32.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=32.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=32.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=32.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=32.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=32.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=32.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=32.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=32.0, B=-0.072892306371455, ws=0.013953605120761 integral=-4.163336342344337E-17, at i=2, j=0 i=0, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=16.0, B=-0.072892306371455, ws=0.013953605120761 integral=-2.0816681711721685E-17, at i=3, j=0 finished k=1 computing local stiffness matrix for triangle 2 nodes i1=0, i2=3, i3=4 coord (0.5,0.5) (0.0,1.0) (0.5,1.0) cm matrix of am cm = ym solve for cm 6.0 0.0 -14.0 0.0 0.0 8.0 cm matrix of am cm = ym solve for cm 1.0 -6.0 0.0 8.0 0.0 -0.0 cm matrix of am cm = ym solve for cm 10.0 -18.0 -18.0 8.0 16.0 8.0 tri_int1a running with np=3, nv=9 i=0, xs=0.405795297023964, ys=0.6061702691195765, fs=8.0 ph1=0.45315539364192814, ph2=-0.11741319744964662, ph3=-0.022785734101998578 gs1=16.0, gs2=16.0, gs3=32.0 i=1, xs=0.446914865440212, ys=0.606170269119576, fs=8.0 ph1=0.45315539364192947, ph2=-0.0836260170297296, ph3=-0.0836260170297285 gs1=16.0, gs2=16.0, gs3=32.0 i=2, xs=0.4880344338564595, ys=0.6061702691195765, fs=8.0 ph1=0.45315539364192814, ph2=-0.022785734101997468, ph3=-0.11741319744964729 gs1=16.0, gs2=16.0, gs3=32.0 i=3, xs=0.23801046613994947, ys=0.7952665677796326, fs=8.0 ph1=-0.07414063829087958, ph2=0.02512905909755142, ph3=-0.05769517994728268 gs1=16.0, gs2=16.0, gs3=32.0 i=4, xs=0.3523667161101835, ys=0.7952665677796329, fs=8.0 ph1=-0.07414063829088136, ph2=-0.12090187568290367, ph3=-0.120901875682903 gs1=16.0, gs2=16.0, gs3=32.0 i=5, xs=0.4667229660804175, ys=0.795266567779633, fs=8.0 ph1=-0.07414063829088136, ph2=-0.05769517994728468, ph3=0.025129059097552364 gs1=16.0, gs2=16.0, gs3=32.0 i=6, xs=0.09565280716116498, ys=0.955706020243648, fs=8.0 ph1=-0.07289230637145572, ph2=0.49927883317549826, ph3=-0.08161582159044656 gs1=16.0, gs2=16.0, gs3=32.0 i=7, xs=0.272146989878176, ys=0.955706020243648, fs=8.0 ph1=-0.07289230637145572, ph2=-0.04037006647103969, ph3=-0.04037006647103869 gs1=16.0, gs2=16.0, gs3=32.0 i=8, xs=0.448641172595187, ys=0.955706020243648, fs=8.0 ph1=-0.07289230637145572, ph2=-0.08161582159044722, ph3=0.4992788331754987 gs1=16.0, gs2=16.0, gs3=32.0 i=0, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 integral=-3.608224830031759E-16, at i=0, j=0 i=0, A=8.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=8.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=8.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=8.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=8.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=8.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=8.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=8.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=8.0, B=-0.07289230637145572, ws=0.013953605120761 intg fg =-1.8041124150158794E-16, at i=0 i=0, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 integral=-3.608224830031759E-16, at i=3, j=0 i=0, A=32.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=32.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=32.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=32.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=32.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=32.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=32.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=32.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=32.0, B=-0.07289230637145572, ws=0.013953605120761 integral=-7.216449660063518E-16, at i=4, j=0 finished k=2 computing local stiffness matrix for triangle 3 nodes i1=0, i2=4, i3=5 coord (0.5,0.5) (0.5,1.0) (1.0,1.0) cm matrix of am cm = ym solve for cm 6.0 0.0 -14.0 0.0 0.0 8.0 cm matrix of am cm = ym solve for cm 0.0 2.0 -2.0 8.0 -16.0 8.0 cm matrix of am cm = ym solve for cm 3.0 -10.0 0.0 8.0 0.0 -0.0 tri_int1a running with np=3, nv=9 i=0, xs=0.5119655661435405, ys=0.6061702691195765, fs=8.0 ph1=0.45315539364192814, ph2=-0.11741319744964684, ph3=-0.022785734101997246 gs1=16.0, gs2=32.0, gs3=16.0 i=1, xs=0.553085134559788, ys=0.606170269119576, fs=8.0 ph1=0.45315539364192947, ph2=-0.08362601702972983, ph3=-0.08362601702972983 gs1=16.0, gs2=32.0, gs3=16.0 i=2, xs=0.594204702976036, ys=0.6061702691195765, fs=8.0 ph1=0.45315539364192814, ph2=-0.022785734101998134, ph3=-0.11741319744964684 gs1=16.0, gs2=32.0, gs3=16.0 i=3, xs=0.533277033919582, ys=0.7952665677796326, fs=8.0 ph1=-0.07414063829087958, ph2=0.025129059097551476, ph3=-0.05769517994728357 gs1=16.0, gs2=32.0, gs3=16.0 i=4, xs=0.6476332838898164, ys=0.7952665677796329, fs=8.0 ph1=-0.07414063829088136, ph2=-0.12090187568290389, ph3=-0.12090187568290389 gs1=16.0, gs2=32.0, gs3=16.0 i=5, xs=0.7619895338600505, ys=0.795266567779633, fs=8.0 ph1=-0.07414063829088136, ph2=-0.05769517994728446, ph3=0.025129059097551476 gs1=16.0, gs2=32.0, gs3=16.0 i=6, xs=0.551358827404813, ys=0.955706020243648, fs=8.0 ph1=-0.07289230637145572, ph2=0.4992788331754978, ph3=-0.08161582159044745 gs1=16.0, gs2=32.0, gs3=16.0 i=7, xs=0.727853010121824, ys=0.955706020243648, fs=8.0 ph1=-0.07289230637145572, ph2=-0.04037006647104047, ph3=-0.04037006647103958 gs1=16.0, gs2=32.0, gs3=16.0 i=8, xs=0.904347192838835, ys=0.955706020243648, fs=8.0 ph1=-0.07289230637145572, ph2=-0.08161582159044745, ph3=0.4992788331754987 gs1=16.0, gs2=32.0, gs3=16.0 i=0, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 integral=-3.608224830031759E-16, at i=0, j=0 i=0, A=8.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=8.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=8.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=8.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=8.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=8.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=8.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=8.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=8.0, B=-0.07289230637145572, ws=0.013953605120761 intg fg =-1.8041124150158794E-16, at i=0 i=0, A=32.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=32.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=32.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=32.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=32.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=32.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=32.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=32.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=32.0, B=-0.07289230637145572, ws=0.013953605120761 integral=-7.216449660063518E-16, at i=4, j=0 i=0, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 integral=-3.608224830031759E-16, at i=5, j=0 finished k=3 computing local stiffness matrix for triangle 4 nodes i1=0, i2=5, i3=6 coord (0.5,0.5) (1.0,1.0) (1.0,0.5) cm matrix of am cm = ym solve for cm 6.0 -14.0 0.0 8.0 0.0 -0.0 cm matrix of am cm = ym solve for cm 3.0 0.0 -10.0 0.0 0.0 8.0 cm matrix of am cm = ym solve for cm 0.0 -2.0 2.0 8.0 -16.0 8.0 tri_int1a running with np=3, nv=9 i=0, xs=0.6061702691195765, ys=0.594204702976036, fs=8.0 ph1=0.45315539364192814, ph2=-0.11741319744964684, ph3=-0.022785734101998134 gs1=16.0, gs2=16.0, gs3=32.0 i=1, xs=0.606170269119576, ys=0.553085134559788, fs=8.0 ph1=0.45315539364192947, ph2=-0.08362601702972983, ph3=-0.08362601702972983 gs1=16.0, gs2=16.0, gs3=32.0 i=2, xs=0.6061702691195765, ys=0.5119655661435405, fs=8.0 ph1=0.45315539364192814, ph2=-0.022785734101997246, ph3=-0.11741319744964684 gs1=16.0, gs2=16.0, gs3=32.0 i=3, xs=0.7952665677796326, ys=0.7619895338600505, fs=8.0 ph1=-0.07414063829087958, ph2=0.025129059097551476, ph3=-0.05769517994728268 gs1=16.0, gs2=16.0, gs3=32.0 i=4, xs=0.7952665677796329, ys=0.6476332838898164, fs=8.0 ph1=-0.07414063829088136, ph2=-0.12090187568290389, ph3=-0.12090187568290389 gs1=16.0, gs2=16.0, gs3=32.0 i=5, xs=0.795266567779633, ys=0.5332770339195825, fs=8.0 ph1=-0.07414063829088136, ph2=-0.0576951799472849, ph3=0.025129059097551032 gs1=16.0, gs2=16.0, gs3=32.0 i=6, xs=0.955706020243648, ys=0.904347192838835, fs=8.0 ph1=-0.07289230637145572, ph2=0.4992788331754987, ph3=-0.08161582159044745 gs1=16.0, gs2=16.0, gs3=32.0 i=7, xs=0.955706020243648, ys=0.727853010121824, fs=8.0 ph1=-0.07289230637145572, ph2=-0.04037006647103958, ph3=-0.04037006647104047 gs1=16.0, gs2=16.0, gs3=32.0 i=8, xs=0.955706020243648, ys=0.551358827404813, fs=8.0 ph1=-0.07289230637145572, ph2=-0.08161582159044745, ph3=0.4992788331754978 gs1=16.0, gs2=16.0, gs3=32.0 i=0, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 integral=-3.608224830031759E-16, at i=0, j=0 i=0, A=8.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=8.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=8.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=8.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=8.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=8.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=8.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=8.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=8.0, B=-0.07289230637145572, ws=0.013953605120761 intg fg =-1.8041124150158794E-16, at i=0 i=0, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 integral=-3.608224830031759E-16, at i=5, j=0 i=0, A=32.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=32.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=32.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=32.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=32.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=32.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=32.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=32.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=32.0, B=-0.07289230637145572, ws=0.013953605120761 integral=-7.216449660063518E-16, at i=6, j=0 finished k=4 computing local stiffness matrix for triangle 5 nodes i1=0, i2=6, i3=7 coord (0.5,0.5) (1.0,0.5) (1.0,0.0) cm matrix of am cm = ym solve for cm 6.0 -14.0 0.0 8.0 0.0 -0.0 cm matrix of am cm = ym solve for cm 10.0 -18.0 -18.0 8.0 16.0 8.0 cm matrix of am cm = ym solve for cm 1.0 0.0 -6.0 0.0 0.0 8.0 tri_int1a running with np=3, nv=9 i=0, xs=0.6061702691195765, ys=0.4880344338564595, fs=8.0 ph1=0.45315539364192814, ph2=-0.11741319744964662, ph3=-0.022785734101997468 gs1=16.0, gs2=32.0, gs3=16.0 i=1, xs=0.606170269119576, ys=0.446914865440212, fs=8.0 ph1=0.45315539364192947, ph2=-0.08362601702972872, ph3=-0.0836260170297296 gs1=16.0, gs2=32.0, gs3=16.0 i=2, xs=0.6061702691195765, ys=0.405795297023964, fs=8.0 ph1=0.45315539364192814, ph2=-0.022785734101997912, ph3=-0.11741319744964662 gs1=16.0, gs2=32.0, gs3=16.0 i=3, xs=0.7952665677796326, ys=0.466722966080418, fs=8.0 ph1=-0.07414063829087958, ph2=0.02512905909755192, ph3=-0.057695179947284014 gs1=16.0, gs2=32.0, gs3=16.0 i=4, xs=0.7952665677796329, ys=0.3523667161101835, fs=8.0 ph1=-0.07414063829088136, ph2=-0.12090187568290323, ph3=-0.12090187568290367 gs1=16.0, gs2=32.0, gs3=16.0 i=5, xs=0.795266567779633, ys=0.23801046613994947, fs=8.0 ph1=-0.07414063829088136, ph2=-0.057695179947283404, ph3=0.02512905909755142 gs1=16.0, gs2=32.0, gs3=16.0 i=6, xs=0.955706020243648, ys=0.448641172595187, fs=8.0 ph1=-0.07289230637145572, ph2=0.49927883317549804, ph3=-0.08161582159044722 gs1=16.0, gs2=32.0, gs3=16.0 i=7, xs=0.955706020243648, ys=0.272146989878176, fs=8.0 ph1=-0.07289230637145572, ph2=-0.04037006647103947, ph3=-0.04037006647103969 gs1=16.0, gs2=32.0, gs3=16.0 i=8, xs=0.955706020243648, ys=0.09565280716116498, fs=8.0 ph1=-0.07289230637145572, ph2=-0.08161582159044657, ph3=0.49927883317549826 gs1=16.0, gs2=32.0, gs3=16.0 i=0, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 integral=-3.608224830031759E-16, at i=0, j=0 i=0, A=8.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=8.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=8.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=8.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=8.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=8.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=8.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=8.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=8.0, B=-0.07289230637145572, ws=0.013953605120761 intg fg =-1.8041124150158794E-16, at i=0 i=0, A=32.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=32.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=32.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=32.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=32.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=32.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=32.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=32.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=32.0, B=-0.07289230637145572, ws=0.013953605120761 integral=-7.216449660063518E-16, at i=6, j=0 i=0, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192947, ws=0.00775855332238375 i=2, A=16.0, B=0.45315539364192814, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829087958, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088136, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088136, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145572, ws=0.02232576819321775 i=8, A=16.0, B=-0.07289230637145572, ws=0.013953605120761 integral=-3.608224830031759E-16, at i=7, j=0 finished k=5 computing local stiffness matrix for triangle 6 nodes i1=0, i2=7, i3=8 coord (0.5,0.5) (1.0,0.0) (0.5,0.0) cm matrix of am cm = ym solve for cm 0.0 0.0 -2.0 0.0 0.0 8.0 cm matrix of am cm = ym solve for cm 3.0 -10.0 0.0 8.0 0.0 -0.0 cm matrix of am cm = ym solve for cm 6.0 -14.0 -14.0 8.0 16.0 8.0 tri_int1a running with np=3, nv=9 i=0, xs=0.594204702976036, ys=0.3938297308804235, fs=8.0 ph1=0.4531553936419275, ph2=-0.11741319744964684, ph3=-0.022785734101997468 gs1=16.0, gs2=16.0, gs3=32.0 i=1, xs=0.553085134559788, ys=0.39382973088042394, fs=8.0 ph1=0.45315539364192925, ph2=-0.08362601702972983, ph3=-0.0836260170297296 gs1=16.0, gs2=16.0, gs3=32.0 i=2, xs=0.5119655661435405, ys=0.3938297308804235, fs=8.0 ph1=0.4531553936419275, ph2=-0.022785734101997246, ph3=-0.11741319744964573 gs1=16.0, gs2=16.0, gs3=32.0 i=3, xs=0.7619895338600505, ys=0.20473343222036747, fs=8.0 ph1=-0.07414063829088052, ph2=0.025129059097551476, ph3=-0.05769517994728468 gs1=16.0, gs2=16.0, gs3=32.0 i=4, xs=0.6476332838898164, ys=0.204733432220367, fs=8.0 ph1=-0.07414063829088113, ph2=-0.12090187568290389, ph3=-0.12090187568290345 gs1=16.0, gs2=16.0, gs3=32.0 i=5, xs=0.5332770339195825, ys=0.204733432220367, fs=8.0 ph1=-0.07414063829088113, ph2=-0.0576951799472849, ph3=0.025129059097552364 gs1=16.0, gs2=16.0, gs3=32.0 i=6, xs=0.904347192838835, ys=0.04429397975635198, fs=8.0 ph1=-0.07289230637145501, ph2=0.4992788331754987, ph3=-0.08161582159044727 gs1=16.0, gs2=16.0, gs3=32.0 i=7, xs=0.727853010121824, ys=0.04429397975635199, fs=8.0 ph1=-0.07289230637145502, ph2=-0.04037006647103958, ph3=-0.04037006647104016 gs1=16.0, gs2=16.0, gs3=32.0 i=8, xs=0.551358827404813, ys=0.044293979756351964, fs=8.0 ph1=-0.072892306371455, ph2=-0.08161582159044745, ph3=0.4992788331754983 gs1=16.0, gs2=16.0, gs3=32.0 i=0, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=16.0, B=-0.072892306371455, ws=0.013953605120761 integral=-2.0816681711721685E-17, at i=0, j=0 i=0, A=8.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=8.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=8.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=8.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=8.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=8.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=8.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=8.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=8.0, B=-0.072892306371455, ws=0.013953605120761 intg fg =-1.0408340855860843E-17, at i=0 i=0, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=16.0, B=-0.072892306371455, ws=0.013953605120761 integral=-2.0816681711721685E-17, at i=7, j=0 i=0, A=32.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=32.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=32.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=32.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=32.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=32.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=32.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=32.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=32.0, B=-0.072892306371455, ws=0.013953605120761 integral=-4.163336342344337E-17, at i=8, j=0 finished k=6 computing local stiffness matrix for triangle 7 nodes i1=0, i2=8, i3=1 coord (0.5,0.5) (0.5,0.0) (0.0,0.0) cm matrix of am cm = ym solve for cm 0.0 0.0 -2.0 0.0 0.0 8.0 cm matrix of am cm = ym solve for cm 0.0 -2.0 2.0 8.0 -16.0 8.0 cm matrix of am cm = ym solve for cm 1.0 -6.0 0.0 8.0 0.0 -0.0 tri_int1a running with np=3, nv=9 i=0, xs=0.4880344338564595, ys=0.3938297308804235, fs=8.0 ph1=0.4531553936419275, ph2=-0.1174131974496464, ph3=-0.022785734101997468 gs1=16.0, gs2=32.0, gs3=16.0 i=1, xs=0.446914865440212, ys=0.39382973088042394, fs=8.0 ph1=0.45315539364192925, ph2=-0.08362601702972983, ph3=-0.0836260170297296 gs1=16.0, gs2=32.0, gs3=16.0 i=2, xs=0.405795297023964, ys=0.3938297308804235, fs=8.0 ph1=0.4531553936419275, ph2=-0.022785734101997246, ph3=-0.11741319744964662 gs1=16.0, gs2=32.0, gs3=16.0 i=3, xs=0.466722966080418, ys=0.20473343222036747, fs=8.0 ph1=-0.07414063829088052, ph2=0.025129059097551476, ph3=-0.057695179947284014 gs1=16.0, gs2=32.0, gs3=16.0 i=4, xs=0.3523667161101835, ys=0.204733432220367, fs=8.0 ph1=-0.07414063829088113, ph2=-0.12090187568290378, ph3=-0.12090187568290367 gs1=16.0, gs2=32.0, gs3=16.0 i=5, xs=0.23801046613994947, ys=0.204733432220367, fs=8.0 ph1=-0.07414063829088113, ph2=-0.057695179947284625, ph3=0.02512905909755142 gs1=16.0, gs2=32.0, gs3=16.0 i=6, xs=0.448641172595187, ys=0.04429397975635198, fs=8.0 ph1=-0.07289230637145501, ph2=0.4992788331754982, ph3=-0.08161582159044722 gs1=16.0, gs2=32.0, gs3=16.0 i=7, xs=0.272146989878176, ys=0.04429397975635199, fs=8.0 ph1=-0.07289230637145502, ph2=-0.040370066471039745, ph3=-0.04037006647103969 gs1=16.0, gs2=32.0, gs3=16.0 i=8, xs=0.09565280716116498, ys=0.044293979756351964, fs=8.0 ph1=-0.072892306371455, ph2=-0.08161582159044706, ph3=0.49927883317549826 gs1=16.0, gs2=32.0, gs3=16.0 i=0, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=16.0, B=-0.072892306371455, ws=0.013953605120761 integral=-2.0816681711721685E-17, at i=0, j=0 i=0, A=8.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=8.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=8.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=8.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=8.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=8.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=8.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=8.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=8.0, B=-0.072892306371455, ws=0.013953605120761 intg fg =-1.0408340855860843E-17, at i=0 i=0, A=32.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=32.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=32.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=32.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=32.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=32.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=32.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=32.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=32.0, B=-0.072892306371455, ws=0.013953605120761 integral=-4.163336342344337E-17, at i=8, j=0 i=0, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=1, A=16.0, B=0.45315539364192925, ws=0.00775855332238375 i=2, A=16.0, B=0.4531553936419275, ws=0.004849095826489875 i=3, A=16.0, B=-0.07414063829088052, ws=0.01591952127497125 i=4, A=16.0, B=-0.07414063829088113, ws=0.025471234039954 i=5, A=16.0, B=-0.07414063829088113, ws=0.01591952127497125 i=6, A=16.0, B=-0.07289230637145501, ws=0.013953605120761 i=7, A=16.0, B=-0.07289230637145502, ws=0.02232576819321775 i=8, A=16.0, B=-0.072892306371455, ws=0.013953605120761 integral=-2.0816681711721685E-17, at i=1, j=0 finished k=7 tri 0 has vertices 0, 1, 2 tri 1 has vertices 0, 2, 3 tri 2 has vertices 0, 3, 4 tri 3 has vertices 0, 4, 5 tri 4 has vertices 0, 5, 6 tri 5 has vertices 0, 6, 7 tri 6 has vertices 0, 7, 8 tri 7 has vertices 0, 8, 1 vertices, coordinates and analytic values coordinate 0 at 0.5, 0.5, uana=12.0 coordinate 1 at 0.0, 0.0, uana=6.0 coordinate 2 at 0.0, 0.5, uana=9.25 coordinate 3 at 0.0, 1.0, uana=14.0 coordinate 4 at 0.5, 1.0, uana=17.25 coordinate 5 at 1.0, 1.0, uana=21.0 coordinate 6 at 1.0, 0.5, uana=15.25 coordinate 7 at 1.0, 0.0, uana=11.0 coordinate 8 at 0.5, 0.0, uana=8.25 k computed stiffness matrix, if debug K(0,0)=-1.5265566588595902E-15 K(0,1)=-4.163336342344337E-17 K(0,2)=-8.326672684688674E-17 K(0,3)=-3.8163916471489756E-16 K(0,4)=-1.4432899320127035E-15 K(0,5)=-7.216449660063518E-16 K(0,6)=-1.4432899320127035E-15 K(0,7)=-3.8163916471489756E-16 K(0,8)=-8.326672684688674E-17 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(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(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(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(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)=1.0 K(5,6)=0.0 K(5,7)=0.0 K(5,8)=0.0 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)=1.0 K(6,7)=0.0 K(6,8)=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)=1.0 K(7,8)=0.0 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 f computed forcing function and boundary, if debug F(0)=-7.632783294297951E-16, ug[0]=0.0 F(1)=6.0, ug[1]=0.0 F(2)=9.25, ug[2]=0.0 F(3)=14.0, ug[3]=0.0 F(4)=17.25, ug[4]=0.0 F(5)=21.0, ug[5]=0.0 F(6)=15.25, ug[6]=0.0 F(7)=11.0, ug[7]=0.0 F(8)=8.25, ug[8]=0.0 redundant row (singular) 0 ug computed Galerkin, Ua analytic, error ug[0]=0.0, Ua=12.0, err=-12.0 ug[1]=6.0, Ua=6.0, err=0.0 ug[2]=9.25, Ua=9.25, err=0.0 ug[3]=14.0, Ua=14.0, err=0.0 ug[4]=17.25, Ua=17.25, err=0.0 ug[5]=21.0, Ua=21.0, err=0.0 ug[6]=15.25, Ua=15.25, err=0.0 ug[7]=11.0, Ua=11.0, err=0.0 ug[8]=8.25, Ua=8.25, err=0.0 maxerr=12.0, avgerr=1.3333333333333333