@@ -1,282 +1,282 | |||||
1 | #include <stdio.h> |
|
1 | #include <stdio.h> | |
2 | #include <math.h> |
|
2 | #include <math.h> | |
3 | #include <string.h> |
|
3 | #include <string.h> | |
4 | #include <stdlib.h> |
|
4 | #include <stdlib.h> | |
5 | #include <fcntl.h> |
|
5 | #include <fcntl.h> | |
6 | #include <time.h> |
|
6 | #include <time.h> | |
7 | #include "complex.h" |
|
7 | #include "complex.h" | |
8 |
|
8 | |||
9 | #define NFREQ 512 |
|
9 | #define NFREQ 512 | |
10 | #define LNES 12 |
|
10 | #define LNES 12 | |
11 | #define LTES 16 |
|
11 | #define LTES 16 | |
12 | #define NANG 19 |
|
12 | #define NANG 19 | |
13 |
|
13 | |||
14 | #define MAX(x,y) ((x)>(y) ? (x) : (y)) |
|
14 | #define MAX(x,y) ((x)>(y) ? (x) : (y)) | |
15 | #define MIN(x,y) ((x)<(y) ? (x) : (y)) |
|
15 | #define MIN(x,y) ((x)<(y) ? (x) : (y)) | |
16 |
|
16 | |||
17 | fcomplex ****exlib; |
|
17 | fcomplex ****exlib; | |
18 | int first=1; |
|
18 | int first=1; | |
19 |
|
19 | |||
20 | void swab_test(void); |
|
20 | void swab_test(void); | |
21 | //float swab(signed char *a); |
|
21 | //float swab(signed char *a); | |
22 | void read_exlib(char *s,int nfreq, int lnes, int lTes, int nang, |
|
22 | void read_exlib(char *s,int nfreq, int lnes, int lTes, int nang, | |
23 | fcomplex ****exlib); |
|
23 | fcomplex ****exlib); | |
24 |
|
24 | |||
25 | void initialize(void); |
|
25 | void initialize(void); | |
26 | //void collision_(float *np, float *tp, float *fp, float *ap, float *yr, float *yi)// step 1 |
|
26 | //void collision_(float *np, float *tp, float *fp, float *ap, float *yr, float *yi)// step 1 | |
27 | void collision_(float *n, float *t, float *f, float *a, fcomplex *y); |
|
27 | void collision_(float *n, float *t, float *f, float *a, fcomplex *y); | |
28 |
|
28 | |||
29 | /* |
|
29 | /* | |
30 | void main(void){ |
|
30 | void main(void){ | |
31 | float n=4.9e11,t=1100.0,f=300.0,a=1.0; |
|
31 | float n=4.9e11,t=1100.0,f=300.0,a=1.0; | |
32 | int i; |
|
32 | int i; | |
33 | fcomplex j; |
|
33 | fcomplex j; | |
34 | for(i=0;i<10;i++){ |
|
34 | for(i=0;i<10;i++){ | |
35 | a+=0.05; |
|
35 | a+=0.05; | |
36 | collision_(&n,&t,&f,&a,&j); |
|
36 | collision_(&n,&t,&f,&a,&j); | |
37 | printf("%f %f %f \n",a,j.r,j.i); |
|
37 | printf("%f %f %f \n",a,j.r,j.i); | |
38 | } |
|
38 | } | |
39 | } |
|
39 | } | |
40 | */ |
|
40 | */ | |
41 |
|
41 | |||
42 | //void collision_(float *np, float *tp, float *fp, float *ap, float *yr, float *yi){ // step 2 |
|
42 | //void collision_(float *np, float *tp, float *fp, float *ap, float *yr, float *yi){ // step 2 | |
43 | void collision_(float *np, float *tp, float *fp, float *ap, fcomplex *y){ |
|
43 | void collision_(float *np, float *tp, float *fp, float *ap, fcomplex *y){ | |
44 | float ne[]={5.0e10, 1.0e11, 2.0e11, 3.5e11, 5.0e11, 7.5e11, |
|
44 | float ne[]={5.0e10, 1.0e11, 2.0e11, 3.5e11, 5.0e11, 7.5e11, | |
45 | 1.0e12, 1.5e12, 2.0e12, 2.5e12, 3.0e12, 3.5e12}; |
|
45 | 1.0e12, 1.5e12, 2.0e12, 2.5e12, 3.0e12, 3.5e12}; | |
46 | float te[]={600.0, 700.0, 800.0, 900.0, 1000.0, 1150.0, 1300.0, 1500.0, |
|
46 | float te[]={600.0, 700.0, 800.0, 900.0, 1000.0, 1150.0, 1300.0, 1500.0, | |
47 | 1750.0, 2000.0, 2250.0, 2500.0, 2750.0, 3000.0, 3500.0, 4000.0}; |
|
47 | 1750.0, 2000.0, 2250.0, 2500.0, 2750.0, 3000.0, 3500.0, 4000.0}; | |
48 | float alpha[]={.125, .1875, .25, .3125, .375, .5, .625, .75, 1.0, 1.25, |
|
48 | float alpha[]={.125, .1875, .25, .3125, .375, .5, .625, .75, 1.0, 1.25, | |
49 | 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 6.0}; |
|
49 | 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 6.0}; | |
50 | float freq[2],fmax=6237.79; |
|
50 | float freq[2],fmax=6237.79; | |
51 | float n,t,f,a; |
|
51 | float n,t,f,a; | |
52 | float v,w1,w2,w3,w4,weight; |
|
52 | float v,w1,w2,w3,w4,weight; | |
53 | fcomplex g; |
|
53 | fcomplex g; | |
54 | int i,j,k,l,m,i1,j1,k1,l1; |
|
54 | int i,j,k,l,m,i1,j1,k1,l1; | |
55 |
|
55 | |||
56 | if(first==1){ |
|
56 | if(first==1){ | |
57 | //printf("initi"); |
|
57 | //printf("initi"); | |
58 | initialize(); |
|
58 | initialize(); | |
59 | first=0; |
|
59 | first=0; | |
60 | } |
|
60 | } | |
61 | //printf(" inputs in collision_ \n y.r %f y.i %f \n *************\n",(*y).r,(*y).i); |
|
61 | //printf(" inputs in collision_ \n y.r %f y.i %f \n *************\n",(*y).r,(*y).i); | |
62 |
|
62 | |||
63 | //printf("After:print\n"); |
|
63 | //printf("After:print\n"); | |
64 | a=*ap; t=*tp; f=*fp; n=*np; |
|
64 | a=*ap; t=*tp; f=*fp; n=*np; | |
65 |
|
65 | |||
66 | a=MAX(MIN(a,6.0),0.125); |
|
66 | a=MAX(MIN(a,6.0),0.125); | |
67 | t=MAX(MIN(t,4000.0),600.0); |
|
67 | t=MAX(MIN(t,4000.0),600.0); | |
68 | f=MIN(f,fmax-1.0); |
|
68 | f=MIN(f,fmax-1.0); | |
69 | n=MAX(MIN(n,3.5e12),5.0e10); |
|
69 | n=MAX(MIN(n,3.5e12),5.0e10); | |
70 |
|
70 | |||
71 | i1=((float)(NFREQ-1))*f/fmax; |
|
71 | i1=((float)(NFREQ-1))*f/fmax; | |
72 | i1=MIN(i1,NFREQ-1); |
|
72 | i1=MIN(i1,NFREQ-1); | |
73 | freq[0]=i1*fmax/(float)(NFREQ-1); |
|
73 | freq[0]=i1*fmax/(float)(NFREQ-1); | |
74 | freq[1]=(i1+1)*fmax/(float)(NFREQ-1); |
|
74 | freq[1]=(i1+1)*fmax/(float)(NFREQ-1); | |
75 |
|
75 | |||
76 | j1=0; |
|
76 | j1=0; | |
77 | while(n>ne[j1+1]) j1++; |
|
77 | while(n>ne[j1+1]) j1++; | |
78 | j1=MIN(j1,LNES-1); |
|
78 | j1=MIN(j1,LNES-1); | |
79 |
|
79 | |||
80 | k1=0; |
|
80 | k1=0; | |
81 | while(t>te[k1+1]) k1++; |
|
81 | while(t>te[k1+1]) k1++; | |
82 | k1=MIN(k1,LTES); |
|
82 | k1=MIN(k1,LTES); | |
83 |
|
83 | |||
84 | l1=0; |
|
84 | l1=0; | |
85 | while(a>alpha[l1+1]) l1++; |
|
85 | while(a>alpha[l1+1]) l1++; | |
86 | l1=MIN(l1,NANG); |
|
86 | l1=MIN(l1,NANG); | |
87 |
|
87 | |||
88 | // printf("i,j,k,l: %i %i %i %i \n",i1,j1,k1,l1); |
|
88 | // printf("i,j,k,l: %i %i %i %i \n",i1,j1,k1,l1); | |
89 |
|
89 | |||
90 | v=(fmax/(float)(NFREQ-1))* |
|
90 | v=(fmax/(float)(NFREQ-1))* | |
91 | (ne[j1+1]-ne[j1])*(te[k1+1]-te[k1])*(alpha[l1+1]-alpha[l1]); |
|
91 | (ne[j1+1]-ne[j1])*(te[k1+1]-te[k1])*(alpha[l1+1]-alpha[l1]); | |
92 | g=Complex(0.0,0.0); |
|
92 | g=Complex(0.0,0.0); | |
93 |
|
93 | |||
94 | for(i=0;i<2;i++){ |
|
94 | for(i=0;i<2;i++){ | |
95 | if(i==0) |
|
95 | if(i==0) | |
96 | w1=freq[1]-f; |
|
96 | w1=freq[1]-f; | |
97 | else |
|
97 | else | |
98 | w1=f-freq[0]; |
|
98 | w1=f-freq[0]; | |
99 | for(j=0;j<2;j++){ |
|
99 | for(j=0;j<2;j++){ | |
100 | if(j==0) |
|
100 | if(j==0) | |
101 | w2=ne[j1+1]-n; |
|
101 | w2=ne[j1+1]-n; | |
102 | else |
|
102 | else | |
103 | w2=n-ne[j1]; |
|
103 | w2=n-ne[j1]; | |
104 | for(k=0;k<2;k++){ |
|
104 | for(k=0;k<2;k++){ | |
105 | if(k==0) |
|
105 | if(k==0) | |
106 | w3=te[k1+1]-t; |
|
106 | w3=te[k1+1]-t; | |
107 | else |
|
107 | else | |
108 | w3=t-te[k1]; |
|
108 | w3=t-te[k1]; | |
109 | for(l=0;l<2;l++){ |
|
109 | for(l=0;l<2;l++){ | |
110 | if(l==0) |
|
110 | if(l==0) | |
111 | w4=alpha[l1+1]-a; |
|
111 | w4=alpha[l1+1]-a; | |
112 | else |
|
112 | else | |
113 | w4=a-alpha[l1]; |
|
113 | w4=a-alpha[l1]; | |
114 | weight=w1*w2*w3*w4; |
|
114 | weight=w1*w2*w3*w4; | |
115 | g.r+=exlib[i1+i][j1+j][k1+k][l1+l].r*weight; |
|
115 | g.r+=exlib[i1+i][j1+j][k1+k][l1+l].r*weight; | |
116 | g.i+=exlib[i1+i][j1+j][k1+k][l1+l].i*weight; |
|
116 | g.i+=exlib[i1+i][j1+j][k1+k][l1+l].i*weight; | |
117 | // printf(" in collision i1+i %d j1+j %d k1+k %d l1+l %d exlib.r %f exlib.i %f \n",i1+i,j1+j,k1+k,l1+l,exlib[i1+i][j1+j][k1+k][l1+l].r,exlib[i1+i][j1+j][k1+k][l1+l].i); |
|
117 | // printf(" in collision i1+i %d j1+j %d k1+k %d l1+l %d exlib.r %f exlib.i %f \n",i1+i,j1+j,k1+k,l1+l,exlib[i1+i][j1+j][k1+k][l1+l].r,exlib[i1+i][j1+j][k1+k][l1+l].i); | |
118 | // printf(" in collision g.r %f g.i %f \n",g.r,g.i); |
|
118 | // printf(" in collision g.r %f g.i %f \n",g.r,g.i); | |
119 | } |
|
119 | } | |
120 | } |
|
120 | } | |
121 | } |
|
121 | } | |
122 | } |
|
122 | } | |
123 | g.r/=v; g.i/=v; |
|
123 | g.r/=v; g.i/=v; | |
124 |
|
124 | |||
125 | *y=g; |
|
125 | *y=g; | |
126 | //*yr=g.r |
|
126 | //*yr=g.r | |
127 | //*yi=g.i |
|
127 | //*yi=g.i | |
128 | //printf(" outputs in collision_ \n g.r %f g.i %f \n *************\n",g.r,g.i); |
|
128 | //printf(" outputs in collision_ \n g.r %f g.i %f \n *************\n",g.r,g.i); | |
129 | //exit(-2); |
|
129 | //exit(-2); | |
130 | //getchar(); |
|
130 | //getchar(); | |
131 | return; |
|
131 | return; | |
132 | } |
|
132 | } | |
133 |
|
133 | |||
134 | void initialize(void){ |
|
134 | void initialize(void){ | |
135 | // char s[]="/usr/local/lib/faraday/jlib26feb2002"; |
|
135 | // char s[]="/usr/local/lib/faraday/jlib26feb2002"; | |
136 | // char s[]="/usr/local/lib/faraday/jlib26feb2002"; |
|
136 | // char s[]="/usr/local/lib/faraday/jlib26feb2002"; | |
137 | // FILE *s=fopen("/usr/local/lib/faraday/jlib26feb2002","r"); |
|
137 | // FILE *s=fopen("/usr/local/lib/faraday/jlib26feb2002","r"); | |
138 | int nfreq=NFREQ, lnes=LNES, lTes=LTES, nang=NANG; |
|
138 | int nfreq=NFREQ, lnes=LNES, lTes=LTES, nang=NANG; | |
139 | int i,j,k; |
|
139 | int i,j,k; | |
140 | int crr,crr2; |
|
140 | int crr,crr2; | |
141 | int the_number=1025; |
|
141 | int the_number=1025; | |
142 |
|
142 | |||
143 | char the_path[1025]; |
|
143 | char the_path[1025]; | |
144 |
|
144 | |||
145 | get_path_(the_path,&the_number); |
|
145 | get_path_reader_(the_path,&the_number); | |
146 | //printf("C the_number: %i",the_number); |
|
146 | //printf("C the_number: %i",the_number); | |
147 | //printf("After\n"); |
|
147 | //printf("After\n"); | |
148 | //printf("The Path: %s", the_path); |
|
148 | //printf("The Path: %s", the_path); | |
149 | //printf("END\n"); |
|
149 | //printf("END\n"); | |
150 | //printf("%s%n",the_path,&crr); |
|
150 | //printf("%s%n",the_path,&crr); | |
151 | //printf("=%d",crr); |
|
151 | //printf("=%d",crr); | |
152 | //printf("END2\n"); |
|
152 | //printf("END2\n"); | |
153 | char the_path_true[the_number]; |
|
153 | char the_path_true[the_number]; | |
154 | memcpy(the_path_true,the_path,the_number); |
|
154 | memcpy(the_path_true,the_path,the_number); | |
155 | the_path_true[the_number] = '\0'; |
|
155 | the_path_true[the_number] = '\0'; | |
156 | //printf("%s%n",the_path_true,&crr2); |
|
156 | //printf("%s%n",the_path_true,&crr2); | |
157 | //printf("=%d",crr2); |
|
157 | //printf("=%d",crr2); | |
158 | //printf("END3\n"); |
|
158 | //printf("END3\n"); | |
159 | //printf("PATH: %s",s); |
|
159 | //printf("PATH: %s",s); | |
160 | //fprintf(s); |
|
160 | //fprintf(s); | |
161 | exlib=malloc(nfreq*8); |
|
161 | exlib=malloc(nfreq*8); | |
162 | for(i=0;i<nfreq;i++){ |
|
162 | for(i=0;i<nfreq;i++){ | |
163 | exlib[i]=malloc(lnes*8); |
|
163 | exlib[i]=malloc(lnes*8); | |
164 | for(j=0;j<lnes;j++){ |
|
164 | for(j=0;j<lnes;j++){ | |
165 | exlib[i][j]=malloc(lTes*8); |
|
165 | exlib[i][j]=malloc(lTes*8); | |
166 | for(k=0;k<lTes;k++){ |
|
166 | for(k=0;k<lTes;k++){ | |
167 | exlib[i][j][k]=(fcomplex *)malloc(nang*sizeof(fcomplex)); |
|
167 | exlib[i][j][k]=(fcomplex *)malloc(nang*sizeof(fcomplex)); | |
168 | } |
|
168 | } | |
169 | } |
|
169 | } | |
170 | } |
|
170 | } | |
171 |
|
171 | |||
172 | read_exlib(the_path_true,nfreq,lnes,lTes,nang,exlib); |
|
172 | read_exlib(the_path_true,nfreq,lnes,lTes,nang,exlib); | |
173 | } |
|
173 | } | |
174 |
|
174 | |||
175 | float swab(signed char *a) |
|
175 | float swab(signed char *a) | |
176 | // table is in big endian format |
|
176 | // table is in big endian format | |
177 | { |
|
177 | { | |
178 | float mant; |
|
178 | float mant; | |
179 | unsigned char *amant; |
|
179 | unsigned char *amant; | |
180 | printf(" in swab 0 *** a[0] %c a[1] %c a[2] %c a[3] %c\n",a[0],a[1],a[2],a[3]); |
|
180 | printf(" in swab 0 *** a[0] %c a[1] %c a[2] %c a[3] %c\n",a[0],a[1],a[2],a[3]); | |
181 | amant=(unsigned char *) &mant; |
|
181 | amant=(unsigned char *) &mant; | |
182 | amant[3]=a[0]; amant[2]=a[1]; amant[1]=a[2]; amant[0]=a[3]; |
|
182 | amant[3]=a[0]; amant[2]=a[1]; amant[1]=a[2]; amant[0]=a[3]; | |
183 | //printf(" in swab 1 ***\n"); |
|
183 | //printf(" in swab 1 ***\n"); | |
184 | return(mant); |
|
184 | return(mant); | |
185 | } |
|
185 | } | |
186 |
|
186 | |||
187 | void swab_test(void) |
|
187 | void swab_test(void) | |
188 | // table is in big endian format |
|
188 | // table is in big endian format | |
189 | { |
|
189 | { | |
190 | printf("*********in swab stest ********\n"); |
|
190 | printf("*********in swab stest ********\n"); | |
191 | getchar(); |
|
191 | getchar(); | |
192 | } |
|
192 | } | |
193 |
|
193 | |||
194 |
|
194 | |||
195 | void read_exlib(char *the_path_true,int nfreq, int lnes, int lTes, int nang, |
|
195 | void read_exlib(char *the_path_true,int nfreq, int lnes, int lTes, int nang, | |
196 | fcomplex ****exlib) |
|
196 | fcomplex ****exlib) | |
197 | { |
|
197 | { | |
198 | FILE *fp; |
|
198 | FILE *fp; | |
199 | int i,j,k,l; |
|
199 | int i,j,k,l; | |
200 | float a; |
|
200 | float a; | |
201 | float t1, t2; |
|
201 | float t1, t2; | |
202 | signed char *t1a,*t2b; |
|
202 | signed char *t1a,*t2b; | |
203 | float mant; |
|
203 | float mant; | |
204 | unsigned char *amant; |
|
204 | unsigned char *amant; | |
205 | printf("Start reading exlib file \n"); |
|
205 | printf("Start reading exlib file \n"); | |
206 | if( (fp=fopen(the_path_true,"r")) == NULL){ |
|
206 | if( (fp=fopen(the_path_true,"r")) == NULL){ | |
207 | printf("Je Library file opening error\n"); |
|
207 | printf("Je Library file opening error\n"); | |
208 | //printf("%s\n", strerror(errno)); |
|
208 | //printf("%s\n", strerror(errno)); | |
209 | exit(-2); |
|
209 | exit(-2); | |
210 | } |
|
210 | } | |
211 | for(l=0;l<nang;l++) |
|
211 | for(l=0;l<nang;l++) | |
212 | for(k=0;k<lTes;k++) |
|
212 | for(k=0;k<lTes;k++) | |
213 | for(j=0;j<lnes;j++) |
|
213 | for(j=0;j<lnes;j++) | |
214 | for(i=0;i<nfreq;i++){ |
|
214 | for(i=0;i<nfreq;i++){ | |
215 | //printf(" in read_exlib 0***\n"); |
|
215 | //printf(" in read_exlib 0***\n"); | |
216 | fread(&t1,sizeof(float),1,fp); |
|
216 | fread(&t1,sizeof(float),1,fp); | |
217 | fread(&t2,sizeof(float),1,fp); |
|
217 | fread(&t2,sizeof(float),1,fp); | |
218 | //printf(" in read_exlib 1***\n"); |
|
218 | //printf(" in read_exlib 1***\n"); | |
219 | //printf(" i %d j %d k %d l %d t1 %f t2 %f \n",i,j,k,l,t1,t2); |
|
219 | //printf(" i %d j %d k %d l %d t1 %f t2 %f \n",i,j,k,l,t1,t2); | |
220 | t1a=(signed char *)&t1; |
|
220 | t1a=(signed char *)&t1; | |
221 | t2b=(signed char *)&t2; |
|
221 | t2b=(signed char *)&t2; | |
222 | //printf(" in read_exlib 1.5 t1a %d t2b %d ***\n",t1a,t2b); |
|
222 | //printf(" in read_exlib 1.5 t1a %d t2b %d ***\n",t1a,t2b); | |
223 | ////t1=swab(t1a); |
|
223 | ////t1=swab(t1a); | |
224 | //printf(" t1 -- in readexlib 1.8 0 *** t1a[0] %c t1a[1] %c t1a[2] %c t1a[3] %c\n",t1a[0],t1a[1],t1a[2],t1a[3]); |
|
224 | //printf(" t1 -- in readexlib 1.8 0 *** t1a[0] %c t1a[1] %c t1a[2] %c t1a[3] %c\n",t1a[0],t1a[1],t1a[2],t1a[3]); | |
225 | amant=(unsigned char *) &mant; |
|
225 | amant=(unsigned char *) &mant; | |
226 | amant[3]=t1a[0]; amant[2]=t1a[1]; amant[1]=t1a[2]; amant[0]=t1a[3]; |
|
226 | amant[3]=t1a[0]; amant[2]=t1a[1]; amant[1]=t1a[2]; amant[0]=t1a[3]; | |
227 | t1=mant; |
|
227 | t1=mant; | |
228 |
|
228 | |||
229 | ////t2=swab(t2b); |
|
229 | ////t2=swab(t2b); | |
230 | //printf(" t2 -- in readexlib 1.8 0 *** t2b[0] %c t2b[1] %c t2b[2] %c t2b[3] %c\n",t2b[0],t2b[1],t2b[2],t2b[3]); |
|
230 | //printf(" t2 -- in readexlib 1.8 0 *** t2b[0] %c t2b[1] %c t2b[2] %c t2b[3] %c\n",t2b[0],t2b[1],t2b[2],t2b[3]); | |
231 | amant=(unsigned char *) &mant; |
|
231 | amant=(unsigned char *) &mant; | |
232 | amant[3]=t2b[0]; amant[2]=t2b[1]; amant[1]=t2b[2]; amant[0]=t2b[3]; |
|
232 | amant[3]=t2b[0]; amant[2]=t2b[1]; amant[1]=t2b[2]; amant[0]=t2b[3]; | |
233 | t2=mant; |
|
233 | t2=mant; | |
234 |
|
234 | |||
235 | ////t1=swab((signed char *)&t1); |
|
235 | ////t1=swab((signed char *)&t1); | |
236 | ////t2=swab((signed char *)&t2); |
|
236 | ////t2=swab((signed char *)&t2); | |
237 | //printf(" in read_exlib 2 i %d j %d k %d l %d t1 %f t2 %f \n",i,j,k,l,t1,t2); |
|
237 | //printf(" in read_exlib 2 i %d j %d k %d l %d t1 %f t2 %f \n",i,j,k,l,t1,t2); | |
238 | exlib[i][j][k][l]=Complex(t1,t2); |
|
238 | exlib[i][j][k][l]=Complex(t1,t2); | |
239 | //printf(" in read_exlib 3 i %d j %d k %d l %d t1 %f t2 %f \n",i,j,k,l,exlib[i][j][k][l].r,exlib[i][j][k][l].i); |
|
239 | //printf(" in read_exlib 3 i %d j %d k %d l %d t1 %f t2 %f \n",i,j,k,l,exlib[i][j][k][l].r,exlib[i][j][k][l].i); | |
240 | //getchar(); |
|
240 | //getchar(); | |
241 | } |
|
241 | } | |
242 | // getchar(); |
|
242 | // getchar(); | |
243 |
|
243 | |||
244 | fclose(fp); |
|
244 | fclose(fp); | |
245 | printf("Done reading exlib file \n"); |
|
245 | printf("Done reading exlib file \n"); | |
246 | } |
|
246 | } | |
247 | /* |
|
247 | /* | |
248 | void read_exlib(char *s,int nfreq, int lnes, int lTes, int nang, |
|
248 | void read_exlib(char *s,int nfreq, int lnes, int lTes, int nang, | |
249 | fcomplex ****exlib) |
|
249 | fcomplex ****exlib) | |
250 | { |
|
250 | { | |
251 | FILE *fp; |
|
251 | FILE *fp; | |
252 | int i,j,k,l; |
|
252 | int i,j,k,l; | |
253 | float a; |
|
253 | float a; | |
254 | float t1, t2; |
|
254 | float t1, t2; | |
255 |
|
255 | |||
256 | printf("Start reading exlib fileee \n"); |
|
256 | printf("Start reading exlib fileee \n"); | |
257 |
|
257 | |||
258 | if( (fp=fopen(s,"r")) == NULL){ |
|
258 | if( (fp=fopen(s,"r")) == NULL){ | |
259 | printf("Je Library file opening error\n"); |
|
259 | printf("Je Library file opening error\n"); | |
260 | exit(-2); |
|
260 | exit(-2); | |
261 | } |
|
261 | } | |
262 |
|
262 | |||
263 | for(l=0;l<nang;l++) |
|
263 | for(l=0;l<nang;l++) | |
264 | for(k=0;k<lTes;k++) |
|
264 | for(k=0;k<lTes;k++) | |
265 | for(j=0;j<lnes;j++) |
|
265 | for(j=0;j<lnes;j++) | |
266 | for(i=0;i<nfreq;i++){ |
|
266 | for(i=0;i<nfreq;i++){ | |
267 | printf(" in read_exlib 0***\n"); |
|
267 | printf(" in read_exlib 0***\n"); | |
268 | fread(&t1,sizeof(float),1,fp); |
|
268 | fread(&t1,sizeof(float),1,fp); | |
269 | fread(&t2,sizeof(float),1,fp); |
|
269 | fread(&t2,sizeof(float),1,fp); | |
270 | printf(" in read_exlib 1***\n"); |
|
270 | printf(" in read_exlib 1***\n"); | |
271 | printf(" i %d j %d k %d l %d t1 %f t2 %f \n",i,j,k,l,t1,t2); |
|
271 | printf(" i %d j %d k %d l %d t1 %f t2 %f \n",i,j,k,l,t1,t2); | |
272 | t1=swab((char *)&t1); t2=swab((char *)&t2); |
|
272 | t1=swab((char *)&t1); t2=swab((char *)&t2); | |
273 | printf(" in read_exlib 2 i %d j %d k %d l %d t1 %f t2 %f \n",i,j,k,l,t1,t2); |
|
273 | printf(" in read_exlib 2 i %d j %d k %d l %d t1 %f t2 %f \n",i,j,k,l,t1,t2); | |
274 | exlib[i][j][k][l]=Complex(t1,t2); |
|
274 | exlib[i][j][k][l]=Complex(t1,t2); | |
275 | printf(" in read_exlib 3 i %d j %d k %d l %d t1 %f t2 %f \n",i,j,k,l,exlib[i][j][k][l].r,exlib[i][j][k][l].i); |
|
275 | printf(" in read_exlib 3 i %d j %d k %d l %d t1 %f t2 %f \n",i,j,k,l,exlib[i][j][k][l].r,exlib[i][j][k][l].i); | |
276 | getchar(); |
|
276 | getchar(); | |
277 | } |
|
277 | } | |
278 |
|
278 | |||
279 | fclose(fp); |
|
279 | fclose(fp); | |
280 | printf("Done reading exlib file \n"); |
|
280 | printf("Done reading exlib file \n"); | |
281 | } |
|
281 | } | |
282 | */ |
|
282 | */ |
@@ -1,406 +1,407 | |||||
1 | # Copyright (c) 2012-2020 Jicamarca Radio Observatory |
|
1 | # Copyright (c) 2012-2020 Jicamarca Radio Observatory | |
2 | # All rights reserved. |
|
2 | # All rights reserved. | |
3 | # |
|
3 | # | |
4 | # Distributed under the terms of the BSD 3-clause license. |
|
4 | # Distributed under the terms of the BSD 3-clause license. | |
5 | """schainpy is an open source library to read, write and process radar data |
|
5 | """schainpy is an open source library to read, write and process radar data | |
6 |
|
6 | |||
7 | Signal Chain is a radar data processing library wich includes modules to read, |
|
7 | Signal Chain is a radar data processing library wich includes modules to read, | |
8 | and write different files formats, besides modules to process and visualize the |
|
8 | and write different files formats, besides modules to process and visualize the | |
9 | data. |
|
9 | data. | |
10 | """ |
|
10 | """ | |
11 |
|
11 | |||
12 | import os |
|
12 | import os | |
13 | from setuptools import setup, Extension |
|
13 | from setuptools import setup, Extension | |
14 | from setuptools.command.build_ext import build_ext as _build_ext |
|
14 | from setuptools.command.build_ext import build_ext as _build_ext | |
15 | from schainpy import __version__ |
|
15 | from schainpy import __version__ | |
16 |
|
16 | |||
17 | DOCLINES = __doc__.split("\n") |
|
17 | DOCLINES = __doc__.split("\n") | |
18 |
|
18 | |||
19 | class build_ext(_build_ext): |
|
19 | class build_ext(_build_ext): | |
20 | def finalize_options(self): |
|
20 | def finalize_options(self): | |
21 | _build_ext.finalize_options(self) |
|
21 | _build_ext.finalize_options(self) | |
22 | # Prevent numpy from thinking it is still in its setup process: |
|
22 | # Prevent numpy from thinking it is still in its setup process: | |
23 | __builtins__.__NUMPY_SETUP__ = False |
|
23 | __builtins__.__NUMPY_SETUP__ = False | |
24 | import numpy |
|
24 | import numpy | |
25 | self.include_dirs.append(numpy.get_include()) |
|
25 | self.include_dirs.append(numpy.get_include()) | |
26 |
|
26 | |||
27 | setup( |
|
27 | setup( | |
28 | name = "schainpy", |
|
28 | name = "schainpy", | |
29 | version = __version__, |
|
29 | version = __version__, | |
30 | description = DOCLINES[0], |
|
30 | description = DOCLINES[0], | |
31 | long_description = "\n".join(DOCLINES[2:]), |
|
31 | long_description = "\n".join(DOCLINES[2:]), | |
32 | url = "https://github.com/JRO-Peru/schainpy", |
|
32 | url = "https://github.com/JRO-Peru/schainpy", | |
33 | author = "Jicamarca Radio Observatory", |
|
33 | author = "Jicamarca Radio Observatory", | |
34 | author_email = "jro-developers@jro.igp.gob.pe", |
|
34 | author_email = "jro-developers@jro.igp.gob.pe", | |
35 | license="BSD-3-Clause", |
|
35 | license="BSD-3-Clause", | |
36 | classifiers=[ |
|
36 | classifiers=[ | |
37 | "Development Status :: 4 - Beta", |
|
37 | "Development Status :: 4 - Beta", | |
38 | "Environment :: Console", |
|
38 | "Environment :: Console", | |
39 | "Intended Audience :: Science/Research", |
|
39 | "Intended Audience :: Science/Research", | |
40 | "License :: OSI Approved :: BSD License", |
|
40 | "License :: OSI Approved :: BSD License", | |
41 | "Operating System :: MacOS :: MacOS X", |
|
41 | "Operating System :: MacOS :: MacOS X", | |
42 | "Operating System :: POSIX :: Linux", |
|
42 | "Operating System :: POSIX :: Linux", | |
43 | "Programming Language :: Python :: 2", |
|
43 | "Programming Language :: Python :: 2", | |
44 | "Programming Language :: Python :: 2.7", |
|
44 | "Programming Language :: Python :: 2.7", | |
45 | "Programming Language :: Python :: 3", |
|
45 | "Programming Language :: Python :: 3", | |
46 | "Programming Language :: Python :: 3.5", |
|
46 | "Programming Language :: Python :: 3.5", | |
47 | "Programming Language :: Python :: 3.6", |
|
47 | "Programming Language :: Python :: 3.6", | |
48 | "Programming Language :: Python :: 3.7", |
|
48 | "Programming Language :: Python :: 3.7", | |
49 | "Programming Language :: Python :: 3.8", |
|
49 | "Programming Language :: Python :: 3.8", | |
50 | "Topic :: Scientific/Engineering", |
|
50 | "Topic :: Scientific/Engineering", | |
51 | ], |
|
51 | ], | |
52 | packages = { |
|
52 | packages = { | |
53 | 'schainpy', |
|
53 | 'schainpy', | |
54 | 'schainpy.model', |
|
54 | 'schainpy.model', | |
55 | 'schainpy.model.data', |
|
55 | 'schainpy.model.data', | |
56 | 'schainpy.model.graphics', |
|
56 | 'schainpy.model.graphics', | |
57 | 'schainpy.model.io', |
|
57 | 'schainpy.model.io', | |
58 | 'schainpy.model.proc', |
|
58 | 'schainpy.model.proc', | |
59 | 'schainpy.model.utils', |
|
59 | 'schainpy.model.utils', | |
60 | 'schainpy.utils', |
|
60 | 'schainpy.utils', | |
61 | 'schainpy.gui', |
|
61 | 'schainpy.gui', | |
62 | 'schainpy.cli', |
|
62 | 'schainpy.cli', | |
63 | }, |
|
63 | }, | |
64 | package_data = {'': ['schain.conf.template'], |
|
64 | package_data = {'': ['schain.conf.template'], | |
65 | 'schainpy.files': ['*.oga'] |
|
65 | 'schainpy.files': ['*.oga'] | |
66 | }, |
|
66 | }, | |
67 | include_package_data = False, |
|
67 | include_package_data = False, | |
68 | scripts = ['schainpy/gui/schainGUI'], |
|
68 | scripts = ['schainpy/gui/schainGUI'], | |
69 | entry_points = { |
|
69 | entry_points = { | |
70 | 'console_scripts': [ |
|
70 | 'console_scripts': [ | |
71 | 'schain = schainpy.cli.cli:main', |
|
71 | 'schain = schainpy.cli.cli:main', | |
72 | ], |
|
72 | ], | |
73 | }, |
|
73 | }, | |
74 | cmdclass = {'build_ext': build_ext}, |
|
74 | cmdclass = {'build_ext': build_ext}, | |
75 | ext_modules=[ |
|
75 | ext_modules=[ | |
76 | Extension("schainpy.model.data._noise", ["schainc/_noise.c"]), |
|
76 | Extension("schainpy.model.data._noise", ["schainc/_noise.c"]), | |
77 | Extension("schainpy.model.data._HS_algorithm", ["schainc/_HS_algorithm.c"]), |
|
77 | Extension("schainpy.model.data._HS_algorithm", ["schainc/_HS_algorithm.c"]), | |
78 | ], |
|
78 | ], | |
79 | setup_requires = ["numpy"], |
|
79 | setup_requires = ["numpy"], | |
80 | install_requires = [ |
|
80 | install_requires = [ | |
81 | "scipy", |
|
81 | "scipy", | |
82 | "h5py", |
|
82 | "h5py", | |
83 | "matplotlib", |
|
83 | "matplotlib", | |
84 | "pyzmq", |
|
84 | "pyzmq", | |
85 | "fuzzywuzzy", |
|
85 | "fuzzywuzzy", | |
86 | "click", |
|
86 | "click", | |
87 | ], |
|
87 | ], | |
88 | ) |
|
88 | ) | |
89 |
|
89 | |||
90 | main_path = os.getcwd() |
|
90 | main_path = os.getcwd() | |
91 | child_path = '/schainf/Ffiles/bfmodel/' |
|
91 | child_path = '/schainf/Ffiles/bfmodel/' | |
92 | child_path_r = '/schainf/Ffiles/jlib26feb2001' |
|
92 | child_path_r = '/schainf/Ffiles/jlib26feb2001' | |
93 | #print("main_path: ", main_path) |
|
93 | #print("main_path: ", main_path) | |
94 | my_str = ' '+"ppath = "+"\""+main_path+"\"\n" |
|
94 | my_str = ' '+"ppath = "+"\""+main_path+"\"\n" | |
95 | my_str_2 = ' '+"cpath = "+"\""+child_path+"\"" |
|
95 | my_str_2 = ' '+"cpath = "+"\""+child_path+"\"" | |
96 | my_str_2_r = ' '+"cpath = "+"\""+child_path_r+"\"" |
|
96 | my_str_2_r = ' '+"cpath = "+"\""+child_path_r+"\"" | |
97 |
|
97 | |||
98 | #''' |
|
98 | #''' | |
99 | with open('./schainf/Ffiles/get_path_1.f', 'r') as path1: data1 = path1.read() |
|
99 | with open('./schainf/Ffiles/get_path_1.f', 'r') as path1: data1 = path1.read() | |
100 | with open('./schainf/Ffiles/get_path_2.f', 'r') as path2: data2 = path2.read() |
|
100 | with open('./schainf/Ffiles/get_path_2.f', 'r') as path2: data2 = path2.read() | |
101 | with open('./schainf/Ffiles/get_path.f', 'w') as final: final.write(data1+my_str+my_str_2+data2) |
|
101 | with open('./schainf/Ffiles/get_path.f', 'w') as final: final.write(data1+my_str+my_str_2+data2) | |
102 |
|
102 | |||
103 | with open('./schainf/Ffiles/get_path_1_reader.f', 'r') as p1: data1_r = p1.read() |
|
103 | with open('./schainf/Ffiles/get_path_1_reader.f', 'r') as p1: data1_r = p1.read() | |
104 | with open('./schainf/Ffiles/get_path_2_reader.f', 'r') as p2: data2_r = p2.read() |
|
104 | with open('./schainf/Ffiles/get_path_2_reader.f', 'r') as p2: data2_r = p2.read() | |
105 | with open('./schainf/Ffiles/get_path_reader.f', 'w') as final_r: final_r.write(data1_r+my_str+my_str_2_r+data2_r) |
|
105 | with open('./schainf/Ffiles/get_path_reader.f', 'w') as final_r: final_r.write(data1_r+my_str+my_str_2_r+data2_r) | |
106 | #''' |
|
106 | #''' | |
107 |
|
107 | |||
108 | from numpy.distutils.core import Extension, setup |
|
108 | from numpy.distutils.core import Extension, setup | |
109 |
|
109 | |||
110 | setup(name='schainpy', |
|
110 | setup(name='schainpy', | |
111 | ext_modules = [ |
|
111 | ext_modules = [ | |
112 | Extension("schainpy.model.proc.mkfact_short_2020_2", |
|
112 | Extension("schainpy.model.proc.mkfact_short_2020_2", | |
113 | sources=[ |
|
113 | sources=[ | |
114 | "schainf/Ffiles/mkfact_short_2020_2.pyf", |
|
114 | "schainf/Ffiles/mkfact_short_2020_2.pyf", | |
115 | "schainf/Ffiles/lmdif1.f", |
|
115 | "schainf/Ffiles/lmdif1.f", | |
116 | "schainf/Ffiles/mkfact.f", |
|
116 | "schainf/Ffiles/mkfact.f", | |
117 | "schainf/Ffiles/r1mach.f", |
|
117 | "schainf/Ffiles/r1mach.f", | |
118 | "schainf/Ffiles/bfield2.f", |
|
118 | "schainf/Ffiles/bfield2.f", | |
119 | "schainf/Ffiles/get_path.f"], |
|
119 | "schainf/Ffiles/get_path.f"], | |
120 | extra_f77_compile_args=["-fallow-argument-mismatch"]), |
|
120 | extra_f77_compile_args=["-fallow-argument-mismatch"]), | |
121 | Extension("schainpy.model.proc.fitacf_guess", |
|
121 | Extension("schainpy.model.proc.fitacf_guess", | |
122 | sources=[ |
|
122 | sources=[ | |
123 | "schainf/Ffiles/fitacf_guess.pyf", |
|
123 | "schainf/Ffiles/fitacf_guess.pyf", | |
124 | "schainf/Ffiles/fitacf_guess.f",], |
|
124 | "schainf/Ffiles/fitacf_guess.f",], | |
125 | extra_f77_compile_args=["-fallow-argument-mismatch"]), |
|
125 | extra_f77_compile_args=["-fallow-argument-mismatch"]), | |
126 | Extension("schainpy.model.proc.fitacf_acf2", |
|
126 | Extension("schainpy.model.proc.fitacf_acf2", | |
127 | sources = [ |
|
127 | sources = [ | |
128 | "schainf/Ffiles/fitacf_acf2.pyf", |
|
128 | "schainf/Ffiles/fitacf_acf2.pyf", | |
129 | "schainf/Ffiles/full_profile_profile.f", |
|
129 | "schainf/Ffiles/full_profile_profile.f", | |
130 | "schainf/Ffiles/fitacf.f", |
|
130 | "schainf/Ffiles/fitacf.f", | |
131 | "schainf/Ffiles/r1mach.f", |
|
131 | "schainf/Ffiles/r1mach.f", | |
132 | "schainf/Ffiles/lmdif1.f", |
|
132 | "schainf/Ffiles/lmdif1.f", | |
133 | "schainf/Ffiles/lagp.f", |
|
133 | "schainf/Ffiles/lagp.f", | |
134 | "schainf/Ffiles/reader.c", |
|
134 | "schainf/Ffiles/reader.c", | |
135 | "schainf/Ffiles/cbesi.f", |
|
135 | "schainf/Ffiles/cbesi.f", | |
136 | "schainf/Ffiles/i1mach.f", |
|
136 | "schainf/Ffiles/i1mach.f", | |
137 | "schainf/Ffiles/zeta.f", |
|
137 | "schainf/Ffiles/zeta.f", | |
138 | "schainf/Ffiles/qc25f.f", |
|
138 | "schainf/Ffiles/qc25f.f", | |
139 | "schainf/Ffiles/qwgtf.f", |
|
139 | "schainf/Ffiles/qwgtf.f", | |
140 | "schainf/Ffiles/qcheb.f", |
|
140 | "schainf/Ffiles/qcheb.f", | |
141 | "schainf/Ffiles/sgtsl.f", |
|
141 | "schainf/Ffiles/sgtsl.f", | |
142 | "schainf/Ffiles/qk15w.f", |
|
142 | "schainf/Ffiles/qk15w.f", | |
143 | "schainf/Ffiles/complex.c", |
|
143 | "schainf/Ffiles/complex.c", | |
144 | "schainf/Ffiles/cbinu.f", |
|
144 | "schainf/Ffiles/cbinu.f", | |
145 | "schainf/Ffiles/cseri.f", |
|
145 | "schainf/Ffiles/cseri.f", | |
146 | "schainf/Ffiles/cwrsk.f", |
|
146 | "schainf/Ffiles/cwrsk.f", | |
147 | "schainf/Ffiles/crati.f", |
|
147 | "schainf/Ffiles/crati.f", | |
148 | "schainf/Ffiles/casyi.f", |
|
148 | "schainf/Ffiles/casyi.f", | |
149 | "schainf/Ffiles/cbuni.f", |
|
149 | "schainf/Ffiles/cbuni.f", | |
150 | "schainf/Ffiles/cuni2.f", |
|
150 | "schainf/Ffiles/cuni2.f", | |
151 | "schainf/Ffiles/gamln.f", |
|
151 | "schainf/Ffiles/gamln.f", | |
152 | "schainf/Ffiles/cuchk.f", |
|
152 | "schainf/Ffiles/cuchk.f", | |
153 | "schainf/Ffiles/cbknu.f", |
|
153 | "schainf/Ffiles/cbknu.f", | |
154 | "schainf/Ffiles/cshch.f", |
|
154 | "schainf/Ffiles/cshch.f", | |
155 | "schainf/Ffiles/ckscl.f", |
|
155 | "schainf/Ffiles/ckscl.f", | |
156 | "schainf/Ffiles/cuoik.f", |
|
156 | "schainf/Ffiles/cuoik.f", | |
157 | "schainf/Ffiles/cunik.f", |
|
157 | "schainf/Ffiles/cunik.f", | |
158 | "schainf/Ffiles/cuni1.f", |
|
158 | "schainf/Ffiles/cuni1.f", | |
159 | "schainf/Ffiles/cairy.f", |
|
159 | "schainf/Ffiles/cairy.f", | |
160 | "schainf/Ffiles/cmlri.f", |
|
160 | "schainf/Ffiles/cmlri.f", | |
161 | "schainf/Ffiles/cunhj.f", |
|
161 | "schainf/Ffiles/cunhj.f", | |
162 | "schainf/Ffiles/cacai.f", |
|
162 | "schainf/Ffiles/cacai.f", | |
163 | "schainf/Ffiles/csisl.f", |
|
163 | "schainf/Ffiles/csisl.f", | |
164 | "schainf/Ffiles/caxpy.f", |
|
164 | "schainf/Ffiles/caxpy.f", | |
165 | "schainf/Ffiles/cs1s2.f", |
|
165 | "schainf/Ffiles/cs1s2.f", | |
166 | "schainf/Ffiles/scabs1.f", |
|
166 | "schainf/Ffiles/scabs1.f", | |
167 | "schainf/Ffiles/cdotu.f", |
|
167 | "schainf/Ffiles/cdotu.f", | |
168 | "schainf/Ffiles/rs.f", |
|
168 | "schainf/Ffiles/rs.f", | |
169 | "schainf/Ffiles/sppfa.f", |
|
169 | "schainf/Ffiles/sppfa.f", | |
170 | "schainf/Ffiles/sdot.f", |
|
170 | "schainf/Ffiles/sdot.f", | |
171 | "schainf/Ffiles/tred2.f", |
|
171 | "schainf/Ffiles/tred2.f", | |
172 | "schainf/Ffiles/tql2.f", |
|
172 | "schainf/Ffiles/tql2.f", | |
173 | "schainf/Ffiles/sppdi.f", |
|
173 | "schainf/Ffiles/sppdi.f", | |
174 | "schainf/Ffiles/saxpy.f", |
|
174 | "schainf/Ffiles/saxpy.f", | |
175 | "schainf/Ffiles/sscal.f", |
|
175 | "schainf/Ffiles/sscal.f", | |
176 | "schainf/Ffiles/pythag.f", |
|
176 | "schainf/Ffiles/pythag.f", | |
177 | "schainf/Ffiles/tql1.f", |
|
177 | "schainf/Ffiles/tql1.f", | |
178 | "schainf/Ffiles/get_path.f", |
|
178 | "schainf/Ffiles/get_path.f", | |
179 | "schainf/Ffiles/tred1.f"], |
|
179 | "schainf/Ffiles/tred1.f"], | |
180 | extra_f77_compile_args=["-fallow-argument-mismatch"]), |
|
180 | extra_f77_compile_args=["-fallow-argument-mismatch"]), | |
181 | Extension("schainpy.model.proc.fitacf_fit_short", |
|
181 | Extension("schainpy.model.proc.fitacf_fit_short", | |
182 | sources = [ |
|
182 | sources = [ | |
183 | "schainf/Ffiles/fitacf_fit_short.pyf", |
|
183 | "schainf/Ffiles/fitacf_fit_short.pyf", | |
184 | "schainf/Ffiles/cairy.f", |
|
184 | "schainf/Ffiles/cairy.f", | |
185 | "schainf/Ffiles/casyi.f", |
|
185 | "schainf/Ffiles/casyi.f", | |
186 | "schainf/Ffiles/cbesi.f", |
|
186 | "schainf/Ffiles/cbesi.f", | |
187 | "schainf/Ffiles/cbinu.f", |
|
187 | "schainf/Ffiles/cbinu.f", | |
188 | "schainf/Ffiles/cbknu.f", |
|
188 | "schainf/Ffiles/cbknu.f", | |
189 | "schainf/Ffiles/cbuni.f", |
|
189 | "schainf/Ffiles/cbuni.f", | |
190 | "schainf/Ffiles/ckscl.f", |
|
190 | "schainf/Ffiles/ckscl.f", | |
191 | "schainf/Ffiles/crati.f", |
|
191 | "schainf/Ffiles/crati.f", | |
192 | "schainf/Ffiles/cacai.f", |
|
192 | "schainf/Ffiles/cacai.f", | |
193 | "schainf/Ffiles/cmlri.f", |
|
193 | "schainf/Ffiles/cmlri.f", | |
194 | "schainf/Ffiles/cs1s2.f", |
|
194 | "schainf/Ffiles/cs1s2.f", | |
195 | "schainf/Ffiles/cseri.f", |
|
195 | "schainf/Ffiles/cseri.f", | |
196 | "schainf/Ffiles/cshch.f", |
|
196 | "schainf/Ffiles/cshch.f", | |
197 | "schainf/Ffiles/cuchk.f", |
|
197 | "schainf/Ffiles/cuchk.f", | |
198 | "schainf/Ffiles/cunhj.f", |
|
198 | "schainf/Ffiles/cunhj.f", | |
199 | "schainf/Ffiles/cuni1.f", |
|
199 | "schainf/Ffiles/cuni1.f", | |
200 | "schainf/Ffiles/complex.c", |
|
200 | "schainf/Ffiles/complex.c", | |
201 | "schainf/Ffiles/cuni2.f", |
|
201 | "schainf/Ffiles/cuni2.f", | |
202 | "schainf/Ffiles/cunik.f", |
|
202 | "schainf/Ffiles/cunik.f", | |
203 | "schainf/Ffiles/cuoik.f", |
|
203 | "schainf/Ffiles/cuoik.f", | |
204 | "schainf/Ffiles/cwrsk.f", |
|
204 | "schainf/Ffiles/cwrsk.f", | |
205 | "schainf/Ffiles/fitacf_fit_short.f", |
|
205 | "schainf/Ffiles/fitacf_fit_short.f", | |
206 | "schainf/Ffiles/gamln.f", |
|
206 | "schainf/Ffiles/gamln.f", | |
207 | "schainf/Ffiles/i1mach.f", |
|
207 | "schainf/Ffiles/i1mach.f", | |
208 | "schainf/Ffiles/lmdif1.f", |
|
208 | "schainf/Ffiles/lmdif1.f", | |
209 | "schainf/Ffiles/pythag.f", |
|
209 | "schainf/Ffiles/pythag.f", | |
210 | "schainf/Ffiles/qc25f.f", |
|
210 | "schainf/Ffiles/qc25f.f", | |
211 | "schainf/Ffiles/qcheb.f", |
|
211 | "schainf/Ffiles/qcheb.f", | |
212 | "schainf/Ffiles/qk15w.f", |
|
212 | "schainf/Ffiles/qk15w.f", | |
213 | "schainf/Ffiles/qwgtf.f", |
|
213 | "schainf/Ffiles/qwgtf.f", | |
214 | "schainf/Ffiles/r1mach.f", |
|
214 | "schainf/Ffiles/r1mach.f", | |
215 | "schainf/Ffiles/reader.c", |
|
215 | "schainf/Ffiles/reader.c", | |
216 | "schainf/Ffiles/rs.f", |
|
216 | "schainf/Ffiles/rs.f", | |
217 | "schainf/Ffiles/saxpy.f", |
|
217 | "schainf/Ffiles/saxpy.f", | |
218 | "schainf/Ffiles/sdot.f", |
|
218 | "schainf/Ffiles/sdot.f", | |
219 | "schainf/Ffiles/sgtsl.f", |
|
219 | "schainf/Ffiles/sgtsl.f", | |
220 | "schainf/Ffiles/sppdi.f", |
|
220 | "schainf/Ffiles/sppdi.f", | |
221 | "schainf/Ffiles/sppfa.f", |
|
221 | "schainf/Ffiles/sppfa.f", | |
222 | "schainf/Ffiles/sscal.f", |
|
222 | "schainf/Ffiles/sscal.f", | |
223 | "schainf/Ffiles/tql1.f", |
|
223 | "schainf/Ffiles/tql1.f", | |
224 | "schainf/Ffiles/tql2.f", |
|
224 | "schainf/Ffiles/tql2.f", | |
225 | "schainf/Ffiles/tred1.f", |
|
225 | "schainf/Ffiles/tred1.f", | |
226 | "schainf/Ffiles/tred2.f", |
|
226 | "schainf/Ffiles/tred2.f", | |
227 | "schainf/Ffiles/get_path.f", |
|
227 | "schainf/Ffiles/get_path.f", | |
|
228 | "schainf/Ffiles/get_path_reader.f", | |||
228 | "schainf/Ffiles/zeta.f"], |
|
229 | "schainf/Ffiles/zeta.f"], | |
229 | extra_f77_compile_args=["-fallow-argument-mismatch"]), |
|
230 | extra_f77_compile_args=["-fallow-argument-mismatch"]), | |
230 | Extension("schainpy.model.proc.full_profile_profile", |
|
231 | Extension("schainpy.model.proc.full_profile_profile", | |
231 | sources = [ |
|
232 | sources = [ | |
232 | "schainf/Ffiles/full_profile_profile.pyf", |
|
233 | "schainf/Ffiles/full_profile_profile.pyf", | |
233 | "schainf/Ffiles/full_profile_profile.f", |
|
234 | "schainf/Ffiles/full_profile_profile.f", | |
234 | "schainf/Ffiles/fitacf.f", |
|
235 | "schainf/Ffiles/fitacf.f", | |
235 | "schainf/Ffiles/r1mach.f", |
|
236 | "schainf/Ffiles/r1mach.f", | |
236 | "schainf/Ffiles/lmdif1.f", |
|
237 | "schainf/Ffiles/lmdif1.f", | |
237 | "schainf/Ffiles/reader.c", |
|
238 | "schainf/Ffiles/reader.c", | |
238 | "schainf/Ffiles/cbesi.f", |
|
239 | "schainf/Ffiles/cbesi.f", | |
239 | "schainf/Ffiles/lagp.f", |
|
240 | "schainf/Ffiles/lagp.f", | |
240 | "schainf/Ffiles/i1mach.f", |
|
241 | "schainf/Ffiles/i1mach.f", | |
241 | "schainf/Ffiles/zeta.f", |
|
242 | "schainf/Ffiles/zeta.f", | |
242 | "schainf/Ffiles/qc25f.f", |
|
243 | "schainf/Ffiles/qc25f.f", | |
243 | "schainf/Ffiles/qwgtf.f", |
|
244 | "schainf/Ffiles/qwgtf.f", | |
244 | "schainf/Ffiles/qcheb.f", |
|
245 | "schainf/Ffiles/qcheb.f", | |
245 | "schainf/Ffiles/sgtsl.f", |
|
246 | "schainf/Ffiles/sgtsl.f", | |
246 | "schainf/Ffiles/qk15w.f", |
|
247 | "schainf/Ffiles/qk15w.f", | |
247 | "schainf/Ffiles/cbinu.f", |
|
248 | "schainf/Ffiles/cbinu.f", | |
248 | "schainf/Ffiles/complex.c", |
|
249 | "schainf/Ffiles/complex.c", | |
249 | "schainf/Ffiles/cseri.f", |
|
250 | "schainf/Ffiles/cseri.f", | |
250 | "schainf/Ffiles/cwrsk.f", |
|
251 | "schainf/Ffiles/cwrsk.f", | |
251 | "schainf/Ffiles/crati.f", |
|
252 | "schainf/Ffiles/crati.f", | |
252 | "schainf/Ffiles/casyi.f", |
|
253 | "schainf/Ffiles/casyi.f", | |
253 | "schainf/Ffiles/cbuni.f", |
|
254 | "schainf/Ffiles/cbuni.f", | |
254 | "schainf/Ffiles/cuni2.f", |
|
255 | "schainf/Ffiles/cuni2.f", | |
255 | "schainf/Ffiles/gamln.f", |
|
256 | "schainf/Ffiles/gamln.f", | |
256 | "schainf/Ffiles/cuchk.f", |
|
257 | "schainf/Ffiles/cuchk.f", | |
257 | "schainf/Ffiles/cbknu.f", |
|
258 | "schainf/Ffiles/cbknu.f", | |
258 | "schainf/Ffiles/cshch.f", |
|
259 | "schainf/Ffiles/cshch.f", | |
259 | "schainf/Ffiles/ckscl.f", |
|
260 | "schainf/Ffiles/ckscl.f", | |
260 | "schainf/Ffiles/cuoik.f", |
|
261 | "schainf/Ffiles/cuoik.f", | |
261 | "schainf/Ffiles/cunik.f", |
|
262 | "schainf/Ffiles/cunik.f", | |
262 | "schainf/Ffiles/cuni1.f", |
|
263 | "schainf/Ffiles/cuni1.f", | |
263 | "schainf/Ffiles/cairy.f", |
|
264 | "schainf/Ffiles/cairy.f", | |
264 | "schainf/Ffiles/cmlri.f", |
|
265 | "schainf/Ffiles/cmlri.f", | |
265 | "schainf/Ffiles/cunhj.f", |
|
266 | "schainf/Ffiles/cunhj.f", | |
266 | "schainf/Ffiles/cacai.f", |
|
267 | "schainf/Ffiles/cacai.f", | |
267 | "schainf/Ffiles/csisl.f", |
|
268 | "schainf/Ffiles/csisl.f", | |
268 | "schainf/Ffiles/caxpy.f", |
|
269 | "schainf/Ffiles/caxpy.f", | |
269 | "schainf/Ffiles/cs1s2.f", |
|
270 | "schainf/Ffiles/cs1s2.f", | |
270 | "schainf/Ffiles/scabs1.f", |
|
271 | "schainf/Ffiles/scabs1.f", | |
271 | "schainf/Ffiles/cdotu.f", |
|
272 | "schainf/Ffiles/cdotu.f", | |
272 | "schainf/Ffiles/rs.f", |
|
273 | "schainf/Ffiles/rs.f", | |
273 | "schainf/Ffiles/sppfa.f", |
|
274 | "schainf/Ffiles/sppfa.f", | |
274 | "schainf/Ffiles/sdot.f", |
|
275 | "schainf/Ffiles/sdot.f", | |
275 | "schainf/Ffiles/tred2.f", |
|
276 | "schainf/Ffiles/tred2.f", | |
276 | "schainf/Ffiles/tql2.f", |
|
277 | "schainf/Ffiles/tql2.f", | |
277 | "schainf/Ffiles/sppdi.f", |
|
278 | "schainf/Ffiles/sppdi.f", | |
278 | "schainf/Ffiles/saxpy.f", |
|
279 | "schainf/Ffiles/saxpy.f", | |
279 | "schainf/Ffiles/sscal.f", |
|
280 | "schainf/Ffiles/sscal.f", | |
280 | "schainf/Ffiles/pythag.f", |
|
281 | "schainf/Ffiles/pythag.f", | |
281 | "schainf/Ffiles/tql1.f", |
|
282 | "schainf/Ffiles/tql1.f", | |
282 | "schainf/Ffiles/get_path.f", |
|
283 | "schainf/Ffiles/get_path.f", | |
283 | "schainf/Ffiles/tred1.f"], |
|
284 | "schainf/Ffiles/tred1.f"], | |
284 | extra_f77_compile_args=["-fallow-argument-mismatch"]) |
|
285 | extra_f77_compile_args=["-fallow-argument-mismatch"]) | |
285 | ] |
|
286 | ] | |
286 | ) |
|
287 | ) | |
287 |
|
288 | |||
288 | ''' |
|
289 | ''' | |
289 | setup(name='schainpy', |
|
290 | setup(name='schainpy', | |
290 | ext_modules = [ |
|
291 | ext_modules = [ | |
291 | Extension("schainpy.model.proc.mkfact_short_2020_2", |
|
292 | Extension("schainpy.model.proc.mkfact_short_2020_2", | |
292 | sources=[ |
|
293 | sources=[ | |
293 | "schainf/mkfact/mkfact_short_2020_2.pyf", |
|
294 | "schainf/mkfact/mkfact_short_2020_2.pyf", | |
294 | "schainf/mkfact/lmdif1.f", |
|
295 | "schainf/mkfact/lmdif1.f", | |
295 | "schainf/mkfact/mkfact.f", |
|
296 | "schainf/mkfact/mkfact.f", | |
296 | "schainf/mkfact/r1mach.f", |
|
297 | "schainf/mkfact/r1mach.f", | |
297 | "schainf/mkfact/bfield2.f"]), |
|
298 | "schainf/mkfact/bfield2.f"]), | |
298 | Extension("schainpy.model.proc.full_profile_profile", |
|
299 | Extension("schainpy.model.proc.full_profile_profile", | |
299 | sources=[ |
|
300 | sources=[ | |
300 | "schainf/full_profile/full_profile_profile.pyf", |
|
301 | "schainf/full_profile/full_profile_profile.pyf", | |
301 | "schainf/full_profile/full_profile_profile.f", |
|
302 | "schainf/full_profile/full_profile_profile.f", | |
302 | "schainf/full_profile/fitacf.f", |
|
303 | "schainf/full_profile/fitacf.f", | |
303 | "schainf/full_profile/r1mach.f", |
|
304 | "schainf/full_profile/r1mach.f", | |
304 | "schainf/full_profile/lmdif1.f", |
|
305 | "schainf/full_profile/lmdif1.f", | |
305 | "schainf/full_profile/lagp.f", |
|
306 | "schainf/full_profile/lagp.f", | |
306 | "schainf/full_profile/reader.c", |
|
307 | "schainf/full_profile/reader.c", | |
307 | "schainf/full_profile/cbesi.f", |
|
308 | "schainf/full_profile/cbesi.f", | |
308 | "schainf/full_profile/i1mach.f", |
|
309 | "schainf/full_profile/i1mach.f", | |
309 | "schainf/full_profile/zeta.f", |
|
310 | "schainf/full_profile/zeta.f", | |
310 | "schainf/full_profile/qc25f.f", |
|
311 | "schainf/full_profile/qc25f.f", | |
311 | "schainf/full_profile/qwgtf.f", |
|
312 | "schainf/full_profile/qwgtf.f", | |
312 | "schainf/full_profile/qcheb.f", |
|
313 | "schainf/full_profile/qcheb.f", | |
313 | "schainf/full_profile/sgtsl.f", |
|
314 | "schainf/full_profile/sgtsl.f", | |
314 | "schainf/full_profile/qk15w.f", |
|
315 | "schainf/full_profile/qk15w.f", | |
315 | "schainf/full_profile/complex.c", |
|
316 | "schainf/full_profile/complex.c", | |
316 | "schainf/full_profile/cbinu.f", |
|
317 | "schainf/full_profile/cbinu.f", | |
317 | "schainf/full_profile/cseri.f", |
|
318 | "schainf/full_profile/cseri.f", | |
318 | "schainf/full_profile/cwrsk.f", |
|
319 | "schainf/full_profile/cwrsk.f", | |
319 | "schainf/full_profile/crati.f", |
|
320 | "schainf/full_profile/crati.f", | |
320 | "schainf/full_profile/casyi.f", |
|
321 | "schainf/full_profile/casyi.f", | |
321 | "schainf/full_profile/cbuni.f", |
|
322 | "schainf/full_profile/cbuni.f", | |
322 | "schainf/full_profile/cuni2.f", |
|
323 | "schainf/full_profile/cuni2.f", | |
323 | "schainf/full_profile/gamln.f", |
|
324 | "schainf/full_profile/gamln.f", | |
324 | "schainf/full_profile/cuchk.f", |
|
325 | "schainf/full_profile/cuchk.f", | |
325 | "schainf/full_profile/cbknu.f", |
|
326 | "schainf/full_profile/cbknu.f", | |
326 | "schainf/full_profile/cshch.f", |
|
327 | "schainf/full_profile/cshch.f", | |
327 | "schainf/full_profile/ckscl.f", |
|
328 | "schainf/full_profile/ckscl.f", | |
328 | "schainf/full_profile/cuoik.f", |
|
329 | "schainf/full_profile/cuoik.f", | |
329 | "schainf/full_profile/cunik.f", |
|
330 | "schainf/full_profile/cunik.f", | |
330 | "schainf/full_profile/cuni1.f", |
|
331 | "schainf/full_profile/cuni1.f", | |
331 | "schainf/full_profile/cairy.f", |
|
332 | "schainf/full_profile/cairy.f", | |
332 | "schainf/full_profile/cmlri.f", |
|
333 | "schainf/full_profile/cmlri.f", | |
333 | "schainf/full_profile/cunhj.f", |
|
334 | "schainf/full_profile/cunhj.f", | |
334 | "schainf/full_profile/cacai.f", |
|
335 | "schainf/full_profile/cacai.f", | |
335 | "schainf/full_profile/csisl.f", |
|
336 | "schainf/full_profile/csisl.f", | |
336 | "schainf/full_profile/caxpy.f", |
|
337 | "schainf/full_profile/caxpy.f", | |
337 | "schainf/full_profile/cs1s2.f", |
|
338 | "schainf/full_profile/cs1s2.f", | |
338 | "schainf/full_profile/scabs1.f", |
|
339 | "schainf/full_profile/scabs1.f", | |
339 | "schainf/full_profile/cdotu.f", |
|
340 | "schainf/full_profile/cdotu.f", | |
340 | "schainf/full_profile/rs.f", |
|
341 | "schainf/full_profile/rs.f", | |
341 | "schainf/full_profile/sppfa.f", |
|
342 | "schainf/full_profile/sppfa.f", | |
342 | "schainf/full_profile/sdot.f", |
|
343 | "schainf/full_profile/sdot.f", | |
343 | "schainf/full_profile/tred2.f", |
|
344 | "schainf/full_profile/tred2.f", | |
344 | "schainf/full_profile/tql2.f", |
|
345 | "schainf/full_profile/tql2.f", | |
345 | "schainf/full_profile/sppdi.f", |
|
346 | "schainf/full_profile/sppdi.f", | |
346 | "schainf/full_profile/saxpy.f", |
|
347 | "schainf/full_profile/saxpy.f", | |
347 | "schainf/full_profile/sscal.f", |
|
348 | "schainf/full_profile/sscal.f", | |
348 | "schainf/full_profile/pythag.f", |
|
349 | "schainf/full_profile/pythag.f", | |
349 | "schainf/full_profile/tql1.f", |
|
350 | "schainf/full_profile/tql1.f", | |
350 | "schainf/full_profile/tred1.f"]), |
|
351 | "schainf/full_profile/tred1.f"]), | |
351 | Extension("schainpy.model.proc.fitacf_acf2", |
|
352 | Extension("schainpy.model.proc.fitacf_acf2", | |
352 | sources = [ |
|
353 | sources = [ | |
353 | "schainf/acf2/fitacf_acf2.pyf", |
|
354 | "schainf/acf2/fitacf_acf2.pyf", | |
354 | "schainf/acf2/full_profile_profile.f", |
|
355 | "schainf/acf2/full_profile_profile.f", | |
355 | "schainf/acf2/fitacf.f", |
|
356 | "schainf/acf2/fitacf.f", | |
356 | "schainf/acf2/r1mach.f", |
|
357 | "schainf/acf2/r1mach.f", | |
357 | "schainf/acf2/lmdif1.f", |
|
358 | "schainf/acf2/lmdif1.f", | |
358 | "schainf/acf2/lagp.f", |
|
359 | "schainf/acf2/lagp.f", | |
359 | "schainf/acf2/reader.c", |
|
360 | "schainf/acf2/reader.c", | |
360 | "schainf/acf2/cbesi.f", |
|
361 | "schainf/acf2/cbesi.f", | |
361 | "schainf/acf2/i1mach.f", |
|
362 | "schainf/acf2/i1mach.f", | |
362 | "schainf/acf2/zeta.f", |
|
363 | "schainf/acf2/zeta.f", | |
363 | "schainf/acf2/qc25f.f", |
|
364 | "schainf/acf2/qc25f.f", | |
364 | "schainf/acf2/qwgtf.f", |
|
365 | "schainf/acf2/qwgtf.f", | |
365 | "schainf/acf2/qcheb.f", |
|
366 | "schainf/acf2/qcheb.f", | |
366 | "schainf/acf2/sgtsl.f", |
|
367 | "schainf/acf2/sgtsl.f", | |
367 | "schainf/acf2/qk15w.f", |
|
368 | "schainf/acf2/qk15w.f", | |
368 | "schainf/acf2/complex.c", |
|
369 | "schainf/acf2/complex.c", | |
369 | "schainf/acf2/cbinu.f", |
|
370 | "schainf/acf2/cbinu.f", | |
370 | "schainf/acf2/cseri.f", |
|
371 | "schainf/acf2/cseri.f", | |
371 | "schainf/acf2/cwrsk.f", |
|
372 | "schainf/acf2/cwrsk.f", | |
372 | "schainf/acf2/crati.f", |
|
373 | "schainf/acf2/crati.f", | |
373 | "schainf/acf2/casyi.f", |
|
374 | "schainf/acf2/casyi.f", | |
374 | "schainf/acf2/cbuni.f", |
|
375 | "schainf/acf2/cbuni.f", | |
375 | "schainf/acf2/cuni2.f", |
|
376 | "schainf/acf2/cuni2.f", | |
376 | "schainf/acf2/gamln.f", |
|
377 | "schainf/acf2/gamln.f", | |
377 | "schainf/acf2/cuchk.f", |
|
378 | "schainf/acf2/cuchk.f", | |
378 | "schainf/acf2/cbknu.f", |
|
379 | "schainf/acf2/cbknu.f", | |
379 | "schainf/acf2/cshch.f", |
|
380 | "schainf/acf2/cshch.f", | |
380 | "schainf/acf2/ckscl.f", |
|
381 | "schainf/acf2/ckscl.f", | |
381 | "schainf/acf2/cuoik.f", |
|
382 | "schainf/acf2/cuoik.f", | |
382 | "schainf/acf2/cunik.f", |
|
383 | "schainf/acf2/cunik.f", | |
383 | "schainf/acf2/cuni1.f", |
|
384 | "schainf/acf2/cuni1.f", | |
384 | "schainf/acf2/cairy.f", |
|
385 | "schainf/acf2/cairy.f", | |
385 | "schainf/acf2/cmlri.f", |
|
386 | "schainf/acf2/cmlri.f", | |
386 | "schainf/acf2/cunhj.f", |
|
387 | "schainf/acf2/cunhj.f", | |
387 | "schainf/acf2/cacai.f", |
|
388 | "schainf/acf2/cacai.f", | |
388 | "schainf/acf2/csisl.f", |
|
389 | "schainf/acf2/csisl.f", | |
389 | "schainf/acf2/caxpy.f", |
|
390 | "schainf/acf2/caxpy.f", | |
390 | "schainf/acf2/cs1s2.f", |
|
391 | "schainf/acf2/cs1s2.f", | |
391 | "schainf/acf2/scabs1.f", |
|
392 | "schainf/acf2/scabs1.f", | |
392 | "schainf/acf2/cdotu.f", |
|
393 | "schainf/acf2/cdotu.f", | |
393 | "schainf/acf2/rs.f", |
|
394 | "schainf/acf2/rs.f", | |
394 | "schainf/acf2/sppfa.f", |
|
395 | "schainf/acf2/sppfa.f", | |
395 | "schainf/acf2/sdot.f", |
|
396 | "schainf/acf2/sdot.f", | |
396 | "schainf/acf2/tred2.f", |
|
397 | "schainf/acf2/tred2.f", | |
397 | "schainf/acf2/tql2.f", |
|
398 | "schainf/acf2/tql2.f", | |
398 | "schainf/acf2/sppdi.f", |
|
399 | "schainf/acf2/sppdi.f", | |
399 | "schainf/acf2/saxpy.f", |
|
400 | "schainf/acf2/saxpy.f", | |
400 | "schainf/acf2/sscal.f", |
|
401 | "schainf/acf2/sscal.f", | |
401 | "schainf/acf2/pythag.f", |
|
402 | "schainf/acf2/pythag.f", | |
402 | "schainf/acf2/tql1.f", |
|
403 | "schainf/acf2/tql1.f", | |
403 | "schainf/acf2/tred1.f"]) |
|
404 | "schainf/acf2/tred1.f"]) | |
404 | ] |
|
405 | ] | |
405 | ) |
|
406 | ) | |
406 | ''' |
|
407 | ''' |
General Comments 0
You need to be logged in to leave comments.
Login now