diff options
| author | Abhinav Dronamraju | 2017-09-29 22:00:40 +0530 |
|---|---|---|
| committer | Abhinav Dronamraju | 2017-09-29 22:00:40 +0530 |
| commit | 9bc7ad78e8d7d7acc4b9387aa592542832e80b31 (patch) | |
| tree | 7fce060665a91de5e5adb12d02003351c3d1fdfc /src/c/signalProcessing/ell1mag/zell1maga.c | |
| parent | 33755eb085a3ca8154cf83773b23fbb8aac4ba3e (diff) | |
| parent | ac0045f12ad3d0938758e9742f4107a334e1afaa (diff) | |
| download | scilab2c-9bc7ad78e8d7d7acc4b9387aa592542832e80b31.tar.gz scilab2c-9bc7ad78e8d7d7acc4b9387aa592542832e80b31.tar.bz2 scilab2c-9bc7ad78e8d7d7acc4b9387aa592542832e80b31.zip | |
NEW FEATURES AND NEW FUNCTIONS
Diffstat (limited to 'src/c/signalProcessing/ell1mag/zell1maga.c')
| -rw-r--r-- | src/c/signalProcessing/ell1mag/zell1maga.c | 49 |
1 files changed, 49 insertions, 0 deletions
diff --git a/src/c/signalProcessing/ell1mag/zell1maga.c b/src/c/signalProcessing/ell1mag/zell1maga.c new file mode 100644 index 00000000..35b49281 --- /dev/null +++ b/src/c/signalProcessing/ell1mag/zell1maga.c @@ -0,0 +1,49 @@ +/* Copyright (C) 2017 - IIT Bombay - FOSSEE + + This file must be used under the terms of the CeCILL. + This source file is licensed as described in the file COPYING, which + you should have received as part of this distribution. The terms + are also available at + http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt + Author: Ankit Raj + Organization: FOSSEE, IIT Bombay + Email: toolbox@scilab.in + */ + +#include<stdio.h> +#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 ) +{ + doubleComplex s[size]; + int i; + for(i=0;i<size;i++) + { + s[i]=zmodsns(z[i],m1); + } + double un[size]; + int j; + for(j=0;j<size;j++) + { + un[j]=1; + } + doubleComplex v; + int k; + double ml; + ml=eps*eps; + doubleComplex tp1,tp2,tp3; + for(k=0;k<size;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); + } +} + |