| # | 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 |