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