Performance Improvement
-- Add a vec_mult_matrix function for much efficient P*U and Q*V
-- Optimize the computation process to reduce I/O overhead