[RELEASE] Cosh is now stable.

1 files changed, 96 insertions, 10 deletions

diff --git a/src/elementaryFunctions/cosh/testCosh.c b/src/elementaryFunctions/cosh/testCosh.c
index ddb6fc79..cb2bb614 100644
--- a/src/elementaryFunctions/cosh/testCosh.c
+++ b/src/elementaryFunctions/cosh/testCosh.c
@@ -5,7 +5,7 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Dec  8 15:05:44 2006 jofret
-** Last update Fri Mar 30 11:42:57 2007 jofret
+** Last update Fri Apr 20 16:41:30 2007 jofret
 **
 ** Copyright INRIA 2006
 */
@@ -42,19 +42,105 @@ void dcoshsTest(void) {
   printf("dcoshs(-PI/6) = %e\n", dcoshs(-PI/6));
 }
 
+void ccoshsTest(void) {
+  floatComplex pi_pi = FloatComplex(PI, PI);
+  floatComplex pi_2_pi_2 = FloatComplex(PI/2, PI/2);
+  floatComplex pi_2_pi_3 = FloatComplex(PI/2, PI/3);
+  floatComplex pi_2_pi_4 = FloatComplex(PI/2, PI/4);
+  floatComplex out;
+
+  printf(">> Float Complex scalar\n");
+  out = ccoshs(pi_pi);
+  printf("ccoshs(PI + I*PI) = %f + I * %f\n", creals(out), cimags(out));
+  out = ccoshs(pi_2_pi_2);
+  printf("ccoshs(PI/2 + I*PI/2) = %f + I * %f\n", creals(out), cimags(out));
+  out = ccoshs(pi_2_pi_3);
+  printf("ccoshs(PI/2 + I*PI/3) = %f + I * %f\n", creals(out), cimags(out));
+  out = ccoshs(pi_2_pi_4);
+  printf("ccoshs(PI/2 + I*PI/4) = %f + I * %f\n", creals(out), cimags(out));
+}
+
+void zcoshsTest(void) {
+  doubleComplex pi_pi = DoubleComplex(PI, PI);
+  doubleComplex pi_2_pi_2 = DoubleComplex(PI/2, PI/2);
+  doubleComplex pi_2_pi_3 = DoubleComplex(PI/2, PI/3);
+  doubleComplex pi_2_pi_4 = DoubleComplex(PI/2, PI/4);
+  doubleComplex out;
+
+  printf(">> Double Complex scalar\n");
+  out = zcoshs(pi_pi);
+  printf("zcoshs(PI + I*PI) = %e + I * %e\n", zreals(out), zimags(out));
+  out = zcoshs(pi_2_pi_2);
+  printf("zcoshs(PI/2 + I*PI/2) = %e + I * %e\n", zreals(out), zimags(out));
+  out = zcoshs(pi_2_pi_3);
+  printf("zcoshs(PI/2 + I*PI/3) = %e + I * %e\n", zreals(out), zimags(out));
+  out = zcoshs(pi_2_pi_4);
+  printf("zcoshs(PI/2 + I*PI/4) = %e + I * %e\n", zreals(out), zimags(out));
+}
+
+void scoshaTest(void) {
+  float out[5];
+  float in[5] = {PI, PI/2, PI/3, PI/4, PI/6};
+  int i = 0;
+
+  printf(">> Float array\n");
+  scosha(in, out, 5);
+  for (i = 0 ; i < 5 ; ++i)
+    printf("scosha(array) = %f\n", out[i]);
+}
+
+void dcoshaTest(void) {
+  double out[5];
+  double in[5] = {PI, PI/2, PI/3, PI/4, PI/6};
+  int i = 0;
+
+  printf(">> Double Array\n");
+  dcosha(in, out, 5);
+  for (i = 0 ; i < 5 ; ++i)
+    printf("scosha(array) = %f\n", out[i]);
+
+}
+
+void ccoshaTest(void) {
+  floatComplex pi_pi = FloatComplex(PI, PI);
+  floatComplex pi_2_pi_2 = FloatComplex(PI/2, PI/2);
+  floatComplex pi_2_pi_3 = FloatComplex(PI/2, PI/3);
+  floatComplex pi_2_pi_4 = FloatComplex(PI/2, PI/4);
+  floatComplex in[4] = {pi_pi, pi_2_pi_2, pi_2_pi_3, pi_2_pi_4};
+  floatComplex out[4];
+  int i = 0;
+
+  ccosha(in, out, 4);
+  printf(">> Float Complex Array\n");
+  for (i = 0 ; i < 4 ; ++i)
+    printf("ccosha(array) = %e + I * %e\n", creals(out[i]), cimags(out[i]));
+}
+
+void zcoshaTest(void) {
+  doubleComplex pi_pi = DoubleComplex(PI, PI);
+  doubleComplex pi_2_pi_2 = DoubleComplex(PI/2, PI/2);
+  doubleComplex pi_2_pi_3 = DoubleComplex(PI/2, PI/3);
+  doubleComplex pi_2_pi_4 = DoubleComplex(PI/2, PI/4);
+  doubleComplex in[4] = {pi_pi, pi_2_pi_2, pi_2_pi_3, pi_2_pi_4 };
+  doubleComplex out[4];
+  int i = 0;
+
+  zcosha(in, out, 4);
+  printf(">> Double Complex Array\n");
+  for (i = 0 ; i < 4 ; ++i)
+    printf("zcosha(array) = %e + I * %e\n", zreals(out[i]), zimags(out[i]));
+}
+
 int testCosh(void) {
   printf("\n>>>> Hyperbolic Cosine Tests\n");
   scoshsTest();
   dcoshsTest();
-  /* FIXME : Implement some test here ... */
-  /*
-    ccoshsTest();
-    zcoshsTest();
-    scoshaTest();
-    dcoshaTest();
-    ccoshaTest();
-    zcoshaTest();
-  */
+  ccoshsTest();
+  zcoshsTest();
+  scoshaTest();
+  dcoshaTest();
+  ccoshaTest();
+  zcoshaTest();
   return 0;
 }