summaryrefslogtreecommitdiff
path: root/2.3-1/src/c/signalProcessing/modk
diff options
context:
space:
mode:
authorSandeep Gupta2017-07-05 12:41:25 +0530
committerSandeep Gupta2017-07-05 12:41:25 +0530
commitdd50e95a8193fb0faa846ccaa971a115ba69e71c (patch)
tree8d66abedd68dc5255ad323b95a16c9b592657f35 /2.3-1/src/c/signalProcessing/modk
parent3308ea7a973e2c1e1c8bea99ac7cc783ce8e8b06 (diff)
downloadScilab2C-dd50e95a8193fb0faa846ccaa971a115ba69e71c.tar.gz
Scilab2C-dd50e95a8193fb0faa846ccaa971a115ba69e71c.tar.bz2
Scilab2C-dd50e95a8193fb0faa846ccaa971a115ba69e71c.zip
LinearAlgebra and MatrixOperation Update
Diffstat (limited to '2.3-1/src/c/signalProcessing/modk')
-rw-r--r--2.3-1/src/c/signalProcessing/modk/dmodka.c97
1 files changed, 0 insertions, 97 deletions
diff --git a/2.3-1/src/c/signalProcessing/modk/dmodka.c b/2.3-1/src/c/signalProcessing/modk/dmodka.c
deleted file mode 100644
index c0630ec0..00000000
--- a/2.3-1/src/c/signalProcessing/modk/dmodka.c
+++ /dev/null
@@ -1,97 +0,0 @@
-#include<stdio.h>
-#include<math.h>
-#include "modk.h"
-double max_calc(double* ptr,int sz)
-{
- int i=0;
- double mx;
- if(ptr[0]<0)
- {
- ptr[0]=-1*ptr[0];
- }
- mx=(ptr[0]);
- //printf("%lf\n",mx);
- for(i=1;i<sz;i++)
- {
- if(ptr[i]<0)
- {
- ptr[i]=-1*ptr[i];
- }
- if(mx<(ptr[i]))
- mx=(ptr[i]);
- }
- return mx;
-}
-void dmodka(double* inp,int size,double* oup)
-{
- double ones[size],PI=M_PI;
- double eps=pow(2,-52);
- int i;
- for(i=0;i<size;i++)
- {
- ones[i]=1;
- }
- double a[size],b[size],c[size],an[size],bn[size],cn[size],kans[size];
-
- int j,kk,l,m;
- for(j=0;j<size;j++)
- {
- a[j]=1;
- }
- for(kk=0;kk<size;kk++)
- {
- b[kk]=sqrt(ones[kk]-inp[kk]);
- }
- for(l=0;l<size;l++)
- {
- c[l]=sqrt(inp[l]);
-
- }
- int x=0;
- //double maxi;
- //maxi=max_calc(c,size);
- //printf("%lf",maxi);
-
- while(max_calc(c,size)>eps)
- {
- int q,w,r;
- for(q=0;q<size;q++)
- {
- an[q]=0.5*(a[q]+b[q]);
- }
- for(w=0;w<size;w++)
- {
- bn[w]=sqrt(a[w]*b[w]);
- }
- for(r=0;r<size;r++)
- {
- cn[r]=0.5*(a[r]-b[r]);
- }
- int x,y,z;
- for(x=0;x<size;x++)
- {
- a[x]=an[x];
- }
- for(y=0;y<size;y++)
- {
- b[y]=bn[y];
- }
- for(z=0;z<size;z++)
- {
- c[z]=cn[z];
- }
- }
- int q,w;
- for(q=0;q<size;q++)
- {
- oup[q]=PI*(ones[q]/(2*a[q]));
- }
-
-}
-/*
-int main()
-{
- double m[3]={0.1,0.2,0.3};
- dka(m,3);
-}
-*/