more parallelized! pulled RHS summation into parallel loop, eliminated SumRHSMatrix function