17 February 2014 8:53:20.520 PM C4LIB_PRB FORTRAN90 version Test the C4LIB library. TEST01 C4_ABS computes the absolute value of a C4. C1=C4_UNIFORM_01 R2=C4_ABS(C1) R3=ABS(C1) --------------------- --------------------- --------------------- 0.449860 -0.126667 0.467352 0.467352 -0.843197 -0.344280 0.910774 0.910774 0.589627 0.260090 0.644443 0.644443 0.391140 0.323400 0.507521 0.507521 -0.139466 -0.156136 0.209354 0.209354 -0.236066 0.077459 0.248450 0.248450 0.018599 -0.633214 0.633487 0.633487 0.892850 0.010314 0.892909 0.892909 -0.560465 0.763795 0.947367 0.947367 0.306357 0.026275 0.307481 0.307481 TEST02 C4_ACOS computes the inverse cosine. C4_COS computes the cosine. C1=C4_UNIFORM_01 C2 = C4_ACOS(C1) C3 = C4_COS(C2) --------------------- --------------------- --------------------- 0.449860 -0.126667 1.109149 0.141009 0.449860 -0.126667 -0.843197 -0.344280 2.416418 0.498201 -0.843197 -0.344280 0.589627 0.260090 0.973444 -0.309595 0.589627 0.260090 0.391140 0.323400 1.192424 -0.341349 0.391140 0.323400 -0.139466 -0.156136 1.709000 0.156993 -0.139466 -0.156136 -0.236066 0.077459 1.808344 -0.079613 -0.236066 0.077459 0.018599 -0.633214 1.555082 0.597163 0.018599 -0.633214 0.892850 0.010314 0.467681 -0.022876 0.892850 0.010314 -0.560465 0.763795 2.013147 -0.767531 -0.560465 0.763795 0.306357 0.026275 1.259556 -0.027598 0.306357 0.026275 TEST03 C4_ACOSH computes the inverse hyperbolic cosine. C4_COSH computes the hyperbolic cosine. C1=C4_UNIFORM_01 C2 = C4_ACOSH(C1) C3 = COSH(C2) --------------------- --------------------- --------------------- 0.449860 -0.126667 -0.141009 1.109149 0.449860 -0.126667 -0.843197 -0.344280 -0.498201 2.416418 -0.843197 -0.344280 0.589627 0.260090 0.309595 0.973444 0.589627 0.260090 0.391140 0.323400 0.341349 1.192424 0.391140 0.323400 -0.139466 -0.156136 -0.156993 1.709000 -0.139466 -0.156136 -0.236066 0.077459 0.079613 1.808344 -0.236066 0.077459 0.018599 -0.633214 -0.597163 1.555082 0.018599 -0.633214 0.892850 0.010314 0.022876 0.467681 0.892850 0.010314 -0.560465 0.763795 0.767531 2.013147 -0.560465 0.763795 0.306357 0.026275 0.027598 1.259556 0.306357 0.026275 TEST04 C4_ADD adds two C4s C1=C4_UNIFORM_01 C2=C4_UNIFORM_01 C3 = C4_ADD(C1,C2) --------------------- --------------------- --------------------- 0.449860 -0.126667 -0.843197 -0.344280 -0.393337 -0.470947 0.589627 0.260090 0.391140 0.323400 0.980767 0.583490 -0.139466 -0.156136 -0.236066 0.077459 -0.375532 -0.078676 0.018599 -0.633214 0.892850 0.010314 0.911449 -0.622901 -0.560465 0.763795 0.306357 0.026275 -0.254108 0.790070 0.500804 -0.779931 0.350471 0.016555 0.851275 -0.763376 0.434989 -0.266623 -0.200947 0.270711 0.234042 0.004088 -0.097460 0.901881 -0.770245 -0.314308 -0.867705 0.587573 -0.889181 0.265655 -0.779871 -0.551160 -1.669052 -0.285505 0.031360 -0.433556 0.249089 0.578743 0.280450 0.145186 TEST05 C4_ARG computes the argument of a C4. C1=C4_UNIFORM_01 R2=C4_ARG(C1) --------------------- --------------------- 0.4499 -0.1267 -0.2745 -0.8432 -0.3443 -2.7539 0.5896 0.2601 0.4154 0.3911 0.3234 0.6909 -0.1395 -0.1561 -2.2999 -0.2361 0.0775 2.8245 0.0186 -0.6332 -1.5414 0.8928 0.0103 0.0116 -0.5605 0.7638 2.2038 0.3064 0.0263 0.0856 TEST06 C4_ASIN computes the inverse sine. C4_SIN computes the sine. C1=C4_UNIFORM_01 C2 = C4_ASIN(C1) C3 = C4_SIN(C2) --------------------- --------------------- --------------------- 0.449860 -0.126667 0.461648 -0.141009 0.449860 -0.126667 -0.843197 -0.344280 -0.845621 -0.498201 -0.843197 -0.344280 0.589627 0.260090 0.597353 0.309595 0.589627 0.260090 0.391140 0.323400 0.378372 0.341349 0.391140 0.323400 -0.139466 -0.156136 -0.138204 -0.156993 -0.139466 -0.156136 -0.236066 0.077459 -0.237548 0.079613 -0.236066 0.077459 0.018599 -0.633214 0.015714 -0.597163 0.018599 -0.633214 0.892850 0.010314 1.103115 0.022876 0.892850 0.010314 -0.560465 0.763795 -0.442351 0.767531 -0.560465 0.763795 0.306357 0.026275 0.311241 0.027598 0.306357 0.026275 TEST07 C4_ASINH computes the inverse hyperbolic sine. C4_SINH computes the hyperbolic sine. C1=C4_UNIFORM_01 C2 = C4_ASINH(C1) C3 = C4_SINH(C2) --------------------- --------------------- --------------------- 0.449860 -0.126667 0.438679 -0.115643 0.449860 -0.126667 -0.843197 -0.344280 -0.788715 -0.262344 -0.843197 -0.344280 0.589627 0.260090 0.572932 0.224433 0.589627 0.260090 0.391140 0.323400 0.399232 0.303890 0.391140 0.323400 -0.139466 -0.156136 -0.140698 -0.155225 -0.139466 -0.156136 -0.236066 0.077459 -0.234582 0.075447 -0.236066 0.077459 0.018599 -0.633214 0.024024 -0.685463 0.018599 -0.633214 0.892850 0.010314 0.803563 0.007693 0.892850 0.010314 -0.560465 0.763795 -0.665100 0.670334 -0.560465 0.763795 0.306357 0.026275 0.301849 0.025125 0.306357 0.026275 TEST08 C4_ATAN computes the inverse tangent. C4_TAN computes the tangent. C1=C4_UNIFORM_01 C2 = C4_ATAN(C1) C3 = C4_TAN(C2) --------------------- --------------------- --------------------- 0.449860 -0.126667 0.427774 -0.105499 0.449860 -0.126667 -0.843197 -0.344280 -0.735020 -0.197908 -0.843197 -0.344280 0.589627 0.260090 0.555251 0.192787 0.589627 0.260090 0.391140 0.323400 0.405767 0.284294 0.391140 0.323400 -0.139466 -0.156136 -0.141920 -0.154298 -0.139466 -0.156136 -0.236066 0.077459 -0.233099 0.073480 -0.236066 0.077459 0.018599 -0.633214 0.031027 -0.746154 0.018599 -0.633214 0.892850 0.010314 0.728880 0.005739 0.892850 0.010314 -0.560465 0.763795 -0.739806 0.556404 -0.560465 0.763795 0.306357 0.026275 0.297455 0.024024 0.306357 0.026275 TEST09 C4_ATANH computes the inverse hyperbolic tangent. C4_TANH computes the hyperbolic tangent. C1=C4_UNIFORM_01 C2 = C4_ATANH(C1) C3 = C4_TANH(C2) --------------------- --------------------- --------------------- 0.449860 -0.126667 0.473511 -0.156723 0.449860 -0.126667 -0.843197 -0.344280 -0.800349 -0.664037 -0.843197 -0.344280 0.589627 0.260090 0.599279 0.363538 0.589627 0.260090 0.391140 0.323400 0.364173 0.358338 0.391140 0.323400 -0.139466 -0.156136 -0.136933 -0.157832 -0.139466 -0.156136 -0.236066 0.077459 -0.239027 0.081817 -0.236066 0.077459 0.018599 -0.633214 0.013276 -0.564596 0.018599 -0.633214 0.892850 0.010314 1.433505 0.050703 0.892850 0.010314 -0.560465 0.763795 -0.339400 0.751900 -0.560465 0.763795 0.306357 0.026275 0.316262 0.028986 0.306357 0.026275 TEST10 C4_CONJ computes the conjugate of a C4. C1=C4_UNIFORM_01 C2=C4_CONJ(C1) C3=CONJG(C1) --------------------- --------------------- --------------------- 0.4499 -0.1267 0.4499 0.1267 0.4499 0.1267 -0.8432 -0.3443 -0.8432 0.3443 -0.8432 0.3443 0.5896 0.2601 0.5896 -0.2601 0.5896 -0.2601 0.3911 0.3234 0.3911 -0.3234 0.3911 -0.3234 -0.1395 -0.1561 -0.1395 0.1561 -0.1395 0.1561 -0.2361 0.0775 -0.2361 -0.0775 -0.2361 -0.0775 0.0186 -0.6332 0.0186 0.6332 0.0186 0.6332 0.8928 0.0103 0.8928 -0.0103 0.8928 -0.0103 -0.5605 0.7638 -0.5605 -0.7638 -0.5605 -0.7638 0.3064 0.0263 0.3064 -0.0263 0.3064 -0.0263 TEST11 C4_COS computes the cosine of a C4. C1=C4_UNIFORM_01 C2=C4_COS(C1) C3=COS(C1) --------------------- --------------------- --------------------- 0.4499 -0.1267 0.9077 0.0552 0.9077 0.0552 -0.8432 -0.3443 0.7049 -0.2622 0.7049 -0.2622 0.5896 0.2601 0.8594 -0.1463 0.8594 -0.1463 0.3911 0.3234 0.9732 -0.1255 0.9732 -0.1255 -0.1395 -0.1561 1.0024 -0.0218 1.0024 -0.0218 -0.2361 0.0775 0.9752 0.0181 0.9752 0.0181 0.0186 -0.6332 1.2071 0.0126 1.2071 0.0126 0.8928 0.0103 0.6272 -0.0080 0.6272 -0.0080 -0.5605 0.7638 1.1063 0.4467 1.1063 0.4467 0.3064 0.0263 0.9538 -0.0079 0.9538 -0.0079 TEST12 C4_COSH computes the hyperbolic cosine of a C4. C1=C4_UNIFORM_01 C2=C4_COSH(C1) --------------------- --------------------- 0.4499 -0.1267 1.0941 -0.0588 -0.8432 -0.3443 1.2963 0.3195 0.5896 0.2601 1.1393 0.1606 0.3911 0.3234 1.0216 0.1275 -0.1395 -0.1561 0.9975 0.0218 -0.2361 0.0775 1.0249 -0.0184 0.0186 -0.6332 0.8063 -0.0110 0.8928 0.0103 1.4257 0.0105 -0.5605 0.7638 0.8386 -0.4083 0.3064 0.0263 1.0469 0.0082 TEST13 C4_CUBE_ROOT computes the principal cube root of a C4. C1=C4_UNIFORM_01 C2=C4_CUBE_ROOT(C1) C3=C2*C2*C2 --------------------- --------------------- --------------------- 0.4499 -0.1267 0.7728 -0.0709 0.4499 -0.1267 -0.8432 -0.3443 0.5888 -0.7700 -0.8432 -0.3443 0.5896 0.2601 0.8555 0.1192 0.5896 0.2601 0.3911 0.3234 0.7766 0.1821 0.3911 0.3234 -0.1395 -0.1561 0.4277 -0.4119 -0.1395 -0.1561 -0.2361 0.0775 0.3700 0.5082 -0.2361 0.0775 0.0186 -0.6332 0.7479 -0.4221 0.0186 -0.6332 0.8928 0.0103 0.9629 0.0037 0.8928 0.0103 -0.5605 0.7638 0.7288 0.6583 -0.5605 0.7638 0.3064 0.0263 0.6747 0.0192 0.3064 0.0263 C1=C4_UNIFORM_01 C2=C1**(1.0/3.0) C3=C2*C2*C2 --------------------- --------------------- --------------------- 0.4499 -0.1267 0.7728 -0.0709 0.4499 -0.1267 -0.8432 -0.3443 0.5888 -0.7700 -0.8432 -0.3443 0.5896 0.2601 0.8555 0.1192 0.5896 0.2601 0.3911 0.3234 0.7766 0.1821 0.3911 0.3234 -0.1395 -0.1561 0.4277 -0.4119 -0.1395 -0.1561 -0.2361 0.0775 0.3700 0.5082 -0.2361 0.0775 0.0186 -0.6332 0.7479 -0.4221 0.0186 -0.6332 0.8928 0.0103 0.9629 0.0037 0.8928 0.0103 -0.5605 0.7638 0.7288 0.6583 -0.5605 0.7638 0.3064 0.0263 0.6747 0.0192 0.3064 0.0263 TEST14 C4_DIV computes C3 = C1 / C2. C1=C4_UNIFORM_01 C2=C4_UNIFORM_01 C3=C4_DIV(C1,C2) --------------------- --------------------- --------------------- 0.4499 -0.1267 -0.8432 -0.3443 -0.4047 0.3155 0.5896 0.2601 0.3911 0.3234 1.2219 -0.3453 -0.1395 -0.1561 -0.2361 0.0775 0.3374 0.7721 0.0186 -0.6332 0.8928 0.0103 0.0126 -0.7094 -0.5605 0.7638 0.3064 0.0263 -1.6038 2.6307 0.5008 -0.7799 0.3505 0.0166 1.3209 -2.2878 0.4350 -0.2666 -0.2009 0.2707 -1.4040 -0.5646 -0.0975 0.9019 -0.7702 -0.3143 -0.3011 -1.0480 -0.8892 0.2657 -0.7799 -0.5512 0.5998 -0.7646 0.0314 -0.4336 0.2491 0.5787 -0.6124 -0.3178 TEST15 C4_EXP computes exp ( Z ). C1=C4_UNIFORM_01 C2=C4_EXP(C1) C3=EXP(C1) --------------------- --------------------- --------------------- 0.4499 -0.1267 1.5555 -0.1981 1.5555 -0.1981 -0.8432 -0.3443 0.4051 -0.1452 0.4051 -0.1452 0.5896 0.2601 1.7427 0.4638 1.7427 0.4638 0.3911 0.3234 1.4020 0.4699 1.4020 0.4699 -0.1395 -0.1561 0.8592 -0.1353 0.8592 -0.1353 -0.2361 0.0775 0.7874 0.0611 0.7874 0.0611 0.0186 -0.6332 0.8213 -0.6028 0.8213 -0.6028 0.8928 0.0103 2.4419 0.0252 2.4419 0.0252 -0.5605 0.7638 0.4123 0.3949 0.4123 0.3949 0.3064 0.0263 1.3580 0.0357 1.3580 0.0357 TEST16 C4_INV computes C2 = 1 / C1. C1=C4_UNIFORM_01 C2=C4_INV(C1) C3=C4_INV(C2) --------------------- --------------------- --------------------- 0.4499 -0.1267 2.0596 0.5799 0.4499 -0.1267 -0.8432 -0.3443 -1.0165 0.4150 -0.8432 -0.3443 0.5896 0.2601 1.4197 -0.6263 0.5896 0.2601 0.3911 0.3234 1.5185 -1.2555 0.3911 0.3234 -0.1395 -0.1561 -3.1820 3.5624 -0.1395 -0.1561 -0.2361 0.0775 -3.8243 -1.2549 -0.2361 0.0775 0.0186 -0.6332 0.0463 1.5779 0.0186 -0.6332 0.8928 0.0103 1.1199 -0.0129 0.8928 0.0103 -0.5605 0.7638 -0.6245 -0.8510 -0.5605 0.7638 0.3064 0.0263 3.2403 -0.2779 0.3064 0.0263 TEST17 C4_LOG computes log ( Z ). C1=C4_UNIFORM_01 C2=C4_LOG(C1) C3=C4_EXP(C2) --------------------- --------------------- --------------------- 0.4499 -0.1267 -0.7607 -0.2745 0.4499 -0.1267 -0.8432 -0.3443 -0.0935 -2.7539 -0.8432 -0.3443 0.5896 0.2601 -0.4394 0.4154 0.5896 0.2601 0.3911 0.3234 -0.6782 0.6909 0.3911 0.3234 -0.1395 -0.1561 -1.5637 -2.2999 -0.1395 -0.1561 -0.2361 0.0775 -1.3925 2.8245 -0.2361 0.0775 0.0186 -0.6332 -0.4565 -1.5414 0.0186 -0.6332 0.8928 0.0103 -0.1133 0.0116 0.8928 0.0103 -0.5605 0.7638 -0.0541 2.2038 -0.5605 0.7638 0.3064 0.0263 -1.1793 0.0856 0.3064 0.0263 C1=C4_UNIFORM_01 C2=LOG(C1) C3=EXP(C2) --------------------- --------------------- --------------------- 0.4499 -0.1267 -0.7607 -0.2745 0.4499 -0.1267 -0.8432 -0.3443 -0.0935 -2.7539 -0.8432 -0.3443 0.5896 0.2601 -0.4394 0.4154 0.5896 0.2601 0.3911 0.3234 -0.6782 0.6909 0.3911 0.3234 -0.1395 -0.1561 -1.5637 -2.2999 -0.1395 -0.1561 -0.2361 0.0775 -1.3925 2.8245 -0.2361 0.0775 0.0186 -0.6332 -0.4565 -1.5414 0.0186 -0.6332 0.8928 0.0103 -0.1133 0.0116 0.8928 0.0103 -0.5605 0.7638 -0.0541 2.2038 -0.5605 0.7638 0.3064 0.0263 -1.1793 0.0856 0.3064 0.0263 TEST18 C4_MAG computes the magnitude of a C4. C1=C4_UNIFORM_01 R2=C4_MAG(C1) --------------------- --------------------- 0.4499 -0.1267 0.4674 -0.8432 -0.3443 0.9108 0.5896 0.2601 0.6444 0.3911 0.3234 0.5075 -0.1395 -0.1561 0.2094 -0.2361 0.0775 0.2484 0.0186 -0.6332 0.6335 0.8928 0.0103 0.8929 -0.5605 0.7638 0.9474 0.3064 0.0263 0.3075 TEST19 C4_MUL computes C3 = C1 * C2. C1=C4_UNIFORM_01 C2=C4_UNIFORM_01 C3=C4_MUL(C1,C2) --------------------- --------------------- --------------------- 0.4499 -0.1267 -0.8432 -0.3443 -0.4229 -0.0481 0.5896 0.2601 0.3911 0.3234 0.1465 0.2924 -0.1395 -0.1561 -0.2361 0.0775 0.0450 0.0261 0.0186 -0.6332 0.8928 0.0103 0.0231 -0.5652 -0.5605 0.7638 0.3064 0.0263 -0.1918 0.2193 0.5008 -0.7799 0.3505 0.0166 0.1884 -0.2651 0.4350 -0.2666 -0.2009 0.2707 -0.0152 0.1713 -0.0975 0.9019 -0.7702 -0.3143 0.3585 -0.6640 -0.8892 0.2657 -0.7799 -0.5512 0.8399 0.2829 0.0314 -0.4336 0.2491 0.5787 0.2587 -0.0898 TEST20 C4_NORMAL_01 generates unit pseudonormal C4s C1=C4_NORMAL_01(SEED) --------------------- 1.67904 -0.472768 -0.566060 -0.231124 1.21293 0.535037 1.26938 1.04954 -1.66609 -1.86523 -2.24246 0.735809 0.396752E-01 -1.35074 0.673068 0.777484E-02 -0.275127 0.374940 2.16400 0.185600 0.297785 -0.463758 2.04454 0.965770E-01 1.39882 -0.857394 -1.24299 1.67452 -0.670837E-01 0.620783 -0.794396 -0.324163 -0.523769 0.156483 -0.350567 -0.247757 0.131700 -1.82076 0.537380 1.24857 TEST21 C4_SIN computes the sine of a C4. C1=C4_UNIFORM_01 C2=C4_SIN(C1) C3=SIN(C1) --------------------- --------------------- --------------------- 0.4499 -0.1267 0.4383 -0.1144 0.4383 -0.1144 -0.8432 -0.3443 -0.7915 -0.2335 -0.7915 -0.2335 0.5896 0.2601 0.5750 0.2186 0.5750 0.2186 0.3911 0.3234 0.4014 0.3042 0.4014 0.3042 -0.1395 -0.1561 -0.1407 -0.1552 -0.1407 -0.1552 -0.2361 0.0775 -0.2346 0.0754 -0.2346 0.0754 0.0186 -0.6332 0.0225 -0.6763 0.0225 -0.6763 0.8928 0.0103 0.7789 0.0065 0.7789 0.0065 -0.5605 0.7638 -0.6943 0.7117 -0.6943 0.7117 0.3064 0.0263 0.3017 0.0251 0.3017 0.0251 TEST22 C4_SINH computes the hyperbolic sine of a C4. C1=C4_UNIFORM_01 C2=C4_SINH(C1) --------------------- --------------------- 0.4499 -0.1267 0.4615 -0.1393 -0.8432 -0.3443 -0.8912 -0.4648 0.5896 0.2601 0.6034 0.3032 0.3911 0.3234 0.3804 0.3424 -0.1395 -0.1561 -0.1382 -0.1570 -0.2361 0.0775 -0.2376 0.0795 0.0186 -0.6332 0.0150 -0.5918 0.8928 0.0103 1.0162 0.0147 -0.5605 0.7638 -0.4263 0.8032 0.3064 0.0263 0.3111 0.0275 TEST23 C4_SQRT computes the principal square root of a C4. C1=C4_UNIFORM_01 C2=C4_SQRT(C1) C3=C2*C2 --------------------- --------------------- --------------------- 0.4499 -0.1267 0.6772 -0.0935 0.4499 -0.1267 -0.8432 -0.3443 0.1838 -0.9365 -0.8432 -0.3443 0.5896 0.2601 0.7855 0.1656 0.5896 0.2601 0.3911 0.3234 0.6703 0.2412 0.3911 0.3234 -0.1395 -0.1561 0.1869 -0.4176 -0.1395 -0.1561 -0.2361 0.0775 0.0787 0.4922 -0.2361 0.0775 0.0186 -0.6332 0.5710 -0.5545 0.0186 -0.6332 0.8928 0.0103 0.9449 0.0055 0.8928 0.0103 -0.5605 0.7638 0.4398 0.8683 -0.5605 0.7638 0.3064 0.0263 0.5540 0.0237 0.3064 0.0263 C1=C4_UNIFORM_01 C2=SQRT(C1) C3=C2*C2 --------------------- --------------------- --------------------- 0.4499 -0.1267 0.6772 -0.0935 0.4499 -0.1267 -0.8432 -0.3443 0.1838 -0.9365 -0.8432 -0.3443 0.5896 0.2601 0.7855 0.1656 0.5896 0.2601 0.3911 0.3234 0.6703 0.2412 0.3911 0.3234 -0.1395 -0.1561 0.1869 -0.4176 -0.1395 -0.1561 -0.2361 0.0775 0.0787 0.4922 -0.2361 0.0775 0.0186 -0.6332 0.5710 -0.5545 0.0186 -0.6332 0.8928 0.0103 0.9449 0.0055 0.8928 0.0103 -0.5605 0.7638 0.4398 0.8683 -0.5605 0.7638 0.3064 0.0263 0.5540 0.0237 0.3064 0.0263 TEST24 C4_SUB subtracts two C4s C1=C4_UNIFORM_01 C2=C4_UNIFORM_01 C3 = C4_SUB(C1,C2) --------------------- --------------------- --------------------- 0.449860 -0.126667 -0.843197 -0.344280 1.293056 0.217613 0.589627 0.260090 0.391140 0.323400 0.198487 -0.063310 -0.139466 -0.156136 -0.236066 0.077459 0.096601 -0.233595 0.018599 -0.633214 0.892850 0.010314 -0.874250 -0.643528 -0.560465 0.763795 0.306357 0.026275 -0.866822 0.737520 0.500804 -0.779931 0.350471 0.016555 0.150333 -0.796486 0.434989 -0.266623 -0.200947 0.270711 0.635936 -0.537333 -0.097460 0.901881 -0.770245 -0.314308 0.672785 1.216188 -0.889181 0.265655 -0.779871 -0.551160 -0.109310 0.816815 0.031360 -0.433556 0.249089 0.578743 -0.217729 -1.012299 TEST25 C4_TAN computes the tangent of a C4. C1=C4_UNIFORM_01 C2=C4_TAN(C1) --------------------- --------------------- 0.4499 -0.1267 0.4735 -0.1548 -0.8432 -0.3443 -0.8781 -0.6579 0.5896 0.2601 0.6081 0.3579 0.3911 0.3234 0.3660 0.3598 -0.1395 -0.1561 -0.1369 -0.1579 -0.2361 0.0775 -0.2390 0.0817 0.0186 -0.6332 0.0128 -0.5604 0.8928 0.0103 1.2415 0.0262 -0.5605 0.7638 -0.3163 0.7710 0.3064 0.0263 0.3161 0.0289 TEST26 C4_TO_CARTESIAN converts C4 to (X,Y). CARTESIAN_TO_C4 converts (X,Y) to C4. C1=C4_UNIFORM_01 (X2,Y2)=C4_TO_CARTESIAN(C1) C3=CARTESIAN_TO_C4(X2,Y2) --------------------- --------------------- --------------------- 0.4499 -0.1267 0.4499 -0.1267 0.4499 -0.1267 -0.8432 -0.3443 -0.8432 -0.3443 -0.8432 -0.3443 0.5896 0.2601 0.5896 0.2601 0.5896 0.2601 0.3911 0.3234 0.3911 0.3234 0.3911 0.3234 -0.1395 -0.1561 -0.1395 -0.1561 -0.1395 -0.1561 -0.2361 0.0775 -0.2361 0.0775 -0.2361 0.0775 0.0186 -0.6332 0.0186 -0.6332 0.0186 -0.6332 0.8928 0.0103 0.8928 0.0103 0.8928 0.0103 -0.5605 0.7638 -0.5605 0.7638 -0.5605 0.7638 0.3064 0.0263 0.3064 0.0263 0.3064 0.0263 TEST27 C4_TO_POLAR converts C4 to (R,T). POLAR_TO_C4 converts (R,T) to C4. C1=C4_UNIFORM_01 (R2,T2)=C4_TO_POLAR(C1) C3=POLAR_TO_C4(R2,T2) --------------------- --------------------- --------------------- 0.4499 -0.1267 0.4674 -0.2745 0.4499 -0.1267 -0.8432 -0.3443 0.9108 -2.7539 -0.8432 -0.3443 0.5896 0.2601 0.6444 0.4154 0.5896 0.2601 0.3911 0.3234 0.5075 0.6909 0.3911 0.3234 -0.1395 -0.1561 0.2094 -2.2999 -0.1395 -0.1561 -0.2361 0.0775 0.2484 2.8245 -0.2361 0.0775 0.0186 -0.6332 0.6335 -1.5414 0.0186 -0.6332 0.8928 0.0103 0.8929 0.0116 0.8928 0.0103 -0.5605 0.7638 0.9474 2.2038 -0.5605 0.7638 0.3064 0.0263 0.3075 0.0856 0.3064 0.0263 TEST28 C4_UNIFORM_01 returns a uniformly random "unit" C4. C1=C4_UNIFORM_01(SEED) --------------------- 0.4499 -0.1267 -0.8432 -0.3443 0.5896 0.2601 0.3911 0.3234 -0.1395 -0.1561 -0.2361 0.0775 0.0186 -0.6332 0.8928 0.0103 -0.5605 0.7638 0.3064 0.0263 TEST29 C4MAT_UNIFORM_01 computes a "random" complex matrix. The matrix: Col: 1 2 3 4 Row --- 1: 0.450 -0.127 -0.236 0.775E-01 0.501 -0.780 -0.770 -0.314 2:-0.843 -0.344 0.186E-01-0.633 0.350 0.166E-01-0.889 0.266 3: 0.590 0.260 0.893 0.103E-01 0.435 -0.267 -0.780 -0.551 4: 0.391 0.323 -0.560 0.764 -0.201 0.271 0.314E-01-0.434 5:-0.139 -0.156 0.306 0.263E-01-0.975E-01 0.902 0.249 0.579 TEST30 C4VEC_INDICATOR sets A = (1-1i,2-2i,...,N-Ni) The "indicator" vector: 1: 1.00000 -1.00000 2: 2.00000 -2.00000 3: 3.00000 -3.00000 4: 4.00000 -4.00000 5: 5.00000 -5.00000 6: 6.00000 -6.00000 7: 7.00000 -7.00000 8: 8.00000 -8.00000 9: 9.00000 -9.00000 10: 10.0000 -10.0000 TEST31 C4VEC_SPIRAL returns N points on a spiral which includes M complete turns. The spiral points: 1: 5.00000 0.00000 2: 4.18579 2.41667 3: 2.33333 4.04145 4: -0.196701E-06 4.50000 5: -2.16667 3.75278 6: -3.60844 2.08333 7: -4.00000 -0.349691E-06 8: -3.31976 -1.91667 9: -1.83333 -3.17543 10: 0.417371E-07 -3.50000 11: 1.66667 -2.88675 12: 2.74241 -1.58333 13: 3.00000 0.524537E-06 TEST32 C4VEC_UNITY returns the N roots of unity The N roots of unity: 1: 1.00000 0.00000 2: 0.809017 0.587785 3: 0.309017 0.951057 4: -0.309017 0.951056 5: -0.809017 0.587785 6: -1.00000 -0.874228E-07 7: -0.809017 -0.587785 8: -0.309017 -0.951057 9: 0.309017 -0.951056 10: 0.809017 -0.587785 C4LIB_PRB Normal end of execution. 17 February 2014 8:53:20.525 PM