Dear Dune, I have a question about BCRS row iterator and column iterator. Is it by any chance the case that the row iterator and the column iterator accesses the matrix entries strictly in ascending local row / column index? If this is the case, I could get rid of some extra temporary storage which is otherwise unnecessary. Regards, Aleksejs