On Thu, 8 May 2008, Bin wrote:
> My code is developed based on libmesh. Now I want to remove the elements
> which have large plastic strain in the process of computing.
> Can anyone tell me what is the easy way to do it?
If you're not using adaptive refinement, I think you can get away with
first calling MeshBase::delete_elem() on each of the elements you want
deleted, then MeshBase::prepare_for_use() to do some internal cleanup
of things like neighbor links, then EquationSystems::reinit() to
properly resize vectors.
If you're using adaptive refinement, things are harder. I've got a
small library patch with changes Victor Velten had to make to get that
working correctly for him, but IIRC he didn't get back to me about a
couple minor problems, and I never got it into SVN. You might want to
check through the mailing list archives on that.