On Mon, 25 Aug 2014, Miguel Angel Salazar de Troya wrote: > I solved this problem, I was calling DofMap::_dof_indices with the wrong > element, it was an inactive element. That's why it didn't have the right > indices. Once I call it when an active element, it works. My bad.
Glad you caught this! Sorry I didn't have time to help debug, but that looked like a nightmare. You say "My bad", but that really does look like the sort of user bug we ought to be catching via an assertion. We can't actually assert(elem->active()) in dof_indices, though, because we call it on parent elements to speed up compute_constraints on Lagrange FE. No, wait... we call it on parent element *side* proxies, which ought to return "true" to elem->active(). I'll put in a PR for that assertion; probably don't want to merge it until MooseBuild is back up and running though. --- Roy ------------------------------------------------------------------------------ Slashdot TV. Video for Nerds. Stuff that matters. http://tv.slashdot.org/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
