09 May 2025 9:47:37.496 PM wathen_matrix_test(): FORTRAN90 version test wathen_matrix(). test01(): Assemble, factor and solve a Wathen system defined by WATHEN_GE(). Elements in X direction NX = 4 Elements in Y direction NY = 4 Number of elements = 16 Number of nodes N = 65 Maximum solution error is 0.410783E-14 test02(): Assemble, factor and solve a Wathen system using WATHEN_GB(). Elements in X direction NX = 4 Elements in Y direction NY = 4 Number of elements = 16 Number of nodes N = 65 Lower bandwidth ML = 16 Upper bandwidth MU = 16 Maximum solution error is 0.310862E-14 test05(): For various problem sizes and storage schemes, measure the storage used for the Wathen system. Predicted Observed GE Band Band Band Sparse NX Elements Nodes storage width width storage storage 1 1 8 64 15 15 176 64 2 4 21 441 21 21 651 221 4 16 65 4225 33 33 3185 817 8 64 225 50625 57 57 19125 3137 16 256 833 693889 105 105 130781 12289 32 1024 3201 10246401 201 201 963501 48641 test06(): For various problem sizes, time the assembly and factorization of a Wathen system using the WATHEN_GE() function. NX Elements Nodes Storage Assembly Factor Error 1 1 8 64 0.20E-05 0.30E-05 0.13E-14 2 4 21 441 0.30E-05 0.90E-05 0.22E-14 4 16 65 4225 0.90E-05 0.16E-03 0.62E-14 8 64 225 50625 0.17E-03 0.52E-02 0.21E-13 16 256 833 693889 0.20E-02 0.22E+00 0.57E-14 32 1024 3201 10246401 0.27E-01 0.14E+02 0.70E-14 test07() For various problem sizes, time the assembly and factorization of a Wathen system using the WATHEN_GB() function. NX Elements Nodes Storage Assembly Factor Error 1 1 8 176 0.20E-05 0.20E-05 0.33E-15 2 4 21 651 0.30E-05 0.80E-05 0.18E-14 4 16 65 3185 0.70E-05 0.40E-04 0.30E-14 8 64 225 19125 0.35E-04 0.36E-03 0.43E-14 16 256 833 130781 0.23E-03 0.37E-02 0.40E-14 32 1024 3201 963501 0.31E-02 0.48E-01 0.10E-13 test08() For various problem sizes, time the assembly and factorization of a Wathen system WATHEN_GE() and WATHEN_GB() NX Elements Nodes Storage Assembly Factor Error WATHEN_GE 1 1 8 64 0.10E-05 0.20E-05 0.67E-15 WATHEN_GB 1 1 8 176 0.10E-05 0.20E-05 0.83E-15 WATHEN_GE 2 4 21 441 0.20E-05 0.80E-05 0.51E-14 WATHEN_GB 2 4 21 651 0.20E-05 0.70E-05 0.19E-14 WATHEN_GE 4 16 65 4225 0.70E-05 0.16E-03 0.26E-14 WATHEN_GB 4 16 65 3185 0.80E-05 0.39E-04 0.13E-14 WATHEN_GE 8 64 225 50625 0.55E-04 0.44E-02 0.98E-14 WATHEN_GB 8 64 225 19125 0.34E-04 0.34E-03 0.42E-14 WATHEN_GE 16 256 833 693889 0.88E-03 0.21E+00 0.14E-13 WATHEN_GB 16 256 833 130781 0.21E-03 0.36E-02 0.57E-14 WATHEN_GE 32 1024 3201 10246401 0.28E-01 0.14E+02 0.80E-14 WATHEN_GB 32 1024 3201 963501 0.20E-02 0.48E-01 0.10E-13 test10() Assemble, factor and solve a Wathen system defined by WATHEN_GE() and CG_GE(). Elements in X direction NX = 1 Elements in Y direction NY = 1 Number of elements = 1 Number of nodes N = 8 Maximum solution error is 0.149880E-14 test11(): Assemble, factor and solve a Wathen system defined by WATHEN_ST() and CG_ST(). Elements in X direction NX = 1 Elements in Y direction NY = 1 Number of elements = 1 Number of nodes N = 8 Number of nonzeros NZ_NUM = 64 Maximum solution error is 0.388578E-15 test115(): Assemble, factor and solve a Wathen system using WATHEN_GB() and CG_GB(). Elements in X direction NX = 4 Elements in Y direction NY = 4 Number of elements = 16 Number of nodes N = 65 Lower bandwidth ML = 16 Upper bandwidth MU = 16 Maximum solution error is 0.294472E-06 wathen_xy_test(): wathen_xy() returns the (X,Y) coordinates of nodes in the Wathen finite element system. k i j x y 1 1 1 0.0000 0.0000 2 2 1 0.1667 0.0000 3 3 1 0.3333 0.0000 4 4 1 0.5000 0.0000 5 5 1 0.6667 0.0000 6 6 1 0.8333 0.0000 7 7 1 1.0000 0.0000 8 1 2 0.0000 0.1667 9 2 2 0.3333 0.1667 10 3 2 0.6667 0.1667 11 4 2 1.0000 0.1667 12 1 3 0.0000 0.3333 13 2 3 0.1667 0.3333 14 3 3 0.3333 0.3333 15 4 3 0.5000 0.3333 16 5 3 0.6667 0.3333 17 6 3 0.8333 0.3333 18 7 3 1.0000 0.3333 19 1 4 0.0000 0.5000 20 2 4 0.3333 0.5000 21 3 4 0.6667 0.5000 22 4 4 1.0000 0.5000 23 1 5 0.0000 0.6667 24 2 5 0.1667 0.6667 25 3 5 0.3333 0.6667 26 4 5 0.5000 0.6667 27 5 5 0.6667 0.6667 28 6 5 0.8333 0.6667 29 7 5 1.0000 0.6667 30 1 6 0.0000 0.8333 31 2 6 0.3333 0.8333 32 3 6 0.6667 0.8333 33 4 6 1.0000 0.8333 34 1 7 0.0000 1.0000 35 2 7 0.1667 1.0000 36 3 7 0.3333 1.0000 37 4 7 0.5000 1.0000 38 5 7 0.6667 1.0000 39 6 7 0.8333 1.0000 40 7 7 1.0000 1.0000 wathen_matrix_test(): Normal end of execution. 09 May 2025 9:48:06.261 PM