| 50 |
|
using namespace std; |
| 51 |
|
namespace OpenMD { |
| 52 |
|
|
| 53 |
< |
SwitchingFunction::SwitchingFunction() : np_(150), haveSpline_(false), |
| 54 |
< |
isCubic_(true), functionType_(cubic) { |
| 53 |
> |
SwitchingFunction::SwitchingFunction() : functionType_(cubic), |
| 54 |
> |
haveSpline_(false), isCubic_(true), |
| 55 |
> |
np_(150) { |
| 56 |
|
switchSpline_ = new CubicSpline(); |
| 57 |
|
} |
| 58 |
|
|
| 160 |
|
} else { |
| 161 |
|
in_switching_region = true; |
| 162 |
|
r = sqrt(r2); |
| 163 |
< |
pair<RealType, RealType> result = switchSpline_->getValueAndDerivativeAt(r); |
| 163 |
< |
sw = result.first; |
| 164 |
< |
dswdr = result.second; |
| 163 |
> |
switchSpline_->getValueAndDerivativeAt(r, sw, dswdr); |
| 164 |
|
} |
| 165 |
|
} |
| 166 |
|
return in_switching_region; |