Moving array functions prototype and add some C99 Compliance

76 files changed, 272 insertions, 256 deletions

diff --git a/src/elementaryFunctions/acos/cacosa.c b/src/elementaryFunctions/acos/cacosa.c
index 8bbfaf2a..dcfa3826 100644
--- a/src/elementaryFunctions/acos/cacosa.c
+++ b/src/elementaryFunctions/acos/cacosa.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 10:23:49 2007 jofret
-** Last update Mon Mar 19 11:29:15 2007 jofret
+** Last update Mon Oct 22 10:00:02 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "acos.h"
 
-void cacosa(floatComplex* x, floatComplex* y, int size) {
+void cacosa(floatComplex* x, int size, floatComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = cacoss(x[i]);
diff --git a/src/elementaryFunctions/acos/dacosa.c b/src/elementaryFunctions/acos/dacosa.c
index c5e291a7..79a516c7 100644
--- a/src/elementaryFunctions/acos/dacosa.c
+++ b/src/elementaryFunctions/acos/dacosa.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 11:29:20 2007 jofret
-** Last update Mon Mar 19 11:29:03 2007 jofret
+** Last update Mon Oct 22 09:59:52 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "acos.h"
 
-void dacosa(double* x, double* y, int size) {
+void dacosa(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = dacoss(x[i]);
diff --git a/src/elementaryFunctions/acos/sacosa.c b/src/elementaryFunctions/acos/sacosa.c
index c905c76c..89209e00 100644
--- a/src/elementaryFunctions/acos/sacosa.c
+++ b/src/elementaryFunctions/acos/sacosa.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 10:25:57 2007 jofret
-** Last update Mon Mar 19 11:28:52 2007 jofret
+** Last update Mon Oct 22 09:58:58 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "acos.h"
 
-void sacosa(float* x, float* y, int size) {
+void sacosa(float* x, int size, float* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = sacoss(x[i]);
diff --git a/src/elementaryFunctions/acos/testAcos.c b/src/elementaryFunctions/acos/testAcos.c
index 3dc5e3c0..2ddd799d 100644
--- a/src/elementaryFunctions/acos/testAcos.c
+++ b/src/elementaryFunctions/acos/testAcos.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 Apr 20 10:24:42 2007 jofret
+** Last update Mon Oct 22 11:55:22 2007 bruno
 **
 ** Copyright INRIA 2006
 */
@@ -84,7 +84,7 @@ void sacosaTest(void) {
   int i = 0;
 
   printf(">> Float array\n");
-  sacosa(in, out, 5);
+  sacosa(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("sacosa(array) = %f\n", out[i]);
 }
@@ -95,37 +95,39 @@ void dacosaTest(void) {
   int i = 0;
 
   printf(">> Double Array\n");
-  dacosa(in, out, 5);
+  dacosa(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("sacosa(array) = %f\n", out[i]);
 
 }
 
 void cacosaTest(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 in[4];
   floatComplex out[4];
   int i = 0;
 
-  cacosa(in, out, 4);
+  in[0] = FloatComplex(PI, PI);
+  in[1] = FloatComplex(PI/2, PI/2);
+  in[2] = FloatComplex(PI/2, PI/3);
+  in[3] = FloatComplex(PI/2, PI/4);
+
+  cacosa(in, 4, out);
   printf(">> Float Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("cacosa(array) = %e + I * %e\n", creals(out[i]), cimags(out[i]));
 }
 
 void zacosaTest(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 in[4];
   doubleComplex out[4];
   int i = 0;
 
-  zacosa(in, out, 4);
+  in[0] = DoubleComplex(PI, PI);
+  in[1] = DoubleComplex(PI/2, PI/2);
+  in[2] = DoubleComplex(PI/2, PI/3);
+  in[3] = DoubleComplex(PI/2, PI/4);
+
+  zacosa(in, 4, out);
   printf(">> Double Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("zacosa(array) = %e + I * %e\n", zreals(out[i]), zimags(out[i]));
diff --git a/src/elementaryFunctions/acos/zacosa.c b/src/elementaryFunctions/acos/zacosa.c
index 9a01f76b..45caa6b4 100644
--- a/src/elementaryFunctions/acos/zacosa.c
+++ b/src/elementaryFunctions/acos/zacosa.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 10:25:14 2007 jofret
-** Last update Mon Mar 19 11:28:38 2007 jofret
+** Last update Mon Oct 22 09:58:47 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "acos.h"
 
-void zacosa(doubleComplex* x, doubleComplex* y, int size) {
+void zacosa(doubleComplex* x, int size, doubleComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = zacoss(x[i]);
diff --git a/src/elementaryFunctions/acosh/cacosha.c b/src/elementaryFunctions/acosh/cacosha.c
index 2191ab31..dd4fa6ba 100644
--- a/src/elementaryFunctions/acosh/cacosha.c
+++ b/src/elementaryFunctions/acosh/cacosha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 10:23:49 2007 jofret
-** Last update Thu Sep  6 10:44:18 2007 bruno
+** Last update Mon Oct 22 09:59:42 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "acosh.h"
 
-void cacosha(floatComplex* x, floatComplex* y, int size) {
+void cacosha(floatComplex* x, int size, floatComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = cacoshs(x[i]);
diff --git a/src/elementaryFunctions/acosh/dacosha.c b/src/elementaryFunctions/acosh/dacosha.c
index 27536bcb..f6ddc6fc 100644
--- a/src/elementaryFunctions/acosh/dacosha.c
+++ b/src/elementaryFunctions/acosh/dacosha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 11:29:20 2007 jofret
-** Last update Thu Sep  6 10:43:55 2007 bruno
+** Last update Mon Oct 22 09:59:29 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "acosh.h"
 
-void dacosha(double* x, double* y, int size) {
+void dacosha(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = dacoshs(x[i]);
diff --git a/src/elementaryFunctions/acosh/sacosha.c b/src/elementaryFunctions/acosh/sacosha.c
index c0c21ad5..450a8a3e 100644
--- a/src/elementaryFunctions/acosh/sacosha.c
+++ b/src/elementaryFunctions/acosh/sacosha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 10:25:57 2007 jofret
-** Last update Thu Sep  6 10:43:26 2007 bruno
+** Last update Mon Oct 22 09:59:20 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "acosh.h"
 
-void sacosha(float* x, float* y, int size) {
+void sacosha(float* x, int size, float* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = sacoshs(x[i]);
diff --git a/src/elementaryFunctions/acosh/testAcosh.c b/src/elementaryFunctions/acosh/testAcosh.c
index 8681dc60..fb49d8de 100644
--- a/src/elementaryFunctions/acosh/testAcosh.c
+++ b/src/elementaryFunctions/acosh/testAcosh.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 Thu Sep  6 10:42:52 2007 bruno
+** Last update Mon Oct 22 11:58:49 2007 bruno
 **
 ** Copyright INRIA 2006
 */
@@ -84,7 +84,7 @@ void sacoshaTest(void) {
   int i = 0;
 
   printf(">> Float array\n");
-  sacosha(in, out, 5);
+  sacosha(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("sacosha(array) = %f\n", out[i]);
 }
@@ -95,37 +95,39 @@ void dacoshaTest(void) {
   int i = 0;
 
   printf(">> Double Array\n");
-  dacosha(in, out, 5);
+  dacosha(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("sacosha(array) = %f\n", out[i]);
 
 }
 
 void cacoshaTest(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 in[4];
   floatComplex out[4];
   int i = 0;
 
-  cacosha(in, out, 4);
+  in[0] = FloatComplex(PI, PI);
+  in[1] = FloatComplex(PI/2, PI/2);
+  in[2] = FloatComplex(PI/2, PI/3);
+  in[3] = FloatComplex(PI/2, PI/4);
+
+  cacosha(in, 4, out);
   printf(">> Float Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("cacosha(array) = %e + I * %e\n", creals(out[i]), cimags(out[i]));
 }
 
 void zacoshaTest(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 in[4];
   doubleComplex out[4];
   int i = 0;
 
-  zacosha(in, out, 4);
+  in[0] = DoubleComplex(PI, PI);
+  in[1] = DoubleComplex(PI/2, PI/2);
+  in[2] = DoubleComplex(PI/2, PI/3);
+  in[3] = DoubleComplex(PI/2, PI/4);
+
+  zacosha(in, 4, out);
   printf(">> Double Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("zacosha(array) = %e + I * %e\n", zreals(out[i]), zimags(out[i]));
diff --git a/src/elementaryFunctions/acosh/zacosha.c b/src/elementaryFunctions/acosh/zacosha.c
index 1fca9a3d..08329151 100644
--- a/src/elementaryFunctions/acosh/zacosha.c
+++ b/src/elementaryFunctions/acosh/zacosha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 10:25:14 2007 jofret
-** Last update Thu Sep  6 10:42:27 2007 bruno
+** Last update Mon Oct 22 09:59:06 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "acosh.h"
 
-void zacosha(doubleComplex* x, doubleComplex* y, int size) {
+void zacosha(doubleComplex* x, int size, doubleComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = zacoshs(x[i]);
diff --git a/src/elementaryFunctions/asin/casina.c b/src/elementaryFunctions/asin/casina.c
index ed14f1f2..d10c5d9d 100644
--- a/src/elementaryFunctions/asin/casina.c
+++ b/src/elementaryFunctions/asin/casina.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 10:23:49 2007 jofret
-** Last update Mon Mar 19 11:30:55 2007 jofret
+** Last update Mon Oct 22 09:58:38 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "asin.h"
 
-void casina(floatComplex* x, floatComplex* y, int size) {
+void casina(floatComplex* x, int size, floatComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = casins(x[i]);
diff --git a/src/elementaryFunctions/asin/dasina.c b/src/elementaryFunctions/asin/dasina.c
index f6f3b1c1..46e231ae 100644
--- a/src/elementaryFunctions/asin/dasina.c
+++ b/src/elementaryFunctions/asin/dasina.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 11:29:20 2007 jofret
-** Last update Mon Mar 19 11:31:13 2007 jofret
+** Last update Mon Oct 22 09:58:27 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "asin.h"
 
-void dasina(double* x, double* y, int size) {
+void dasina(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = dasins(x[i]);
diff --git a/src/elementaryFunctions/asin/sasina.c b/src/elementaryFunctions/asin/sasina.c
index 2264ef34..98e4355f 100644
--- a/src/elementaryFunctions/asin/sasina.c
+++ b/src/elementaryFunctions/asin/sasina.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 10:25:57 2007 jofret
-** Last update Mon Mar 19 11:31:04 2007 jofret
+** Last update Mon Oct 22 09:57:35 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "asin.h"
 
-void sasina(float* x, float* y, int size) {
+void sasina(float* x, int size, float* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = sasins(x[i]);
diff --git a/src/elementaryFunctions/asin/testAsin.c b/src/elementaryFunctions/asin/testAsin.c
index 5ac58411..452077c9 100644
--- a/src/elementaryFunctions/asin/testAsin.c
+++ b/src/elementaryFunctions/asin/testAsin.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 Apr 20 10:32:47 2007 jofret
+** Last update Mon Oct 22 14:09:57 2007 bruno
 **
 ** Copyright INRIA 2006
 */
@@ -84,7 +84,7 @@ void sasinaTest(void) {
   int i = 0;
 
   printf(">> Float array\n");
-  sasina(in, out, 5);
+  sasina(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("sasina(array) = %f\n", out[i]);
 }
@@ -95,37 +95,39 @@ void dasinaTest(void) {
   int i = 0;
 
   printf(">> Double Array\n");
-  dasina(in, out, 5);
+  dasina(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("sasina(array) = %f\n", out[i]);
 
 }
 
 void casinaTest(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 in[4];
   floatComplex out[4];
   int i = 0;
 
-  casina(in, out, 4);
+  in[0] = FloatComplex(PI, PI);
+  in[1] = FloatComplex(PI/2, PI/2);
+  in[2] = FloatComplex(PI/2, PI/3);
+  in[3] = FloatComplex(PI/2, PI/4);
+
+  casina(in, 4, out);
   printf(">> Float Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("casina(array) = %e + I * %e\n", creals(out[i]), cimags(out[i]));
 }
 
 void zasinaTest(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 in[4];
   doubleComplex out[4];
   int i = 0;
 
-  zasina(in, out, 4);
+  in[0] = DoubleComplex(PI, PI);
+  in[1] = DoubleComplex(PI/2, PI/2);
+  in[2] = DoubleComplex(PI/2, PI/3);
+  in[3] = DoubleComplex(PI/2, PI/4);
+
+  zasina(in, 4, out);
   printf(">> Double Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("zasina(array) = %e + I * %e\n", zreals(out[i]), zimags(out[i]));
diff --git a/src/elementaryFunctions/asin/zasina.c b/src/elementaryFunctions/asin/zasina.c
index c11d9b07..cae0a271 100644
--- a/src/elementaryFunctions/asin/zasina.c
+++ b/src/elementaryFunctions/asin/zasina.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 10:25:14 2007 jofret
-** Last update Mon Mar 19 11:31:24 2007 jofret
+** Last update Mon Oct 22 09:57:21 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "asin.h"
 
-void zasina(doubleComplex* x, doubleComplex* y, int size) {
+void zasina(doubleComplex* x, int size, doubleComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = zasins(x[i]);
diff --git a/src/elementaryFunctions/asinh/casinha.c b/src/elementaryFunctions/asinh/casinha.c
index f8a286fa..32fc2aca 100644
--- a/src/elementaryFunctions/asinh/casinha.c
+++ b/src/elementaryFunctions/asinh/casinha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 10:23:49 2007 jofret
-** Last update Thu Sep  6 11:23:09 2007 bruno
+** Last update Mon Oct 22 09:58:15 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "asinh.h"
 
-void casinha(floatComplex* x, floatComplex* y, int size) {
+void casinha(floatComplex* x, int size, floatComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = casinhs(x[i]);
diff --git a/src/elementaryFunctions/asinh/dasinha.c b/src/elementaryFunctions/asinh/dasinha.c
index 8c37ab77..8e2acc75 100644
--- a/src/elementaryFunctions/asinh/dasinha.c
+++ b/src/elementaryFunctions/asinh/dasinha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 11:29:20 2007 jofret
-** Last update Thu Sep  6 11:22:33 2007 bruno
+** Last update Mon Oct 22 09:58:05 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "asinh.h"
 
-void dasinha(double* x, double* y, int size) {
+void dasinha(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = dasinhs(x[i]);
diff --git a/src/elementaryFunctions/asinh/sasinha.c b/src/elementaryFunctions/asinh/sasinha.c
index 65a1e9db..16916ef1 100644
--- a/src/elementaryFunctions/asinh/sasinha.c
+++ b/src/elementaryFunctions/asinh/sasinha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 10:25:57 2007 jofret
-** Last update Thu Sep  6 11:21:13 2007 bruno
+** Last update Mon Oct 22 09:57:55 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "asinh.h"
 
-void sasinha(float* x, float* y, int size) {
+void sasinha(float* x, int size, float* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = sasinhs(x[i]);
diff --git a/src/elementaryFunctions/asinh/testAsinh.c b/src/elementaryFunctions/asinh/testAsinh.c
index daefffc6..5c78625a 100644
--- a/src/elementaryFunctions/asinh/testAsinh.c
+++ b/src/elementaryFunctions/asinh/testAsinh.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 Thu Sep  6 11:20:54 2007 bruno
+** Last update Mon Oct 22 14:17:39 2007 bruno
 **
 ** Copyright INRIA 2006
 */
@@ -84,7 +84,7 @@ void sasinhaTest(void) {
   int i = 0;
 
   printf(">> Float array\n");
-  sasinha(in, out, 5);
+  sasinha(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("sasinha(array) = %f\n", out[i]);
 }
@@ -95,37 +95,39 @@ void dasinhaTest(void) {
   int i = 0;
 
   printf(">> Double Array\n");
-  dasinha(in, out, 5);
+  dasinha(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("sasinha(array) = %f\n", out[i]);
 
 }
 
 void casinhaTest(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 in[4];
   floatComplex out[4];
   int i = 0;
 
-  casinha(in, out, 4);
+  in[0] = FloatComplex(PI, PI);
+  in[1] = FloatComplex(PI/2, PI/2);
+  in[2] = FloatComplex(PI/2, PI/3);
+  in[3] = FloatComplex(PI/2, PI/4);
+
+  casinha(in, 4, out);
   printf(">> Float Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("casinha(array) = %e + I * %e\n", creals(out[i]), cimags(out[i]));
 }
 
 void zasinhaTest(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 in[4];
   doubleComplex out[4];
   int i = 0;
 
-  zasinha(in, out, 4);
+  in[0] = DoubleComplex(PI, PI);
+  in[1] = DoubleComplex(PI/2, PI/2);
+  in[2] = DoubleComplex(PI/2, PI/3);
+  in[3] = DoubleComplex(PI/2, PI/4);
+
+  zasinha(in, 4, out);
   printf(">> Double Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("zasinha(array) = %e + I * %e\n", zreals(out[i]), zimags(out[i]));
diff --git a/src/elementaryFunctions/asinh/zasinha.c b/src/elementaryFunctions/asinh/zasinha.c
index 9ac20b54..49b15bf3 100644
--- a/src/elementaryFunctions/asinh/zasinha.c
+++ b/src/elementaryFunctions/asinh/zasinha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Jan  5 10:25:14 2007 jofret
-** Last update Thu Sep  6 11:20:42 2007 bruno
+** Last update Mon Oct 22 09:57:45 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "asinh.h"
 
-void zasinha(doubleComplex* x, doubleComplex* y, int size) {
+void zasinha(doubleComplex* x, int size, doubleComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = zasinhs(x[i]);
diff --git a/src/elementaryFunctions/atan/catana.c b/src/elementaryFunctions/atan/catana.c
index d1652c22..66d23e62 100644
--- a/src/elementaryFunctions/atan/catana.c
+++ b/src/elementaryFunctions/atan/catana.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:24 2006 jofret
-** Last update Thu Sep  6 11:56:04 2007 bruno
+** Last update Mon Oct 22 09:57:09 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "atan.h"
 
-void catana(floatComplex* x, floatComplex* y, int size) {
+void catana(floatComplex* x, int size, floatComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = catans(x[i]);
diff --git a/src/elementaryFunctions/atan/datana.c b/src/elementaryFunctions/atan/datana.c
index c238ccc0..28a008ac 100644
--- a/src/elementaryFunctions/atan/datana.c
+++ b/src/elementaryFunctions/atan/datana.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:56 2006 jofret
-** Last update Thu Sep  6 11:54:07 2007 bruno
+** Last update Mon Oct 22 09:56:59 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "atan.h"
 
-void datana(double* x, double* y, int size) {
+void datana(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = datans(x[i]);
diff --git a/src/elementaryFunctions/atan/satana.c b/src/elementaryFunctions/atan/satana.c
index 9d9086f3..9e259f2f 100644
--- a/src/elementaryFunctions/atan/satana.c
+++ b/src/elementaryFunctions/atan/satana.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:03:27 2006 jofret
-** Last update Thu Sep  6 11:53:48 2007 bruno
+** Last update Mon Oct 22 09:56:10 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "atan.h"
 
-void satana(float* x, float* y, int size) {
+void satana(float* x, int size, float* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = satans(x[i]);
diff --git a/src/elementaryFunctions/atan/zatana.c b/src/elementaryFunctions/atan/zatana.c
index da5a9826..476d4662 100644
--- a/src/elementaryFunctions/atan/zatana.c
+++ b/src/elementaryFunctions/atan/zatana.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:12:02 2006 jofret
-** Last update Thu Sep  6 11:53:18 2007 bruno
+** Last update Mon Oct 22 09:56:00 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "atan.h"
 
-void zatana(doubleComplex* x, doubleComplex* y, int size) {
+void zatana(doubleComplex* x, int size, doubleComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = zatans(x[i]);
diff --git a/src/elementaryFunctions/atanh/catanha.c b/src/elementaryFunctions/atanh/catanha.c
index 6cae472d..d6b1986e 100644
--- a/src/elementaryFunctions/atanh/catanha.c
+++ b/src/elementaryFunctions/atanh/catanha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:24 2006 jofret
-** Last update Thu Sep  6 15:12:44 2007 bruno
+** Last update Mon Oct 22 09:56:50 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "atanh.h"
 
-void catanha(floatComplex* x, floatComplex* y, int size) {
+void catanha(floatComplex* x, int size, floatComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = catanhs(x[i]);
diff --git a/src/elementaryFunctions/atanh/datanha.c b/src/elementaryFunctions/atanh/datanha.c
index 3c2f62b6..149c6f86 100644
--- a/src/elementaryFunctions/atanh/datanha.c
+++ b/src/elementaryFunctions/atanh/datanha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:56 2006 jofret
-** Last update Thu Sep  6 15:12:16 2007 bruno
+** Last update Mon Oct 22 09:56:40 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "atanh.h"
 
-void datanha(double* x, double* y, int size) {
+void datanha(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = datanhs(x[i]);
diff --git a/src/elementaryFunctions/atanh/satanha.c b/src/elementaryFunctions/atanh/satanha.c
index 7aa3b544..d65f5876 100644
--- a/src/elementaryFunctions/atanh/satanha.c
+++ b/src/elementaryFunctions/atanh/satanha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:03:27 2006 jofret
-** Last update Thu Sep  6 15:10:59 2007 bruno
+** Last update Mon Oct 22 09:56:30 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "atanh.h"
 
-void satanha(float* x, float* y, int size) {
+void satanha(float* x, int size, float* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = satanhs(x[i]);
diff --git a/src/elementaryFunctions/atanh/zatanha.c b/src/elementaryFunctions/atanh/zatanha.c
index bef76801..659fae61 100644
--- a/src/elementaryFunctions/atanh/zatanha.c
+++ b/src/elementaryFunctions/atanh/zatanha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:12:02 2006 jofret
-** Last update Thu Sep  6 15:09:32 2007 bruno
+** Last update Mon Oct 22 09:56:21 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "atanh.h"
 
-void zatanha(doubleComplex* x, doubleComplex* y, int size) {
+void zatanha(doubleComplex* x, int size, doubleComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = zatanhs(x[i]);
diff --git a/src/elementaryFunctions/cos/ccosa.c b/src/elementaryFunctions/cos/ccosa.c
index 4b0ec511..1a626307 100644
--- a/src/elementaryFunctions/cos/ccosa.c
+++ b/src/elementaryFunctions/cos/ccosa.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:24 2006 jofret
-** Last update Mon Mar 19 11:30:33 2007 jofret
+** Last update Mon Oct 22 09:44:15 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "cos.h"
 
-void ccosa(floatComplex* x, floatComplex* y, int size) {
+void ccosa(floatComplex* x, int size, floatComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = ccoss(x[i]);
diff --git a/src/elementaryFunctions/cos/dcosa.c b/src/elementaryFunctions/cos/dcosa.c
index 07becd7e..0d04ddab 100644
--- a/src/elementaryFunctions/cos/dcosa.c
+++ b/src/elementaryFunctions/cos/dcosa.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:56 2006 jofret
-** Last update Mon Mar 19 11:30:22 2007 jofret
+** Last update Mon Oct 22 09:55:47 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "cos.h"
 
-void dcosa(double* x, double* y, int size) {
+void dcosa(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = dcoss(x[i]);
diff --git a/src/elementaryFunctions/cos/scosa.c b/src/elementaryFunctions/cos/scosa.c
index 9559bed9..6abece41 100644
--- a/src/elementaryFunctions/cos/scosa.c
+++ b/src/elementaryFunctions/cos/scosa.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:03:27 2006 jofret
-** Last update Mon Mar 19 11:30:11 2007 jofret
+** Last update Mon Oct 22 09:54:53 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "cos.h"
 
-void scosa(float* x, float* y, int size) {
+void scosa(float* x, int size, float* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = scoss(x[i]);
diff --git a/src/elementaryFunctions/cos/testCos.c b/src/elementaryFunctions/cos/testCos.c
index 3059567e..583908dd 100644
--- a/src/elementaryFunctions/cos/testCos.c
+++ b/src/elementaryFunctions/cos/testCos.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:12:04 2007 jofret
+** Last update Mon Oct 22 11:36:49 2007 bruno
 **
 ** Copyright INRIA 2006
 */
@@ -84,7 +84,7 @@ void scosaTest(void) {
   int i = 0;
 
   printf(">> Float array\n");
-  scosa(in, out, 5);
+  scosa(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("scosa(array) = %f\n", out[i]);
 }
@@ -95,37 +95,39 @@ void dcosaTest(void) {
   int i = 0;
 
   printf(">> Double Array\n");
-  dcosa(in, out, 5);
+  dcosa(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("scosa(array) = %f\n", out[i]);
 
 }
 
 void ccosaTest(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 in[4];
   floatComplex out[4];
   int i = 0;
 
-  ccosa(in, out, 4);
+  in[0] = FloatComplex(PI, PI);
+  in[1] = FloatComplex(PI/2, PI/2);
+  in[2] = FloatComplex(PI/2, PI/3);
+  in[3] = FloatComplex(PI/2, PI/4);
+
+  ccosa(in, 4, out);
   printf(">> Float Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("ccosa(array) = %e + I * %e\n", creals(out[i]), cimags(out[i]));
 }
 
 void zcosaTest(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 in[4];
   doubleComplex out[4];
   int i = 0;
 
-  zcosa(in, out, 4);
+  in[0] = DoubleComplex(PI, PI);
+  in[1] = DoubleComplex(PI/2, PI/2);
+  in[2] = DoubleComplex(PI/2, PI/3);
+  in[3] = DoubleComplex(PI/2, PI/4);
+
+  zcosa(in, 4, out);
   printf(">> Double Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("zcosa(array) = %e + I * %e\n", zreals(out[i]), zimags(out[i]));
diff --git a/src/elementaryFunctions/cos/zcosa.c b/src/elementaryFunctions/cos/zcosa.c
index a8936c98..a2964947 100644
--- a/src/elementaryFunctions/cos/zcosa.c
+++ b/src/elementaryFunctions/cos/zcosa.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:12:02 2006 jofret
-** Last update Mon Mar 19 11:30:44 2007 jofret
+** Last update Mon Oct 22 09:44:40 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "cos.h"
 
-void zcosa(doubleComplex* x, doubleComplex* y, int size) {
+void zcosa(doubleComplex* x, int size, doubleComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = zcoss(x[i]);
diff --git a/src/elementaryFunctions/cosh/ccosha.c b/src/elementaryFunctions/cosh/ccosha.c
index 28731c11..62f548f8 100644
--- a/src/elementaryFunctions/cosh/ccosha.c
+++ b/src/elementaryFunctions/cosh/ccosha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:24 2006 jofret
-** Last update Mon Mar 19 11:32:02 2007 jofret
+** Last update Mon Oct 22 09:55:37 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "cosh.h"
 
-void ccosha(floatComplex* x, floatComplex* y, int size) {
+void ccosha(floatComplex* x, int size, floatComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = ccoshs(x[i]);
diff --git a/src/elementaryFunctions/cosh/dcosha.c b/src/elementaryFunctions/cosh/dcosha.c
index 56cbce38..67c5ac7f 100644
--- a/src/elementaryFunctions/cosh/dcosha.c
+++ b/src/elementaryFunctions/cosh/dcosha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:56 2006 jofret
-** Last update Mon Mar 19 11:31:34 2007 jofret
+** Last update Mon Oct 22 09:55:27 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "cosh.h"
 
-void dcosha(double* x, double* y, int size) {
+void dcosha(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = dcoshs(x[i]);
diff --git a/src/elementaryFunctions/cosh/scosha.c b/src/elementaryFunctions/cosh/scosha.c
index 302f19f2..f63c2c61 100644
--- a/src/elementaryFunctions/cosh/scosha.c
+++ b/src/elementaryFunctions/cosh/scosha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:03:27 2006 jofret
-** Last update Mon Mar 19 11:31:44 2007 jofret
+** Last update Mon Oct 22 09:55:17 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "cosh.h"
 
-void scosha(float* x, float* y, int size) {
+void scosha(float* x, int size, float* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = scoshs(x[i]);
diff --git a/src/elementaryFunctions/cosh/testCosh.c b/src/elementaryFunctions/cosh/testCosh.c
index cb2bb614..35933d3b 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 Apr 20 16:41:30 2007 jofret
+** Last update Mon Oct 22 11:52:54 2007 bruno
 **
 ** Copyright INRIA 2006
 */
@@ -84,7 +84,7 @@ void scoshaTest(void) {
   int i = 0;
 
   printf(">> Float array\n");
-  scosha(in, out, 5);
+  scosha(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("scosha(array) = %f\n", out[i]);
 }
@@ -95,37 +95,39 @@ void dcoshaTest(void) {
   int i = 0;
 
   printf(">> Double Array\n");
-  dcosha(in, out, 5);
+  dcosha(in, 5, out);
   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 in[4];
   floatComplex out[4];
   int i = 0;
 
-  ccosha(in, out, 4);
+  in[0] = FloatComplex(PI, PI);
+  in[1] = FloatComplex(PI/2, PI/2);
+  in[2] = FloatComplex(PI/2, PI/3);
+  in[3] = FloatComplex(PI/2, PI/4);
+
+  ccosha(in, 4, out);
   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 in[4];
   doubleComplex out[4];
   int i = 0;
 
-  zcosha(in, out, 4);
+  in[0] = DoubleComplex(PI, PI);
+  in[1] = DoubleComplex(PI/2, PI/2);
+  in[2] = DoubleComplex(PI/2, PI/3);
+  in[3] = DoubleComplex(PI/2, PI/4);
+
+  zcosha(in, 4, out);
   printf(">> Double Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("zcosha(array) = %e + I * %e\n", zreals(out[i]), zimags(out[i]));
diff --git a/src/elementaryFunctions/cosh/zcosha.c b/src/elementaryFunctions/cosh/zcosha.c
index 73599d1f..5c79933e 100644
--- a/src/elementaryFunctions/cosh/zcosha.c
+++ b/src/elementaryFunctions/cosh/zcosha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:12:02 2006 jofret
-** Last update Mon Mar 19 11:31:53 2007 jofret
+** Last update Mon Oct 22 09:55:03 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "cosh.h"
 
-void zcosha(doubleComplex* x, doubleComplex* y, int size) {
+void zcosha(doubleComplex* x, int size, doubleComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = zcoshs(x[i]);
diff --git a/src/elementaryFunctions/exp/cexpa.c b/src/elementaryFunctions/exp/cexpa.c
index 3dca3bf5..0380e43f 100644
--- a/src/elementaryFunctions/exp/cexpa.c
+++ b/src/elementaryFunctions/exp/cexpa.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:24 2006 jofret
-** Last update Mon Mar 19 11:35:22 2007 jofret
+** Last update Mon Oct 22 09:53:24 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "exp.h"
 
-void cexpa(floatComplex* x, floatComplex* y, int size) {
+void cexpa(floatComplex* x, int size, floatComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = cexps(x[i]);
diff --git a/src/elementaryFunctions/exp/dexpa.c b/src/elementaryFunctions/exp/dexpa.c
index c4a258d0..037d89e5 100644
--- a/src/elementaryFunctions/exp/dexpa.c
+++ b/src/elementaryFunctions/exp/dexpa.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:56 2006 jofret
-** Last update Mon Mar 19 11:34:48 2007 jofret
+** Last update Mon Oct 22 09:53:15 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "exp.h"
 
-void dexpa(double* x, double* y, int size) {
+void dexpa(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = dexps(x[i]);
diff --git a/src/elementaryFunctions/exp/sexpa.c b/src/elementaryFunctions/exp/sexpa.c
index 4b2a2bd8..3c576971 100644
--- a/src/elementaryFunctions/exp/sexpa.c
+++ b/src/elementaryFunctions/exp/sexpa.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:03:27 2006 jofret
-** Last update Mon Mar 19 11:35:01 2007 jofret
+** Last update Mon Oct 22 09:53:06 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "exp.h"
 
-void sexpa(float* x, float* y, int size) {
+void sexpa(float* x, int size, float* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = sexps(x[i]);
diff --git a/src/elementaryFunctions/exp/zexpa.c b/src/elementaryFunctions/exp/zexpa.c
index 54006f66..dda49422 100644
--- a/src/elementaryFunctions/exp/zexpa.c
+++ b/src/elementaryFunctions/exp/zexpa.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:12:02 2006 jofret
-** Last update Mon Mar 19 11:35:12 2007 jofret
+** Last update Mon Oct 22 09:52:57 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "exp.h"
 
-void zexpa(doubleComplex* x, doubleComplex* y, int size) {
+void zexpa(doubleComplex* x, int size, doubleComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = zexps(x[i]);
diff --git a/src/elementaryFunctions/exp10/cexp10a.c b/src/elementaryFunctions/exp10/cexp10a.c
index 932f42ba..a1245ffb 100644
--- a/src/elementaryFunctions/exp10/cexp10a.c
+++ b/src/elementaryFunctions/exp10/cexp10a.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:24 2006 jofret
-** Last update Thu Sep  6 16:58:36 2007 bruno
+** Last update Mon Oct 22 09:54:11 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "exp10.h"
 
-void cexp10a(floatComplex* x, floatComplex* y, int size) {
+void cexp10a(floatComplex* x, int size, floatComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = cexp10s(x[i]);
diff --git a/src/elementaryFunctions/exp10/dexp10a.c b/src/elementaryFunctions/exp10/dexp10a.c
index c243508a..270abc53 100644
--- a/src/elementaryFunctions/exp10/dexp10a.c
+++ b/src/elementaryFunctions/exp10/dexp10a.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:56 2006 jofret
-** Last update Thu Sep  6 16:57:38 2007 bruno
+** Last update Mon Oct 22 09:53:56 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "exp10.h"
 
-void dexp10a(double* x, double* y, int size) {
+void dexp10a(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = dexp10s(x[i]);
diff --git a/src/elementaryFunctions/exp10/sexp10a.c b/src/elementaryFunctions/exp10/sexp10a.c
index f76089f7..259661dd 100644
--- a/src/elementaryFunctions/exp10/sexp10a.c
+++ b/src/elementaryFunctions/exp10/sexp10a.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:03:27 2006 jofret
-** Last update Thu Sep  6 16:57:07 2007 bruno
+** Last update Mon Oct 22 09:53:46 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "exp10.h"
 
-void sexp10a(float* x, float* y, int size) {
+void sexp10a(float* x, int size, float* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = sexp10s(x[i]);
diff --git a/src/elementaryFunctions/exp10/zexp10a.c b/src/elementaryFunctions/exp10/zexp10a.c
index 7f2ae965..1ce8bb1c 100644
--- a/src/elementaryFunctions/exp10/zexp10a.c
+++ b/src/elementaryFunctions/exp10/zexp10a.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:12:02 2006 jofret
-** Last update Thu Sep  6 16:56:32 2007 bruno
+** Last update Mon Oct 22 09:53:37 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "exp10.h"
 
-void zexp10a(doubleComplex* x, doubleComplex* y, int size) {
+void zexp10a(doubleComplex* x, int size, doubleComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = zexp10s(x[i]);
diff --git a/src/elementaryFunctions/log/cloga.c b/src/elementaryFunctions/log/cloga.c
index e2dc87b2..533c0dfd 100644
--- a/src/elementaryFunctions/log/cloga.c
+++ b/src/elementaryFunctions/log/cloga.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Sep  7 10:32:52 2007 bruno
-** Last update Fri Sep  7 10:37:56 2007 bruno
+** Last update Mon Oct 22 09:52:01 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "log.h"
 
-void cloga(floatComplex* in, floatComplex* out, int size) {
+void cloga(floatComplex* in, int size, floatComplex* out) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     out[i] = clogs(in[i]);
diff --git a/src/elementaryFunctions/log/dloga.c b/src/elementaryFunctions/log/dloga.c
index dbc2d49d..35aabcad 100644
--- a/src/elementaryFunctions/log/dloga.c
+++ b/src/elementaryFunctions/log/dloga.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Sep  7 10:31:35 2007 bruno
-** Last update Fri Sep  7 10:32:23 2007 bruno
+** Last update Mon Oct 22 09:51:49 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "log.h"
 
-void dloga(double* x, double* y, int size) {
+void dloga(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = slogs(x[i]);
diff --git a/src/elementaryFunctions/log/sloga.c b/src/elementaryFunctions/log/sloga.c
index 73ca53e3..e6da9104 100644
--- a/src/elementaryFunctions/log/sloga.c
+++ b/src/elementaryFunctions/log/sloga.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Sep  7 10:29:01 2007 bruno
-** Last update Fri Sep  7 10:30:59 2007 bruno
+** Last update Mon Oct 22 09:51:33 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "log.h"
 
-void		sloga(float* in, float* out, int size) {
+void		sloga(float* in, int size, float* out) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     out[i] = slogs(in[i]);
diff --git a/src/elementaryFunctions/log/zloga.c b/src/elementaryFunctions/log/zloga.c
index 79048803..8ea92bcd 100644
--- a/src/elementaryFunctions/log/zloga.c
+++ b/src/elementaryFunctions/log/zloga.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Sep  7 10:38:10 2007 bruno
-** Last update Fri Sep  7 10:38:57 2007 bruno
+** Last update Mon Oct 22 09:51:22 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "log.h"
 
-void zloga(doubleComplex* in, doubleComplex* out, int size) {
+void zloga(doubleComplex* in, int size, doubleComplex* out) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     out[i] = zlogs(in[i]);
diff --git a/src/elementaryFunctions/log10/clog10a.c b/src/elementaryFunctions/log10/clog10a.c
index f38c2651..2228c869 100644
--- a/src/elementaryFunctions/log10/clog10a.c
+++ b/src/elementaryFunctions/log10/clog10a.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Sep  6 16:23:40 2007 bruno
-** Last update Fri Sep  7 10:49:01 2007 bruno
+** Last update Mon Oct 22 09:52:46 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "log10.h"
 
-void		clog10a(floatComplex* in, floatComplex* out, int size) {
+void		clog10a(floatComplex* in, int size, floatComplex* out) {
   int i = 0;
   for (i = 0 ; i < size ; ++i) {
     out[i] = clog10s(in[i]);
diff --git a/src/elementaryFunctions/log10/dlog10a.c b/src/elementaryFunctions/log10/dlog10a.c
index 14ceb067..0eadbc52 100644
--- a/src/elementaryFunctions/log10/dlog10a.c
+++ b/src/elementaryFunctions/log10/dlog10a.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Sep  6 16:23:09 2007 bruno
-** Last update Fri Sep  7 10:43:49 2007 bruno
+** Last update Mon Oct 22 09:52:35 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "log10.h"
 
-void		dlog10a(double* in, double* out, int size) {
+void		dlog10a(double* in, int size, double* out) {
   int i = 0;
   for (i = 0 ; i < size ; ++i) {
     out[i] = dlog10s(in[i]);
diff --git a/src/elementaryFunctions/log10/slog10a.c b/src/elementaryFunctions/log10/slog10a.c
index e161d8d0..7eac52e7 100644
--- a/src/elementaryFunctions/log10/slog10a.c
+++ b/src/elementaryFunctions/log10/slog10a.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Sep  6 16:22:18 2007 bruno
-** Last update Fri Sep  7 10:42:46 2007 bruno
+** Last update Mon Oct 22 09:52:23 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "log10.h"
 
-void		slog10a(float* in, float* out, int size) {
+void		slog10a(float* in, int size, float* out) {
   int i = 0;
   for (i = 0 ; i < size ; ++i) {
     out[i] = slog10s(in[i]);
diff --git a/src/elementaryFunctions/log10/zlog10a.c b/src/elementaryFunctions/log10/zlog10a.c
index c4036c81..789f8f98 100644
--- a/src/elementaryFunctions/log10/zlog10a.c
+++ b/src/elementaryFunctions/log10/zlog10a.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Sep  6 16:21:57 2007 bruno
-** Last update Fri Sep  7 10:50:59 2007 bruno
+** Last update Mon Oct 22 09:52:13 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "log10.h"
 
-void		zlog10a(doubleComplex* in, doubleComplex* out, int size) {
+void		zlog10a(doubleComplex* in, int size, doubleComplex* out) {
   int i = 0;
   for (i = 0 ; i < size ; ++i) {
     out[i] = zlog10s(in[i]);
diff --git a/src/elementaryFunctions/sin/csina.c b/src/elementaryFunctions/sin/csina.c
index 4e400055..6392af41 100644
--- a/src/elementaryFunctions/sin/csina.c
+++ b/src/elementaryFunctions/sin/csina.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Dec  8 13:07:37 2006 jofret
-** Last update Mon Mar 19 11:36:02 2007 jofret
+** Last update Mon Oct 22 09:51:11 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "sin.h"
 
-void csina(floatComplex* x, floatComplex* y, int size) {
+void csina(floatComplex* x, int size, floatComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = csins(x[i]);
diff --git a/src/elementaryFunctions/sin/dsina.c b/src/elementaryFunctions/sin/dsina.c
index 4cd138b1..f86508f3 100644
--- a/src/elementaryFunctions/sin/dsina.c
+++ b/src/elementaryFunctions/sin/dsina.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Dec  8 11:02:19 2006 jofret
-** Last update Mon Mar 19 11:35:52 2007 jofret
+** Last update Mon Oct 22 09:51:01 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "sin.h"
 
-void dsina(double* x, double* y, int size) {
+void dsina(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = dsins(x[i]);
diff --git a/src/elementaryFunctions/sin/ssina.c b/src/elementaryFunctions/sin/ssina.c
index 585c5136..5d4c43be 100644
--- a/src/elementaryFunctions/sin/ssina.c
+++ b/src/elementaryFunctions/sin/ssina.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Dec  8 10:52:14 2006 jofret
-** Last update Mon Mar 19 11:35:43 2007 jofret
+** Last update Mon Oct 22 09:50:09 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "sin.h"
 
-void ssina(float* x, float* y, int size) {
+void ssina(float* x, int size, float* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = ssins(x[i]);
diff --git a/src/elementaryFunctions/sin/testSin.c b/src/elementaryFunctions/sin/testSin.c
index 1ed78db0..b7bc9ca6 100644
--- a/src/elementaryFunctions/sin/testSin.c
+++ b/src/elementaryFunctions/sin/testSin.c
@@ -5,7 +5,7 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Dec  8 15:06:16 2006 jofret
-** Last update Fri Apr 20 11:25:31 2007 jofret
+** Last update Mon Oct 22 12:05:02 2007 bruno
 **
 ** Copyright INRIA 2006
 */
@@ -84,7 +84,7 @@ void ssinaTest(void) {
   int i = 0;
 
   printf(">> Float array\n");
-  ssina(in, out, 5);
+  ssina(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("ssina(array) = %f\n", out[i]);
 }
@@ -95,37 +95,39 @@ void dsinaTest(void) {
   int i = 0;
 
   printf(">> Double Array\n");
-  dsina(in, out, 5);
+  dsina(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("ssina(array) = %f\n", out[i]);
 
 }
 
 void csinaTest(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 in[4];
   floatComplex out[4];
   int i = 0;
 
-  csina(in, out, 4);
+  in[0] = FloatComplex(PI, PI);
+  in[1] = FloatComplex(PI/2, PI/2);
+  in[2] = FloatComplex(PI/2, PI/3);
+  in[3] = FloatComplex(PI/2, PI/4);
+
+  csina(in, 4, out);
   printf(">> Float Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("csina(array) = %e + I * %e\n", creals(out[i]), cimags(out[i]));
 }
 
 void zsinaTest(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 in[4];
   doubleComplex out[4];
   int i = 0;
 
-  zsina(in, out, 4);
+  in[0] = DoubleComplex(PI, PI);
+  in[1] = DoubleComplex(PI/2, PI/2);
+  in[2] = DoubleComplex(PI/2, PI/3);
+  in[3] = DoubleComplex(PI/2, PI/4);
+
+  zsina(in, 4, out);
   printf(">> Double Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("zsina(array) = %e + I * %e\n", zreals(out[i]), zimags(out[i]));
diff --git a/src/elementaryFunctions/sin/zsina.c b/src/elementaryFunctions/sin/zsina.c
index c50efcd9..e5b25272 100644
--- a/src/elementaryFunctions/sin/zsina.c
+++ b/src/elementaryFunctions/sin/zsina.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Dec  8 14:02:04 2006 jofret
-** Last update Mon Mar 19 11:36:29 2007 jofret
+** Last update Mon Oct 22 09:49:58 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "sin.h"
 
-void zsina(doubleComplex* x, doubleComplex* y, int size) {
+void zsina(doubleComplex* x, int size, doubleComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = zsins(x[i]);
diff --git a/src/elementaryFunctions/sinh/csinha.c b/src/elementaryFunctions/sinh/csinha.c
index b63ac1af..fca522cc 100644
--- a/src/elementaryFunctions/sinh/csinha.c
+++ b/src/elementaryFunctions/sinh/csinha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Dec  8 13:07:37 2006 jofret
-** Last update Mon Mar 19 11:32:14 2007 jofret
+** Last update Mon Oct 22 09:50:51 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "sinh.h"
 
-void csinha(floatComplex* x, floatComplex* y, int size) {
+void csinha(floatComplex* x, int size, floatComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = csinhs(x[i]);
diff --git a/src/elementaryFunctions/sinh/dsinha.c b/src/elementaryFunctions/sinh/dsinha.c
index c6d7ede6..93dcaadf 100644
--- a/src/elementaryFunctions/sinh/dsinha.c
+++ b/src/elementaryFunctions/sinh/dsinha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Dec  8 11:02:19 2006 jofret
-** Last update Mon Mar 19 11:32:24 2007 jofret
+** Last update Mon Oct 22 09:50:41 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "sinh.h"
 
-void dsinha(double* x, double* y, int size) {
+void dsinha(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = dsinhs(x[i]);
diff --git a/src/elementaryFunctions/sinh/ssinha.c b/src/elementaryFunctions/sinh/ssinha.c
index 43846371..f4ca7ea5 100644
--- a/src/elementaryFunctions/sinh/ssinha.c
+++ b/src/elementaryFunctions/sinh/ssinha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Dec  8 10:52:14 2006 jofret
-** Last update Mon Mar 19 11:32:34 2007 jofret
+** Last update Mon Oct 22 09:50:30 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "sinh.h"
 
-void ssinha(float* x, float* y, int size) {
+void ssinha(float* x, int size, float* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = ssinhs(x[i]);
diff --git a/src/elementaryFunctions/sinh/testSinh.c b/src/elementaryFunctions/sinh/testSinh.c
index a63d08ff..7fe00e6d 100644
--- a/src/elementaryFunctions/sinh/testSinh.c
+++ b/src/elementaryFunctions/sinh/testSinh.c
@@ -5,7 +5,7 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Dec  8 15:06:16 2006 jofret
-** Last update Fri Apr 20 16:44:27 2007 jofret
+** Last update Mon Oct 22 13:39:10 2007 bruno
 **
 ** Copyright INRIA 2006
 */
@@ -84,7 +84,7 @@ void ssinhaTest(void) {
   int i = 0;
 
   printf(">> Float array\n");
-  ssinha(in, out, 5);
+  ssinha(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("ssinha(array) = %f\n", out[i]);
 }
@@ -95,37 +95,39 @@ void dsinhaTest(void) {
   int i = 0;
 
   printf(">> Double Array\n");
-  dsinha(in, out, 5);
+  dsinha(in, 5, out);
   for (i = 0 ; i < 5 ; ++i)
     printf("ssinha(array) = %f\n", out[i]);
 
 }
 
 void csinhaTest(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 in[4];
   floatComplex out[4];
   int i = 0;
 
-  csinha(in, out, 4);
+  in[0] = FloatComplex(PI, PI);
+  in[1] = FloatComplex(PI/2, PI/2);
+  in[2] = FloatComplex(PI/2, PI/3);
+  in[3] = FloatComplex(PI/2, PI/4);
+
+  csinha(in, 4, out);
   printf(">> Float Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("csinha(array) = %e + I * %e\n", creals(out[i]), cimags(out[i]));
 }
 
 void zsinhaTest(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 in[4];
   doubleComplex out[4];
   int i = 0;
 
-  zsinha(in, out, 4);
+  in[0] = DoubleComplex(PI, PI);
+  in[1] = DoubleComplex(PI/2, PI/2);
+  in[2] = DoubleComplex(PI/2, PI/3);
+  in[3] = DoubleComplex(PI/2, PI/4);
+
+  zsinha(in, 4, out);
   printf(">> Double Complex Array\n");
   for (i = 0 ; i < 4 ; ++i)
     printf("zsinha(array) = %e + I * %e\n", zreals(out[i]), zimags(out[i]));
diff --git a/src/elementaryFunctions/sinh/zsinha.c b/src/elementaryFunctions/sinh/zsinha.c
index c4c36d67..e3725dbb 100644
--- a/src/elementaryFunctions/sinh/zsinha.c
+++ b/src/elementaryFunctions/sinh/zsinha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Dec  8 14:02:04 2006 jofret
-** Last update Mon Mar 19 11:32:43 2007 jofret
+** Last update Mon Oct 22 09:50:21 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "sinh.h"
 
-void zsinha(doubleComplex* x, doubleComplex* y, int size) {
+void zsinha(doubleComplex* x, int size, doubleComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = zsinhs(x[i]);
diff --git a/src/elementaryFunctions/sqrt/csqrta.c b/src/elementaryFunctions/sqrt/csqrta.c
index e2b1b795..62851ff1 100644
--- a/src/elementaryFunctions/sqrt/csqrta.c
+++ b/src/elementaryFunctions/sqrt/csqrta.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Sep  7 11:40:55 2007 bruno
-** Last update Fri Sep  7 11:41:27 2007 bruno
+** Last update Mon Oct 22 09:49:44 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "sqrt.h"
 
-void		csqrta(floatComplex* in, floatComplex* out, int size) {
+void		csqrta(floatComplex* in, int size, floatComplex* out) {
   int i = 0;
   for (i = 0 ; i < size ; ++i) {
     out[i] = csqrts(in[i]);
diff --git a/src/elementaryFunctions/sqrt/dsqrta.c b/src/elementaryFunctions/sqrt/dsqrta.c
index 0206002e..56d6ea1f 100644
--- a/src/elementaryFunctions/sqrt/dsqrta.c
+++ b/src/elementaryFunctions/sqrt/dsqrta.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Sep  7 11:40:21 2007 bruno
-** Last update Fri Sep  7 11:40:42 2007 bruno
+** Last update Mon Oct 22 09:49:34 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "sqrt.h"
 
-void		dsqrta(double* in, double* out, int size) {
+void		dsqrta(double* in, int size, double* out) {
   int i = 0;
   for (i = 0 ; i < size ; ++i) {
     out[i] = dsqrts(in[i]);
diff --git a/src/elementaryFunctions/sqrt/ssqrta.c b/src/elementaryFunctions/sqrt/ssqrta.c
index bd585ec5..6ffa5c26 100644
--- a/src/elementaryFunctions/sqrt/ssqrta.c
+++ b/src/elementaryFunctions/sqrt/ssqrta.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Sep  7 11:06:03 2007 bruno
-** Last update Fri Sep  7 11:39:54 2007 bruno
+** Last update Mon Oct 22 09:49:22 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "sqrt.h"
 
-void		ssqrta(float* in, float* out, int size) {
+void		ssqrta(float* in, int size, float* out) {
   int i = 0;
   for (i = 0 ; i < size ; ++i) {
     out[i] = ssqrts(in[i]);
diff --git a/src/elementaryFunctions/sqrt/zsqrta.c b/src/elementaryFunctions/sqrt/zsqrta.c
index d1553237..92e6fb1c 100644
--- a/src/elementaryFunctions/sqrt/zsqrta.c
+++ b/src/elementaryFunctions/sqrt/zsqrta.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Fri Sep  7 11:41:40 2007 bruno
-** Last update Fri Sep  7 11:42:01 2007 bruno
+** Last update Mon Oct 22 09:49:12 2007 bruno
 **
 ** Copyright INRIA 2007
 */
 
 #include "sqrt.h"
 
-void		zsqrta(doubleComplex* in, doubleComplex* out, int size) {
+void		zsqrta(doubleComplex* in, int size, doubleComplex* out) {
   int i = 0;
   for (i = 0 ; i < size ; ++i) {
     out[i] = zsqrts(in[i]);
diff --git a/src/elementaryFunctions/tan/ctana.c b/src/elementaryFunctions/tan/ctana.c
index 54f589dd..339dda27 100644
--- a/src/elementaryFunctions/tan/ctana.c
+++ b/src/elementaryFunctions/tan/ctana.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:24 2006 jofret
-** Last update Mon Mar 19 11:32:52 2007 jofret
+** Last update Mon Oct 22 09:49:00 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "tan.h"
 
-void ctana(floatComplex* x, floatComplex* y, int size) {
+void ctana(floatComplex* x, int size, floatComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = ctans(x[i]);
diff --git a/src/elementaryFunctions/tan/dtana.c b/src/elementaryFunctions/tan/dtana.c
index 3c3efe7e..d74dd3bf 100644
--- a/src/elementaryFunctions/tan/dtana.c
+++ b/src/elementaryFunctions/tan/dtana.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:56 2006 jofret
-** Last update Mon Mar 19 11:33:02 2007 jofret
+** Last update Mon Oct 22 09:48:46 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "tan.h"
 
-void dtana(double* x, double* y, int size) {
+void dtana(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = dtans(x[i]);
diff --git a/src/elementaryFunctions/tan/stana.c b/src/elementaryFunctions/tan/stana.c
index abb2cd86..95cc287e 100644
--- a/src/elementaryFunctions/tan/stana.c
+++ b/src/elementaryFunctions/tan/stana.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:03:27 2006 jofret
-** Last update Mon Mar 19 11:33:11 2007 jofret
+** Last update Mon Oct 22 09:46:19 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "tan.h"
 
-void stana(float* x, float* y, int size) {
+void stana(float* x, int size, float* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = stans(x[i]);
diff --git a/src/elementaryFunctions/tan/ztana.c b/src/elementaryFunctions/tan/ztana.c
index b0d2477e..ef296d6a 100644
--- a/src/elementaryFunctions/tan/ztana.c
+++ b/src/elementaryFunctions/tan/ztana.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:12:02 2006 jofret
-** Last update Mon Mar 19 11:33:44 2007 jofret
+** Last update Mon Oct 22 09:46:07 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "tan.h"
 
-void ztana(doubleComplex* x, doubleComplex* y, int size) {
+void ztana(doubleComplex* x, int size, doubleComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = ztans(x[i]);
diff --git a/src/elementaryFunctions/tanh/ctanha.c b/src/elementaryFunctions/tanh/ctanha.c
index 8fd0d72d..cd6f7acb 100644
--- a/src/elementaryFunctions/tanh/ctanha.c
+++ b/src/elementaryFunctions/tanh/ctanha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:24 2006 jofret
-** Last update Mon Mar 19 11:33:54 2007 jofret
+** Last update Mon Oct 22 09:47:06 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "tanh.h"
 
-void ctanha(floatComplex* x, floatComplex* y, int size) {
+void ctanha(floatComplex* x, int size, floatComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = ctanhs(x[i]);
diff --git a/src/elementaryFunctions/tanh/dtanha.c b/src/elementaryFunctions/tanh/dtanha.c
index 9cf0e8de..71330f9c 100644
--- a/src/elementaryFunctions/tanh/dtanha.c
+++ b/src/elementaryFunctions/tanh/dtanha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 14:54:56 2006 jofret
-** Last update Mon Mar 19 11:34:04 2007 jofret
+** Last update Mon Oct 22 09:46:56 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "tanh.h"
 
-void dtanha(double* x, double* y, int size) {
+void dtanha(double* x, int size, double* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = dtanhs(x[i]);
diff --git a/src/elementaryFunctions/tanh/stanha.c b/src/elementaryFunctions/tanh/stanha.c
index f93e1bca..5b41f68b 100644
--- a/src/elementaryFunctions/tanh/stanha.c
+++ b/src/elementaryFunctions/tanh/stanha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:03:27 2006 jofret
-** Last update Mon Mar 19 11:34:25 2007 jofret
+** Last update Mon Oct 22 09:46:45 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "tanh.h"
 
-void stanha(float* x, float* y, int size) {
+void stanha(float* x, int size, float* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = stanhs(x[i]);
diff --git a/src/elementaryFunctions/tanh/ztanha.c b/src/elementaryFunctions/tanh/ztanha.c
index b1039b9f..a27ff436 100644
--- a/src/elementaryFunctions/tanh/ztanha.c
+++ b/src/elementaryFunctions/tanh/ztanha.c
@@ -5,14 +5,14 @@
 ** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
 **
 ** Started on  Thu Dec  7 16:12:02 2006 jofret
-** Last update Mon Mar 19 11:34:36 2007 jofret
+** Last update Mon Oct 22 09:46:32 2007 bruno
 **
 ** Copyright INRIA 2006
 */
 
 #include "tanh.h"
 
-void ztanha(doubleComplex* x, doubleComplex* y, int size) {
+void ztanha(doubleComplex* x, int size, doubleComplex* y) {
   int i = 0;
   for (i = 0; i < size; ++i) {
     y[i] = ztanhs(x[i]);