matrix_tc.cpp