libMesh::OldSolutionBase< Output, point_output > Class Template Reference

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 >:

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 ()
template<>
void	get_shape_outputs (FEAbstract &fe)
template<>
void	get_shape_outputs (FEAbstract &fe)
template<>
void	get_shape_outputs (FEAbstract &fe)
template<>
void	get_shape_outputs (FEAbstract &fe)
template<>
void	get_shape_outputs (FEAbstract &fe)
template<>
void	get_shape_outputs (FEAbstract &fe)
template<>
void	get_shape_outputs (FEAbstract &fe)
template<>
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 524 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 527 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 529 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 531 of file generic_projector.h.

References libMesh::make_range().

                                                        :
    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);
  }

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 555 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 596 of file generic_projector.h.

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().

  {
    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());
      }
  }

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 629 of file generic_projector.h.

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().

  {
    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;
  }

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]

template<>

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

inline

Definition at line 944 of file generic_projector.h.

References libMesh::FEAbstract::request_phi().

{
  fe.request_phi();
}

get_shape_outputs() [3/9]

template<>

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

inline

Definition at line 952 of file generic_projector.h.

References libMesh::FEAbstract::request_dphi().

{
  fe.request_dphi();
}

get_shape_outputs() [4/9]

template<>

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

inline

Definition at line 959 of file generic_projector.h.

References libMesh::FEAbstract::request_phi().

{
  fe.request_phi();
}

get_shape_outputs() [5/9]

template<>

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

inline

Definition at line 967 of file generic_projector.h.

References libMesh::FEAbstract::request_dphi().

{
  fe.request_dphi();
}

get_shape_outputs() [6/9]

template<>

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

inline

Definition at line 976 of file generic_projector.h.

References libMesh::FEAbstract::request_phi().

{
  fe.request_phi();
}

get_shape_outputs() [7/9]

template<>

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

inline

Definition at line 984 of file generic_projector.h.

References libMesh::FEAbstract::request_dphi().

{
  fe.request_dphi();
}

get_shape_outputs() [8/9]

template<>

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

inline

Definition at line 991 of file generic_projector.h.

References libMesh::FEAbstract::request_phi().

{
  fe.request_phi();
}

get_shape_outputs() [9/9]

template<>

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

inline

Definition at line 999 of file generic_projector.h.

References libMesh::FEAbstract::request_dphi().

{
  fe.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 567 of file generic_projector.h.

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

  {
    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);
            }
      }
  }

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 593 of file generic_projector.h.

{ 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 700 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 697 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 699 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 702 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 698 of file generic_projector.h.

---

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

• include/systems/generic_projector.h