Dear Josh,
> But I think it should instead be > > for (unsigned int i=0; i<my_nonzeros; ++i) > values[i] = in_values[i]; > > That is, values[i] and in_values[i] should be walked through in tandem, > rather than just assigning the diagonal entry in_values[0] throughout. You are absolutely right, this was a bug in the implementation. I've prepared a fix that I will check in as soon as the deal.II svn permits me to do so. Can you please try again and report if there are any further issues? Thanks for pointing this out! Best, Martin _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
