| # | Line 177 | Line 177 | template<typename T> void Integrator<T>::integrate(voi | |
|---|---|---|
| 177 | // initialize the forces before the first step | |
| 178 | ||
| 179 | calcForce(1, 1); | |
| 180 | + | |
| 181 | + | if (nConstrained){ |
| 182 | + | preMove(); |
| 183 | + | constrainA(); |
| 184 | + | calcForce(1, 1); |
| 185 | + | constrainB(); |
| 186 | + | } |
| 187 | ||
| 188 | if (info->setTemp){ | |
| 189 | thermalize(); | |
| – | Removed lines |
| + | Added lines |
| < | Changed lines |
| > | Changed lines |