The primitive `grb::set( Matrix<D>&, const Matrix<T> )` contains a bug related to the cast of the values from T to D.