I think matrix_sparce function : transpose_in_place() is wrong.

[xinlnix]

template<class T> matrix_sparse<T> matrix_sparse<T>::transpose_in_place() {
     orientation = other_orientation(orientation);
     std::swap(number_rows, number_columns);
     return *this;
  }

Why transpose in place just swap number rows and columns and change orientation?. I can't find where it has been transpose and thank for your reply in advance.

[xinlnix]

And in practice, CCS format calculate of ILU0 is wrong, hoping your reply.

[c-f-h]

template<class T> matrix_sparse<T> matrix_sparse<T>::transpose_in_place() {
     orientation = other_orientation(orientation);
     std::swap(number_rows, number_columns);
     return *this;
  }

Why transpose in place just swap number rows and columns and change orientation?. I can't find where it has been transpose and thank for your reply in advance.

Turning a row-major into a column-major matrix in effect transposes it. I'm pretty sure it's correct. Have you tried it?

[c-f-h]

And in practice, CCS format calculate of ILU0 is wrong, hoping your reply.

Can you give a reproducible test case where it's wrong?