Public Member Functions | List of all members
oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT > Class Template Reference

Axisymmetric poro elasticity upgraded to become projectable. More...

#include <axisym_poroelasticity_elements.h>

+ Inheritance diagram for oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >:

Public Member Functions

 ProjectableAxisymmetricPoroelasticityElement ()
 
Vector< std::pair< Data *, unsigned > > data_values_of_field (const unsigned &fld)
 
unsigned nfields_for_projection ()
 
unsigned nhistory_values_for_projection (const unsigned &fld)
 
unsigned nhistory_values_for_coordinate_projection ()
 
double jacobian_and_shape_of_field (const unsigned &fld, const Vector< double > &s, Shape &psi)
 
double get_field (const unsigned &t, const unsigned &fld, const Vector< double > &s)
 
unsigned nvalue_of_field (const unsigned &fld)
 Return number of values in field fld. More...
 
int local_equation (const unsigned &fld, const unsigned &j)
 Return local equation number of value j in field fld. More...
 
void output (std::ostream &outfile, const unsigned &nplot)
 Output FE representation of soln as in underlying element. More...
 
void residual_for_projection (Vector< double > &residuals, DenseMatrix< double > &jacobian, const unsigned &flag)
 
- Public Member Functions inherited from oomph::ProjectableElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >
 ProjectableElement ()
 
double zeta_nodal (const unsigned &n, const unsigned &k, const unsigned &i) const
 
void enable_projection ()
 
void disable_projection ()
 
void set_project_coordinates ()
 Project (history values of) coordintes. More...
 
void set_project_values ()
 Project (history values of) values. More...
 
void set_project_lagrangian ()
 Project (current and only values of) lagrangian coordinates. More...
 
unsigned & projected_field ()
 Field that is currently being projected. More...
 
unsigned & time_level_for_projection ()
 Which history value are we projecting? More...
 
unsigned & projected_coordinate ()
 
unsigned & projected_lagrangian_coordinate ()
 
- Public Member Functions inherited from oomph::ProjectableElementBase
 ProjectableElementBase ()
 
virtual ~ProjectableElementBase ()
 Virtual destructor. More...
 
- Public Member Functions inherited from oomph::ElementWithExternalElement
 ElementWithExternalElement ()
 
virtual ~ElementWithExternalElement ()
 The destructor, clean up any allocated memory. More...
 
 ElementWithExternalElement (const ElementWithExternalElement &)=delete
 Broken copy constructor. More...
 
void operator= (const ElementWithExternalElement &)=delete
 Broken assignment operator. More...
 
bool storage_has_been_allocated () const
 Helper function to check if storage has actually been allocated. More...
 
FiniteElement *& external_element_pt (const unsigned &interaction_index, const unsigned &ipt)
 
FiniteElement *const & external_element_pt (const unsigned &interaction_index, const unsigned &ipt) const
 Access function to source element, const version. More...
 
Vector< double > & external_element_local_coord (const unsigned &interaction_index, const unsigned &ipt)
 
Vector< double > const & external_element_local_coord (const unsigned &interaction_index, const unsigned &ipt) const
 Access function to get source element's coords, const version. More...
 
void output_external_elements (std::ostream &outfile, const unsigned &interaction_index)
 
void initialise_external_element_storage ()
 
void flush_all_external_element_storage ()
 Flush the storage for external elements. More...
 
void set_ninteraction (const unsigned &n_interaction)
 
unsigned ninteraction () const
 Return the number of interactions in the element. More...
 
virtual void identify_all_field_data_for_external_interaction (Vector< std::set< FiniteElement * >> const &external_elements_pt, std::set< std::pair< Data *, unsigned >> &paired_interaction_data)
 
virtual void identify_all_geometric_data_for_external_interaction (Vector< std::set< FiniteElement * >> const &external_elements_pt, std::set< Data * > &external_geometric_data_pt)
 
unsigned nexternal_interaction_field_data () const
 
Vector< Data * > external_interaction_field_data_pt () const
 
unsigned nexternal_interaction_geometric_data () const
 
bool add_external_geometric_data ()
 Are we including external geometric data in the element's Jacobian. More...
 
bool add_external_interaction_data ()
 Are we including external data in the element's Jacobian. More...
 
Vector< Data * > external_interaction_geometric_data_pt () const
 
void ignore_external_geometric_data ()
 
void ignore_external_interaction_data ()
 
void include_external_geometric_data ()
 
void include_external_interaction_data ()
 
bool external_geometric_data_is_included () const
 
- Public Member Functions inherited from oomph::FiniteElement
void set_dimension (const unsigned &dim)
 
void set_nodal_dimension (const unsigned &nodal_dim)
 
void set_nnodal_position_type (const unsigned &nposition_type)
 Set the number of types required to interpolate the coordinate. More...
 
void set_n_node (const unsigned &n)
 
int nodal_local_eqn (const unsigned &n, const unsigned &i) const
 
double dJ_eulerian_at_knot (const unsigned &ipt, Shape &psi, DenseMatrix< double > &djacobian_dX) const
 
 FiniteElement ()
 Constructor. More...
 
virtual ~FiniteElement ()
 
 FiniteElement (const FiniteElement &)=delete
 Broken copy constructor. More...
 
virtual bool local_coord_is_valid (const Vector< double > &s)
 Broken assignment operator. More...
 
virtual void move_local_coord_back_into_element (Vector< double > &s) const
 
void get_centre_of_gravity_and_max_radius_in_terms_of_zeta (Vector< double > &cog, double &max_radius) const
 
virtual void local_coordinate_of_node (const unsigned &j, Vector< double > &s) const
 
virtual void local_fraction_of_node (const unsigned &j, Vector< double > &s_fraction)
 
virtual double local_one_d_fraction_of_node (const unsigned &n1d, const unsigned &i)
 
virtual void set_macro_elem_pt (MacroElement *macro_elem_pt)
 
MacroElement * macro_elem_pt ()
 Access function to pointer to macro element. More...
 
void get_x (const Vector< double > &s, Vector< double > &x) const
 
void get_x (const unsigned &t, const Vector< double > &s, Vector< double > &x)
 
virtual void get_x_from_macro_element (const Vector< double > &s, Vector< double > &x) const
 
virtual void get_x_from_macro_element (const unsigned &t, const Vector< double > &s, Vector< double > &x)
 
virtual void set_integration_scheme (Integral *const &integral_pt)
 Set the spatial integration scheme. More...
 
Integral *const & integral_pt () const
 Return the pointer to the integration scheme (const version) More...
 
virtual void shape (const Vector< double > &s, Shape &psi) const =0
 
virtual void shape_at_knot (const unsigned &ipt, Shape &psi) const
 
virtual void dshape_local (const Vector< double > &s, Shape &psi, DShape &dpsids) const
 
virtual void dshape_local_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsids) const
 
virtual void d2shape_local (const Vector< double > &s, Shape &psi, DShape &dpsids, DShape &d2psids) const
 
virtual void d2shape_local_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsids, DShape &d2psids) const
 
virtual double J_eulerian (const Vector< double > &s) const
 
virtual double J_eulerian_at_knot (const unsigned &ipt) const
 
void check_J_eulerian_at_knots (bool &passed) const
 
void check_jacobian (const double &jacobian) const
 
double dshape_eulerian (const Vector< double > &s, Shape &psi, DShape &dpsidx) const
 
virtual double dshape_eulerian_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsidx) const
 
virtual double dshape_eulerian_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsi, DenseMatrix< double > &djacobian_dX, RankFourTensor< double > &d_dpsidx_dX) const
 
double d2shape_eulerian (const Vector< double > &s, Shape &psi, DShape &dpsidx, DShape &d2psidx) const
 
virtual double d2shape_eulerian_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsidx, DShape &d2psidx) const
 
virtual void assign_nodal_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void describe_nodal_local_dofs (std::ostream &out, const std::string &current_string) const
 
virtual void assign_all_generic_local_eqn_numbers (const bool &store_local_dof_pt)
 
Node *& node_pt (const unsigned &n)
 Return a pointer to the local node n. More...
 
Node *const & node_pt (const unsigned &n) const
 Return a pointer to the local node n (const version) More...
 
unsigned nnode () const
 Return the number of nodes. More...
 
virtual unsigned nnode_1d () const
 
double raw_nodal_position (const unsigned &n, const unsigned &i) const
 
double raw_nodal_position (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double raw_dnodal_position_dt (const unsigned &n, const unsigned &i) const
 
double raw_dnodal_position_dt (const unsigned &n, const unsigned &j, const unsigned &i) const
 
double raw_nodal_position_gen (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_nodal_position_gen (const unsigned &t, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_dnodal_position_gen_dt (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double raw_dnodal_position_gen_dt (const unsigned &j, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double nodal_position (const unsigned &n, const unsigned &i) const
 
double nodal_position (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double dnodal_position_dt (const unsigned &n, const unsigned &i) const
 Return the i-th component of nodal velocity: dx/dt at local node n. More...
 
double dnodal_position_dt (const unsigned &n, const unsigned &j, const unsigned &i) const
 
double nodal_position_gen (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double nodal_position_gen (const unsigned &t, const unsigned &n, const unsigned &k, const unsigned &i) const
 
double dnodal_position_gen_dt (const unsigned &n, const unsigned &k, const unsigned &i) const
 
double dnodal_position_gen_dt (const unsigned &j, const unsigned &n, const unsigned &k, const unsigned &i) const
 
virtual void get_dresidual_dnodal_coordinates (RankThreeTensor< double > &dresidual_dnodal_coordinates)
 
virtual void disable_ALE ()
 
virtual void enable_ALE ()
 
virtual unsigned required_nvalue (const unsigned &n) const
 
unsigned nnodal_position_type () const
 
bool has_hanging_nodes () const
 
unsigned nodal_dimension () const
 Return the required Eulerian dimension of the nodes in this element. More...
 
virtual unsigned nvertex_node () const
 
virtual Node * vertex_node_pt (const unsigned &j) const
 
virtual Node * construct_node (const unsigned &n)
 
virtual Node * construct_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 
virtual Node * construct_boundary_node (const unsigned &n)
 
virtual Node * construct_boundary_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 
int get_node_number (Node *const &node_pt) const
 
virtual Node * get_node_at_local_coordinate (const Vector< double > &s) const
 
double raw_nodal_value (const unsigned &n, const unsigned &i) const
 
double raw_nodal_value (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double nodal_value (const unsigned &n, const unsigned &i) const
 
double nodal_value (const unsigned &t, const unsigned &n, const unsigned &i) const
 
unsigned dim () const
 
virtual ElementGeometry::ElementGeometry element_geometry () const
 Return the geometry type of the element (either Q or T usually). More...
 
virtual double interpolated_x (const Vector< double > &s, const unsigned &i) const
 Return FE interpolated coordinate x[i] at local coordinate s. More...
 
virtual double interpolated_x (const unsigned &t, const Vector< double > &s, const unsigned &i) const
 
virtual void interpolated_x (const Vector< double > &s, Vector< double > &x) const
 Return FE interpolated position x[] at local coordinate s as Vector. More...
 
virtual void interpolated_x (const unsigned &t, const Vector< double > &s, Vector< double > &x) const
 
virtual double interpolated_dxdt (const Vector< double > &s, const unsigned &i, const unsigned &t)
 
virtual void interpolated_dxdt (const Vector< double > &s, const unsigned &t, Vector< double > &dxdt)
 
unsigned ngeom_data () const
 
Data * geom_data_pt (const unsigned &j)
 
void position (const Vector< double > &zeta, Vector< double > &r) const
 
void position (const unsigned &t, const Vector< double > &zeta, Vector< double > &r) const
 
void dposition_dt (const Vector< double > &zeta, const unsigned &t, Vector< double > &drdt)
 
void interpolated_zeta (const Vector< double > &s, Vector< double > &zeta) const
 
void locate_zeta (const Vector< double > &zeta, GeomObject *&geom_object_pt, Vector< double > &s, const bool &use_coordinate_as_initial_guess=false)
 
virtual void node_update ()
 
virtual void identify_field_data_for_interactions (std::set< std::pair< Data *, unsigned >> &paired_field_data)
 
virtual void identify_geometric_data (std::set< Data * > &geometric_data_pt)
 
virtual double s_min () const
 Min value of local coordinate. More...
 
virtual double s_max () const
 Max. value of local coordinate. More...
 
double size () const
 
virtual double compute_physical_size () const
 
virtual void point_output_data (const Vector< double > &s, Vector< double > &data)
 
void point_output (std::ostream &outfile, const Vector< double > &s)
 
virtual unsigned nplot_points_paraview (const unsigned &nplot) const
 
virtual unsigned nsub_elements_paraview (const unsigned &nplot) const
 
void output_paraview (std::ofstream &file_out, const unsigned &nplot) const
 
virtual void write_paraview_output_offset_information (std::ofstream &file_out, const unsigned &nplot, unsigned &counter) const
 
virtual void write_paraview_type (std::ofstream &file_out, const unsigned &nplot) const
 
virtual void write_paraview_offsets (std::ofstream &file_out, const unsigned &nplot, unsigned &offset_sum) const
 
virtual unsigned nscalar_paraview () const
 
virtual void scalar_value_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot) const
 
virtual void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt) const
 
virtual void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt) const
 
virtual std::string scalar_name_paraview (const unsigned &i) const
 
virtual void output (std::ostream &outfile)
 
virtual void output (const unsigned &t, std::ostream &outfile, const unsigned &n_plot) const
 
virtual void output (FILE *file_pt)
 
virtual void output (FILE *file_pt, const unsigned &n_plot)
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 Output an exact solution over the element. More...
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 Output a time-dependent exact solution over the element. More...
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, const SolutionFunctorBase &exact_soln) const
 Output a time-dependent exact solution over the element. More...
 
virtual void get_s_plot (const unsigned &i, const unsigned &nplot, Vector< double > &s, const bool &shifted_to_interior=false) const
 
virtual std::string tecplot_zone_string (const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (std::ostream &outfile, const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (FILE *file_pt, const unsigned &nplot) const
 
virtual unsigned nplot_points (const unsigned &nplot) const
 
virtual void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 Calculate the norm of the error and that of the exact solution. More...
 
virtual void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 Calculate the norm of the error and that of the exact solution. More...
 
virtual void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_abs_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error)
 
void integrate_fct (FiniteElement::SteadyExactSolutionFctPt integrand_fct_pt, Vector< double > &integral)
 Integrate Vector-valued function over element. More...
 
void integrate_fct (FiniteElement::UnsteadyExactSolutionFctPt integrand_fct_pt, const double &time, Vector< double > &integral)
 Integrate Vector-valued time-dep function over element. More...
 
virtual void build_face_element (const int &face_index, FaceElement *face_element_pt)
 
virtual unsigned self_test ()
 
virtual unsigned get_bulk_node_number (const int &face_index, const unsigned &i) const
 
virtual int face_outer_unit_normal_sign (const int &face_index) const
 Get the sign of the outer unit normal on the face given by face_index. More...
 
virtual unsigned nnode_on_face () const
 
void face_node_number_error_check (const unsigned &i) const
 Range check for face node numbers. More...
 
virtual CoordinateMappingFctPt face_to_bulk_coordinate_fct_pt (const int &face_index) const
 
virtual BulkCoordinateDerivativesFctPt bulk_coordinate_derivatives_fct_pt (const int &face_index) const
 
- Public Member Functions inherited from oomph::GeneralisedElement
 GeneralisedElement ()
 Constructor: Initialise all pointers and all values to zero. More...
 
virtual ~GeneralisedElement ()
 Virtual destructor to clean up any memory allocated by the object. More...
 
 GeneralisedElement (const GeneralisedElement &)=delete
 Broken copy constructor. More...
 
void operator= (const GeneralisedElement &)=delete
 Broken assignment operator. More...
 
Data *& internal_data_pt (const unsigned &i)
 Return a pointer to i-th internal data object. More...
 
Data *const & internal_data_pt (const unsigned &i) const
 Return a pointer to i-th internal data object (const version) More...
 
Data *& external_data_pt (const unsigned &i)
 Return a pointer to i-th external data object. More...
 
Data *const & external_data_pt (const unsigned &i) const
 Return a pointer to i-th external data object (const version) More...
 
unsigned long eqn_number (const unsigned &ieqn_local) const
 
int local_eqn_number (const unsigned long &ieqn_global) const
 
unsigned add_external_data (Data *const &data_pt, const bool &fd=true)
 
bool external_data_fd (const unsigned &i) const
 
void exclude_external_data_fd (const unsigned &i)
 
void include_external_data_fd (const unsigned &i)
 
void flush_external_data ()
 Flush all external data. More...
 
void flush_external_data (Data *const &data_pt)
 Flush the object addressed by data_pt from the external data array. More...
 
unsigned ninternal_data () const
 Return the number of internal data objects. More...
 
unsigned nexternal_data () const
 Return the number of external data objects. More...
 
unsigned ndof () const
 Return the number of equations/dofs in the element. More...
 
void dof_vector (const unsigned &t, Vector< double > &dof)
 Return the vector of dof values at time level t. More...
 
void dof_pt_vector (Vector< double * > &dof_pt)
 Return the vector of pointers to dof values. More...
 
void set_internal_data_time_stepper (const unsigned &i, TimeStepper *const &time_stepper_pt, const bool &preserve_existing_data)
 
void assign_internal_eqn_numbers (unsigned long &global_number, Vector< double * > &Dof_pt)
 
void describe_dofs (std::ostream &out, const std::string &current_string) const
 
void add_internal_value_pt_to_map (std::map< unsigned, double * > &map_of_value_pt)
 
virtual void assign_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void complete_setup_of_dependencies ()
 
virtual void get_residuals (Vector< double > &residuals)
 
virtual void get_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
virtual void get_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &mass_matrix)
 
virtual void get_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
virtual void get_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
virtual void get_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
virtual void get_djacobian_and_dmass_matrix_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam)
 
virtual void get_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
virtual void get_inner_products (Vector< std::pair< unsigned, unsigned >> const &history_index, Vector< double > &inner_product)
 
virtual void get_inner_product_vectors (Vector< unsigned > const &history_index, Vector< Vector< double >> &inner_product_vector)
 
virtual void compute_norm (Vector< double > &norm)
 
virtual void compute_norm (double &norm)
 
virtual unsigned ndof_types () const
 
virtual void get_dof_numbers_for_unknowns (std::list< std::pair< unsigned long, unsigned >> &dof_lookup_list) const
 
- Public Member Functions inherited from oomph::GeomObject
 GeomObject ()
 Default constructor. More...
 
 GeomObject (const unsigned &ndim)
 
 GeomObject (const unsigned &nlagrangian, const unsigned &ndim)
 
 GeomObject (const unsigned &nlagrangian, const unsigned &ndim, TimeStepper *time_stepper_pt)
 
 GeomObject (const GeomObject &dummy)=delete
 Broken copy constructor. More...
 
void operator= (const GeomObject &)=delete
 Broken assignment operator. More...
 
virtual ~GeomObject ()
 (Empty) destructor More...
 
unsigned nlagrangian () const
 Access function to # of Lagrangian coordinates. More...
 
unsigned ndim () const
 Access function to # of Eulerian coordinates. More...
 
void set_nlagrangian_and_ndim (const unsigned &n_lagrangian, const unsigned &n_dim)
 Set # of Lagrangian and Eulerian coordinates. More...
 
TimeStepper *& time_stepper_pt ()
 
TimeStepper * time_stepper_pt () const
 
virtual void position (const double &t, const Vector< double > &zeta, Vector< double > &r) const
 
virtual void dposition (const Vector< double > &zeta, DenseMatrix< double > &drdzeta) const
 
virtual void d2position (const Vector< double > &zeta, RankThreeTensor< double > &ddrdzeta) const
 
virtual void d2position (const Vector< double > &zeta, Vector< double > &r, DenseMatrix< double > &drdzeta, RankThreeTensor< double > &ddrdzeta) const
 

Additional Inherited Members

- Public Types inherited from oomph::FiniteElement
typedef void(* SteadyExactSolutionFctPt) (const Vector< double > &, Vector< double > &)
 
typedef void(* UnsteadyExactSolutionFctPt) (const double &, const Vector< double > &, Vector< double > &)
 
- Static Public Attributes inherited from oomph::FiniteElement
static double Tolerance_for_singular_jacobian = 1.0e-16
 Tolerance below which the jacobian is considered singular. More...
 
static bool Accept_negative_jacobian = false
 
static bool Suppress_output_while_checking_for_inverted_elements
 
- Static Public Attributes inherited from oomph::GeneralisedElement
static bool Suppress_warning_about_repeated_internal_data
 
static bool Suppress_warning_about_repeated_external_data = true
 
static double Default_fd_jacobian_step = 1.0e-8
 
- Protected Types inherited from oomph::ProjectableElementBase
enum  Projection_Type { Coordinate , Lagrangian , Value }
 
- Protected Member Functions inherited from oomph::ProjectableElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Overloaded version of fill_in_contribution_to_residuals. More...
 
void describe_local_dofs (std::ostream &out, const std::string &current_string) const
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 Overloaded version of fill_in_contribution_to_jacobian. More...
 
- Protected Member Functions inherited from oomph::ElementWithExternalElement
void assign_internal_and_external_local_eqn_numbers (const bool &store_local_dof_pt)
 
void assign_external_interaction_data_local_eqn_numbers (const bool &store_local_dof_pt)
 
void fill_in_jacobian_from_external_interaction_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_external_interaction_by_fd (DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_external_interaction_field_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_external_interaction_field_by_fd (DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_external_interaction_geometric_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_external_interaction_geometric_by_fd (DenseMatrix< double > &jacobian)
 
virtual void update_before_external_interaction_field_fd ()
 
virtual void reset_after_external_interaction_field_fd ()
 
virtual void update_in_external_interaction_field_fd (const unsigned &i)
 
virtual void reset_in_external_interaction_field_fd (const unsigned &i)
 
virtual void update_before_external_interaction_geometric_fd ()
 
virtual void reset_after_external_interaction_geometric_fd ()
 
virtual void update_in_external_interaction_geometric_fd (const unsigned &i)
 
virtual void reset_in_external_interaction_geometric_fd (const unsigned &i)
 
- Protected Member Functions inherited from oomph::FiniteElement
virtual void assemble_local_to_eulerian_jacobian (const DShape &dpsids, DenseMatrix< double > &jacobian) const
 
virtual void assemble_local_to_eulerian_jacobian2 (const DShape &d2psids, DenseMatrix< double > &jacobian2) const
 
virtual void assemble_eulerian_base_vectors (const DShape &dpsids, DenseMatrix< double > &interpolated_G) const
 
template<unsigned DIM>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
virtual double invert_jacobian_mapping (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
virtual double local_to_eulerian_mapping (const DShape &dpsids, DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
double local_to_eulerian_mapping (const DShape &dpsids, DenseMatrix< double > &inverse_jacobian) const
 
virtual double local_to_eulerian_mapping_diagonal (const DShape &dpsids, DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
virtual void dJ_eulerian_dnodal_coordinates (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<unsigned DIM>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
virtual void d_dshape_eulerian_dnodal_coordinates (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<unsigned DIM>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
virtual void transform_derivatives (const DenseMatrix< double > &inverse_jacobian, DShape &dbasis) const
 
void transform_derivatives_diagonal (const DenseMatrix< double > &inverse_jacobian, DShape &dbasis) const
 
virtual void transform_second_derivatives (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<unsigned DIM>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<unsigned DIM>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
virtual void fill_in_jacobian_from_nodal_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_jacobian_from_nodal_by_fd (DenseMatrix< double > &jacobian)
 
virtual void update_before_nodal_fd ()
 
virtual void reset_after_nodal_fd ()
 
virtual void update_in_nodal_fd (const unsigned &i)
 
virtual void reset_in_nodal_fd (const unsigned &i)
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 Zero-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 One-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 Two-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
- Protected Member Functions inherited from oomph::GeneralisedElement
unsigned add_internal_data (Data *const &data_pt, const bool &fd=true)
 
bool internal_data_fd (const unsigned &i) const
 
void exclude_internal_data_fd (const unsigned &i)
 
void include_internal_data_fd (const unsigned &i)
 
void clear_global_eqn_numbers ()
 
void add_global_eqn_numbers (std::deque< unsigned long > const &global_eqn_numbers, std::deque< double * > const &global_dof_pt)
 
virtual void assign_additional_local_eqn_numbers ()
 
int internal_local_eqn (const unsigned &i, const unsigned &j) const
 
int external_local_eqn (const unsigned &i, const unsigned &j)
 
void fill_in_jacobian_from_internal_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_internal_by_fd (DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_external_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_external_by_fd (DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
virtual void update_before_internal_fd ()
 
virtual void reset_after_internal_fd ()
 
virtual void update_in_internal_fd (const unsigned &i)
 
virtual void reset_in_internal_fd (const unsigned &i)
 
virtual void update_before_external_fd ()
 
virtual void reset_after_external_fd ()
 
virtual void update_in_external_fd (const unsigned &i)
 
virtual void reset_in_external_fd (const unsigned &i)
 
virtual void fill_in_contribution_to_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &mass_matrix)
 
virtual void fill_in_contribution_to_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
virtual void fill_in_contribution_to_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
virtual void fill_in_contribution_to_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
virtual void fill_in_contribution_to_djacobian_and_dmass_matrix_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam)
 
virtual void fill_in_contribution_to_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
virtual void fill_in_contribution_to_inner_products (Vector< std::pair< unsigned, unsigned >> const &history_index, Vector< double > &inner_product)
 
virtual void fill_in_contribution_to_inner_product_vectors (Vector< unsigned > const &history_index, Vector< Vector< double >> &inner_product_vector)
 
- Protected Attributes inherited from oomph::ProjectableElementBase
unsigned Projected_field
 Field that is currently being projected. More...
 
unsigned Time_level_for_projection
 Time level we are projecting (0=current values; >0: history values) More...
 
unsigned Projected_coordinate
 
unsigned Projected_lagrangian
 
Projection_Type Projection_type
 
bool Do_projection
 
unsigned Backup_ninteraction
 
bool Backup_external_geometric_data
 
bool Backup_external_interaction_data
 
- Protected Attributes inherited from oomph::ElementWithExternalElement
bool Add_external_interaction_data
 Boolean flag to indicate whether to include the external data. More...
 
bool Add_external_geometric_data
 Boolean flag to indicate whether to include the external geometric data. More...
 
Data ** External_interaction_field_data_pt
 
Data ** External_interaction_geometric_data_pt
 
- Protected Attributes inherited from oomph::FiniteElement
MacroElement * Macro_elem_pt
 Pointer to the element's macro element (NULL by default) More...
 
- Protected Attributes inherited from oomph::GeomObject
unsigned NLagrangian
 Number of Lagrangian (intrinsic) coordinates. More...
 
unsigned Ndim
 Number of Eulerian coordinates. More...
 
TimeStepper * Geom_object_time_stepper_pt
 
- Static Protected Attributes inherited from oomph::FiniteElement
static const unsigned Default_Initial_Nvalue = 0
 Default value for the number of values at a node. More...
 
static const double Node_location_tolerance = 1.0e-14
 
static const unsigned N2deriv [] = {0, 1, 3, 6}
 
- Static Protected Attributes inherited from oomph::GeneralisedElement
static DenseMatrix< double > Dummy_matrix
 
static std::deque< double * > Dof_pt_deque
 

Detailed Description

template<class AXISYMMETRIC_POROELASTICITY_ELEMENT>
class oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >

Axisymmetric poro elasticity upgraded to become projectable.

///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////////////

Constructor & Destructor Documentation

◆ ProjectableAxisymmetricPoroelasticityElement()

template<class AXISYMMETRIC_POROELASTICITY_ELEMENT >
oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::ProjectableAxisymmetricPoroelasticityElement ( )
inline

Constructor [this was only required explicitly from gcc 4.5.2 onwards...]

1403 {}

Member Function Documentation

◆ data_values_of_field()

template<class AXISYMMETRIC_POROELASTICITY_ELEMENT >
Vector<std::pair<Data*, unsigned> > oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::data_values_of_field ( const unsigned &  fld)
inlinevirtual

Specify the values associated with field fld. The information is returned in a vector of pairs which comprise the Data object and the value within it, that correspond to field fld.

Implements oomph::ProjectableElementBase.

1409  {
1410 #ifdef PARANOID
1411  if (fld > 3)
1412  {
1413  std::stringstream error_stream;
1414  error_stream
1415  << "AxisymmetricPoroelasticity elements only store 4 fields so fld = "
1416  << fld << " is illegal \n";
1417  throw OomphLibError(
1418  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1419  }
1420 #endif
1421 
1422  // Create the vector
1423  Vector<std::pair<Data*, unsigned>> data_values;
1424 
1425 
1426  // Pressure
1427  //---------
1428  if (fld == 2)
1429  {
1430  Data* dat_pt = this->p_data_pt();
1431  unsigned nvalue = dat_pt->nvalue();
1432  for (unsigned i = 0; i < nvalue; i++)
1433  {
1434  data_values.push_back(std::make_pair(dat_pt, i));
1435  }
1436  }
1437  // Darcy flux
1438  //-----------
1439  else if (fld == 3)
1440  {
1441  Vector<Data*> edge_dat_pt = this->q_edge_data_pt();
1442  unsigned n = edge_dat_pt.size();
1443  for (unsigned j = 0; j < n; j++)
1444  {
1445  unsigned nvalue = this->nedge_flux_interpolation_point();
1446  for (unsigned i = 0; i < nvalue; i++)
1447  {
1448  data_values.push_back(
1449  std::make_pair(edge_dat_pt[j], this->q_edge_index(i)));
1450  }
1451  }
1452  if (this->nq_basis_internal() > 0)
1453  {
1454  Data* int_dat_pt = this->q_internal_data_pt();
1455  unsigned nvalue = int_dat_pt->nvalue();
1456  for (unsigned i = 0; i < nvalue; i++)
1457  {
1458  data_values.push_back(std::make_pair(int_dat_pt, i));
1459  }
1460  }
1461  }
1462  // Solid displacements
1463  else
1464  {
1465  // Loop over all nodes
1466  unsigned nnod = this->nnode();
1467  for (unsigned j = 0; j < nnod; j++)
1468  {
1469  // Add the data value associated with the displacement components
1470  // (stored first)
1471  data_values.push_back(std::make_pair(
1472  this->node_pt(j), this->u_index_axisym_poroelasticity(fld)));
1473  }
1474  }
1475 
1476  // Return the vector
1477  return data_values;
1478  }
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9
n
const unsigned n
Definition: CG3DPackingUnitTest.cpp:11
oomph::FiniteElement::node_pt
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175
oomph::FiniteElement::nnode
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210
OOMPH_EXCEPTION_LOCATION
#define OOMPH_EXCEPTION_LOCATION
Definition: oomph_definitions.h:61
OOMPH_CURRENT_FUNCTION
#define OOMPH_CURRENT_FUNCTION
Definition: oomph_definitions.h:86
j
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References i, j, n, oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), oomph::Data::nvalue(), OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

◆ get_field()

template<class AXISYMMETRIC_POROELASTICITY_ELEMENT >
double oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::get_field ( const unsigned &  t,
const unsigned &  fld,
const Vector< double > &  s 
)
inlinevirtual

Return interpolated field fld at local coordinate s, at time level t (t=0: present; t>0: history values)

Implements oomph::ProjectableElementBase.

1613  {
1614 #ifdef PARANOID
1615  if (fld > 3)
1616  {
1617  std::stringstream error_stream;
1618  error_stream
1619  << "AxisymmetricPoroelasticity elements only store 4 fields so fld = "
1620  << fld << " is illegal\n";
1621  throw OomphLibError(
1622  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1623  }
1624 #endif
1625 
1626  double return_value = 0.0;
1627 
1628  // Pressure
1629  if (fld == 2)
1630  {
1631  // No time-dependence in here
1632 #ifdef PARANOID
1633  if (t != 0)
1634  {
1635  throw OomphLibError(
1636  "Pressure in AxisymmetricPoroelasticity has no time-dep.!",
1637  OOMPH_CURRENT_FUNCTION,
1638  OOMPH_EXCEPTION_LOCATION);
1639  }
1640 #endif
1641  this->interpolated_p(s, return_value);
1642  }
1643  // Darcy flux -- doesn't really work as it's a vector field
1644  else if (fld == 3)
1645  {
1646  throw OomphLibError(
1647  "Don't call this function for AxisymmetricPoroelasticity!",
1648  OOMPH_CURRENT_FUNCTION,
1649  OOMPH_EXCEPTION_LOCATION);
1650  }
1651  // Displacement components
1652  else
1653  {
1654  return_value = this->interpolated_u(t, s, fld);
1655  }
1656  return return_value;
1657  }
s
RealScalar s
Definition: level1_cplx_impl.h:130
plotPSD.t
t
Definition: plotPSD.py:36

References OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, s, and plotPSD::t.

Referenced by oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::residual_for_projection().

◆ jacobian_and_shape_of_field()

template<class AXISYMMETRIC_POROELASTICITY_ELEMENT >
double oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::jacobian_and_shape_of_field ( const unsigned &  fld,
const Vector< double > &  s,
Shape &  psi 
)
inlinevirtual

Return Jacobian of mapping and shape functions of field fld at local coordinate s

Implements oomph::ProjectableElementBase.

1529  {
1530 #ifdef PARANOID
1531  if (fld > 3)
1532  {
1533  std::stringstream error_stream;
1534  error_stream
1535  << "AxisymmetricPoroelasticity elements only store 4 fields so fld = "
1536  << fld << " is illegal.\n";
1537  throw OomphLibError(
1538  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1539  }
1540 #endif
1541 
1542 
1543  // Get the number of geometric nodes, total number of basis functions,
1544  // and number of edges basis functions
1545  const unsigned n_dim = this->dim();
1546  const unsigned n_node = this->nnode();
1547  const unsigned n_q_basis = this->nq_basis();
1548  const unsigned n_p_basis = this->np_basis();
1549 
1550  // Storage for the geometric and computational bases and the test
1551  // functions
1552  Shape psi_geom(n_node), q_basis(n_q_basis, n_dim),
1553  q_test(n_q_basis, n_dim), p_basis(n_p_basis), p_test(n_p_basis),
1554  div_q_basis_ds(n_q_basis), div_q_test_ds(n_q_basis);
1555  DShape dpsidx_geom(n_node, n_dim);
1556  Shape u_basis(n_node);
1557  Shape u_test(n_node);
1558  DShape du_basis_dx(n_node, n_dim);
1559  DShape du_test_dx(n_node, n_dim);
1560  double J = this->shape_basis_test_local(s,
1561  psi_geom,
1562  dpsidx_geom,
1563  u_basis,
1564  u_test,
1565  du_basis_dx,
1566  du_test_dx,
1567  q_basis,
1568  q_test,
1569  p_basis,
1570  p_test,
1571  div_q_basis_ds,
1572  div_q_test_ds);
1573  // Pressure basis functions
1574  if (fld == 2)
1575  {
1576  unsigned n = p_basis.nindex1();
1577  for (unsigned i = 0; i < n; i++)
1578  {
1579  psi[i] = p_basis[i];
1580  }
1581  }
1582  // Flux basis functions
1583  else if (fld == 3)
1584  {
1585  unsigned n = q_basis.nindex1();
1586  unsigned m = q_basis.nindex2();
1587  for (unsigned i = 0; i < n; i++)
1588  {
1589  for (unsigned j = 0; j < m; j++)
1590  {
1591  psi(i, j) = q_basis(i, j);
1592  }
1593  }
1594  }
1595  // Displacement components
1596  else
1597  {
1598  for (unsigned j = 0; j < n_node; j++)
1599  {
1600  psi[j] = u_basis[j];
1601  }
1602  }
1603 
1604  return J;
1605  }
J
JacobiRotation< float > J
Definition: Jacobi_makeJacobi.cpp:3
oomph::FiniteElement::dim
unsigned dim() const
Definition: elements.h:2611
m
int * m
Definition: level2_cplx_impl.h:294

References oomph::FiniteElement::dim(), i, J, j, m, n, oomph::FiniteElement::nnode(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, and s.

Referenced by oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::residual_for_projection().

◆ local_equation()

template<class AXISYMMETRIC_POROELASTICITY_ELEMENT >
int oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::local_equation ( const unsigned &  fld,
const unsigned &  j 
)
inlinevirtual

Return local equation number of value j in field fld.

Implements oomph::ProjectableElementBase.

1698  {
1699 #ifdef PARANOID
1700  if (fld > 3)
1701  {
1702  std::stringstream error_stream;
1703  error_stream
1704  << "AxisymmetricPoroelasticity elements only store 4 fields so fld = "
1705  << fld << " is illegal\n";
1706  throw OomphLibError(
1707  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1708  }
1709 #endif
1710 
1711  int return_value = 0;
1712 
1713  // Pressure
1714  if (fld == 2)
1715  {
1716  return_value = this->p_local_eqn(j);
1717  }
1718  // Darcy flux
1719  else if (fld == 3)
1720  {
1721  unsigned nedge = this->nq_basis_edge();
1722  if (j < nedge)
1723  {
1724  return_value = this->q_edge_local_eqn(j);
1725  }
1726  else
1727  {
1728  return_value = this->q_internal_local_eqn(j - nedge);
1729  }
1730  }
1731  // Displacement
1732  else
1733  {
1734  return_value =
1735  this->nodal_local_eqn(j, this->u_index_axisym_poroelasticity(fld));
1736  }
1737 
1738  return return_value;
1739  }
oomph::FiniteElement::nodal_local_eqn
int nodal_local_eqn(const unsigned &n, const unsigned &i) const
Definition: elements.h:1432

References j, oomph::FiniteElement::nodal_local_eqn(), OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

Referenced by oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::residual_for_projection().

◆ nfields_for_projection()

template<class AXISYMMETRIC_POROELASTICITY_ELEMENT >
unsigned oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::nfields_for_projection ( )
inlinevirtual

Number of fields to be projected: 4 (two displacement components, pressure, Darcy flux)

Implements oomph::ProjectableElementBase.

1483  {
1484  return 4;
1485  }

◆ nhistory_values_for_coordinate_projection()

template<class AXISYMMETRIC_POROELASTICITY_ELEMENT >
unsigned oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::nhistory_values_for_coordinate_projection ( )
inlinevirtual

Number of positional history values (Note: count includes current value!)

Implements oomph::ProjectableElementBase.

1520  {
1521  return this->node_pt(0)->position_time_stepper_pt()->ntstorage();
1522  }
oomph::Node::position_time_stepper_pt
TimeStepper *& position_time_stepper_pt()
Return a pointer to the position timestepper.
Definition: nodes.h:1022
oomph::TimeStepper::ntstorage
unsigned ntstorage() const
Definition: timesteppers.h:601

References oomph::FiniteElement::node_pt(), oomph::TimeStepper::ntstorage(), and oomph::Node::position_time_stepper_pt().

◆ nhistory_values_for_projection()

template<class AXISYMMETRIC_POROELASTICITY_ELEMENT >
unsigned oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::nhistory_values_for_projection ( const unsigned &  fld)
inlinevirtual

Number of history values to be stored for fld-th field. (Note: count includes current value!)

Implements oomph::ProjectableElementBase.

1490  {
1491 #ifdef PARANOID
1492  if (fld > 3)
1493  {
1494  std::stringstream error_stream;
1495  error_stream
1496  << "AxisymmetricPoroelasticity elements only store 4 fields so fld = "
1497  << fld << " is illegal\n";
1498  throw OomphLibError(
1499  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1500  }
1501 #endif
1502 
1503  // Displacements -- infer from first (vertex) node
1504  unsigned return_value = this->node_pt(0)->ntstorage();
1505 
1506  // Pressure: No history values (just present one!)
1507  if (fld == 2) return_value = 1;
1508 
1509  // Flux: infer from first midside node
1510  if (fld == 3)
1511  {
1512  return_value = this->node_pt(3)->ntstorage();
1513  }
1514  return return_value;
1515  }
oomph::Data::ntstorage
unsigned ntstorage() const
Definition: nodes.cc:879

References oomph::FiniteElement::node_pt(), oomph::Data::ntstorage(), OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

◆ nvalue_of_field()

template<class AXISYMMETRIC_POROELASTICITY_ELEMENT >
unsigned oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::nvalue_of_field ( const unsigned &  fld)
inlinevirtual

Return number of values in field fld.

Implements oomph::ProjectableElementBase.

1662  {
1663 #ifdef PARANOID
1664  if (fld > 3)
1665  {
1666  std::stringstream error_stream;
1667  error_stream
1668  << "AxisymmetricPoroelasticity elements only store 4 fields so fld = "
1669  << fld << " is illegal\n";
1670  throw OomphLibError(
1671  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1672  }
1673 #endif
1674 
1675  unsigned return_value = 0;
1676  // Pressure
1677  if (fld == 2)
1678  {
1679  return_value = this->np_basis();
1680  }
1681  // Darcy flux
1682  else if (fld == 3)
1683  {
1684  return_value = this->nq_basis();
1685  }
1686  // Displacements
1687  else
1688  {
1689  return_value = this->nnode();
1690  }
1691 
1692  return return_value;
1693  }

References oomph::FiniteElement::nnode(), OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

Referenced by oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::residual_for_projection().

◆ output()

template<class AXISYMMETRIC_POROELASTICITY_ELEMENT >
void oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::output ( std::ostream &  outfile,
const unsigned &  nplot 
)
inlinevirtual

Output FE representation of soln as in underlying element.

Reimplemented from oomph::FiniteElement.

1744  {
1745  AXISYMMETRIC_POROELASTICITY_ELEMENT::output(outfile, nplot);
1746  }
output
void output(std::ostream &outfile, const unsigned &nplot)
Overload output function.
Definition: overloaded_element_body.h:490

References output().

◆ residual_for_projection()

template<class AXISYMMETRIC_POROELASTICITY_ELEMENT >
void oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::residual_for_projection ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian,
const unsigned &  flag 
)
inlinevirtual

Residual for the projection step. Flag indicates if we want the Jacobian (1) or not (0). Virtual so it can be overloaded if necessary

Reimplemented from oomph::ProjectableElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >.

1755  {
1756  // Am I a solid element
1757  SolidFiniteElement* solid_el_pt = dynamic_cast<SolidFiniteElement*>(this);
1758 
1759  unsigned n_dim = this->dim();
1760 
1761  // Allocate storage for local coordinates
1762  Vector<double> s(n_dim);
1763 
1764  // Current field
1765  unsigned fld = this->Projected_field;
1766 
1767  // Number of nodes
1768  const unsigned n_node = this->nnode();
1769 
1770  // Number of positional dofs
1771  const unsigned n_position_type = this->nnodal_position_type();
1772 
1773  // Number of dof for current field
1774  const unsigned n_value = nvalue_of_field(fld);
1775 
1776  // Set the value of n_intpt
1777  const unsigned n_intpt = this->integral_pt()->nweight();
1778 
1779  // Loop over the integration points
1780  for (unsigned ipt = 0; ipt < n_intpt; ipt++)
1781  {
1782  // Get the local coordinates of Gauss point
1783  for (unsigned i = 0; i < n_dim; i++)
1784  s[i] = this->integral_pt()->knot(ipt, i);
1785 
1786  // Get the integral weight
1787  double w = this->integral_pt()->weight(ipt);
1788 
1789  // Find same point in base mesh using external storage
1790  FiniteElement* other_el_pt = 0;
1791  other_el_pt = this->external_element_pt(0, ipt);
1792  Vector<double> other_s(n_dim);
1793  other_s = this->external_element_local_coord(0, ipt);
1794  ProjectableElement<AXISYMMETRIC_POROELASTICITY_ELEMENT>* cast_el_pt =
1795  dynamic_cast<
1796  ProjectableElement<AXISYMMETRIC_POROELASTICITY_ELEMENT>*>(
1797  other_el_pt);
1798 
1799  // Storage for the local equation and local unknown
1800  int local_eqn = 0, local_unknown = 0;
1801 
1802  // Now set up the three different projection problems
1803  switch (Projection_type)
1804  {
1805  case Lagrangian:
1806  {
1807  // If we don't have a solid element there is a problem
1808  if (solid_el_pt == 0)
1809  {
1810  throw OomphLibError("Trying to project Lagrangian coordinates in "
1811  "non-SolidElement\n",
1812  OOMPH_CURRENT_FUNCTION,
1813  OOMPH_EXCEPTION_LOCATION);
1814  }
1815 
1816  // Find the position shape functions
1817  Shape psi(n_node, n_position_type);
1818  // Get the position shape functions
1819  this->shape(s, psi);
1820  // Get the jacobian
1821  double J = this->J_eulerian(s);
1822 
1823  // Premultiply the weights and the Jacobian
1824  double W = w * J;
1825 
1826  // Get the value of the current position of the 0th coordinate
1827  // in the current element
1828  // at the current time level (which is the unkown)
1829  double interpolated_xi_proj = this->interpolated_x(s, 0);
1830 
1831  // Get the Lagrangian position in the other element
1832  double interpolated_xi_bar =
1833  dynamic_cast<SolidFiniteElement*>(cast_el_pt)
1834  ->interpolated_xi(other_s, Projected_lagrangian);
1835 
1836  // Now loop over the nodes and position dofs
1837  for (unsigned l = 0; l < n_node; ++l)
1838  {
1839  // Loop over position unknowns
1840  for (unsigned k = 0; k < n_position_type; ++k)
1841  {
1842  // The local equation is going to be the positional local
1843  // equation
1844  local_eqn = solid_el_pt->position_local_eqn(l, k, 0);
1845 
1846  // Now assemble residuals
1847  if (local_eqn >= 0)
1848  {
1849  // calculate residuals
1850  residuals[local_eqn] +=
1851  (interpolated_xi_proj - interpolated_xi_bar) * psi(l, k) *
1852  W;
1853 
1854  // Calculate the jacobian
1855  if (flag == 1)
1856  {
1857  for (unsigned l2 = 0; l2 < n_node; l2++)
1858  {
1859  // Loop over position dofs
1860  for (unsigned k2 = 0; k2 < n_position_type; k2++)
1861  {
1862  local_unknown =
1863  solid_el_pt->position_local_eqn(l2, k2, 0);
1864  if (local_unknown >= 0)
1865  {
1866  jacobian(local_eqn, local_unknown) +=
1867  psi(l2, k2) * psi(l, k) * W;
1868  }
1869  }
1870  }
1871  } // end of jacobian
1872  }
1873  }
1874  }
1875  } // End of Lagrangian coordinate case
1876 
1877  break;
1878 
1879  // Now the coordinate history case
1880  case Coordinate:
1881  {
1882  // Find the position shape functions
1883  Shape psi(n_node, n_position_type);
1884  // Get the position shape functions
1885  this->shape(s, psi);
1886  // Get the jacobian
1887  double J = this->J_eulerian(s);
1888 
1889  // Premultiply the weights and the Jacobian
1890  double W = w * J;
1891 
1892  // Get the value of the current position in the current element
1893  // at the current time level (which is the unkown)
1894  double interpolated_x_proj = 0.0;
1895  // If we are a solid element read it out directly from the data
1896  if (solid_el_pt != 0)
1897  {
1898  interpolated_x_proj =
1899  this->interpolated_x(s, Projected_coordinate);
1900  }
1901  // Otherwise it's the field value at the current time
1902  else
1903  {
1904  interpolated_x_proj = this->get_field(0, fld, s);
1905  }
1906 
1907  // Get the position in the other element
1908  double interpolated_x_bar = cast_el_pt->interpolated_x(
1909  Time_level_for_projection, other_s, Projected_coordinate);
1910 
1911  // Now loop over the nodes and position dofs
1912  for (unsigned l = 0; l < n_node; ++l)
1913  {
1914  // Loop over position unknowns
1915  for (unsigned k = 0; k < n_position_type; ++k)
1916  {
1917  // If I'm a solid element
1918  if (solid_el_pt != 0)
1919  {
1920  // The local equation is going to be the positional local
1921  // equation
1922  local_eqn =
1923  solid_el_pt->position_local_eqn(l, k, Projected_coordinate);
1924  }
1925  // Otherwise just pick the local unknown of a field to
1926  // project into
1927  else
1928  {
1929  // Complain if using generalised position types
1930  // but this is not a solid element, because the storage
1931  // is then not clear
1932  if (n_position_type != 1)
1933  {
1934  throw OomphLibError("Trying to project generalised "
1935  "positions not in SolidElement\n",
1936  OOMPH_CURRENT_FUNCTION,
1937  OOMPH_EXCEPTION_LOCATION);
1938  }
1939  local_eqn = local_equation(fld, l);
1940  }
1941 
1942  // Now assemble residuals
1943  if (local_eqn >= 0)
1944  {
1945  // calculate residuals
1946  residuals[local_eqn] +=
1947  (interpolated_x_proj - interpolated_x_bar) * psi(l, k) * W;
1948 
1949  // Calculate the jacobian
1950  if (flag == 1)
1951  {
1952  for (unsigned l2 = 0; l2 < n_node; l2++)
1953  {
1954  // Loop over position dofs
1955  for (unsigned k2 = 0; k2 < n_position_type; k2++)
1956  {
1957  // If I'm a solid element
1958  if (solid_el_pt != 0)
1959  {
1960  local_unknown = solid_el_pt->position_local_eqn(
1961  l2, k2, Projected_coordinate);
1962  }
1963  else
1964  {
1965  local_unknown = local_equation(fld, l2);
1966  }
1967 
1968  if (local_unknown >= 0)
1969  {
1970  jacobian(local_eqn, local_unknown) +=
1971  psi(l2, k2) * psi(l, k) * W;
1972  }
1973  }
1974  }
1975  } // end of jacobian
1976  }
1977  }
1978  }
1979  } // End of coordinate case
1980  break;
1981 
1982  // Now the value case
1983  case Value:
1984  {
1985  // Pressure or displacements -- "normal" procedure
1986  if (fld <= 2)
1987  {
1988  // Field shape function
1989  Shape psi(n_value);
1990 
1991  // Calculate jacobian and shape functions for that field
1992  double J = jacobian_and_shape_of_field(fld, s, psi);
1993 
1994  // Premultiply the weights and the Jacobian
1995  double W = w * J;
1996 
1997  // Value of field in current element at current time level
1998  //(the unknown)
1999  unsigned now = 0;
2000  double interpolated_value_proj = this->get_field(now, fld, s);
2001 
2002  // Value of the interpolation of element located in base mesh
2003  double interpolated_value_bar =
2004  cast_el_pt->get_field(Time_level_for_projection, fld, other_s);
2005 
2006  // Loop over dofs of field fld
2007  for (unsigned l = 0; l < n_value; l++)
2008  {
2009  local_eqn = local_equation(fld, l);
2010  if (local_eqn >= 0)
2011  {
2012  // calculate residuals
2013  residuals[local_eqn] +=
2014  (interpolated_value_proj - interpolated_value_bar) *
2015  psi[l] * W;
2016 
2017  // Calculate the jacobian
2018  if (flag == 1)
2019  {
2020  for (unsigned l2 = 0; l2 < n_value; l2++)
2021  {
2022  local_unknown = local_equation(fld, l2);
2023  if (local_unknown >= 0)
2024  {
2025  jacobian(local_eqn, local_unknown) +=
2026  psi[l2] * psi[l] * W;
2027  }
2028  }
2029  } // end of jacobian
2030  }
2031  }
2032  }
2033  // Flux -- need inner product
2034  else if (fld == 3)
2035  {
2036  // Field shape function
2037  Shape psi(n_value, n_dim);
2038 
2039  // Calculate jacobian and shape functions for that field
2040  double J = jacobian_and_shape_of_field(fld, s, psi);
2041 
2042  // Premultiply the weights and the Jacobian
2043  double W = w * J;
2044 
2045  // Value of flux in current element at current time level
2046  //(the unknown)
2047  unsigned now = 0;
2048  Vector<double> q_proj(n_dim);
2049  this->interpolated_q(now, s, q_proj);
2050 
2051  // Value of the interpolation of element located in base mesh
2052  Vector<double> q_bar(n_dim);
2053  cast_el_pt->interpolated_q(
2054  Time_level_for_projection, other_s, q_bar);
2055 
2056  // Loop over dofs of field fld
2057  for (unsigned l = 0; l < n_value; l++)
2058  {
2059  local_eqn = local_equation(fld, l);
2060  if (local_eqn >= 0)
2061  {
2062  // Loop over spatial dimension for dot product
2063  for (unsigned i = 0; i < n_dim; i++)
2064  {
2065  // Add to residuals
2066  residuals[local_eqn] +=
2067  (q_proj[i] - q_bar[i]) * psi(l, i) * W;
2068 
2069  // Calculate the jacobian
2070  if (flag == 1)
2071  {
2072  for (unsigned l2 = 0; l2 < n_value; l2++)
2073  {
2074  local_unknown = local_equation(fld, l2);
2075  if (local_unknown >= 0)
2076  {
2077  jacobian(local_eqn, local_unknown) +=
2078  psi(l2, i) * psi(l, i) * W;
2079  }
2080  }
2081  } // end of jacobian
2082  }
2083  }
2084  }
2085  }
2086  else
2087  {
2088  throw OomphLibError(
2089  "Wrong field specified. This should never happen\n",
2090  OOMPH_CURRENT_FUNCTION,
2091  OOMPH_EXCEPTION_LOCATION);
2092  }
2093 
2094 
2095  break;
2096 
2097  default:
2098  throw OomphLibError("Wrong type specificied in Projection_type. "
2099  "This should never happen\n",
2100  OOMPH_CURRENT_FUNCTION,
2101  OOMPH_EXCEPTION_LOCATION);
2102  }
2103  } // End of the switch statement
2104 
2105  } // End of loop over ipt
2106 
2107  } // End of residual_for_projection function
w
RowVector3d w
Definition: Matrix_resize_int.cpp:3
oomph::ElementWithExternalElement::external_element_local_coord
Vector< double > & external_element_local_coord(const unsigned &interaction_index, const unsigned &ipt)
Definition: element_with_external_element.h:136
oomph::ElementWithExternalElement::external_element_pt
FiniteElement *& external_element_pt(const unsigned &interaction_index, const unsigned &ipt)
Definition: element_with_external_element.h:107
oomph::FiniteElement::nnodal_position_type
unsigned nnodal_position_type() const
Definition: elements.h:2463
oomph::FiniteElement::interpolated_x
virtual double interpolated_x(const Vector< double > &s, const unsigned &i) const
Return FE interpolated coordinate x[i] at local coordinate s.
Definition: elements.cc:3962
oomph::FiniteElement::shape
virtual void shape(const Vector< double > &s, Shape &psi) const =0
oomph::FiniteElement::integral_pt
Integral *const & integral_pt() const
Return the pointer to the integration scheme (const version)
Definition: elements.h:1963
oomph::FiniteElement::FiniteElement
FiniteElement()
Constructor.
Definition: elements.h:1782
oomph::FiniteElement::J_eulerian
virtual double J_eulerian(const Vector< double > &s) const
Definition: elements.cc:4103
oomph::Integral::knot
virtual double knot(const unsigned &i, const unsigned &j) const =0
Return local coordinate s[j] of i-th integration point.
oomph::Integral::nweight
virtual unsigned nweight() const =0
Return the number of integration points of the scheme.
oomph::Integral::weight
virtual double weight(const unsigned &i) const =0
Return weight of i-th integration point.
oomph::ProjectableAxisymmetricPoroelasticityElement::local_equation
int local_equation(const unsigned &fld, const unsigned &j)
Return local equation number of value j in field fld.
Definition: axisym_poroelasticity_elements.h:1697
oomph::ProjectableAxisymmetricPoroelasticityElement::jacobian_and_shape_of_field
double jacobian_and_shape_of_field(const unsigned &fld, const Vector< double > &s, Shape &psi)
Definition: axisym_poroelasticity_elements.h:1526
oomph::ProjectableAxisymmetricPoroelasticityElement::nvalue_of_field
unsigned nvalue_of_field(const unsigned &fld)
Return number of values in field fld.
Definition: axisym_poroelasticity_elements.h:1661
oomph::ProjectableAxisymmetricPoroelasticityElement::get_field
double get_field(const unsigned &t, const unsigned &fld, const Vector< double > &s)
Definition: axisym_poroelasticity_elements.h:1610
oomph::ProjectableElementBase::Projected_lagrangian
unsigned Projected_lagrangian
Definition: projection.h:78
oomph::ProjectableElementBase::Projection_type
Projection_Type Projection_type
Definition: projection.h:83
oomph::ProjectableElementBase::Time_level_for_projection
unsigned Time_level_for_projection
Time level we are projecting (0=current values; >0: history values)
Definition: projection.h:70
oomph::ProjectableElementBase::Value
@ Value
Definition: projection.h:63
oomph::ProjectableElementBase::Lagrangian
@ Lagrangian
Definition: projection.h:62
oomph::ProjectableElementBase::Coordinate
@ Coordinate
Definition: projection.h:61
oomph::ProjectableElementBase::Projected_field
unsigned Projected_field
Field that is currently being projected.
Definition: projection.h:67
oomph::ProjectableElementBase::Projected_coordinate
unsigned Projected_coordinate
Definition: projection.h:74
k
char char char int int * k
Definition: level2_impl.h:374
oomph::QuadTreeNames::W
@ W
Definition: quadtree.h:63

References oomph::ProjectableElementBase::Coordinate, oomph::FiniteElement::dim(), oomph::ElementWithExternalElement::external_element_local_coord(), oomph::ElementWithExternalElement::external_element_pt(), oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::get_field(), oomph::ProjectableElementBase::get_field(), i, oomph::FiniteElement::integral_pt(), oomph::FiniteElement::interpolated_x(), J, oomph::FiniteElement::J_eulerian(), oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::jacobian_and_shape_of_field(), k, oomph::Integral::knot(), oomph::ProjectableElementBase::Lagrangian, oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::local_equation(), oomph::FiniteElement::nnodal_position_type(), oomph::FiniteElement::nnode(), oomph::ProjectableAxisymmetricPoroelasticityElement< AXISYMMETRIC_POROELASTICITY_ELEMENT >::nvalue_of_field(), oomph::Integral::nweight(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, oomph::SolidFiniteElement::position_local_eqn(), oomph::ProjectableElementBase::Projected_coordinate, oomph::ProjectableElementBase::Projected_field, oomph::ProjectableElementBase::Projected_lagrangian, oomph::ProjectableElementBase::Projection_type, s, oomph::FiniteElement::shape(), oomph::ProjectableElementBase::Time_level_for_projection, oomph::ProjectableElementBase::Value, w, oomph::QuadTreeNames::W, and oomph::Integral::weight().

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