only read part of the matrix for parallel computation (actual parallel computation still missing)