| 58 |
|
#include "primitives/Inversion.hpp" |
| 59 |
|
#include "nonbonded/NonBondedInteraction.hpp" |
| 60 |
|
#include "perturbations/UniformField.hpp" |
| 61 |
+ |
#include "perturbations/UniformGradient.hpp" |
| 62 |
|
#include "parallel/ForceMatrixDecomposition.hpp" |
| 63 |
|
|
| 64 |
|
#include <cstdio> |
| 450 |
|
UniformField* eField = new UniformField(info_); |
| 451 |
|
perturbations_.push_back(eField); |
| 452 |
|
} |
| 453 |
< |
|
| 453 |
> |
if (info_->getSimParams()->haveUniformGradientStrength() || |
| 454 |
> |
info_->getSimParams()->haveUniformGradientDirection1() || |
| 455 |
> |
info_->getSimParams()->haveUniformGradientDirection2() ) { |
| 456 |
> |
UniformGradient* eGrad = new UniformGradient(info_); |
| 457 |
> |
perturbations_.push_back(eGrad); |
| 458 |
> |
} |
| 459 |
> |
|
| 460 |
|
usePeriodicBoundaryConditions_ = info_->getSimParams()->getUsePeriodicBoundaryConditions(); |
| 461 |
|
|
| 462 |
|
fDecomp_->distributeInitialData(); |
