Hi Johannes, thanks for this interesting bug report. I tested your program with the most recent svn version and the bug is no longer present. I have no idea what kind of bug this must have been.
Timo On Fri, Oct 7, 2011 at 8:25 AM, Johannes Reinhardt <[email protected]> wrote: > Hello, > > I found that in deal.ii version 6.3.1 make_sparsity_patterns produces > wrong sparsity patterns in some cases. > > I have not checked whether this problem is present in version 7 or in > trunk, but I didn't find it explicitly mentioned in the changelogs. > > Attached find a short piece of code that shows the behaviour. With first > order DGQ elements the resulting sparsity pattern consists of 4x4 blocks > on the diagonal for a up to six times globally refined unit square. When > refining 7 times, the resulting sparsity pattern is only the main > diagonal, which is wrong. With Q elements this does not happen. > > I have not investigated the problem further. It is not a urgent problem > for me, I will just assemble the Sparsity Pattern differently. > > I am using using deal.ii as it is available in Ubuntu 11.04, on a 64bit > machine running Ubuntu 11.04. > > And many thanks to all people involved with this project. It is quite > convenient to work with deal.ii, so well documented, so well structured. > > Greetings > > Johannes > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii > > -- Timo Heister http://www.math.tamu.edu/~heister/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
