10 September 2021 11:35:49.935 AM WATHEN_TEST FORTRAN90 version Test the WATHEN library. 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.111022E-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.832841E-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.30E-05 0.20E-05 0.56E-15 2 4 21 441 0.10E-05 0.10E-04 0.44E-14 4 16 65 4225 0.80E-05 0.20E-03 0.18E-14 8 64 225 50625 0.15E-03 0.72E-02 0.68E-14 16 256 833 693889 0.16E-02 0.32E+00 0.23E-14 32 1024 3201 10246401 0.24E-01 0.19E+02 0.52E-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.20E-05 0.70E-05 0.52E-14 4 16 65 3185 0.11E-04 0.47E-04 0.14E-13 8 64 225 19125 0.42E-04 0.35E-03 0.39E-14 16 256 833 130781 0.34E-03 0.35E-02 0.21E-13 32 1024 3201 963501 0.32E-02 0.50E-01 0.73E-14 TEST08 For various problem sizes, time the assembly and factorization of a Wathen system WATHEN_GE/WATHEN_GB NX Elements Nodes Storage Assembly Factor Error WATHEN_GE 1 1 8 64 0.10E-05 0.30E-05 0.21E-14 WATHEN_GB 1 1 8 176 0.10E-05 0.20E-05 0.12E-14 WATHEN_GE 2 4 21 441 0.20E-05 0.90E-05 0.34E-14 WATHEN_GB 2 4 21 651 0.20E-05 0.70E-05 0.18E-14 WATHEN_GE 4 16 65 4225 0.70E-05 0.21E-03 0.31E-14 WATHEN_GB 4 16 65 3185 0.90E-05 0.42E-04 0.13E-14 WATHEN_GE 8 64 225 50625 0.64E-04 0.64E-02 0.39E-14 WATHEN_GB 8 64 225 19125 0.39E-04 0.37E-03 0.72E-14 WATHEN_GE 16 256 833 693889 0.91E-03 0.32E+00 0.56E-14 WATHEN_GB 16 256 833 130781 0.26E-03 0.36E-02 0.92E-14 WATHEN_GE 32 1024 3201 10246401 0.24E-01 0.20E+02 0.70E-14 WATHEN_GB 32 1024 3201 963501 0.21E-02 0.54E-01 0.22E-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.999201E-15 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.666134E-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.175069E-07 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_TEST Normal end of execution. 10 September 2021 11:36:30.496 AM