36 # define N_SEQ_BORON 5
58 static bool lgFirst=
true,
193 for( j=0; j < N_SEQ_BORON; j++ )
204 stat[0] = t10->
Lo->
g;
205 stat[1] = t10->
Hi->
g;
206 stat[2] = t20->
Hi->
g;
207 stat[3] = t30->
Hi->
g;
208 stat[4] = t41->
Hi->
g;
209 ASSERT( stat[0]>0. && stat[1]>0. &&stat[2]>0. &&stat[3]>0. &&stat[4]>0.);
210 ASSERT( fabs(t10->
Lo->
g/2.-1.) < FLT_EPSILON);
211 ASSERT( fabs(t10->
Hi->
g/4.-1.) < FLT_EPSILON);
212 ASSERT( fabs(t20->
Lo->
g/2.-1.) < FLT_EPSILON);
213 ASSERT( fabs(t20->
Hi->
g/2.-1.) < FLT_EPSILON);
214 ASSERT( fabs(t30->
Lo->
g/2.-1.) < FLT_EPSILON);
215 ASSERT( fabs(t30->
Hi->
g/4.-1.) < FLT_EPSILON);
216 ASSERT( fabs(t21->
Lo->
g/4.-1.) < FLT_EPSILON);
217 ASSERT( fabs(t21->
Hi->
g/2.-1.) < FLT_EPSILON);
218 ASSERT( fabs(t31->
Lo->
g/4.-1.) < FLT_EPSILON);
219 ASSERT( fabs(t31->
Hi->
g/4.-1.) < FLT_EPSILON);
220 ASSERT( fabs(t41->
Lo->
g/4.-1.) < FLT_EPSILON);
221 ASSERT( fabs(t41->
Hi->
g/6.-1.) < FLT_EPSILON);
229 ASSERT( excit[1]>0. &&excit[2]>0. &&excit[3]>0. &&excit[4]>0.);
234 col_str[1][0] = t10->
Coll.
cs;
238 AulPump[0][1] += pump_rate;
242 col_str[2][0] = t20->
Coll.
cs;
247 col_str[3][0] = t30->
Coll.
cs;
250 AulEscp[4][0] = 1e-8;
252 col_str[4][0] = cs40;
257 col_str[2][1] = t21->
Coll.
cs;
262 col_str[3][1] = t31->
Coll.
cs;
267 col_str[4][1] = t41->
Coll.
cs;
270 AulEscp[3][2] = 1e-8;
272 col_str[3][2] = cs32;
275 AulEscp[4][2] = 1e-8;
277 col_str[4][2] = cs42;
280 AulEscp[4][3] = 1e-8;
282 col_str[4][3] = cs43;
379 for( i=2; i < N_SEQ_BORON; i++ )
391 for( i=0; i< N_SEQ_BORON; ++i )
397 t10->
Lo->
Pop = pops[0];
398 t20->
Lo->
Pop = pops[0];
399 t30->
Lo->
Pop = pops[0];
400 t21->
Lo->
Pop = pops[1];
401 t31->
Lo->
Pop = pops[1];
402 t41->
Lo->
Pop = pops[1];
411 t10->
Hi->
Pop = pops[1];
412 t20->
Hi->
Pop = pops[2];
413 t30->
Hi->
Pop = pops[3];
414 t21->
Hi->
Pop = pops[2];
415 t31->
Hi->
Pop = pops[3];
416 t41->
Hi->
Pop = pops[4];