16 March 2011 11:11:17.913 AM WALSH_PRB FORTRAN77 version Test the WALSH library. TEST01 FWT computes a fast Walsh transform. I X(I) Y=FWT(X)/N Z=FWT(Y)/N 1 0.2184 0.4099 0.2184 2 0.9563 0.0169 0.9563 3 0.8295 0.1552 0.8295 4 0.5617 0.0594 0.5617 5 0.4153 -0.0325 0.4153 6 0.0661 -0.0088 0.0661 7 0.2576 -0.0391 0.2576 8 0.1100 0.0263 0.1100 9 0.0438 -0.1347 0.0438 10 0.6340 0.0342 0.6340 11 0.0617 -0.1346 0.0617 12 0.4495 -0.0163 0.4495 13 0.4013 -0.1492 0.4013 14 0.7547 0.0283 0.7547 15 0.7973 0.0352 0.7973 16 0.0018 -0.0318 0.0018 I X(I) Y=FWT(X)/N Z=FWT(Y)/N 1 1.0000 8.5000 1.0000 2 2.0000 -4.0000 2.0000 3 3.0000 0.0000 3.0000 4 4.0000 -2.0000 4.0000 5 5.0000 0.0000 5.0000 6 6.0000 0.0000 6.0000 7 7.0000 0.0000 7.0000 8 8.0000 -1.0000 8.0000 9 9.0000 0.0000 9.0000 10 10.0000 0.0000 10.0000 11 11.0000 0.0000 11.0000 12 12.0000 0.0000 12.0000 13 13.0000 0.0000 13.0000 14 14.0000 0.0000 14.0000 15 15.0000 0.0000 15.0000 16 16.0000 -0.5000 16.0000 TEST02 WALSH computes a fast Walsh transform. I X(I) Y=FWT(X)/N Z=FWT(Y)/N 1 0.2184 0.4099 0.2184 2 0.9563 0.0169 0.9563 3 0.8295 0.1552 0.8295 4 0.5617 0.0594 0.5617 5 0.4153 -0.0325 0.4153 6 0.0661 -0.0088 0.0661 7 0.2576 -0.0391 0.2576 8 0.1100 0.0263 0.1100 9 0.0438 -0.1347 0.0438 10 0.6340 0.0342 0.6340 11 0.0617 -0.1346 0.0617 12 0.4495 -0.0163 0.4495 13 0.4013 -0.1492 0.4013 14 0.7547 0.0283 0.7547 15 0.7973 0.0352 0.7973 16 0.0018 -0.0318 0.0018 I X(I) Y=FWT(X)/N Z=FWT(Y)/N 1 1.0000 8.5000 1.0000 2 2.0000 -4.0000 2.0000 3 3.0000 0.0000 3.0000 4 4.0000 -2.0000 4.0000 5 5.0000 0.0000 5.0000 6 6.0000 -0.0000 6.0000 7 7.0000 0.0000 7.0000 8 8.0000 -1.0000 8.0000 9 9.0000 0.0000 9.0000 10 10.0000 -0.0000 10.0000 11 11.0000 0.0000 11.0000 12 12.0000 -0.0000 12.0000 13 13.0000 0.0000 13.0000 14 14.0000 -0.0000 14.0000 15 15.0000 0.0000 15.0000 16 16.0000 -0.5000 16.0000 TEST03 HAAR computes a Haar transform. HNORM normalizes the transformed data. HAARIN computes an inverse Haar transform. I X(I) Y=HAAR(X) Z=HNORM(Y) W=HAARIN(Z) 1 0.2184 6.5591 0.4099 0.2184 2 0.9563 0.2707 0.0169 0.9563 3 0.8295 1.7170 0.2146 0.8295 4 0.5617 -0.7660 -0.0958 0.5617 5 0.4153 -0.2165 -0.0541 0.4153 6 0.0661 0.1139 0.0285 0.0661 7 0.2576 0.1665 0.0416 0.2576 8 0.1100 0.3569 0.0892 0.1100 9 0.0438 -0.7379 -0.3689 0.0438 10 0.6340 0.2678 0.1339 0.6340 11 0.0617 0.3492 0.1746 0.0617 12 0.4495 0.1476 0.0738 0.4495 13 0.4013 -0.5901 -0.2951 0.4013 14 0.7547 -0.3878 -0.1939 0.7547 15 0.7973 -0.3534 -0.1767 0.7973 16 0.0018 0.7954 0.3977 0.0018 I X(I) Y=HAAR(X) Z=HNORM(Y) W=HAARIN(Z) 1 1.0000 136.0000 8.5000 1.0000 2 2.0000 -64.0000 -4.0000 2.0000 3 3.0000 -16.0000 -2.0000 3.0000 4 4.0000 -16.0000 -2.0000 4.0000 5 5.0000 -4.0000 -1.0000 5.0000 6 6.0000 -4.0000 -1.0000 6.0000 7 7.0000 -4.0000 -1.0000 7.0000 8 8.0000 -4.0000 -1.0000 8.0000 9 9.0000 -1.0000 -0.5000 9.0000 10 10.0000 -1.0000 -0.5000 10.0000 11 11.0000 -1.0000 -0.5000 11.0000 12 12.0000 -1.0000 -0.5000 12.0000 13 13.0000 -1.0000 -0.5000 13.0000 14 14.0000 -1.0000 -0.5000 14.0000 15 15.0000 -1.0000 -0.5000 15.0000 16 16.0000 -1.0000 -0.5000 16.0000 TEST04 FFWT computes a fast Walsh transform. I X(I) Y=FFWT(X)/N Z=FFWT(Y)/N 1 0.2184 0.4099 0.2184 2 0.9563 0.0169 0.9563 3 0.8295 0.1552 0.8295 4 0.5617 0.0594 0.5617 5 0.4153 -0.0325 0.4153 6 0.0661 -0.0088 0.0661 7 0.2576 -0.0391 0.2576 8 0.1100 0.0263 0.1100 9 0.0438 -0.1347 0.0438 10 0.6340 0.0342 0.6340 11 0.0617 -0.1346 0.0617 12 0.4495 -0.0163 0.4495 13 0.4013 -0.1492 0.4013 14 0.7547 0.0283 0.7547 15 0.7973 0.0352 0.7973 16 0.0018 -0.0318 0.0018 I X(I) Y=FFWT(X)/N Z=FFWT(Y)/N 1 1.0000 8.5000 1.0000 2 2.0000 -4.0000 2.0000 3 3.0000 0.0000 3.0000 4 4.0000 -2.0000 4.0000 5 5.0000 0.0000 5.0000 6 6.0000 0.0000 6.0000 7 7.0000 0.0000 7.0000 8 8.0000 -1.0000 8.0000 9 9.0000 0.0000 9.0000 10 10.0000 0.0000 10.0000 11 11.0000 0.0000 11.0000 12 12.0000 0.0000 12.0000 13 13.0000 0.0000 13.0000 14 14.0000 0.0000 14.0000 15 15.0000 0.0000 15.0000 16 16.0000 -0.5000 16.0000 TEST05 FRT computes a fast Walsh transform that is indifferent to circular shifts. I X(I) Y=FRT(X)/N 1 -7.8477 6.5591 2 -0.4395 -0.5415 3 0.8153 -3.4340 4 0.6162 -1.0992 5 0.4153 3.9301 6 0.3237 -0.1053 7 0.3829 1.0678 8 0.1100 -6.6224 9 0.0438 -1.5242 10 0.6340 1.2439 11 0.0617 3.1388 12 0.4495 -0.6587 13 0.4013 5.4658 14 0.7547 23.2869 15 0.7973 -1.7701 16 0.0018 7.4955 I X(I) Y=FRT(X)/N 1 7.6615 6.5591 2 4.6155 -0.1090 3 -0.6463 -1.4108 4 4.6106 0.0245 5 0.5617 12.4756 6 0.4814 2.3904 7 -0.6984 -1.4147 8 0.2576 1.0470 9 0.1100 -3.6708 10 0.0438 0.4937 11 0.6340 1.1440 12 0.0617 -1.1762 13 0.4495 5.1452 14 0.4013 -2.2151 15 0.7547 -0.5013 16 0.7973 -1.4813 I X(I) Y=FRT(X)/N 1 4.2715 6.5591 2 -8.2478 -4.9660 3 -0.0834 1.2217 4 0.1741 -1.3939 5 0.2184 -0.0670 6 1.7858 -0.6607 7 -0.2536 1.4501 8 0.5617 -1.2590 9 0.4153 -0.7474 10 0.0661 4.3795 11 0.2576 -3.5079 12 0.1100 -0.6045 13 0.0438 2.9883 14 0.6340 -7.1295 15 0.0617 2.4680 16 0.4495 0.9682 WALSH_PRB Normal end of execution. 16 March 2011 11:11:17.916 AM