libMesh
Loading...
Searching...
No Matches
Public Member Functions | Public Attributes | Protected Member Functions | List of all members
libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::ProjectEdges Struct Reference

#include <generic_projector.h>

Inheritance diagram for libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::ProjectEdges:
[legend]

Public Member Functions

 ProjectEdges (GenericProjector &p)
 
 ProjectEdges (ProjectEdges &p_e, Threads::split)
 
void operator() (const node_range &range)
 
template<typename InsertInput , typename std::enable_if< std::is_same< typename ProjectionAction::InsertInput, InsertInput >::value, int >::type = 0>
void insert_id (dof_id_type id, const InsertInput &val, processor_id_type pid)
 
template<typename InsertInput , typename std::enable_if< !std::is_same< typename ProjectionAction::InsertInput, InsertInput >::value, int >::type = 0>
void insert_id (dof_id_type id, const InsertInput &val, processor_id_type pid)
 
template<typename InsertInput , typename std::enable_if< std::is_same< typename ProjectionAction::InsertInput, InsertInput >::value, int >::type = 0>
void insert_ids (const std::vector< dof_id_type > &ids, const std::vector< InsertInput > &vals, processor_id_type pid)
 
template<typename InsertInput , typename std::enable_if< !std::is_same< typename ProjectionAction::InsertInput, InsertInput >::value, int >::type = 0>
void insert_ids (const std::vector< dof_id_type > &ids, const std::vector< InsertInput > &vals, processor_id_type pid)
 
void find_dofs_to_send (const Node &node, const Elem &elem, unsigned short node_num, const var_set &vars)
 
void join (const SubFunctor &other)
 

Public Attributes

ProjectionAction action
 
FFunctor f
 
std::unique_ptr< GFunctor > g
 
const Systemsystem
 
FEMContext context
 
std::vector< FEContinuityconts
 
std::vector< FEFieldTypefield_types
 
GenericProjectorprojector
 
std::unordered_map< dof_id_type, std::pair< typename FFunctor::ValuePushType, processor_id_type > > new_ids_to_push
 
std::unordered_map< dof_id_type, typename FFunctor::ValuePushType > new_ids_to_save
 

Protected Member Functions

void construct_projection (const std::vector< dof_id_type > &dof_indices_var, const std::vector< unsigned int > &involved_dofs, unsigned int var_component, const Node *node, const FEGenericBase< typename FFunctor::RealType > &fe)
 

Detailed Description

template<typename FFunctor, typename GFunctor, typename FValue, typename ProjectionAction>
struct libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::ProjectEdges

Definition at line 327 of file generic_projector.h.

Constructor & Destructor Documentation

◆ ProjectEdges() [1/2]

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::ProjectEdges::ProjectEdges ( GenericProjector p)
inline

Definition at line 328 of file generic_projector.h.

◆ ProjectEdges() [2/2]

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::ProjectEdges::ProjectEdges ( ProjectEdges p_e,
Threads::split   
)
inline

Definition at line 330 of file generic_projector.h.

330: SubProjector(p_e.projector) {}

Member Function Documentation

◆ construct_projection()

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
void libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubProjector::construct_projection ( const std::vector< dof_id_type > &  dof_indices_var,
const std::vector< unsigned int > &  involved_dofs,
unsigned int  var_component,
const Node node,
const FEGenericBase< typename FFunctor::RealType > &  fe 
)
protectedinherited

Definition at line 3003 of file generic_projector.h.

3009{
3010 const auto & JxW = fe.get_JxW();
3011 const auto & phi = fe.get_phi();
3012 const std::vector<std::vector<typename FEGenericBase<typename FFunctor::RealType>::OutputGradient>> * dphi = nullptr;
3013 const std::vector<Point> & xyz_values = fe.get_xyz();
3014 const FEContinuity cont = fe.get_continuity();
3015 const std::unordered_map<dof_id_type, typename FFunctor::ValuePushType> & ids_to_save =
3016 this->projector.ids_to_save;
3017
3018 if (cont == C_ONE)
3019 dphi = &(fe.get_dphi());
3020
3021 const unsigned int n_involved_dofs =
3022 cast_int<unsigned int>(involved_dofs.size());
3023
3024 std::vector<dof_id_type> free_dof_ids;
3025 DenseVector<typename FFunctor::ValuePushType> Uinvolved(n_involved_dofs);
3026 std::vector<char> dof_is_fixed(n_involved_dofs, false); // bools
3027
3028 for (auto i : make_range(n_involved_dofs))
3029 {
3030 const dof_id_type id = dof_indices_var[involved_dofs[i]];
3031 auto iter = ids_to_save.find(id);
3032 if (iter == ids_to_save.end())
3033 free_dof_ids.push_back(id);
3034 else
3035 {
3036 dof_is_fixed[i] = true;
3037 Uinvolved(i) = iter->second;
3038 }
3039 }
3040
3041 const unsigned int free_dofs = free_dof_ids.size();
3042
3043 // There may be nothing to project
3044 if (!free_dofs)
3045 return;
3046
3047 // The element matrix and RHS for projections.
3048 // Note that Ke is always real-valued, whereas
3049 // Fe may be complex valued if complex number
3050 // support is enabled
3051 DenseMatrix<Real> Ke(free_dofs, free_dofs);
3052 DenseVector<typename FFunctor::ValuePushType> Fe(free_dofs);
3053 // The new degree of freedom coefficients to solve for
3054 DenseVector<typename FFunctor::ValuePushType> Ufree(free_dofs);
3055
3056 const unsigned int n_qp =
3057 cast_int<unsigned int>(xyz_values.size());
3058
3059 // Loop over the quadrature points
3060 for (unsigned int qp=0; qp<n_qp; qp++)
3061 {
3062 // solution at the quadrature point
3063 FValue fineval = f.eval_at_point(context,
3064 var_component,
3065 xyz_values[qp],
3066 system.time,
3067 false);
3068 // solution grad at the quadrature point
3069 typename GFunctor::FunctorValue finegrad;
3070 if (cont == C_ONE)
3071 finegrad = g->eval_at_point(context,
3072 var_component,
3073 xyz_values[qp],
3074 system.time,
3075 false);
3076
3077 // Form edge projection matrix
3078 for (unsigned int sidei=0, freei=0;
3079 sidei != n_involved_dofs; ++sidei)
3080 {
3081 unsigned int i = involved_dofs[sidei];
3082 // fixed DoFs aren't test functions
3083 if (dof_is_fixed[sidei])
3084 continue;
3085 for (unsigned int sidej=0, freej=0;
3086 sidej != n_involved_dofs; ++sidej)
3087 {
3088 unsigned int j = involved_dofs[sidej];
3089 if (dof_is_fixed[sidej])
3090 Fe(freei) -= phi[i][qp] * phi[j][qp] *
3091 JxW[qp] * Uinvolved(sidej);
3092 else
3093 Ke(freei,freej) += phi[i][qp] *
3094 phi[j][qp] * JxW[qp];
3095 if (cont == C_ONE)
3096 {
3097 if (dof_is_fixed[sidej])
3098 Fe(freei) -= ( TensorTools::inner_product((*dphi)[i][qp],
3099 (*dphi)[j][qp]) ) *
3100 JxW[qp] * Uinvolved(sidej);
3101 else
3102 Ke(freei,freej) += ( TensorTools::inner_product((*dphi)[i][qp],
3103 (*dphi)[j][qp]) )
3104 * JxW[qp];
3105 }
3106 if (!dof_is_fixed[sidej])
3107 freej++;
3108 }
3109 Fe(freei) += phi[i][qp] * fineval * JxW[qp];
3110 if (cont == C_ONE)
3111 Fe(freei) += (TensorTools::inner_product(finegrad,
3112 (*dphi)[i][qp]) ) *
3113 JxW[qp];
3114 freei++;
3115 }
3116 }
3117
3118 Ke.cholesky_solve(Fe, Ufree);
3119
3120 // Transfer new edge solutions to element
3121 const processor_id_type pid = node ?
3122 node->processor_id() : DofObject::invalid_processor_id;
3123 insert_ids(free_dof_ids, Ufree.get_values(), pid);
3124}
static constexpr processor_id_type invalid_processor_id
An invalid processor_id to distinguish DoFs that have not been assigned to a processor.
Definition dof_object.h:484
std::unordered_map< dof_id_type, typename FFunctor::ValuePushType > ids_to_save
Real time
For time-dependent problems, this is the time t at the beginning of the current timestep.
Definition system.h:1677
std::enable_if< ScalarTraits< T >::value &&ScalarTraits< T2 >::value, typenameCompareTypes< T, T2 >::supertype >::type inner_product(const T &a, const T2 &b)
template class LIBMESH_EXPORT DenseMatrix< Real >
uint8_t dof_id_type
Definition id_types.h:67
uint8_t processor_id_type
Definition id_types.h:104
IntRange< T > make_range(T beg, T end)
The 2-parameter make_range() helper function returns an IntRange<T> when both input parameters are of...
Definition int_range.h:176
void insert_ids(const std::vector< dof_id_type > &ids, const std::vector< InsertInput > &vals, processor_id_type pid)

References libMesh::C_ONE, libMesh::DenseMatrix< T >::cholesky_solve(), libMesh::FEAbstract::get_continuity(), libMesh::FEGenericBase< OutputType >::get_dphi(), libMesh::FEAbstract::get_JxW(), libMesh::FEGenericBase< OutputType >::get_phi(), libMesh::DenseVector< T >::get_values(), libMesh::FEAbstract::get_xyz(), libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::ids_to_save, libMesh::TensorTools::inner_product(), libMesh::DofObject::invalid_processor_id, libMesh::make_range(), libMesh::DofObject::processor_id(), libMesh::DenseVector< T >::size(), libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::system, and libMesh::System::time.

◆ find_dofs_to_send()

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
void libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::find_dofs_to_send ( const Node node,
const Elem elem,
unsigned short  node_num,
const var_set vars 
)
inherited

Definition at line 2844 of file generic_projector.h.

2846{
2847 libmesh_assert (&node == elem.node_ptr(node_num));
2848
2849 // Any ghosted node in our range might have an owner who needs our
2850 // data
2851 const processor_id_type owner = node.processor_id();
2852 if (owner != system.processor_id())
2853 {
2854 const MeshBase & mesh = system.get_mesh();
2855 const DofMap & dof_map = system.get_dof_map();
2856
2857 // But let's check and see if we can be certain the owner can
2858 // compute any or all of its own dof coefficients on that node.
2859 std::vector<dof_id_type> node_dof_ids, patch_dof_ids;
2860 for (const auto & var : vars)
2861 {
2862 const Variable & variable = system.variable(var);
2863
2864 if (!variable.active_on_subdomain(elem.subdomain_id()))
2865 continue;
2866
2867 dof_map.dof_indices(elem, node_num, node_dof_ids, var);
2868 }
2869 libmesh_assert(std::is_sorted(node_dof_ids.begin(),
2870 node_dof_ids.end()));
2871 const std::vector<dof_id_type> & patch =
2872 (*this->projector.nodes_to_elem)[node.id()];
2873 for (const auto & elem_id : patch)
2874 {
2875 const Elem & patch_elem = mesh.elem_ref(elem_id);
2876 if (!patch_elem.active() || owner != patch_elem.processor_id())
2877 continue;
2878 dof_map.dof_indices(&patch_elem, patch_dof_ids);
2879 std::sort(patch_dof_ids.begin(), patch_dof_ids.end());
2880
2881 // Remove any node_dof_ids that we see can be calculated on
2882 // this element
2883 std::vector<dof_id_type> diff_ids(node_dof_ids.size());
2884 auto it = std::set_difference(node_dof_ids.begin(), node_dof_ids.end(),
2885 patch_dof_ids.begin(), patch_dof_ids.end(), diff_ids.begin());
2886 diff_ids.resize(it-diff_ids.begin());
2887 node_dof_ids.swap(diff_ids);
2888 if (node_dof_ids.empty())
2889 break;
2890 }
2891
2892 // Give ids_to_push default invalid pid: not yet computed
2893 for (auto id : node_dof_ids)
2894 new_ids_to_push[id].second = DofObject::invalid_processor_id;
2895 }
2896}
processor_id_type processor_id() const
const Variable & variable(unsigned int var) const
Return a constant reference to Variable var.
Definition system.C:2704
const DofMap & get_dof_map() const
Definition system.h:2417
const MeshBase & get_mesh() const
Definition system.h:2401
MeshBase & mesh
libmesh_assert(ctx)
std::unordered_map< dof_id_type, std::pair< typename FFunctor::ValuePushType, processor_id_type > > new_ids_to_push

References libMesh::Elem::active(), libMesh::Variable::active_on_subdomain(), libMesh::DofMap::dof_indices(), libMesh::System::get_dof_map(), libMesh::System::get_mesh(), libMesh::DofObject::id(), libMesh::DofObject::invalid_processor_id, libMesh::libmesh_assert(), mesh, libMesh::Elem::node_ptr(), libMesh::DofObject::processor_id(), libMesh::ParallelObject::processor_id(), libMesh::Elem::subdomain_id(), libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::system, and libMesh::System::variable().

◆ insert_id() [1/2]

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
template<typename InsertInput , typename std::enable_if< std::is_same< typename ProjectionAction::InsertInput, InsertInput >::value, int >::type = 0>
void libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::insert_id ( dof_id_type  id,
const InsertInput &  val,
processor_id_type  pid 
)

Definition at line 262 of file generic_projector.h.

1380{
1381 libmesh_error_msg("ID insertion should only occur when the provided input matches that "
1382 "expected by the ProjectionAction");
1383}

◆ insert_id() [2/2]

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
template<typename InsertInput , typename std::enable_if< !std::is_same< typename ProjectionAction::InsertInput, InsertInput >::value, int >::type = 0>
void libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::insert_id ( dof_id_type  id,
const InsertInput &  val,
processor_id_type  pid 
)

◆ insert_ids() [1/2]

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
template<typename InsertInput , typename std::enable_if< std::is_same< typename ProjectionAction::InsertInput, InsertInput >::value, int >::type = 0>
void libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::insert_ids ( const std::vector< dof_id_type > &  ids,
const std::vector< InsertInput > &  vals,
processor_id_type  pid 
)

Definition at line 263 of file generic_projector.h.

1424{
1425 libmesh_error_msg("ID insertion should only occur when the provided input matches that "
1426 "expected by the ProjectionAction");
1427}

◆ insert_ids() [2/2]

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
template<typename InsertInput , typename std::enable_if< !std::is_same< typename ProjectionAction::InsertInput, InsertInput >::value, int >::type = 0>
void libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::insert_ids ( const std::vector< dof_id_type > &  ids,
const std::vector< InsertInput > &  vals,
processor_id_type  pid 
)

◆ join()

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
void libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::join ( const SubFunctor other)
inherited

Definition at line 1470 of file generic_projector.h.

1472{
1473 new_ids_to_push.insert(other.new_ids_to_push.begin(), other.new_ids_to_push.end());
1474 new_ids_to_save.insert(other.new_ids_to_save.begin(), other.new_ids_to_save.end());
1475}
std::unordered_map< dof_id_type, typename FFunctor::ValuePushType > new_ids_to_save

References libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::new_ids_to_push, and libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::new_ids_to_save.

Referenced by libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SortAndCopy::join().

◆ operator()()

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
void libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::ProjectEdges::operator() ( const node_range range)

Definition at line 2347 of file generic_projector.h.

2349{
2350 LOG_SCOPE ("project_edges","GenericProjector");
2351
2352 const unsigned int sys_num = system.number();
2353
2354 for (const auto & e_pair : range)
2355 {
2356 const Elem & elem = *std::get<0>(e_pair.second);
2357
2358 // If this is an unchanged element then we already copied all
2359 // its dofs
2360#ifdef LIBMESH_ENABLE_AMR
2361 if (f.is_grid_projection() &&
2362 (elem.refinement_flag() != Elem::JUST_REFINED &&
2363 elem.refinement_flag() != Elem::JUST_COARSENED &&
2364 elem.p_refinement_flag() != Elem::JUST_REFINED &&
2365 elem.p_refinement_flag() != Elem::JUST_COARSENED))
2366 continue;
2367#endif // LIBMESH_ENABLE_AMR
2368
2369 const Node & edge_node = *e_pair.first;
2370 const int dim = elem.dim();
2371 const var_set & edge_vars = std::get<2>(e_pair.second);
2372
2373 const unsigned short edge_num = std::get<1>(e_pair.second);
2374 const unsigned short node_num = elem.n_vertices() + edge_num;
2375 context.edge = edge_num;
2377
2378 this->find_dofs_to_send(edge_node, elem, node_num, edge_vars);
2379
2380 // Look at all the variables we're supposed to interpolate from
2381 // this element on this edge
2382 for (const auto & var : edge_vars)
2383 {
2384 const Variable & variable = system.variable(var);
2385 const FEType & base_fe_type = variable.type();
2386 const unsigned int var_component =
2388
2389 if (base_fe_type.family == SCALAR)
2390 continue;
2391
2392 FEType fe_type = base_fe_type;
2393
2394 // This may be a p refined element
2395 fe_type.order = fe_type.order + elem.p_level();
2396
2397 // If this is a Lagrange element with DoFs on edges then by
2398 // convention we interpolate at the node rather than project
2399 // along the edge.
2400 if (fe_type.family == LAGRANGE || fe_type.family == LAGRANGE_VEC)
2401 {
2402 if (fe_type.order > 1)
2403 {
2404 const processor_id_type pid =
2405 edge_node.processor_id();
2406 FValue fval = f.eval_at_point
2407 (context, var_component, edge_node, system.time,
2408 false);
2409 if (fe_type.family == LAGRANGE_VEC)
2410 {
2411 // We will have a number of nodal value DoFs equal to the elem dim
2412 for (auto i : make_range(elem.dim()))
2413 {
2414 const dof_id_type dof_id =
2415 edge_node.dof_number(sys_num, var, i);
2416
2417 // Need this conversion so that this method
2418 // will compile for TYPE_SCALAR instantiations
2419 const auto insert_val =
2420 raw_value<typename ProjectionAction::InsertInput>(fval, i);
2421
2422 insert_id(dof_id, insert_val, pid);
2423 }
2424 }
2425 else // We are LAGRANGE
2426 {
2427 const dof_id_type dof_id =
2428 edge_node.dof_number(sys_num, var, 0);
2429 insert_id(dof_id, fval, pid);
2430 }
2431 }
2432 continue;
2433 }
2434
2435#ifdef LIBMESH_ENABLE_AMR
2436 // If this is a low order monomial element which has merely
2437 // been h refined then we already copied all its dofs
2438 if (fe_type.family == MONOMIAL &&
2439 fe_type.order == CONSTANT &&
2440 elem.refinement_flag() != Elem::JUST_COARSENED &&
2441 elem.p_refinement_flag() != Elem::JUST_COARSENED)
2442 continue;
2443#endif // LIBMESH_ENABLE_AMR
2444
2445 // FIXME: Need to generalize this to vector-valued elements. [PB]
2446 FEGenericBase<typename FFunctor::RealType> * fe = nullptr;
2447 context.get_element_fe( var, fe, dim );
2448 FEGenericBase<typename FFunctor::RealType> * edge_fe = nullptr;
2449 context.get_edge_fe( var, edge_fe );
2450
2451 // If we're JUST_COARSENED we'll need a custom
2452 // evaluation, not just the standard edge FE
2453 const FEGenericBase<typename FFunctor::RealType> & proj_fe =
2454#ifdef LIBMESH_ENABLE_AMR
2455 (elem.refinement_flag() == Elem::JUST_COARSENED) ?
2456 *fe :
2457#endif
2458 *edge_fe;
2459
2460#ifdef LIBMESH_ENABLE_AMR
2461 if (elem.refinement_flag() == Elem::JUST_COARSENED)
2462 {
2463 std::vector<Point> fine_points;
2464
2465 std::unique_ptr<FEGenericBase<typename FFunctor::RealType>> fine_fe
2467
2468 std::unique_ptr<QBase> qrule
2469 (base_fe_type.default_quadrature_rule(1));
2470 fine_fe->attach_quadrature_rule(qrule.get());
2471
2472 const std::vector<Point> & child_xyz =
2473 fine_fe->get_xyz();
2474
2475 for (unsigned int c = 0, nc = elem.n_children();
2476 c != nc; ++c)
2477 {
2478 if (!elem.is_child_on_edge(c, context.edge))
2479 continue;
2480
2481 fine_fe->edge_reinit(elem.child_ptr(c), context.edge);
2482 fine_points.insert(fine_points.end(),
2483 child_xyz.begin(),
2484 child_xyz.end());
2485 }
2486
2487 std::vector<Point> fine_qp;
2488 FEMap::inverse_map (dim, &elem, fine_points, fine_qp);
2489
2490 context.elem_fe_reinit(&fine_qp);
2491 }
2492 else
2493#endif // LIBMESH_ENABLE_AMR
2495
2496 const std::vector<dof_id_type> & dof_indices =
2498
2499 std::vector<unsigned int> edge_dofs;
2500 FEInterface::dofs_on_edge(&elem, dim, base_fe_type,
2501 context.edge, edge_dofs);
2502
2504 (dof_indices, edge_dofs, var_component,
2505 &edge_node, proj_fe);
2506 }
2507 }
2508}
unsigned int dim
const std::vector< dof_id_type > & get_dof_indices() const
Accessor for element dof indices.
@ JUST_COARSENED
Definition elem.h:1450
static std::unique_ptr< FEGenericBase > build(const unsigned int dim, const FEType &type)
Builds a specific finite element type.
static void dofs_on_edge(const Elem *const elem, const unsigned int dim, const FEType &fe_t, unsigned int e, std::vector< unsigned int > &di, const bool add_p_level=true)
Fills the vector di with the local degree of freedom indices associated with edge e of element elem A...
unsigned char edge
Current edge for edge_* to examine.
void get_edge_fe(unsigned int var, FEGenericBase< OutputShape > *&fe) const
Accessor for edge (3D only!) finite element object for variable var.
virtual void pre_fe_reinit(const System &, const Elem *e)
Reinitializes local data vectors/matrices on the current geometric element.
virtual void elem_fe_reinit(const std::vector< Point > *const pts=nullptr)
Reinitializes interior FE objects on the current geometric element.
void get_element_fe(unsigned int var, FEGenericBase< OutputShape > *&fe) const
Accessor for interior finite element object for variable var for the largest dimension in the mesh.
virtual void edge_fe_reinit()
Reinitializes edge FE objects on the current geometric element.
static Point inverse_map(const unsigned int dim, const Elem *elem, const Point &p, const Real tolerance=TOLERANCE, const bool secure=true, const bool extra_checks=true)
Definition fe_map.C:1512
std::set< unsigned int > var_set
unsigned int variable_scalar_number(std::string_view var, unsigned int component) const
Definition system.h:2474
unsigned int number() const
Definition system.h:2393
const FEType & type() const
Definition variable.h:144
void insert_id(dof_id_type id, const InsertInput &val, processor_id_type pid)
void find_dofs_to_send(const Node &node, const Elem &elem, unsigned short node_num, const var_set &vars)
void construct_projection(const std::vector< dof_id_type > &dof_indices_var, const std::vector< unsigned int > &involved_dofs, unsigned int var_component, const Node *node, const FEGenericBase< typename FFunctor::RealType > &fe)

References libMesh::Elem::child_ptr(), libMesh::CONSTANT, libMesh::FEType::default_quadrature_rule(), dim, libMesh::Elem::dim(), libMesh::DofObject::dof_number(), libMesh::FEInterface::dofs_on_edge(), libMesh::FEType::family, libMesh::FEMap::inverse_map(), libMesh::Elem::is_child_on_edge(), libMesh::Elem::JUST_COARSENED, libMesh::Elem::JUST_REFINED, libMesh::LAGRANGE, libMesh::LAGRANGE_VEC, libMesh::make_range(), libMesh::MONOMIAL, libMesh::Elem::n_children(), libMesh::Elem::n_vertices(), libMesh::System::number(), libMesh::FEType::order, libMesh::Elem::p_level(), libMesh::Elem::p_refinement_flag(), libMesh::DofObject::processor_id(), libMesh::Elem::refinement_flag(), libMesh::SCALAR, libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::system, libMesh::System::time, libMesh::Variable::type(), libMesh::System::variable(), and libMesh::System::variable_scalar_number().

Member Data Documentation

◆ action

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
ProjectionAction libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::action

Definition at line 256 of file generic_projector.h.

◆ context

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
FEMContext libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::context

Definition at line 259 of file generic_projector.h.

◆ conts

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
std::vector<FEContinuity> libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::conts

Definition at line 260 of file generic_projector.h.

◆ f

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
FFunctor libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::f

Definition at line 257 of file generic_projector.h.

◆ field_types

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
std::vector<FEFieldType> libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::field_types

Definition at line 240 of file generic_projector.h.

◆ g

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
std::unique_ptr<GFunctor> libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubProjector::g

Definition at line 267 of file generic_projector.h.

◆ new_ids_to_push

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
std::unordered_map<dof_id_type, std::pair<typename FFunctor::ValuePushType, processor_id_type> > libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::new_ids_to_push
inherited

◆ new_ids_to_save

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
std::unordered_map<dof_id_type, typename FFunctor::ValuePushType> libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::new_ids_to_save
inherited

◆ projector

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
GenericProjector& libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::projector
inherited

◆ system

template<typename FFunctor , typename GFunctor , typename FValue , typename ProjectionAction >
const System& libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::system

Definition at line 258 of file generic_projector.h.


The documentation for this struct was generated from the following file: