January 26 2013 2:06:07.861 PM 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 "baffle_nodes.txt". Element file is "baffle_elements.txt". Number of nodes = 512 First 10 nodes Row 1 2 Col 1 0.00000 1.00000 2 0.00000 1.37500 3 0.00000 1.75000 4 0.00000 2.12500 5 0.00000 2.50000 6 0.00000 2.87500 7 0.00000 3.25000 8 0.00000 3.62500 9 0.00000 4.00000 10 0.00000 4.37500 Element order = 3 Number of elements = 874 First 10 elements Row 1 2 3 Col 1 317 306 298 2 28 3 18 3 4 3 28 4 28 33 4 5 72 69 57 6 49 33 28 7 114 106 122 8 84 80 67 9 88 69 72 10 147 133 131 Quadrature order = 3 Number of nonzero coefficients NZ_NUM = 3308 Diagonal adjacency vector: 1: 1 2: 6 3: 10 4: 15 5: 21 6: 25 7: 31 8: 35 9: 41 10: 46 11: 51 12: 56 13: 60 14: 65 15: 70 16: 76 17: 80 18: 86 19: 93 20: 100 21: 107 22: 114 23: 120 24: 126 25: 131 26: 136 27: 145 28: 152 29: 159 30: 167 31: 174 32: 180 33: 187 34: 194 35: 201 36: 208 37: 215 38: 222 39: 229 40: 236 41: 243 42: 248 43: 253 44: 261 45: 268 46: 273 47: 279 48: 286 49: 294 50: 302 51: 309 52: 316 53: 323 54: 330 55: 338 56: 343 57: 350 58: 355 59: 360 60: 367 61: 375 62: 382 63: 390 64: 396 65: 402 66: 409 67: 416 68: 423 69: 430 70: 437 71: 445 72: 452 73: 457 74: 464 75: 471 76: 478 77: 485 78: 491 79: 497 80: 503 81: 510 82: 516 83: 521 84: 527 85: 533 86: 538 87: 543 88: 549 89: 556 90: 562 91: 566 92: 572 93: 580 94: 588 95: 596 96: 602 97: 608 98: 614 99: 620 100: 626 101: 632 102: 637 103: 642 104: 648 105: 656 106: 662 107: 670 108: 676 109: 682 110: 688 111: 695 112: 701 113: 709 114: 716 115: 723 116: 730 117: 737 118: 744 119: 751 120: 757 121: 763 122: 769 123: 777 124: 783 125: 789 126: 796 127: 803 128: 810 129: 817 130: 824 131: 829 132: 833 133: 840 134: 847 135: 854 136: 862 137: 869 138: 876 139: 883 140: 890 141: 896 142: 903 143: 909 144: 917 145: 923 146: 932 147: 938 148: 944 149: 950 150: 957 151: 963 152: 970 153: 977 154: 984 155: 992 156: 999 157: 1005 158: 1011 159: 1017 160: 1025 161: 1032 162: 1039 163: 1046 164: 1053 165: 1059 166: 1063 167: 1069 168: 1075 169: 1081 170: 1087 171: 1093 172: 1099 173: 1105 174: 1112 175: 1120 176: 1129 177: 1136 178: 1143 179: 1150 180: 1155 181: 1161 182: 1167 183: 1172 184: 1177 185: 1183 186: 1189 187: 1196 188: 1203 189: 1209 190: 1216 191: 1222 192: 1228 193: 1235 194: 1240 195: 1246 196: 1253 197: 1261 198: 1269 199: 1277 200: 1284 201: 1291 202: 1297 203: 1304 204: 1311 205: 1317 206: 1324 207: 1331 208: 1337 209: 1344 210: 1350 211: 1355 212: 1362 213: 1370 214: 1376 215: 1383 216: 1390 217: 1397 218: 1404 219: 1411 220: 1418 221: 1425 222: 1432 223: 1439 224: 1445 225: 1453 226: 1459 227: 1466 228: 1471 229: 1476 230: 1484 231: 1491 232: 1498 233: 1505 234: 1512 235: 1519 236: 1526 237: 1533 238: 1541 239: 1547 240: 1553 241: 1558 242: 1565 243: 1573 244: 1578 245: 1584 246: 1590 247: 1596 248: 1602 249: 1608 250: 1614 251: 1620 252: 1626 253: 1632 254: 1640 255: 1647 256: 1650 257: 1655 258: 1663 259: 1670 260: 1677 261: 1683 262: 1689 263: 1695 264: 1701 265: 1707 266: 1713 267: 1719 268: 1725 269: 1731 270: 1737 271: 1743 272: 1748 273: 1754 274: 1760 275: 1768 276: 1776 277: 1783 278: 1790 279: 1797 280: 1804 281: 1811 282: 1818 283: 1825 284: 1833 285: 1840 286: 1845 287: 1850 288: 1857 289: 1863 290: 1870 291: 1877 292: 1884 293: 1891 294: 1898 295: 1905 296: 1912 297: 1919 298: 1926 299: 1933 300: 1939 301: 1945 302: 1951 303: 1958 304: 1964 305: 1971 306: 1978 307: 1985 308: 1992 309: 1999 310: 2006 311: 2013 312: 2020 313: 2028 314: 2036 315: 2043 316: 2051 317: 2059 318: 2064 319: 2069 320: 2076 321: 2082 322: 2088 323: 2094 324: 2101 325: 2107 326: 2113 327: 2118 328: 2124 329: 2130 330: 2135 331: 2140 332: 2147 333: 2154 334: 2160 335: 2164 336: 2171 337: 2178 338: 2186 339: 2194 340: 2201 341: 2209 342: 2215 343: 2221 344: 2227 345: 2233 346: 2239 347: 2246 348: 2254 349: 2260 350: 2265 351: 2271 352: 2278 353: 2286 354: 2293 355: 2299 356: 2305 357: 2311 358: 2317 359: 2325 360: 2331 361: 2338 362: 2345 363: 2352 364: 2359 365: 2365 366: 2370 367: 2376 368: 2384 369: 2390 370: 2398 371: 2406 372: 2412 373: 2419 374: 2427 375: 2434 376: 2441 377: 2447 378: 2455 379: 2462 380: 2470 381: 2476 382: 2480 383: 2485 384: 2492 385: 2499 386: 2506 387: 2513 388: 2519 389: 2525 390: 2532 391: 2539 392: 2545 393: 2553 394: 2559 395: 2565 396: 2572 397: 2579 398: 2586 399: 2593 400: 2600 401: 2608 402: 2614 403: 2621 404: 2626 405: 2633 406: 2641 407: 2648 408: 2654 409: 2661 410: 2668 411: 2672 412: 2677 413: 2683 414: 2689 415: 2695 416: 2701 417: 2707 418: 2713 419: 2720 420: 2729 421: 2736 422: 2742 423: 2747 424: 2753 425: 2760 426: 2766 427: 2772 428: 2777 429: 2783 430: 2789 431: 2794 432: 2800 433: 2806 434: 2812 435: 2818 436: 2826 437: 2831 438: 2839 439: 2844 440: 2849 441: 2856 442: 2864 443: 2871 444: 2878 445: 2885 446: 2892 447: 2899 448: 2906 449: 2913 450: 2920 451: 2927 452: 2934 453: 2940 454: 2949 455: 2956 456: 2962 457: 2967 458: 2973 459: 2979 460: 2986 461: 2993 462: 3000 463: 3008 464: 3015 465: 3023 466: 3029 467: 3036 468: 3045 469: 3051 470: 3057 471: 3062 472: 3067 473: 3073 474: 3081 475: 3087 476: 3094 477: 3101 478: 3108 479: 3114 480: 3121 481: 3128 482: 3135 483: 3141 484: 3150 485: 3155 486: 3164 487: 3171 488: 3177 489: 3182 490: 3188 491: 3196 492: 3203 493: 3209 494: 3216 495: 3223 496: 3229 497: 3233 498: 3238 499: 3244 500: 3248 501: 3254 502: 3258 503: 3264 504: 3269 505: 3274 506: 3279 507: 3283 508: 3288 509: 3293 510: 3298 511: 3303 512: 3308 Part of Finite Element matrix A: Col: 1 2 3 4 5 Row --- 1 0.925899 -0.405114 0.00000 0.00000 0.00000 2 -0.405114 2.06613 -0.370937 0.00000 0.00000 3 0.00000 -0.370937 1.86995 -0.528324 0.00000 4 0.00000 0.00000 -0.528324 1.84140 -0.412654 5 0.00000 0.00000 0.00000 -0.412654 2.05364 6 0.00000 0.00000 0.00000 0.00000 -0.382925 Col: 6 7 8 9 10 Row --- 5 -0.382925 0.00000 0.00000 0.00000 0.00000 6 1.75240 -0.389630 0.00000 0.00000 0.00000 7 -0.389630 2.04239 -0.428216 0.00000 0.00000 8 0.00000 -0.428216 1.74296 -0.343709 0.00000 9 0.00000 0.00000 -0.343709 1.90310 -0.474999 10 0.00000 0.00000 0.00000 -0.474999 1.88456 Part of right hand side vector F: 1 0.81684830E-01 2 0.71605730E-01 3 0.14574345 4 0.21434168 5 0.74742670E-01 6 0.17805744 7 0.65253598E-01 8 0.14463460 9 0.85275214E-01 10 0.76656028E-01 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.00000 0.00000 0.00000 1.00000 Col: 6 7 8 9 10 Row --- 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.00000 1.00000 0.00000 0.00000 9 0.00000 0.00000 0.00000 1.00000 0.00000 10 0.00000 0.00000 0.00000 0.00000 1.00000 Part of F after adjustment for Dirichlet condition: 1 0.0000000 2 0.0000000 3 0.0000000 4 0.0000000 5 0.0000000 6 0.0000000 7 0.0000000 8 0.0000000 9 0.0000000 10 0.0000000 Step Residual 1 1.94589 2 1.02565 3 0.411712 4 0.271570 5 0.143458 6 0.102652 7 0.490419E-01 8 0.317833E-01 9 0.164217E-01 10 0.860515E-02 11 0.392365E-02 12 0.229199E-02 13 0.156186E-02 14 0.742249E-03 15 0.399329E-03 16 0.205418E-03 17 0.124346E-03 18 0.676339E-04 19 0.322410E-04 20 0.196990E-04 21 0.134579E-04 22 0.703567E-05 23 0.352624E-05 24 0.196816E-05 25 0.107468E-05 26 0.719786E-06 27 0.424296E-06 28 0.244519E-06 29 0.135357E-06 30 0.810820E-07 31 0.470307E-07 Number of iterations was 31 Estimated error is 0.289709E-07 Part of the solution vector U: 1 0.0000000 2 0.0000000 3 0.0000000 4 0.0000000 5 0.0000000 6 0.0000000 7 0.0000000 8 0.0000000 9 0.0000000 10 0.0000000 FEM2D_POISSON_CG: Wrote an ASCII file "baffle_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 2:06:07.899 PM