Hi Timo Thanks, that solved the problem.
Gennadiy On 26 August 2011 18:12, Timo Heister <[email protected]>wrote: > > 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
