so now there is no while loop in the compute_modified_matrices, but we still have no coupling terms: everything is done in one loop on the IGA side