Hi! I'm trying to implement XFEM method for elasticity problem with voids. I use FE_Collection: fe_collection.push_back (FESystem<dim> (FE_Q<dim>(1), 1, FE_Nothing<dim>(), 1)); fe_collection.push_back (FESystem<dim> (FE_Q<dim>(1), 1, FE_Q<dim>(1), 1));
Triangulation remains always unchanged, I do not cut voids in the mesh. Instead I use enriching function V(x) which is 0 inside voids and 1 elsewhere. So I get a singular matrix with zero degrees of freedom inside voids. For this system to be solved I should eliminate rows corresponding to zero degrees of freedom from my system matrix. I need your advise how to implement it in dealii in an optimal way. I do not know beforehand how many rows are to be removed so initially I allocate my sparsity pattern as: DynamicSparsityPattern dsp(dof_handler.n_dofs(), dof_handler.n_dofs()); Or maybe there are some other techniques to deal with this kind of problem that I'm not aware of. I would appreciate any suggestion on how it could be solved. Alex -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
