Function double_cross_matrix 

pub fn double_cross_matrix(a: &Vec3) -> Mat3

Double cross product a × (a × b) expressed as a matrix times b: the matrix is -[a]ײ = a*aᵀ - (a·a)*I.