254 |
|
|
255 |
|
free(working_line); |
256 |
|
return(ntokens); |
257 |
+ |
} |
258 |
+ |
|
259 |
+ |
|
260 |
+ |
double SHAPE::getSigmaAt(double costheta, double phi) { |
261 |
+ |
|
262 |
+ |
vector<SHFunc*>::iterator sigmaIter; |
263 |
+ |
double sigma; |
264 |
+ |
|
265 |
+ |
sigma = 0.0; |
266 |
+ |
|
267 |
+ |
for(sigmaIter = sigmaFuncs.begin(); sigmaIter != sigmaFuncs.end(); |
268 |
+ |
++sigmaIter) |
269 |
+ |
sigma += (*sigmaIter)->getValueAt(costheta, phi); |
270 |
+ |
|
271 |
+ |
return sigma; |
272 |
|
} |
273 |
+ |
|
274 |
+ |
double SHAPE::getSAt(double costheta, double phi) { |
275 |
+ |
|
276 |
+ |
vector<SHFunc*>::iterator sIter; |
277 |
+ |
double s; |
278 |
+ |
|
279 |
+ |
s = 0.0; |
280 |
+ |
|
281 |
+ |
for(sIter = sFuncs.begin(); sIter != sFuncs.end(); ++sIter) |
282 |
+ |
s += (*sIter)->getValueAt(costheta, phi); |
283 |
+ |
|
284 |
+ |
return s; |
285 |
+ |
} |
286 |
+ |
|
287 |
+ |
double SHAPE::getEpsAt(double costheta, double phi) { |
288 |
+ |
|
289 |
+ |
vector<SHFunc*>::iterator epsIter; |
290 |
+ |
double eps; |
291 |
+ |
|
292 |
+ |
eps = 0.0; |
293 |
+ |
|
294 |
+ |
for(epsIter = epsFuncs.begin(); epsIter != epsFuncs.end(); ++epsIter) |
295 |
+ |
eps += (*epsIter)->getValueAt(costheta, phi); |
296 |
+ |
|
297 |
+ |
return eps; |
298 |
+ |
} |