The OldSolutionBase input functor abstract base class is the root of the OldSolutionValue and OldSolutionCoefs classes which allow a GenericProjector to read old solution values or solution interpolation coefficients for a just-refined-and-coarsened mesh. More...

#include <generic_projector.h>

Inheritance diagram for libMesh::OldSolutionBase< Output, point_output >:

[legend]

Public Types
typedef TensorTools::MakeReal< Output >::type	RealType

Public Member Functions
	OldSolutionBase (const libMesh::System &sys_in, const std::vector< unsigned int > *vars)

	OldSolutionBase (const OldSolutionBase &in)

void	init_context (FEMContext &c)

bool	is_grid_projection ()

void	get_shape_outputs (FEAbstract &fe)

void	get_shape_outputs (FEAbstract &fe)

void	get_shape_outputs (FEAbstract &fe)

void	get_shape_outputs (FEAbstract &fe)

void	get_shape_outputs (FEAbstract &fe)

void	get_shape_outputs (FEAbstract &fe)

void	get_shape_outputs (FEAbstract &fe)

void	get_shape_outputs (FEAbstract &fe)

Static Public Member Functions
static void	get_shape_outputs (FEAbstract &fe)

Protected Types
typedef TensorTools::MakeBaseNumber< Output >::type	DofValueType

Protected Member Functions
void	check_old_context (const FEMContext &c)

bool	check_old_context (const FEMContext &c, const Point &p)

Protected Attributes
const Elem *	last_elem

const System &	sys

FEMContext	old_context

std::vector< unsigned int >	component_to_var

Static Protected Attributes
static const Real	out_of_elem_tol = 10 * TOLERANCE

Detailed Description

template<typename Output, void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>
class libMesh::OldSolutionBase< Output, point_output >

The OldSolutionBase input functor abstract base class is the root of the OldSolutionValue and OldSolutionCoefs classes which allow a GenericProjector to read old solution values or solution interpolation coefficients for a just-refined-and-coarsened mesh.

Author: Roy H. Stogner

Date: 2016

Definition at line 537 of file generic_projector.h.

Member Typedef Documentation

◆ DofValueType

template<typename Output , void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>

typedef TensorTools::MakeBaseNumber<Output>::type libMesh::OldSolutionBase< Output, point_output >::DofValueType

protected

Definition at line 540 of file generic_projector.h.

◆ RealType

template<typename Output , void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>

typedef TensorTools::MakeReal<Output>::type libMesh::OldSolutionBase< Output, point_output >::RealType

Definition at line 542 of file generic_projector.h.

Constructor & Destructor Documentation

◆ OldSolutionBase() [1/2]

template<typename Output , void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>

libMesh::OldSolutionBase< Output, point_output >::OldSolutionBase	(	const libMesh::System &	sys_in,
		const std::vector< unsigned int > *	vars
	)

inline

Definition at line 544 of file generic_projector.h.

                                                        :
    last_elem(nullptr),
    sys(sys_in),
    old_context(sys_in, vars, /* allocate_local_matrices= */ false)
  {
    // We'll be queried for components but we'll typically be looking
    // up data by variables, and those indices don't always match
    auto make_lookup = [this](unsigned int v)
      {
        const unsigned int vcomp = sys.variable_scalar_number(v,0);
        if (vcomp >= component_to_var.size())
          component_to_var.resize(vcomp+1, static_cast<unsigned int>(-1));
        component_to_var[vcomp] = v;
      };
 
    if (vars)
      for (auto v : *vars)
        make_lookup(v);
    else
      for (auto v : make_range(sys.n_vars()))
        make_lookup(v);
  }

References libMesh::make_range().

◆ OldSolutionBase() [2/2]

template<typename Output , void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>

libMesh::OldSolutionBase< Output, point_output >::OldSolutionBase ( const OldSolutionBase< Output, point_output > & in )

inline

Definition at line 568 of file generic_projector.h.

                                              :
    last_elem(nullptr),
    sys(in.sys),
    old_context(sys, in.old_context.active_vars(), /* allocate_local_matrices= */ false),
    component_to_var(in.component_to_var)
  {
  }

Member Function Documentation

◆ check_old_context() [1/2]

template<typename Output , void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>

void libMesh::OldSolutionBase< Output, point_output >::check_old_context ( const FEMContext & c )

inlineprotected

Definition at line 609 of file generic_projector.h.

  {
    LOG_SCOPE ("check_old_context(c)", "OldSolutionBase");
    const Elem & elem = c.get_elem();
    if (last_elem != &elem)
      {
        if (elem.refinement_flag() == Elem::JUST_REFINED)
          {
            old_context.pre_fe_reinit(sys, elem.parent());
          }
        else if (elem.refinement_flag() == Elem::JUST_COARSENED)
          {
            libmesh_error();
          }
        else
          {
            if (!elem.get_old_dof_object())
              {
                libmesh_error();
              }
 
            old_context.pre_fe_reinit(sys, &elem);
          }
 
        last_elem = &elem;
      }
    else
      {
        libmesh_assert(old_context.has_elem());
      }
  }

References libMesh::FEMContext::get_elem(), libMesh::DofObject::get_old_dof_object(), libMesh::Elem::JUST_COARSENED, libMesh::Elem::JUST_REFINED, libMesh::libmesh_assert(), libMesh::Elem::parent(), and libMesh::Elem::refinement_flag().

◆ check_old_context() [2/2]

template<typename Output , void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>

bool libMesh::OldSolutionBase< Output, point_output >::check_old_context	(	const FEMContext &	c,
		const Point &	p
	)

inlineprotected

Definition at line 642 of file generic_projector.h.

  {
    LOG_SCOPE ("check_old_context(c,p)", "OldSolutionBase");
    const Elem & elem = c.get_elem();
    if (last_elem != &elem)
      {
        if (elem.refinement_flag() == Elem::JUST_REFINED)
          {
            old_context.pre_fe_reinit(sys, elem.parent());
          }
        else if (elem.refinement_flag() == Elem::JUST_COARSENED)
          {
            // Find the child with this point.  Use out_of_elem_tol
            // (in physical space, which may correspond to a large
            // tolerance in master space!) to allow for out-of-element
            // finite differencing of mixed gradient terms.  Pray we
            // have no quadrature locations which are within 1e-5 of
            // the element subdivision boundary but are not exactly on
            // that boundary.
            const Real master_tol = out_of_elem_tol / elem.hmax() * 2;
 
            for (auto & child : elem.child_ref_range())
              if (child.close_to_point(p, master_tol))
                {
                  old_context.pre_fe_reinit(sys, &child);
                  break;
                }
 
            libmesh_assert
              (old_context.get_elem().close_to_point(p, master_tol));
          }
        else
          {
            if (!elem.get_old_dof_object())
              return false;
 
            old_context.pre_fe_reinit(sys, &elem);
          }
 
        last_elem = &elem;
      }
    else
      {
        libmesh_assert(old_context.has_elem());
 
        const Real master_tol = out_of_elem_tol / elem.hmax() * 2;
 
        if (!old_context.get_elem().close_to_point(p, master_tol))
          {
            libmesh_assert_equal_to
              (elem.refinement_flag(), Elem::JUST_COARSENED);
 
            for (auto & child : elem.child_ref_range())
              if (child.close_to_point(p, master_tol))
                {
                  old_context.pre_fe_reinit(sys, &child);
                  break;
                }
 
            libmesh_assert
              (old_context.get_elem().close_to_point(p, master_tol));
          }
      }
 
    return true;
  }

References libMesh::Elem::child_ref_range(), libMesh::FEMContext::get_elem(), libMesh::DofObject::get_old_dof_object(), libMesh::Elem::hmax(), libMesh::Elem::JUST_COARSENED, libMesh::Elem::JUST_REFINED, libMesh::libmesh_assert(), libMesh::Elem::parent(), libMesh::Real, and libMesh::Elem::refinement_flag().

◆ get_shape_outputs() [1/9]

template<typename Output , void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>

static void libMesh::OldSolutionBase< Output, point_output >::get_shape_outputs ( FEAbstract & fe )

static

◆ get_shape_outputs() [2/9]

void libMesh::OldSolutionBase< Number,&FEMContext::point_value >::get_shape_outputs ( FEAbstract & fe )

inline

Definition at line 957 of file generic_projector.h.

{
  fe.request_phi();
}

References libMesh::FEAbstract::request_phi().

◆ get_shape_outputs() [3/9]

void libMesh::OldSolutionBase< Gradient,&FEMContext::point_gradient >::get_shape_outputs ( FEAbstract & fe )

inline

Definition at line 965 of file generic_projector.h.

{
  fe.request_dphi();
}

References libMesh::FEAbstract::request_dphi().

◆ get_shape_outputs() [4/9]

void libMesh::OldSolutionBase< Gradient,&FEMContext::point_value >::get_shape_outputs ( FEAbstract & fe )

inline

Definition at line 972 of file generic_projector.h.

{
  fe.request_phi();
}

References libMesh::FEAbstract::request_phi().

◆ get_shape_outputs() [5/9]

void libMesh::OldSolutionBase< Tensor,&FEMContext::point_gradient >::get_shape_outputs ( FEAbstract & fe )

inline

Definition at line 980 of file generic_projector.h.

{
  fe.request_dphi();
}

References libMesh::FEAbstract::request_dphi().

◆ get_shape_outputs() [6/9]

void libMesh::OldSolutionBase< Real,&FEMContext::point_value >::get_shape_outputs ( FEAbstract & fe )

inline

Definition at line 989 of file generic_projector.h.

{
  fe.request_phi();
}

References libMesh::FEAbstract::request_phi().

◆ get_shape_outputs() [7/9]

void libMesh::OldSolutionBase< RealGradient,&FEMContext::point_gradient >::get_shape_outputs ( FEAbstract & fe )

inline

Definition at line 997 of file generic_projector.h.

{
  fe.request_dphi();
}

References libMesh::FEAbstract::request_dphi().

◆ get_shape_outputs() [8/9]

void libMesh::OldSolutionBase< RealGradient,&FEMContext::point_value >::get_shape_outputs ( FEAbstract & fe )

inline

Definition at line 1004 of file generic_projector.h.

{
  fe.request_phi();
}

References libMesh::FEAbstract::request_phi().

◆ get_shape_outputs() [9/9]

void libMesh::OldSolutionBase< RealTensor,&FEMContext::point_gradient >::get_shape_outputs ( FEAbstract & fe )

inline

Definition at line 1012 of file generic_projector.h.

{
  fe.request_dphi();
}

References libMesh::FEAbstract::request_dphi().

◆ init_context()

template<typename Output , void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>

void libMesh::OldSolutionBase< Output, point_output >::init_context ( FEMContext & c )

inline

Definition at line 580 of file generic_projector.h.

  {
    c.set_algebraic_type(FEMContext::DOFS_ONLY);
 
    const std::set<unsigned char> & elem_dims =
      old_context.elem_dimensions();
 
    // Loop over variables and dimensions, to prerequest
    for (const auto & dim : elem_dims)
      {
        FEAbstract * fe = nullptr;
        if (old_context.active_vars())
          for (auto var : *old_context.active_vars())
            {
              old_context.get_element_fe(var, fe, dim);
              get_shape_outputs(*fe);
            }
        else
          for (auto var : make_range(sys.n_vars()))
            {
              old_context.get_element_fe(var, fe, dim);
              get_shape_outputs(*fe);
            }
      }
  }

References dim, libMesh::FEMContext::DOFS_ONLY, libMesh::make_range(), and libMesh::FEMContext::set_algebraic_type().

◆ is_grid_projection()

template<typename Output , void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>

bool libMesh::OldSolutionBase< Output, point_output >::is_grid_projection ( )

inline

Definition at line 606 of file generic_projector.h.

606{ return true; }

Member Data Documentation

◆ component_to_var

template<typename Output , void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>

std::vector<unsigned int> libMesh::OldSolutionBase< Output, point_output >::component_to_var

protected

Definition at line 713 of file generic_projector.h.

◆ last_elem

template<typename Output , void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>

const Elem* libMesh::OldSolutionBase< Output, point_output >::last_elem

protected

Definition at line 710 of file generic_projector.h.

◆ old_context

template<typename Output , void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>

FEMContext libMesh::OldSolutionBase< Output, point_output >::old_context

protected

Definition at line 712 of file generic_projector.h.

◆ out_of_elem_tol

template<typename Output , void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>

const Real libMesh::OldSolutionBase< Output, point_output >::out_of_elem_tol = 10 * TOLERANCE

staticprotected

Definition at line 715 of file generic_projector.h.

◆ sys

template<typename Output , void(FEMContext::*)(unsigned int, const Point &, Output &, const Real) const point_output>

const System& libMesh::OldSolutionBase< Output, point_output >::sys

protected

Definition at line 711 of file generic_projector.h.

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

include/systems/generic_projector.h

Public Types

Public Member Functions

Static Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Member Typedef Documentation

◆ DofValueType

◆ RealType

Constructor & Destructor Documentation

◆ OldSolutionBase() [1/2]

◆ OldSolutionBase() [2/2]

Member Function Documentation

◆ check_old_context() [1/2]

◆ check_old_context() [2/2]

◆ get_shape_outputs() [1/9]

◆ get_shape_outputs() [2/9]

◆ get_shape_outputs() [3/9]

◆ get_shape_outputs() [4/9]

◆ get_shape_outputs() [5/9]

◆ get_shape_outputs() [6/9]

◆ get_shape_outputs() [7/9]

◆ get_shape_outputs() [8/9]

◆ get_shape_outputs() [9/9]

◆ init_context()

◆ is_grid_projection()

Member Data Documentation

◆ component_to_var

◆ last_elem

◆ old_context

◆ out_of_elem_tol

◆ sys