| 35 |
|
* |
| 36 |
|
* [1] Meineke, et al., J. Comp. Chem. 26, 252-271 (2005). |
| 37 |
|
* [2] Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006). |
| 38 |
< |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008). |
| 38 |
> |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008). |
| 39 |
|
* [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). |
| 40 |
|
* [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). |
| 41 |
|
*/ |
| 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; |