From 1dc99d862fe3178a3a6e72e2c1338000545ec859 Mon Sep 17 00:00:00 2001
From: Ankitr19
Date: Fri, 7 Jul 2017 12:47:56 +0530
Subject: Demos and RPi compiler compatibility

---
 src/c/signalProcessing/ell1mag/zell1maga.c | 11 ++++++++++-
 1 file changed, 10 insertions(+), 1 deletion(-)

(limited to 'src/c/signalProcessing/ell1mag/zell1maga.c')

diff --git a/src/c/signalProcessing/ell1mag/zell1maga.c b/src/c/signalProcessing/ell1mag/zell1maga.c
index 6e7a6f9..35b4928 100644
--- a/src/c/signalProcessing/ell1mag/zell1maga.c
+++ b/src/c/signalProcessing/ell1mag/zell1maga.c
@@ -14,6 +14,9 @@
 #include "modsn.h"
 #include "ell1mag.h"
 #include "doubleComplex.h"
+#include "multiplication.h"
+#include "addition.h"
+#include "division.h"
 
 void zell1maga(double eps,double m1,doubleComplex* z,int size,double* oup )
 {
@@ -31,9 +34,15 @@ void zell1maga(double eps,double m1,doubleComplex* z,int size,double* oup )
 	}
 	doubleComplex v;
 	int k;
+	double ml;
+	ml=eps*eps;
+	doubleComplex tp1,tp2,tp3;
 	for(k=0;k<size;k++)
 	{
-	v=un[k]/(un[k]+(eps*eps*s[k]*s[k]));
+	tp1=zmuls(s[k],s[k]);
+	tp2=DoubleComplex(ml,0);
+	tp3=DoubleComplex(un[k],0);
+	v=zrdivs(tp3,zadds(tp3,(zmuls(tp2,tp1))));
 	oup[k]=zreals(v);
 	}
 }
-- 
cgit