6 March 2014 9:39:11.739 AM HAAR_PRB FORTRAN77 version Test the HAAR library. TEST01 HAAR_1D computes the Haar transform of a vector. i U(i) H(U)(i) Hinv(H(U))(i) 1 0.2184 1.6398 0.2184 2 0.9563 0.0677 0.9563 3 0.8295 0.6070 0.8295 4 0.5617 -0.2708 0.5617 5 0.4153 -0.1082 0.4153 6 0.0661 0.0569 0.0661 7 0.2576 0.0833 0.2576 8 0.1100 0.1784 0.1100 9 0.0438 -0.5218 0.0438 10 0.6340 0.1894 0.6340 11 0.0617 0.2469 0.0617 12 0.4495 0.1044 0.4495 13 0.4013 -0.4173 0.4013 14 0.7547 -0.2742 0.7547 15 0.7973 -0.2499 0.7973 16 0.0018 0.5625 0.0018 i U(i) H(U)(i) Hinv(H(U))(i) 1 1.0000 2.8284 1.0000 2 1.0000 0.0000 1.0000 3 1.0000 0.0000 1.0000 4 1.0000 0.0000 1.0000 5 1.0000 0.0000 1.0000 6 1.0000 0.0000 1.0000 7 1.0000 0.0000 1.0000 8 1.0000 0.0000 1.0000 i U(i) H(U)(i) Hinv(H(U))(i) 1 1.0000 34.0000 1.0000 2 2.0000 -16.0000 2.0000 3 3.0000 -5.6569 3.0000 4 4.0000 -5.6569 4.0000 5 5.0000 -2.0000 5.0000 6 6.0000 -2.0000 6.0000 7 7.0000 -2.0000 7.0000 8 8.0000 -2.0000 8.0000 9 9.0000 -0.7071 9.0000 10 10.0000 -0.7071 10.0000 11 11.0000 -0.7071 11.0000 12 12.0000 -0.7071 12.0000 13 13.0000 -0.7071 13.0000 14 14.0000 -0.7071 14.0000 15 15.0000 -0.7071 15.0000 16 16.0000 -0.7071 16.0000 i U(i) H(U)(i) Hinv(H(U))(i) 1 25.0000 21.2132 25.0000 2 16.0000 16.9706 16.0000 3 9.0000 14.0000 9.0000 4 4.0000 -2.0000 4.0000 5 1.0000 6.3640 1.0000 6 0.0000 3.5355 0.0000 7 1.0000 0.7071 1.0000 8 4.0000 -2.1213 4.0000 For N = 99, ||u-haar_1d_inverse(haar_1d(u))|| = 0.333568E-14 TEST02 HAAR_2D computes the Haar transform of an array. HAAR_2D_INVERSE inverts the transform. Input array U: Col 1 2 3 4 Row 1: 0.218418 0.897504 0.861216 0.825003 2: 0.956318 0.350752 0.453794 0.824660 3: 0.829509 0.945448E-01 0.911977 0.618618E-01 4: 0.561695 0.136169E-01 0.597917 0.710781 5: 0.415307 0.859097 0.188955 0.882833E-01 6: 0.661187E-01 0.840847 0.761492 0.777994 7: 0.257578 0.123104 0.396988 0.745303 8: 0.109957 0.751236E-02 0.185314 0.308675 9: 0.438290E-01 0.260303 0.574366 0.899373 10: 0.633966 0.912484 0.367027 0.763537 11: 0.617272E-01 0.113664 0.617205 0.761731 12: 0.449539 0.351629 0.361529 0.406970 13: 0.401306 0.822887 0.212930 0.938749 14: 0.754673 0.267132 0.714471 0.562088 15: 0.797287 0.692066 0.117707 0.178200E-01 16: 0.183837E-02 0.561662 0.299329 0.501103 Transformed array V: Col 1 2 3 4 Row 1: 3.81800 -0.386034 -0.107788 -0.277843 2: 0.752059E-02 -0.138549 0.188370 0.283179 3: 0.536878 -0.974688E-01 0.547781 0.197526 4: -0.148802E-01 -0.505449 -0.150938E-01 -0.590401E-01 5: 0.401441 0.603724E-01 -0.479618 -0.378976 6: 0.465916 0.375694 -0.514575 0.196521 7: 0.332723 0.104288 -0.191260 -0.187931 8: 0.421356 -0.324784 0.184047 -0.166718 9: 0.765858E-01 -0.211748 -0.642326 0.203540 10: 0.490840E-02 0.241689 0.934430E-01 0.481489 11: -0.316363 0.576181 0.165469 0.585867E-01 12: 0.322269 -0.136150 0.160147E-01 -0.112477 13: -0.317895 -0.560556 0.310220E-01 0.357514E-01 14: -0.542323E-02 -0.437067 -0.749236E-01 -0.495424E-01 15: 0.274031E-01 0.115707 -0.454561 -0.439101 16: 0.922587E-01 0.562418 0.332522 0.150830 Recovered array W: Col 1 2 3 4 Row 1: 0.218418 0.897504 0.861216 0.825003 2: 0.956318 0.350752 0.453794 0.824660 3: 0.829509 0.945448E-01 0.911977 0.618618E-01 4: 0.561695 0.136169E-01 0.597917 0.710781 5: 0.415307 0.859097 0.188955 0.882833E-01 6: 0.661187E-01 0.840847 0.761492 0.777994 7: 0.257578 0.123104 0.396988 0.745303 8: 0.109957 0.751236E-02 0.185314 0.308675 9: 0.438290E-01 0.260303 0.574366 0.899373 10: 0.633966 0.912484 0.367027 0.763537 11: 0.617272E-01 0.113664 0.617205 0.761731 12: 0.449539 0.351629 0.361529 0.406970 13: 0.401306 0.822887 0.212930 0.938749 14: 0.754673 0.267132 0.714471 0.562088 15: 0.797287 0.692066 0.117707 0.178200E-01 16: 0.183837E-02 0.561662 0.299329 0.501103 M = 37, N = 53, ||haar_2d_inverse(haar_2d(u))-u|| = 0.254887E-13 HAAR_PRB Normal end of execution. 6 March 2014 9:39:11.740 AM