Adding scilab algorithm

Update License

1 files changed, 26 insertions, 12 deletions

diff --git a/src/elementaryFunctions/tan/ctans.c b/src/elementaryFunctions/tan/ctans.c
index ee3396f1..9a6bd8c1 100644
--- a/src/elementaryFunctions/tan/ctans.c
+++ b/src/elementaryFunctions/tan/ctans.c
@@ -10,22 +10,36 @@
  *
  */
 
+#include <math.h>
+#include "lapack.h"
 #include "tan.h"
-/*
-  #include "log.h"
-*/
 #include "sqrt.h"
-/*
-  #include "lapack.h"
-*/
+#include "cos.h"
+#include "sinh.h"
+#include "sin.h"
+#include "log.h"
+#include "abs.h"
+
+#define localSign(x) x >= 0 ? 1.0f : -1.0f
 
 floatComplex		ctans(floatComplex z) {
-  /*
+  float Temp = 0.0f;
+  float Lim = 1.0f + slogs(2.0f / ssqrts((float) getRelativeMachinePrecision()));
+  float RealIn = creals(z);
+  float ImagIn = cimags(z);
+  float RealOut = 0.0f;
+  float ImagOut = 0.0f;
 
-    double out_real = 0;
-    double out_imag = 0;
-    double L = 1 + dlogs(2 / dsqrts(dlamch('e', 1)));
-  */
+  Temp = powf(scoss(RealIn), 2.0f) + powf(ssinhs(ImagIn), 2.0f);
+  RealOut = 0.5f * ssins(2 * RealIn) / Temp;
+  if(sabss(ImagIn) < Lim)
+    {
+      ImagOut = 0.5f * ssinhs(2.0f * ImagIn) / Temp;
+    }
+  else
+    {
+      ImagOut = localSign(ImagIn);
+    }
 
-  return z;
+  return FloatComplex(RealOut, ImagOut);
 }