pub fn angles(&self) -> f64 {
        let mut energy = 0.0;
        for molecule in self.system.molecules() {
            for angle in molecule.angles() {
                let (i, j, k) = (angle.i(), angle.j(), angle.k());
                let theta = self.system.angle(i, j, k);
                energy += self.angle(theta, i, j, k);
            }
        }
        return energy;
    }
/// Compute the energy of all the angles in the system
Compute the energy of all the angles in the system
Compute the energy of all the angles in the system
pub fn angles(&self) -> f64 {
        let mut energy = 0.0;
        for molecule in self.system.molecules() {
            for angle in molecule.angles() {
                let (i, j, k) = (angle.i(), angle.j(), angle.k());
                let theta = self.system.angle(i, j, k);
                energy += self.angle(theta, i, j, k);
            }
        }
        return energy;
    }
pub fn angles ( & self ) -> f64 { let mut energy = 0.0 ; for molecule in self . system . molecules ( ) { for angle in molecule . angles ( ) { let ( i , j , k ) = ( angle . i ( ) , angle . j ( ) , angle . k ( ) ) ; let theta = self . system . angle ( i , j , k ) ; energy += self . angle ( theta , i , j , k ) ; } } return energy ; }
Compute the energy of all the angles in the system
Compute the energy of all the angles in the system
&self