diff options
| author | Abhinav Dronamraju | 2017-07-10 22:08:27 +0530 |
|---|---|---|
| committer | Abhinav Dronamraju | 2017-07-10 22:08:27 +0530 |
| commit | f80b5800ddb8417530f68049f6895beccb89cadf (patch) | |
| tree | 564a23878bc400ecc0873bb4b1f5b172af22cfc6 /2.3-1/src/c/signalProcessing/ell1mag | |
| parent | c1874d367d68cb47b82cc7a1173caaf38e9e3d68 (diff) | |
| parent | 2e9eadc75e51483c41a30544496b4e05c0cf7c41 (diff) | |
| download | Scilab2C-f80b5800ddb8417530f68049f6895beccb89cadf.tar.gz Scilab2C-f80b5800ddb8417530f68049f6895beccb89cadf.tar.bz2 Scilab2C-f80b5800ddb8417530f68049f6895beccb89cadf.zip | |
Pulled from upstream master
Diffstat (limited to '2.3-1/src/c/signalProcessing/ell1mag')
| -rw-r--r-- | 2.3-1/src/c/signalProcessing/ell1mag/dell1maga.c | 2 | ||||
| -rw-r--r-- | 2.3-1/src/c/signalProcessing/ell1mag/zell1maga.c | 11 |
2 files changed, 11 insertions, 2 deletions
diff --git a/2.3-1/src/c/signalProcessing/ell1mag/dell1maga.c b/2.3-1/src/c/signalProcessing/ell1mag/dell1maga.c index 9af0c8e8..58ef3035 100644 --- a/2.3-1/src/c/signalProcessing/ell1mag/dell1maga.c +++ b/2.3-1/src/c/signalProcessing/ell1mag/dell1maga.c @@ -21,7 +21,7 @@ void dell1maga(double eps,double m1,double* z,int size,double* oup ) int i; for(i=0;i<size;i++) { - s[i]=zmodsns(z[i],m1); + s[i]=dmodsns(z[i],m1); } double un[size]; int j; diff --git a/2.3-1/src/c/signalProcessing/ell1mag/zell1maga.c b/2.3-1/src/c/signalProcessing/ell1mag/zell1maga.c index 6e7a6f93..35b49281 100644 --- a/2.3-1/src/c/signalProcessing/ell1mag/zell1maga.c +++ b/2.3-1/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); } } |