real(kind=dp), dimension(3,current_cell%max_num_ions_in_species,
current_cell%num_species,num_ionic_constraints) :: ionic_constraints

The constraints on the ionic motion as an array of coefficients in the sum of coordinates which specify the linear constraint. See Section 2.1.1 for details.