Adding sqrt functions.

Need more precision and think about type returned.

6 files changed, 116 insertions, 0 deletions

diff --git a/src/elementaryFunctions/sqrt/csqrta.c b/src/elementaryFunctions/sqrt/csqrta.c
index e69de29b..e2b1b795 100644
--- a/src/elementaryFunctions/sqrt/csqrta.c
+++ b/src/elementaryFunctions/sqrt/csqrta.c
@@ -0,0 +1,20 @@
+/*
+**  -*- C -*-
+**
+** csqrta.c
+** 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
+**
+** Copyright INRIA 2007
+*/
+
+#include "sqrt.h"
+
+void		csqrta(floatComplex* in, floatComplex* out, int size) {
+  int i = 0;
+  for (i = 0 ; i < size ; ++i) {
+    out[i] = csqrts(in[i]);
+  }
+}
diff --git a/src/elementaryFunctions/sqrt/csqrts.c b/src/elementaryFunctions/sqrt/csqrts.c
index e69de29b..ed1f8345 100644
--- a/src/elementaryFunctions/sqrt/csqrts.c
+++ b/src/elementaryFunctions/sqrt/csqrts.c
@@ -0,0 +1,18 @@
+/*
+**  -*- C -*-
+**
+** csqrts.c
+** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
+**
+** Started on  Fri Sep  7 11:04:08 2007 bruno
+** Last update Fri Sep  7 11:04:56 2007 bruno
+**
+** Copyright INRIA 2007
+*/
+
+#include "sqrt.h"
+
+floatComplex	csqrts(floatComplex in) {
+  /* FIXME : Dummy ... */
+  return FloatComplex(0, 0);
+}
diff --git a/src/elementaryFunctions/sqrt/dsqrta.c b/src/elementaryFunctions/sqrt/dsqrta.c
index e69de29b..0206002e 100644
--- a/src/elementaryFunctions/sqrt/dsqrta.c
+++ b/src/elementaryFunctions/sqrt/dsqrta.c
@@ -0,0 +1,20 @@
+/*
+**  -*- C -*-
+**
+** dsqrta.c
+** 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
+**
+** Copyright INRIA 2007
+*/
+
+#include "sqrt.h"
+
+void		dsqrta(double* in, double* out, int size) {
+  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 e69de29b..bd585ec5 100644
--- a/src/elementaryFunctions/sqrt/ssqrta.c
+++ b/src/elementaryFunctions/sqrt/ssqrta.c
@@ -0,0 +1,20 @@
+/*
+**  -*- C -*-
+**
+** ssqrta.c
+** 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
+**
+** Copyright INRIA 2007
+*/
+
+#include "sqrt.h"
+
+void		ssqrta(float* in, float* out, int size) {
+  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 e69de29b..d1553237 100644
--- a/src/elementaryFunctions/sqrt/zsqrta.c
+++ b/src/elementaryFunctions/sqrt/zsqrta.c
@@ -0,0 +1,20 @@
+/*
+**  -*- C -*-
+**
+** zsqrta.c
+** 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
+**
+** Copyright INRIA 2007
+*/
+
+#include "sqrt.h"
+
+void		zsqrta(doubleComplex* in, doubleComplex* out, int size) {
+  int i = 0;
+  for (i = 0 ; i < size ; ++i) {
+    out[i] = zsqrts(in[i]);
+  }
+}
diff --git a/src/elementaryFunctions/sqrt/zsqrts.c b/src/elementaryFunctions/sqrt/zsqrts.c
index e69de29b..9d01b432 100644
--- a/src/elementaryFunctions/sqrt/zsqrts.c
+++ b/src/elementaryFunctions/sqrt/zsqrts.c
@@ -0,0 +1,18 @@
+/*
+**  -*- C -*-
+**
+** zsqrts.c
+** Made by  Bruno JOFRET <bruno.jofret@inria.fr>
+**
+** Started on  Fri Sep  7 11:05:10 2007 bruno
+** Last update Fri Sep  7 11:05:48 2007 bruno
+**
+** Copyright INRIA 2007
+*/
+
+#include "sqrt.h"
+
+doubleComplex	zsqrts(doubleComplex in) {
+  /* FIXME : Dummy ... */
+  return DoubleComplex(0, 0);
+}