matrix_compute
