Add rare_matrix.cpp