RESULTS FOR TUTORIAL 8 HERE FOR f(x) = 1 NO GUARANTEE FOR THE CORRECTNESS (ERRORS CAN ALSO HAPPEN TO ME) ---------------- FIRST MESH in meshdemo.cc ------------------------------------ pattern = [** ** ] [*** ** ] [ ** * ] [* ** ** ] [** *** **] [ ** ** *] [ * ** ] [ ** ***] [ ** **] stiffness matrix = [ 1.0416666666667 -0.48958333333333 0 -0.48958333333333 0.020833333333333 0 0 0 0 ] [ -0.48958333333333 2.0625 -0.48958333333333 0 -0.97916666666667 0.020833333333333 0 0 0 ] [ 0 -0.48958333333333 1.0208333333333 0 0 -0.48958333333333 0 0 0 ] [ -0.48958333333333 0 0 2.0625 -0.97916666666667 0 -0.48958333333333 0.020833333333333 0 ] [ 0.020833333333333 -0.97916666666667 0 -0.97916666666667 4.125 -0.97916666666667 0 -0.97916666666667 0.020833333333333 ] [ 0 0.020833333333333 -0.48958333333333 0 -0.97916666666667 2.0625 0 0 -0.48958333333333 ] [ 0 0 0 -0.48958333333333 0 0 1.0208333333333 -0.48958333333333 0 ] [ 0 0 0 0.020833333333333 -0.97916666666667 0 -0.48958333333333 2.0625 -0.48958333333333 ] [ 0 0 0 0 0.020833333333333 -0.48958333333333 0 -0.48958333333333 1.0416666666667 ] load vector = [ 0.083333333333333 0.125 0.041666666666667 0.125 0.25 0.125 0.041666666666667 0.125 0.083333333333333 ] ---------------- SECOND MESH in meshdemo.cc (the first one refined) ----------- pattern = [* *** ] [ * * *** ] [ * * * ] [ * * *** ] [ * * * * *** ] [ * ** * * ] [ * * * ] [ * * * **] [ * ** *] [** * * * ] [* * ** * ] [* * *** * * ] [ ** * ** ] [ * * * * ** * ] [ * * **** * ] [ * * * ** ] [ ** ** * * * ] [ * * ** * ] [ * * *** * * ] [ ** ** * ** ] [ * * * * ** *] [ * * **** *] [ * * * ** ] [ ** ** * ] [ ** ** *] stiffness matrix = [ 1.0104166666667 0 0 0 0 0 0 0 0 -0.49739583333333 -0.49739583333333 0.0052083333333333 0 0 0 0 0 0 0 0 0 0 0 0 0 ] [ 0 2.015625 0 0 0 0 0 0 0 -0.49739583333333 0 0 -0.49739583333333 -0.99479166666667 0.0052083333333333 0 0 0 0 0 0 0 0 0 0 ] [ 0 0 1.0052083333333 0 0 0 0 0 0 0 0 0 -0.49739583333333 0 0 -0.49739583333333 0 0 0 0 0 0 0 0 0 ] [ 0 0 0 2.015625 0 0 0 0 0 0 -0.49739583333333 0 0 0 0 0 -0.99479166666667 -0.49739583333333 0.0052083333333333 0 0 0 0 0 0 ] [ 0 0 0 0 4.03125 0 0 0 0 0 0 0.0052083333333333 0 -0.99479166666667 0 0 -0.99479166666667 0 0 -0.99479166666667 -0.99479166666667 0.0052083333333333 0 0 0 ] [ 0 0 0 0 0 2.015625 0 0 0 0 0 0 0 0 0.0052083333333333 -0.49739583333333 0 0 0 -0.99479166666667 0 0 -0.49739583333333 0 0 ] [ 0 0 0 0 0 0 1.0052083333333 0 0 0 0 0 0 0 0 0 0 -0.49739583333333 0 0 0 0 0 -0.49739583333333 0 ] [ 0 0 0 0 0 0 0 2.015625 0 0 0 0 0 0 0 0 0 0 0.0052083333333333 0 -0.99479166666667 0 0 -0.49739583333333 -0.49739583333333 ] [ 0 0 0 0 0 0 0 0 1.0104166666667 0 0 0 0 0 0 0 0 0 0 0 0 0.0052083333333333 -0.49739583333333 0 -0.49739583333333 ] [ -0.49739583333333 -0.49739583333333 0 0 0 0 0 0 0 2.015625 0 -0.99479166666667 0 0.0052083333333333 0 0 0 0 0 0 0 0 0 0 0 ] [ -0.49739583333333 0 0 -0.49739583333333 0 0 0 0 0 0 2.015625 -0.99479166666667 0 0 0 0 0.0052083333333333 0 0 0 0 0 0 0 0 ] [ 0.0052083333333333 0 0 0 0.0052083333333333 0 0 0 0 -0.99479166666667 -0.99479166666667 4.03125 0 -0.99479166666667 0 0 -0.99479166666667 0 0 0 0 0 0 0 0 ] [ 0 -0.49739583333333 -0.49739583333333 0 0 0 0 0 0 0 0 0 2.015625 0 -0.99479166666667 0.0052083333333333 0 0 0 0 0 0 0 0 0 ] [ 0 -0.99479166666667 0 0 -0.99479166666667 0 0 0 0 0.0052083333333333 0 -0.99479166666667 0 4.03125 -0.99479166666667 0 0 0 0 0.0052083333333333 0 0 0 0 0 ] [ 0 0.0052083333333333 0 0 0 0.0052083333333333 0 0 0 0 0 0 -0.99479166666667 -0.99479166666667 4.03125 -0.99479166666667 0 0 0 -0.99479166666667 0 0 0 0 0 ] [ 0 0 -0.49739583333333 0 0 -0.49739583333333 0 0 0 0 0 0 0.0052083333333333 0 -0.99479166666667 2.015625 0 0 0 0 0 0 0 0 0 ] [ 0 0 0 -0.99479166666667 -0.99479166666667 0 0 0 0 0 0.0052083333333333 -0.99479166666667 0 0 0 0 4.03125 0 -0.99479166666667 0 0.0052083333333333 0 0 0 0 ] [ 0 0 0 -0.49739583333333 0 0 -0.49739583333333 0 0 0 0 0 0 0 0 0 0 2.015625 -0.99479166666667 0 0 0 0 0.0052083333333333 0 ] [ 0 0 0 0.0052083333333333 0 0 0 0.0052083333333333 0 0 0 0 0 0 0 0 -0.99479166666667 -0.99479166666667 4.03125 0 -0.99479166666667 0 0 -0.99479166666667 0 ] [ 0 0 0 0 -0.99479166666667 -0.99479166666667 0 0 0 0 0 0 0 0.0052083333333333 -0.99479166666667 0 0 0 0 4.03125 0 -0.99479166666667 0.0052083333333333 0 0 ] [ 0 0 0 0 -0.99479166666667 0 0 -0.99479166666667 0 0 0 0 0 0 0 0 0.0052083333333333 0 -0.99479166666667 0 4.03125 -0.99479166666667 0 0 0.0052083333333333 ] [ 0 0 0 0 0.0052083333333333 0 0 0 0.0052083333333333 0 0 0 0 0 0 0 0 0 0 -0.99479166666667 -0.99479166666667 4.03125 -0.99479166666667 0 -0.99479166666667 ] [ 0 0 0 0 0 -0.49739583333333 0 0 -0.49739583333333 0 0 0 0 0 0 0 0 0 0 0.0052083333333333 0 -0.99479166666667 2.015625 0 0 ] [ 0 0 0 0 0 0 -0.49739583333333 -0.49739583333333 0 0 0 0 0 0 0 0 0 0.0052083333333333 -0.99479166666667 0 0 0 0 2.015625 0 ] [ 0 0 0 0 0 0 0 -0.49739583333333 -0.49739583333333 0 0 0 0 0 0 0 0 0 0 0 0.0052083333333333 -0.99479166666667 0 0 2.015625 ] load vector = [ 0.020833333333333 0.03125 0.010416666666667 0.03125 0.0625 0.03125 0.010416666666667 0.03125 0.020833333333333 0.03125 0.03125 0.0625 0.03125 0.0625 0.0625 0.03125 0.0625 0.03125 0.0625 0.0625 0.0625 0.0625 0.03125 0.03125 0.03125 ]