> for (unsigned int i = 0; i < system_matrix.m(); ++i) {
> SparsityPattern::row_iterator rowBegin =
> sparsity_pattern.row_begin(i);
> SparsityPattern::row_iterator rowEnd = sparsity_pattern.row_end(i);
>
> for ( ; rowBegin != rowEnd; ++rowBegin) {
> j = system_matrix.get_sparsity_pattern().column_number(i,
> *rowBegin);
column_number wants an index between 0 and row_length(i). But
*rowBegin is already the column index you are looking for. So replace
the line with
j = *rowBegin;
--
Timo Heister
http://www.math.tamu.edu/~heister/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
