Function orthonormal_subspace_basis 

pub fn orthonormal_subspace_basis<F>(vs: &[Vector], f: F)
where
    F: FnMut(Vector) -> bool,

Computes an orthonormal basis for the subspace orthogonal to the given vectors. Calls the callback f with each basis vector.

For 3D: given 1 vector, produces 2 orthonormal vectors perpendicular to it.