hp_singular.C
Go to the documentation of this file.
1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2014 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner
3 
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either
7 // version 2.1 of the License, or (at your option) any later version.
8 
9 // This library is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // Lesser General Public License for more details.
13 
14 // You should have received a copy of the GNU Lesser General Public
15 // License along with this library; if not, write to the Free Software
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 
18 // C++ includes
19 
20 // Local Includes
21 #include "libmesh/elem.h"
22 #include "libmesh/hp_singular.h"
24 #include "libmesh/mesh_base.h"
25 #include "libmesh/system.h"
26 
27 #ifdef LIBMESH_ENABLE_AMR
28 
29 namespace libMesh
30 {
31 
32 //-----------------------------------------------------------------
33 // HPSingularity implementations
34 
35 
37 {
38  START_LOG("select_refinement()", "HPSingularity");
39 
40  // The current mesh
41  MeshBase& mesh = system.get_mesh();
42 
45  const MeshBase::element_iterator elem_end =
47 
48  for (; elem_it != elem_end; ++elem_it)
49  {
50  Elem* elem = *elem_it;
51 
52  // We're only checking elements that are already flagged for h
53  // refinement
54  if (elem->refinement_flag() != Elem::REFINE)
55  continue;
56 
59 
60  for (std::list<Point>::iterator ppoint =
61  singular_points.begin();
62  ppoint != singular_points.end(); ++ppoint)
63  {
64  if (elem->contains_point(*ppoint))
65  {
68  break;
69  }
70  }
71  }
72 
73  STOP_LOG("select_refinement()", "HPSingularity");
74 }
75 
76 } // namespace libMesh
77 
78 #endif // #ifdef LIBMESH_ENABLE_AMR

Site Created By: libMesh Developers
Last modified: February 07 2014 16:57:05 UTC

Hosted By:
SourceForge.net Logo