diff options
Diffstat (limited to '2.3-1/src/c/signalProcessing/transforms/dct/zdcta.c')
| -rw-r--r-- | 2.3-1/src/c/signalProcessing/transforms/dct/zdcta.c | 230 |
1 files changed, 0 insertions, 230 deletions
diff --git a/2.3-1/src/c/signalProcessing/transforms/dct/zdcta.c b/2.3-1/src/c/signalProcessing/transforms/dct/zdcta.c deleted file mode 100644 index 3ae2e333..00000000 --- a/2.3-1/src/c/signalProcessing/transforms/dct/zdcta.c +++ /dev/null @@ -1,230 +0,0 @@ -/* Copyright (C) 2016 - 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: Ukasha Noor - Organization: FOSSEE, IIT Bombay - Email: toolbox@scilab.in -*/ - - -#include "dct.h" -#include "addition.h" -#include "types.h" -#include "doubleComplex.h" -#include "multiplication.h" -/*#include "matrixMultiplication"*/ -/*#include <fftw3.h>*/ -#include <math.h> - -void zdcta(doubleComplex *in,int row,int col,int sign,doubleComplex *out) -{ - int i,j,k,u,v; - int n; - int x,y; - double res,ress,vv,ff; - double re,z,q,m; - doubleComplex accu = DoubleComplex(0, 0); - doubleComplex temp,mm,aa,bb,cc; - if(sign==-1) - { - if(row==1) - { - n=col; - for(u=0;u<row;u++) - { - for(v=0;v<col;v++) - { - x=v*row+u; - out[x]=DoubleComplex(0,0); - for(i=0;i<row;i++) - { - for(j=0;j<col;j++) - { - y=row*j+i; - vv = cos(((M_PI)*(y+1-1./2.)*(x))/n); - aa = DoubleComplex(vv,0); - temp=zmuls(in[y],aa); - out[x]=zadds(out[x],temp); - } - } - if(x==0) - { - vv = 1./(sqrt(n)); - aa = DoubleComplex(vv,0); - out[x]=zmuls(out[x],aa); - } - else - { - double res=2./n; - res = sqrt(res); - aa = DoubleComplex(res,0); - out[x]=zmuls(out[x],aa); - } - } - } - } - else - { - n=col*row; - for(u=0;u<row;u++) - { - for(v=0;v<col;v++) - { - x=v*row+u; - out[x]=DoubleComplex(0,0); - for(i=0;i<row;i++) - { - temp=DoubleComplex(0,0); - mm=DoubleComplex(0,0); - for(j=0;j<col;j++) - { - y=j*row+i; - z=(double )(((double )j+1.0/2.0)*(double )v); - q=(double )(M_PI/(double )col); - vv = cos(q*z); - aa = DoubleComplex(vv,0); - mm=zmuls(in[y],aa); - temp=zadds(temp,mm); - } - z=(double )(((double )i+1.0/2.0)*(double )u); - q=(double )(M_PI/(double )row); - ff = cos(q*z); - bb = DoubleComplex(ff,0); - temp=zmuls(temp,bb); - out[x]=zadds(out[x],temp); - } - if(u==0) - { - vv = 1./sqrt((double )row); - aa = DoubleComplex(vv,0); - out[x]=zmuls(out[x],aa); - if(v==0) - { - vv = 1./sqrt((double )col); - aa = DoubleComplex(vv,0); - out[x]=zmuls(out[x],aa); - } - else - { - vv = sqrt(2./col); - aa = DoubleComplex(vv,0); - out[x]=zmuls(out[x],aa); - } - } - else - { - vv = sqrt(2./row); - aa = DoubleComplex(vv,0); - out[x]=zmuls(out[x],aa); - if(v==0) - { - vv = 1./sqrt((double )col); - aa = DoubleComplex(vv,0); - out[x]=zmuls(out[x],aa); - } - else - { - vv = sqrt(2./col); - aa = DoubleComplex(vv,0); - out[x]=zmuls(out[x],aa); - } - } - } - } - } - } - else if(sign==1) - { - n=col; - if(row==1) - { - res=1./sqrt(n); - ress=sqrt(2./n); - for(u=0;u<row;u++) - { - for(v=0;v<col;v++) - { - x=v*row+u; - out[x]=DoubleComplex(0,0); - for(i=0;i<row;i++) - { - for(j=0;j<col;j++) - { - y=row*j+i; - if(y==0) - { - q=res*(cos(((M_PI)*(j)*(v+1./2.))/n)); - aa = DoubleComplex(q,0); - out[x]=zadds(out[x],zmuls(in[y],aa)); - } - else - { - q=ress*(cos(((M_PI)*(j)*(v+1./2.))/n)); - aa = DoubleComplex(q,0); - out[x]=zadds(out[x],zmuls(in[y],aa)); - } - } - } - } - - } - } - else - { - for(u=0;u<row;u++) - { - for(v=0;v<col;v++) - { - x=v*row+u; - out[x]=DoubleComplex(0,0); - for(i=0;i<row;i++) - { - re=0; - mm=DoubleComplex(0,0); - temp=DoubleComplex(0,0); - for(j=0;j<col;j++) - { - y=row*j+i; - mm=in[j*row+i]; - z=(double )(((double )v+1.0/2.0)*(double )j); - q=(double )(M_PI/(double )col); - vv = cos(q*z); - aa = DoubleComplex(vv,0); - mm=zmuls(mm,aa); - if(j==0) - { - vv = 1./sqrt((double )col); - aa = DoubleComplex(vv,0); - temp=zadds(temp,zmuls(mm,aa)); - } - else - { - vv = sqrt(2./col); - aa = DoubleComplex(vv,0); - temp=zadds(temp,zmuls(mm,aa)); - } - } - z=(double )(((double )u+1.0/2.0)*(double )i); - q=(double )(M_PI/(double )row); - if(i==0) - { - vv = (cos(z*q))*(1./sqrt(row)); - aa = DoubleComplex(vv,0); - out[x]=zadds(out[x],zmuls(temp,aa)); - } - else - { - vv = (cos(z*q))*sqrt(2./row); - aa = DoubleComplex(vv,0); - out[x]=zadds(out[x],zmuls(temp,aa)); - } - } - } - } - } - } -} |