リビジョン | 323ab84bc8b3e944b3a6887298843e805a5063e0 (tree) |
---|---|
日時 | 2008-10-01 19:41:42 |
作者 | iselllo |
コミッター | iselllo |
Just added some trivial changes to the code (mainly removed some print
statements).
@@ -78,7 +78,7 @@ | ||
78 | 78 | |
79 | 79 | def fitting_stat(x,y): |
80 | 80 | slope, intercept, r, prob2, see = stats.linregress(x,y) |
81 | - #print "len(x) is, ", len(x) | |
81 | + | |
82 | 82 | if (len(x)>2): |
83 | 83 | |
84 | 84 | see=see*s.sqrt(len(x)/(len(x)-2.)) |
@@ -128,7 +128,7 @@ | ||
128 | 128 | |
129 | 129 | #Specify the bin structure you want to use |
130 | 130 | |
131 | -linear =0 #linear ==1---> use a linear bin structure and solve smoluchowski equation in standard form | |
131 | +linear =1 #linear ==1---> use a linear bin structure and solve smoluchowski equation in standard form | |
132 | 132 | #linear !1---> use a non-linear (log-spaced) bin structure and solve smoluchowski equation |
133 | 133 | # using the splitting operator |
134 | 134 |
@@ -155,7 +155,6 @@ | ||
155 | 155 | |
156 | 156 | |
157 | 157 | |
158 | -print "vol_grid is, ", vol_grid | |
159 | 158 | |
160 | 159 | if (linear !=1): #calculate the splitting operator on the non-uniform grid |
161 | 160 | f_garrick=mycount_garrick(vol_grid) #I calculate the splitting operator on the grid |
@@ -164,15 +163,13 @@ | ||
164 | 163 | y0=s.zeros(n_bin) |
165 | 164 | y0[0]=initial_density #initial state (monodisperse aerosol) |
166 | 165 | |
167 | -print "the initial condition is, ", y0 | |
166 | + | |
168 | 167 | |
169 | 168 | #Generate the kernel matrix |
170 | 169 | |
171 | 170 | kernel=Brow_ker_cont_optim(vol_grid) |
172 | 171 | |
173 | -print "kernel is, ", kernel | |
174 | 172 | |
175 | -print "time is, ", t | |
176 | 173 | |
177 | 174 | if (linear==1): |
178 | 175 | solution = si.odeint(coupling_optim, y0, \ |