pub fn calculate_u_2( d: &[Vec<Vec<Vec<RqVector>>>], h_ij: &[Vec<RqVector>], ep: &EnvironmentParameters, ) -> RqVector
line 20, page 18: calculate u_2 = \sum(D_ijk * h_ij^(k))