Fri Sep 14 08:47:36 2018 WATHEN_TEST Python version: 3.6.5 Test the WATHEN library. WATHEN_TEST01 Python version: 3.6.5 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 Norm of solution error is 6.5081e-15 WATHEN_TEST01: Normal end of execution. WATHEN_TEST02 Python version: 3.6.5 Assemble, factor and solve a Wathen system defined by WATHEN_CSC. Elements in X direction NX = 4 Elements in Y direction NY = 4 Number of elements = 16 Number of nodes N = 65 Norm of solution error is 6.34306e-15 WATHEN_TEST02: Normal end of execution. WATHEN_TEST03 Python version: 3.6.5 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 3.46e-04 1.53e-03 9.86e-16 2 4 21 441 1.09e-03 1.20e-03 4.71e-15 4 16 65 4225 4.09e-03 9.44e-04 6.51e-15 8 64 225 50625 1.76e-02 2.76e-03 7.20e-15 16 256 833 693889 8.07e-02 4.67e-02 5.66e-13 32 1024 3201 10246401 3.29e-01 1.50e+00 4.38e-14 WATHEN_TEST03: Normal end of execution. WATHEN_TEST04 Python version: 3.6.5 For various problem sizes, time the assembly and factorization of a Wathen system using the WATHEN_CSC function. NX Elements Nodes Assembly Factor Error 1 1 8 1.23e-03 3.59e-04 1.33e-15 2 4 21 1.54e-03 4.04e-04 2.52e-15 4 16 65 4.01e-03 8.55e-04 6.34e-15 8 64 225 1.35e-02 2.73e-03 6.87e-15 16 256 833 5.31e-02 1.18e-02 4.27e-13 32 1024 3201 2.03e-01 8.58e-02 3.27e-14 64 4096 12545 3.70e-01 3.01e-01 6.98e-14 WATHEN_TEST04: Normal end of execution. WATHEN_TEST05 Python version: 3.6.5 For various problem sizes, time the assembly and factorization of a Wathen system WATHEN_GE and WATHEN_CSC. NX Elements Nodes Assembly Factor Error WATHEN_GE 1 1 8 3.67e-04 5.54e-04 9.86e-16 WATHEN_CSC 1 1 8 1.11e-03 3.01e-04 9.86e-16 WATHEN_GE 2 4 21 1.09e-03 4.69e-04 4.71e-15 WATHEN_CSC 2 4 21 1.57e-03 3.60e-04 3.41e-15 WATHEN_GE 4 16 65 4.06e-03 8.17e-04 6.51e-15 WATHEN_CSC 4 16 65 4.20e-03 7.36e-04 1.12e-14 WATHEN_GE 8 64 225 1.60e-02 1.12e-03 7.20e-15 WATHEN_CSC 8 64 225 1.37e-02 2.49e-03 6.52e-15 WATHEN_GE 16 256 833 6.89e-02 3.49e-02 5.66e-13 WATHEN_CSC 16 256 833 5.48e-02 1.11e-02 1.52e-13 WATHEN_GE 32 1024 3201 3.27e-01 1.48e+00 4.38e-14 WATHEN_CSC 32 1024 3201 2.11e-01 8.28e-02 3.26e-14 WATHEN_TEST05: Normal end of execution. WATHEN_TEST06 Python version: 3.6.5 Assemble, factor and solve a Wathen system defined by WATHEN_GE and CG_GE. Elements in X direction NX = 2 Elements in Y direction NY = 2 Number of elements = 4 Number of nodes N = 21 Maximum solution error is 2.33828e-05 WATHEN_TEST06: Normal end of execution. WATHEN_TEST07 Python version: 3.6.5 Assemble, factor and solve a Wathen system defined by WATHEN_CSC and CG_CSC. Elements in X direction NX = 2 Elements in Y direction NY = 2 Number of elements = 4 Number of nodes N = 21 Maximum solution error is 1.11836e-05 WATHEN_TEST07: Normal end of execution. WATHEN_TEST08 Python version: 3.6.5 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 = 64 Maximum solution error is 1.66337e-15 WATHEN_TEST08: Normal end of execution. WATHEN_TEST09 Python version: 3.6.5 Display the sparsity of the Wathen matrix. Graphics file saved as "wathen_spy.png" WATHEN_TEST09: Normal end of execution. WATHEN_ORDER_TEST: Python version: 3.6.5 WATHEN_ORDER returns N, the order of a Wathen finite element matrix given NX and NY, the number of rows and columns of nodes in the underlying grid. NX / NY: 1 2 3 4 5 6 1 8 13 18 23 28 33 2 13 21 29 37 45 53 3 18 29 40 51 62 73 4 23 37 51 65 79 93 5 28 45 62 79 96 113 6 33 53 73 93 113 133 7 38 61 84 107 130 153 8 43 69 95 121 147 173 9 48 77 106 135 164 193 10 53 85 117 149 181 213 WATHEN_ORDER_TEST: Normal end of execution. WATHEN_ST_SIZE_TEST_TEST: Python version: 3.6.5 WATHEN_ST_SIZE returns NZ_NUM, the number of nonzeros in a sparse triplet format for a Wathen finite element matrix, given NX and NY, the number of rows and columns of nodes in the underlying grid. NX / NY: 1 2 3 4 5 6 1 64 128 192 256 320 384 2 128 256 384 512 640 768 3 192 384 576 768 960 1152 4 256 512 768 1024 1280 1536 5 320 640 960 1280 1600 1920 6 384 768 1152 1536 1920 2304 7 448 896 1344 1792 2240 2688 8 512 1024 1536 2048 2560 3072 9 576 1152 1728 2304 2880 3456 10 640 1280 1920 2560 3200 3840 WATHEN_ST_SIZE_TEST: Normal end of execution. WATHEN_TEST: Normal end of execution. Fri Sep 14 08:47:42 2018