January 26 2013 11:23:02.353 AM FEM2D_POISSON_CG FORTRAN90 version: A version of FEM2D_POISSON using sparse storage and a conjugate gradient solver. Solution of the Poisson equation in an arbitrary region in 2 dimensions. - DEL H(x,y) DEL U(x,y) + K(x,y) * U(x,y) = F(x,y) in the region U(x,y) = G(x,y) on the boundary. The finite element method is used, with triangular elements, which must be a 3 node linear triangle. Node file is "ell_nodes.txt". Element file is "ell_elements.txt". Number of nodes = 65 First 10 nodes Row 1 2 Col 1 0.00000 0.00000 2 0.00000 0.500000 3 0.500000 0.00000 4 0.00000 1.00000 5 0.500000 0.500000 6 1.00000 0.00000 7 0.00000 1.50000 8 0.500000 1.00000 9 1.00000 0.500000 10 1.50000 0.00000 Element order = 3 Number of elements = 96 First 10 elements Row 1 2 3 Col 1 1 3 2 2 6 5 3 3 4 2 5 4 3 5 2 5 23 22 10 6 21 9 22 7 6 10 9 8 22 9 10 9 19 7 20 10 4 8 7 Quadrature order = 3 Number of nonzero coefficients NZ_NUM = 385 Diagonal adjacency vector: 1: 1 2: 5 3: 11 4: 15 5: 22 6: 28 7: 32 8: 39 9: 46 10: 52 11: 55 12: 59 13: 63 14: 68 15: 74 16: 78 17: 83 18: 88 19: 94 20: 100 21: 107 22: 114 23: 120 24: 126 25: 132 26: 138 27: 146 28: 152 29: 159 30: 167 31: 173 32: 180 33: 187 34: 193 35: 199 36: 205 37: 212 38: 218 39: 226 40: 232 41: 239 42: 246 43: 252 44: 258 45: 265 46: 271 47: 275 48: 282 49: 289 50: 297 51: 301 52: 308 53: 314 54: 321 55: 329 56: 334 57: 339 58: 345 59: 353 60: 358 61: 363 62: 370 63: 375 64: 381 65: 385 Part of Finite Element matrix A: Col: 1 2 3 4 5 Row --- 1 1.02083 -0.489583 -0.489583 0.00000 0.00000 2 -0.489583 2.06250 0.208333E-01 -0.489583 -0.979167 3 -0.489583 0.208333E-01 2.06250 0.00000 -0.979167 4 0.00000 -0.489583 0.00000 2.06250 0.208333E-01 5 0.00000 -0.979167 -0.979167 0.208333E-01 4.12500 6 0.00000 0.00000 -0.489583 0.00000 0.208333E-01 7 0.00000 0.00000 0.00000 -0.489583 0.00000 8 0.00000 0.00000 0.00000 -0.979167 -0.979167 9 0.00000 0.00000 0.00000 0.00000 -0.979167 Col: 6 7 8 9 10 Row --- 3 -0.489583 0.00000 0.00000 0.00000 0.00000 4 0.00000 -0.489583 -0.979167 0.00000 0.00000 5 0.208333E-01 0.00000 -0.979167 -0.979167 0.00000 6 2.06250 0.00000 0.00000 -0.979167 -0.489583 7 0.00000 2.06250 0.208333E-01 0.00000 0.00000 8 0.00000 0.208333E-01 4.12500 0.208333E-01 0.00000 9 -0.979167 0.00000 0.208333E-01 4.12500 0.208333E-01 10 -0.489583 0.00000 0.00000 0.208333E-01 2.06250 Part of right hand side vector F: 1 -.16406250 2 -.46875000 3 -.46875000 4 -.38541667 5 -.85416667 6 -.38541667 7 -.23958333 8 -.66666667 9 -.66666667 10 -.23958333 Part of A after adjustment for Dirichlet condition: Col: 1 2 3 4 5 Row --- 1 1.00000 0.00000 0.00000 0.00000 0.00000 2 0.00000 1.00000 0.00000 0.00000 0.00000 3 0.00000 0.00000 1.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 1.00000 0.00000 5 0.00000 -0.979167 -0.979167 0.208333E-01 4.12500 8 0.00000 0.00000 0.00000 -0.979167 -0.979167 9 0.00000 0.00000 0.00000 0.00000 -0.979167 Col: 6 7 8 9 10 Row --- 5 0.208333E-01 0.00000 -0.979167 -0.979167 0.00000 6 1.00000 0.00000 0.00000 0.00000 0.00000 7 0.00000 1.00000 0.00000 0.00000 0.00000 8 0.00000 0.208333E-01 4.12500 0.208333E-01 0.00000 9 -0.979167 0.00000 0.208333E-01 4.12500 0.208333E-01 10 0.00000 0.00000 0.00000 0.00000 1.00000 Part of F after adjustment for Dirichlet condition: 1 0.0000000 2 0.25000000 3 0.25000000 4 1.0000000 5 -.85416667 6 1.0000000 7 2.2500000 8 -.66666667 9 -.66666667 10 2.2500000 Step Residual 1 39.3312 2 23.1596 3 10.7601 4 6.12728 5 3.47712 6 1.65520 7 0.592650 8 0.162597 9 0.693399E-01 10 0.241027E-01 11 0.609233E-02 12 0.162877E-02 13 0.299049E-03 14 0.670415E-04 15 0.438113E-05 Number of iterations was 15 Estimated error is 0.744849E-07 Part of the solution vector U: 1 0.0000000 2 0.25000000 3 0.25000000 4 1.0000000 5 0.48493558 6 1.0000000 7 2.2500000 8 1.2289069 9 1.2289069 10 2.2500000 FEM2D_POISSON_CG: Wrote an ASCII file "ell_values.txt" of the form U ( X(I), Y(I) ) which can be used for plotting. FEM2D_POISSON_CG: Normal end of execution. January 26 2013 11:23:02.401 AM