* of coordinates working on settle type constraint.
*/
void settle_proj(settledata *settled, ConstraintVariable econq,
- int nsettle, t_iatom iatoms[],
+ int nsettle, const t_iatom iatoms[],
const t_pbc *pbc, /* PBC data pointer, can be NULL */
const rvec x[],
rvec *der, rvec *derp,