On Tue, May 11, 2010 at 3:49 PM, Bangerth, Wolfgang
<[email protected]> wrote:
>
> James,
>
>> http://diku.dk/~avery/aspect/tria_accessor.cc
>> http://diku.dk/~avery/aspect/tria_accessor.h
>
> Outstanding! I've tweaked a couple things but it's now
> in the svn repository. Thanks a lot!
Great! I'm happy to have added my first tiny bit. :)
>> Also, things are apparently not completely resolved: Things now hang
>> in hanging_node_constraints.close() after
>> DoFTools::make_hanging_node_constraints (dof_handler,
>> hanging_node_constraints) -- any ideas about what's going wrong here?
>
> What do you mean with "it hangs"?
Sorry, that wasn't very clear. What happens is: When I do anisotropic
refinement, I run into trouble the next time I update the hanging node
constraints as
hanging_node_constraints.clear ();
DoFTools::make_hanging_node_constraints(dof_handler,
hanging_node_constraints);
hanging_node_constraints.close ();
hanging_node_constraints.condense (sparsity_pattern);
I never get out from hanging_node_constraints.close (), and the memory
consumption explodes. The back trace is
#0 0x00002afc82953879 in std::vector<std::pair<unsigned int, double>,
std::allocator<std::pair<unsigned int, double> > >::_M_insert_aux ()
from /build/deal.II/lib/libdeal_II_2d.so
#1 0x00002afc8436706e in dealii::ConstraintMatrix::close ()
from /build/deal.II/lib/liblac.so
#2 0x00002afc819d4522 in FESpace<3>::update ()
from /work/lib/libspace.so.0
I'll have to compile deal.II with debug symbols to give more
information. I'll also try to reproduce the problem with a small
example program -- I'll do both tomorrow.
Cheers,
--
Med venlig hilsen,
James Avery <[email protected]>
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
