| 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 |
+ |
} |