Function calculate_u_1

pub fn calculate_u_1(
    b: &[Vec<Vec<PolyVector>>],
    c: &[Vec<Vec<Vec<PolyVector>>>],
    t_i: &[Vec<PolyVector>],
    g_ij: &[Vec<PolyVector>],
    ep: &EnvironmentParameters,
) -> PolyVector

line 19, page 18: calculate u_1 = \sum(B_ik * t_i^(k)) + sum(C_ijk * g_ij^(k))