Function geometric_slerp_multiple 

pub fn geometric_slerp_multiple(
    a: Vec3A,
    b: Vec3A,
    indices: &[u32],
    points: &mut [Vec3A],
)

This is an optimization for the case where multiple points require the calculation of varying values of p for the same start and end points.

See the intended use in BaseShape::interpolate_multiple.

Note: a and b should both be normalized for normalized results.