Public Member Functions | Static Private Attributes | List of all members
oomph::TElement< 2, NNODE_1D > Class Template Reference

#include <Telements.h>

+ Inheritance diagram for oomph::TElement< 2, NNODE_1D >:

Public Member Functions

 TElement ()
 Constructor. More...
 
 TElement (const bool &allow_high_order)
 Alternative constructor. More...
 
 TElement (const TElement &)=delete
 Broken copy constructor. More...
 
 ~TElement ()
 Broken assignment operator. More...
 
unsigned nnode_1d () const
 Number of nodes along each element edge. More...
 
unsigned nvertex_node () const
 
unsigned get_bulk_node_number (const int &face_index, const unsigned &i) const
 Public access function for Node_on_face. More...
 
Node * vertex_node_pt (const unsigned &j) const
 Pointer to the j-th vertex node in the element. More...
 
void shape (const Vector< double > &s, Shape &psi) const
 Calculate the geometric shape functions at local coordinate s. More...
 
void dshape_local (const Vector< double > &s, Shape &psi, DShape &dpsids) const
 
void d2shape_local (const Vector< double > &s, Shape &psi, DShape &dpsids, DShape &d2psids) const
 
double invert_jacobian_mapping (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
double s_min () const
 Min. value of local coordinate. More...
 
double s_max () const
 Max. value of local coordinate. More...
 
void local_coordinate_of_node (const unsigned &j, Vector< double > &s) const
 Return local coordinates of node j. More...
 
unsigned nplot_points_paraview (const unsigned &nplot) const
 
unsigned nsub_elements_paraview (const unsigned &nplot) const
 
void write_paraview_output_offset_information (std::ofstream &file_out, const unsigned &nplot, unsigned &counter) const
 
void write_paraview_type (std::ofstream &file_out, const unsigned &nplot) const
 
void write_paraview_offsets (std::ofstream &file_out, const unsigned &nplot, unsigned &offset_sum) const
 
void output (std::ostream &output)
 Output. More...
 
void output (std::ostream &outfile, const unsigned &nplot)
 Output at specified number of plot points. More...
 
void output (FILE *file_pt)
 C-style output. More...
 
void output (FILE *file_pt, const unsigned &n_plot)
 C_style output at n_plot points. More...
 
void get_s_plot (const unsigned &iplot, const unsigned &nplot, Vector< double > &s, const bool &use_equally_spaced_interior_sample_points=false) const
 
std::string tecplot_zone_string (const unsigned &nplot) const
 
void write_tecplot_zone_footer (std::ostream &outfile, const unsigned &nplot) const
 
void write_tecplot_zone_footer (FILE *file_pt, const unsigned &nplot) const
 
unsigned nplot_points (const unsigned &nplot) const
 
void build_face_element (const int &face_index, FaceElement *face_element_pt)
 
- Public Member Functions inherited from oomph::TElementBase
 TElementBase ()
 Empty default constructor. More...
 
 TElementBase (const TElementBase &)=delete
 Broken copy constructor. More...
 
ElementGeometry::ElementGeometry element_geometry () const
 Broken assignment operator. More...
 
bool local_coord_is_valid (const Vector< double > &s)
 Check whether the local coordinates are valid or not. More...
 
void move_local_coord_back_into_element (Vector< double > &s) const
 
- Public Member Functions inherited from oomph::TElementGeometricBase
 TElementGeometricBase ()
 Empty default constructor. More...
 
 TElementGeometricBase (const TElementGeometricBase &)=delete
 Broken copy constructor. More...
 
- 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...
 
void get_centre_of_gravity_and_max_radius_in_terms_of_zeta (Vector< double > &cog, double &max_radius) 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_at_knot (const unsigned &ipt, Shape &psi) const
 
virtual void dshape_local_at_knot (const unsigned &ipt, Shape &psi, DShape &dpsids) 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_local_dofs (std::ostream &out, const std::string &current_string) const
 
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...
 
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 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 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)
 
virtual double zeta_nodal (const unsigned &n, const unsigned &k, const unsigned &i) const
 
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)
 
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)
 
void output_paraview (std::ofstream &file_out, const unsigned &nplot) 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 (const unsigned &t, std::ostream &outfile, const unsigned &n_plot) const
 
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 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 unsigned self_test ()
 
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
 

Static Private Attributes

static const unsigned Node_on_face [3][NNODE_1D]
 Nodal translation scheme for use when generating face elements. More...
 
static TGauss< 2, NNODE_1D > Default_integration_scheme
 

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 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 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)
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
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_internal_and_external_local_eqn_numbers (const bool &store_local_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)
 
virtual void fill_in_contribution_to_residuals (Vector< double > &residuals)
 
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::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<unsigned NNODE_1D>
class oomph::TElement< 2, NNODE_1D >

General TElement class specialised to two spatial dimensions Ordering of nodes as in Zienkiwizc sketches: vertex nodes 0 - 1 - 2 anticlockwise. Midside nodes filled in progressing along the consecutive edges. Central node(s) come(s) last.

Constructor & Destructor Documentation

◆ TElement() [1/3]

template<unsigned NNODE_1D>
oomph::TElement< 2, NNODE_1D >::TElement ( )
inline

Constructor.

1521  {
1522  // Number of nodes
1523  switch (NNODE_1D)
1524  {
1525  case 2:
1526  case 3:
1527  case 4:
1528  break;
1529 
1530  default:
1531  std::string error_message =
1532  "Triangles are currently only implemented for\n";
1533  error_message += "three and six nodes, i.e. NNODE_1D=2 , 3 , 4\n";
1534 
1535  throw OomphLibError(
1536  error_message, OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1537  }
1538 
1539  // Set the number of nodes
1540  unsigned n_node = (NNODE_1D * (NNODE_1D + 1)) / 2;
1541  this->set_n_node(n_node);
1542 
1543  // Set the elemental and nodal dimension
1544  set_dimension(2);
1545 
1546  // Assign default (full) spatial integration scheme
1547  set_integration_scheme(&Default_integration_scheme);
1548  }
oomph::FiniteElement::set_dimension
void set_dimension(const unsigned &dim)
Definition: elements.h:1380
oomph::FiniteElement::set_n_node
void set_n_node(const unsigned &n)
Definition: elements.h:1404
oomph::FiniteElement::set_integration_scheme
virtual void set_integration_scheme(Integral *const &integral_pt)
Set the spatial integration scheme.
Definition: elements.cc:3210
oomph::TElement< 2, NNODE_1D >::Default_integration_scheme
static TGauss< 2, NNODE_1D > Default_integration_scheme
Definition: Telements.h:1515
oomph::Global_string_for_annotation::string
std::string string(const unsigned &i)
Definition: oomph_definitions.cc:286
OOMPH_EXCEPTION_LOCATION
#define OOMPH_EXCEPTION_LOCATION
Definition: oomph_definitions.h:61
OOMPH_CURRENT_FUNCTION
#define OOMPH_CURRENT_FUNCTION
Definition: oomph_definitions.h:86

References OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, and oomph::Global_string_for_annotation::string().

◆ TElement() [2/3]

template<unsigned NNODE_1D>
oomph::TElement< 2, NNODE_1D >::TElement ( const bool &  allow_high_order)
inline

Alternative constructor.

1552  {
1553  // Check if we are overriding the restriction on NNODE_1D
1554  if (!allow_high_order)
1555  {
1556  // Call the default constructor
1557  TElement<2, NNODE_1D>();
1558  }
1559  else
1560  {
1561  // Set the number of nodes
1562  unsigned n_node = (NNODE_1D * (NNODE_1D + 1)) / 2;
1563  this->set_n_node(n_node);
1564 
1565  // Set the elemental and nodal dimension
1566  set_dimension(2);
1567 
1568  // Assign default (full) spatial integration scheme
1569  set_integration_scheme(&Default_integration_scheme);
1570  }
1571  }

◆ TElement() [3/3]

template<unsigned NNODE_1D>
oomph::TElement< 2, NNODE_1D >::TElement ( const TElement< 2, NNODE_1D > &  )
delete

Broken copy constructor.

◆ ~TElement()

template<unsigned NNODE_1D>
oomph::TElement< 2, NNODE_1D >::~TElement ( )
inline

Broken assignment operator.

Destructor

1582 {}

Member Function Documentation

◆ build_face_element()

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::build_face_element ( const int &  face_index,
FaceElement *  face_element_pt 
)
virtual

Build the lower-dimensional FaceElement (an element of type TElement<1,NNODE_1D>). The face index takes three possible values: 0 (Left) s[0] = 0.0 1 (Bottom) s[1] = 0.0 2 (Sloping face) s[0] = 1.0 - s[1]

Function to setup geometrical information for lower-dimensional FaceElements (which are of type TElement<2,NNODE_1D>).

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::SolidTElement< 2, NNODE_1D >.

372  {
373  // Set the nodal dimension from the first node
374  face_element_pt->set_nodal_dimension(nodal_dimension());
375 
376  // Set the pointer to the orginal "bulk" element
377  face_element_pt->bulk_element_pt() = this;
378 
379 #ifdef OOMPH_HAS_MPI
380  // Pass on non-halo proc ID
381  face_element_pt->set_halo(Non_halo_proc_ID);
382 #endif
383 
384  // Calculate the number of nodes in the face element
385  const unsigned n_face_nodes = NNODE_1D;
386 
387  // Resize storage for the number of values originally stored
388  // at the face element's nodes.
389  face_element_pt->nbulk_value_resize(n_face_nodes);
390 
391  // Resize storage for the bulk node numbers corresponding to
392  // the nodes of the face
393  face_element_pt->bulk_node_number_resize(n_face_nodes);
394 
395  // Set the face index in the face element
396  face_element_pt->face_index() = face_index;
397 
398  // So the faces are
399  // 0 : s_0 fixed
400  // 1 : s_1 fixed
401  // 2 : sloping face
402 
403  // Copy nodes across to the face
404  for (unsigned i = 0; i < n_face_nodes; i++)
405  {
406  // The number is just offset by one
407  unsigned bulk_number = Node_on_face[face_index][i];
408  face_element_pt->node_pt(i) = node_pt(bulk_number);
409  face_element_pt->bulk_node_number(i) = bulk_number;
410  // set the number of values originally stored at this node
411  face_element_pt->nbulk_value(i) = required_nvalue(bulk_number);
412  }
413 
414  // Now set up the node pointers and the normal vectors
415  switch (face_index)
416  {
417  //
418  //-----The face s0 is constant
419  case 0:
420 
421  // Set the pointer to the function that determines the bulk
422  // coordinates in the face element
423  face_element_pt->face_to_bulk_coordinate_fct_pt() =
424  &TElement2FaceToBulkCoordinates::face0;
425 
426  // Set the pointer to the function that determines the mapping of
427  // derivatives
428  face_element_pt->bulk_coordinate_derivatives_fct_pt() =
429  &TElement2BulkCoordinateDerivatives::face0;
430 
431  // Outer unit normal is the positive cross product of two in plane
432  // tangent vectors
433  face_element_pt->normal_sign() = +1;
434 
435  break;
436 
437  //-----The face s1 is constant
438  case 1:
439 
440  // Set the pointer to the function that determines the bulk
441  // coordinates in the face element
442  face_element_pt->face_to_bulk_coordinate_fct_pt() =
443  &TElement2FaceToBulkCoordinates::face1;
444 
445  // Set the pointer to the function that determines the mapping of
446  // derivatives
447  face_element_pt->bulk_coordinate_derivatives_fct_pt() =
448  &TElement2BulkCoordinateDerivatives::face1;
449 
450  // Outer unit normal is the positive cross product of two in plane
451  // tangent vectors
452  face_element_pt->normal_sign() = +1;
453 
454  break;
455 
456  //
457  //-----The face s2 is constant
458  case 2:
459 
460  // Set the pointer to the function that determines the bulk
461  // coordinates in the face element
462  face_element_pt->face_to_bulk_coordinate_fct_pt() =
463  &TElement2FaceToBulkCoordinates::face2;
464 
465  // Set the pointer to the function that determines the mapping of
466  // derivatives
467  face_element_pt->bulk_coordinate_derivatives_fct_pt() =
468  &TElement2BulkCoordinateDerivatives::face2;
469 
470  // Outer unit normal is the negative cross product of two in plane
471  // tangent vectors
472  face_element_pt->normal_sign() = -1;
473 
474  break;
475 
476  default:
477 
478  std::ostringstream error_message;
479  error_message << "Face_index should only take "
480  << "the values 0, 1 or 2 not " << face_index << std::endl;
481 
482  throw OomphLibError(error_message.str(),
483  OOMPH_CURRENT_FUNCTION,
484  OOMPH_EXCEPTION_LOCATION);
485  } // end switch
486  }
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9
oomph::FiniteElement::node_pt
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175
oomph::FiniteElement::required_nvalue
virtual unsigned required_nvalue(const unsigned &n) const
Definition: elements.h:2455
oomph::FiniteElement::nodal_dimension
unsigned nodal_dimension() const
Return the required Eulerian dimension of the nodes in this element.
Definition: elements.h:2484
oomph::TElement< 2, NNODE_1D >::Node_on_face
static const unsigned Node_on_face[3][NNODE_1D]
Nodal translation scheme for use when generating face elements.
Definition: Telements.h:1508
oomph::TElement2BulkCoordinateDerivatives::face0
void face0(const Vector< double > &s, DenseMatrix< double > &dsbulk_dsface, unsigned &interior_direction)
Function for the "left" face along which s0 is fixed.
Definition: Telements.cc:320
oomph::TElement2BulkCoordinateDerivatives::face1
void face1(const Vector< double > &s, DenseMatrix< double > &dsbulk_dsface, unsigned &interior_direction)
Function for the "bottom" face along which s1 is fixed.
Definition: Telements.cc:335
oomph::TElement2BulkCoordinateDerivatives::face2
void face2(const Vector< double > &s, DenseMatrix< double > &dsbulk_dsface, unsigned &interior_direction)
Function for the sloping face.
Definition: Telements.cc:349
oomph::TElement2FaceToBulkCoordinates::face1
void face1(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the face s1 = 0.
Definition: Telements.cc:297
oomph::TElement2FaceToBulkCoordinates::face2
void face2(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the face s2 = 0.
Definition: Telements.cc:304
oomph::TElement2FaceToBulkCoordinates::face0
void face0(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the face s0 = 0.
Definition: Telements.cc:290

References oomph::FaceElement::bulk_coordinate_derivatives_fct_pt(), oomph::FaceElement::bulk_element_pt(), oomph::FaceElement::bulk_node_number(), oomph::FaceElement::bulk_node_number_resize(), oomph::TElement2BulkCoordinateDerivatives::face0(), oomph::TElement2FaceToBulkCoordinates::face0(), oomph::TElement2BulkCoordinateDerivatives::face1(), oomph::TElement2FaceToBulkCoordinates::face1(), oomph::TElement2BulkCoordinateDerivatives::face2(), oomph::TElement2FaceToBulkCoordinates::face2(), oomph::FaceElement::face_index(), oomph::FaceElement::face_to_bulk_coordinate_fct_pt(), i, oomph::FaceElement::nbulk_value(), oomph::FaceElement::nbulk_value_resize(), oomph::TElement< DIM, NNODE_1D >::Node_on_face(), oomph::FiniteElement::node_pt(), oomph::FaceElement::normal_sign(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, and oomph::FiniteElement::set_nodal_dimension().

◆ d2shape_local()

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::d2shape_local ( const Vector< double > &  s,
Shape &  psi,
DShape &  dpsids,
DShape &  d2psids 
) const
inlinevirtual

Compute the geometric shape functions, derivatives and second derivatives w.r.t local coordinates at local coordinate s d2psids(i,0) = \( \partial^2 \psi_j / \partial s_0^2 \) d2psids(i,1) = \( \partial^2 \psi_j / \partial s_1^2 \) d2psids(i,2) = \( \partial^2 \psi_j / \partial s_0 \partial s_1 \)

Reimplemented from oomph::FiniteElement.

1647  {
1648  TElementShape<2, NNODE_1D>::d2shape_local(s, psi, dpsids, d2psids);
1649  }
s
RealScalar s
Definition: level1_cplx_impl.h:130

References s.

◆ dshape_local()

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::dshape_local ( const Vector< double > &  s,
Shape &  psi,
DShape &  dpsids 
) const
inlinevirtual

Compute the geometric shape functions and derivatives w.r.t. local coordinates at local coordinate s

Reimplemented from oomph::FiniteElement.

1634  {
1635  TElementShape<2, NNODE_1D>::dshape_local(s, psi, dpsids);
1636  }

References s.

◆ get_bulk_node_number()

template<unsigned NNODE_1D>
unsigned oomph::TElement< 2, NNODE_1D >::get_bulk_node_number ( const int &  face_index,
const unsigned &  i 
) const
inlinevirtual

Public access function for Node_on_face.

Reimplemented from oomph::FiniteElement.

1600  {
1601  return Node_on_face[face_index][i];
1602  }

References i, and oomph::TElement< DIM, NNODE_1D >::Node_on_face().

◆ get_s_plot()

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::get_s_plot ( const unsigned &  iplot,
const unsigned &  nplot,
Vector< double > &  s,
const bool &  use_equally_spaced_interior_sample_points = false 
) const
inlinevirtual

Get vector of local coordinates of plot point i (when plotting nplot points in each "coordinate direction").

Reimplemented from oomph::FiniteElement.

1786  {
1787  if (nplot > 1)
1788  {
1789  unsigned np = 0, i, j;
1790  for (i = 0; i < nplot; ++i)
1791  {
1792  for (j = 0; j < nplot - i; ++j)
1793  {
1794  if (np == iplot)
1795  {
1796  s[0] = double(j) / double(nplot - 1);
1797  s[1] = double(i) / double(nplot - 1);
1798  if (use_equally_spaced_interior_sample_points)
1799  {
1800  double range = 1.0;
1801  double dx_new = range / (double(nplot) + 0.5);
1802  double range_new = double(nplot - 1) * dx_new;
1803  s[0] = 0.5 * dx_new + range_new * s[0] / range;
1804  s[1] = 0.5 * dx_new + range_new * s[1] / range;
1805  }
1806  return;
1807  }
1808  ++np;
1809  }
1810  }
1811  }
1812  else
1813  {
1814  s[0] = 1.0 / 3.0;
1815  s[1] = 1.0 / 3.0;
1816  }
1817  }
j
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References i, j, and s.

◆ invert_jacobian_mapping()

template<unsigned NNODE_1D>
double oomph::TElement< 2, NNODE_1D >::invert_jacobian_mapping ( const DenseMatrix< double > &  jacobian,
DenseMatrix< double > &  inverse_jacobian 
) const
inlinevirtual

Overload the template-free interface for the calculation of inverse jacobian matrix. This is a two dimensional element, so use the 2D version.

Reimplemented from oomph::FiniteElement.

1656  {
1657  return FiniteElement::invert_jacobian<2>(jacobian, inverse_jacobian);
1658  }

◆ local_coordinate_of_node()

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::local_coordinate_of_node ( const unsigned &  j,
Vector< double > &  s 
) const
inlinevirtual

Return local coordinates of node j.

Reimplemented from oomph::FiniteElement.

1675  {
1676  TElementShape<2, NNODE_1D>::local_coordinate_of_node(j, s);
1677  }

References j, and s.

◆ nnode_1d()

template<unsigned NNODE_1D>
unsigned oomph::TElement< 2, NNODE_1D >::nnode_1d ( ) const
inlinevirtual

Number of nodes along each element edge.

Reimplemented from oomph::FiniteElement.

1586  {
1587  return NNODE_1D;
1588  }

◆ nplot_points()

template<unsigned NNODE_1D>
unsigned oomph::TElement< 2, NNODE_1D >::nplot_points ( const unsigned &  nplot) const
inlinevirtual

Return total number of plot points (when plotting nplot points in each "coordinate direction)

Reimplemented from oomph::FiniteElement.

1900  {
1901  unsigned np = 0;
1902  for (unsigned i = 1; i <= nplot; i++)
1903  {
1904  np += i;
1905  }
1906  return np;
1907  }

References i.

◆ nplot_points_paraview()

template<unsigned NNODE_1D>
unsigned oomph::TElement< 2, NNODE_1D >::nplot_points_paraview ( const unsigned &  nplot) const
inlinevirtual

Return the number of actual plot points for paraview plot with parameter nplot.

Reimplemented from oomph::FiniteElement.

1682  {
1683  unsigned node_sum = 0;
1684  for (unsigned i = 1; i <= nplot; i++)
1685  {
1686  node_sum += i;
1687  }
1688  return node_sum;
1689  }

References i.

◆ nsub_elements_paraview()

template<unsigned NNODE_1D>
unsigned oomph::TElement< 2, NNODE_1D >::nsub_elements_paraview ( const unsigned &  nplot) const
inlinevirtual

Return the number of local sub-elements for paraview plot with parameter nplot.

Reimplemented from oomph::FiniteElement.

1694  {
1695  unsigned local_sum = 0;
1696  for (unsigned i = 1; i < nplot; i++)
1697  {
1698  local_sum += 2 * (nplot - i - 1) + 1;
1699  }
1700  return local_sum;
1701  }

References i.

◆ nvertex_node()

template<unsigned NNODE_1D>
unsigned oomph::TElement< 2, NNODE_1D >::nvertex_node ( ) const
inlinevirtual

Number of vertex nodes in the element: One more than spatial dimension

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::TTimeHarmonicFourierDecomposedLinearElasticityElement< NNODE_1D >, and oomph::TPMLFourierDecomposedHelmholtzElement< NNODE_1D >.

1593  {
1594  return 3;
1595  }

◆ output() [1/4]

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::output ( FILE *  file_pt)
virtual

C-style output.

The C-style output function for TElement<2,NNODE_1D>

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::TTimeHarmonicFourierDecomposedLinearElasticityElement< NNODE_1D >, and oomph::TPMLFourierDecomposedHelmholtzElement< NNODE_1D >.

539  {
540  output(file_pt, NNODE_1D);
541  }
oomph::TElement< 2, NNODE_1D >::output
void output(std::ostream &output)
Output.
Definition: Telements.cc:493

References oomph::output().

◆ output() [2/4]

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::output ( FILE *  file_pt,
const unsigned &  n_plot 
)
virtual

C_style output at n_plot points.

The C-style output function for TElement<2,NNODE_1D>

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::TTimeHarmonicFourierDecomposedLinearElasticityElement< NNODE_1D >, and oomph::TPMLFourierDecomposedHelmholtzElement< NNODE_1D >.

549  {
550  // Vector of local coordinates
551  Vector<double> s(2);
552 
553  // Find the dimensions of the nodes
554  unsigned n_dim = this->nodal_dimension();
555 
556  // Tecplot header info
557  fprintf(file_pt, "%s \n", tecplot_zone_string(nplot).c_str());
558 
559  // Loop over plot points
560  unsigned num_plot_points = nplot_points(nplot);
561  for (unsigned iplot = 0; iplot < num_plot_points; iplot++)
562  {
563  // Get local coordinates of plot point
564  get_s_plot(iplot, nplot, s);
565 
566  for (unsigned i = 0; i < n_dim; i++)
567  {
568  fprintf(file_pt, "%g ", interpolated_x(s, i));
569  // outfile << interpolated_x(s,i) << " ";
570  }
571  fprintf(file_pt, "\n");
572  // outfile << std::endl;
573  }
574 
575  // Write tecplot footer (e.g. FE connectivity lists)
576  write_tecplot_zone_footer(file_pt, nplot);
577  }
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::TElement< 2, NNODE_1D >::write_tecplot_zone_footer
void write_tecplot_zone_footer(std::ostream &outfile, const unsigned &nplot) const
Definition: Telements.h:1838
oomph::TElement< 2, NNODE_1D >::tecplot_zone_string
std::string tecplot_zone_string(const unsigned &nplot) const
Definition: Telements.h:1821
oomph::TElement< 2, NNODE_1D >::get_s_plot
void get_s_plot(const unsigned &iplot, const unsigned &nplot, Vector< double > &s, const bool &use_equally_spaced_interior_sample_points=false) const
Definition: Telements.h:1781
oomph::TElement< 2, NNODE_1D >::nplot_points
unsigned nplot_points(const unsigned &nplot) const
Definition: Telements.h:1899

References i, and s.

◆ output() [3/4]

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::output ( std::ostream &  outfile,
const unsigned &  nplot 
)
virtual

Output at specified number of plot points.

The output function for TElement<2,NNODE_1D>

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::TTimeHarmonicFourierDecomposedLinearElasticityElement< NNODE_1D >, and oomph::TPMLFourierDecomposedHelmholtzElement< NNODE_1D >.

506  {
507  // Vector of local coordinates
508  Vector<double> s(2);
509 
510  // Get the dimension of the node
511  unsigned n_dim = this->nodal_dimension();
512 
513  // Tecplot header info
514  outfile << tecplot_zone_string(nplot);
515 
516  // Loop over plot points
517  unsigned num_plot_points = nplot_points(nplot);
518  for (unsigned iplot = 0; iplot < num_plot_points; iplot++)
519  {
520  // Get local coordinates of plot point
521  get_s_plot(iplot, nplot, s);
522 
523  for (unsigned i = 0; i < n_dim; i++)
524  {
525  outfile << interpolated_x(s, i) << " ";
526  }
527  outfile << std::endl;
528  }
529 
530  // Write tecplot footer (e.g. FE connectivity lists)
531  write_tecplot_zone_footer(outfile, nplot);
532  }

References i, and s.

◆ output() [4/4]

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::output ( std::ostream &  output)
virtual

Output.

The output function for TElement<2,NNODE_1D>

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::TTimeHarmonicFourierDecomposedLinearElasticityElement< NNODE_1D >, and oomph::TPMLFourierDecomposedHelmholtzElement< NNODE_1D >.

494  {
495  // QUEHACERES want to perform same output, but at each node
496  output(outfile, NNODE_1D);
497  }

References oomph::output().

◆ s_max()

template<unsigned NNODE_1D>
double oomph::TElement< 2, NNODE_1D >::s_max ( ) const
inlinevirtual

Max. value of local coordinate.

Reimplemented from oomph::FiniteElement.

1668  {
1669  return 1.0;
1670  }

◆ s_min()

template<unsigned NNODE_1D>
double oomph::TElement< 2, NNODE_1D >::s_min ( ) const
inlinevirtual

Min. value of local coordinate.

Reimplemented from oomph::FiniteElement.

1662  {
1663  return 0.0;
1664  }

◆ shape()

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::shape ( const Vector< double > &  s,
Shape &  psi 
) const
inlinevirtual

Calculate the geometric shape functions at local coordinate s.

Implements oomph::FiniteElement.

1625  {
1626  TElementShape<2, NNODE_1D>::shape(s, psi);
1627  }
oomph::OneDimLagrange::shape
void shape(const double &s, double *Psi)
Definition: shape.h:564

References s, and oomph::OneDimLagrange::shape().

◆ tecplot_zone_string()

template<unsigned NNODE_1D>
std::string oomph::TElement< 2, NNODE_1D >::tecplot_zone_string ( const unsigned &  nplot) const
inlinevirtual

Return string for tecplot zone header (when plotting nplot points in each "coordinate direction)

Reimplemented from oomph::FiniteElement.

1822  {
1823  std::ostringstream header;
1824  unsigned nel = 0;
1825  for (unsigned i = 1; i < nplot; i++)
1826  {
1827  nel += 2 * i - 1;
1828  }
1829  header << "ZONE N=" << nplot_points(nplot) << ", E=" << nel
1830  << ", F=FEPOINT, ET=TRIANGLE\n";
1831  return header.str();
1832  }

References i.

◆ vertex_node_pt()

template<unsigned NNODE_1D>
Node* oomph::TElement< 2, NNODE_1D >::vertex_node_pt ( const unsigned &  j) const
inlinevirtual

Pointer to the j-th vertex node in the element.

Reimplemented from oomph::FiniteElement.

Reimplemented in oomph::TTimeHarmonicFourierDecomposedLinearElasticityElement< NNODE_1D >, and oomph::TPMLFourierDecomposedHelmholtzElement< NNODE_1D >.

1606  {
1607  // Vertex nodes come first:
1608 #ifdef PARANOID
1609  if (j > 2)
1610  {
1611  std::ostringstream error_message;
1612  error_message
1613  << "Element only has three vertex nodes; called with node number "
1614  << j << std::endl;
1615  throw OomphLibError(error_message.str(),
1616  OOMPH_CURRENT_FUNCTION,
1617  OOMPH_EXCEPTION_LOCATION);
1618  }
1619 #endif
1620  return node_pt(j);
1621  }

References j, OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

◆ write_paraview_offsets()

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::write_paraview_offsets ( std::ofstream &  file_out,
const unsigned &  nplot,
unsigned &  offset_sum 
) const
inlinevirtual

Return the offsets for the paraview sub-elements. Needs to be implemented for each new geometric element type; see http://www.vtk.org/VTK/img/file-formats.pdf

Reimplemented from oomph::FiniteElement.

1755  {
1756  unsigned local_loop = nsub_elements_paraview(nplot);
1757 
1758  // Loop over all local elements and add its offset to the overall
1759  // offset_sum
1760  for (unsigned i = 0; i < local_loop; i++)
1761  {
1762  offset_sum += 3;
1763  file_out << offset_sum << std::endl;
1764  }
1765  }
oomph::TElement< 2, NNODE_1D >::nsub_elements_paraview
unsigned nsub_elements_paraview(const unsigned &nplot) const
Definition: Telements.h:1693

References i.

◆ write_paraview_output_offset_information()

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::write_paraview_output_offset_information ( std::ofstream &  file_out,
const unsigned &  nplot,
unsigned &  counter 
) const
inlinevirtual

Fill in the offset information for paraview plot. Needs to be implemented for each new geometric element type; see http://www.vtk.org/VTK/img/file-formats.pdf

Reimplemented from oomph::FiniteElement.

1709  {
1710  // Outputs list of connectivity of Paraview elements,
1711  // whilst remembering the overall ordering
1712 
1713  unsigned node_count = 0;
1714  for (unsigned i = 0; i < nplot - 1; i++)
1715  {
1716  for (unsigned j = 0; j < nplot - i - 1; j++)
1717  {
1718  file_out << j + node_count + counter << " "
1719  << j + node_count + 1 + counter << " "
1720  << j + nplot + node_count - i + counter << std::endl;
1721 
1722  if (j < nplot - i - 2)
1723  {
1724  file_out << j + node_count + 1 + counter << " "
1725  << j + nplot + node_count - i + 1 + counter << " "
1726  << j + nplot + node_count - i + counter << std::endl;
1727  }
1728  }
1729  node_count += (nplot - i);
1730  }
1731  counter += nplot_points_paraview(nplot);
1732  }
oomph::TElement< 2, NNODE_1D >::nplot_points_paraview
unsigned nplot_points_paraview(const unsigned &nplot) const
Definition: Telements.h:1681

References i, and j.

◆ write_paraview_type()

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::write_paraview_type ( std::ofstream &  file_out,
const unsigned &  nplot 
) const
inlinevirtual

Return the paraview element type. Needs to be implemented for each new geometric element type; see http://www.vtk.org/VTK/img/file-formats.pdf

Reimplemented from oomph::FiniteElement.

1739  {
1740  unsigned local_loop = nsub_elements_paraview(nplot);
1741 
1742  // Loop over all the local elements and print its paraview type
1743  for (unsigned i = 0; i < local_loop; i++)
1744  {
1745  file_out << "5" << std::endl;
1746  }
1747  }

References i.

◆ write_tecplot_zone_footer() [1/2]

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::write_tecplot_zone_footer ( FILE *  file_pt,
const unsigned &  nplot 
) const
inlinevirtual

Add tecplot zone "footer" to C-style output. (when plotting nplot points in each "coordinate direction). Empty by default – can be used, e.g., to add FE connectivity lists to elements that need it.

Reimplemented from oomph::FiniteElement.

1868  {
1869  // Output node lists for sub elements for Tecplot (node index
1870  // must start at 1)
1871  unsigned nod_count = 1;
1872  for (unsigned i = 0; i < nplot; i++)
1873  {
1874  for (unsigned j = 0; j < nplot - i; j++)
1875  {
1876  if (j < nplot - i - 1)
1877  {
1878  fprintf(file_pt,
1879  "%i %i %i \n",
1880  nod_count,
1881  nod_count + 1,
1882  nod_count + nplot - i);
1883  if (j < nplot - i - 2)
1884  {
1885  fprintf(file_pt,
1886  "%i %i %i \n",
1887  nod_count + 1,
1888  nod_count + nplot - i + 1,
1889  nod_count + nplot - i);
1890  }
1891  }
1892  ++nod_count;
1893  }
1894  }
1895  }

References i, and j.

◆ write_tecplot_zone_footer() [2/2]

template<unsigned NNODE_1D>
void oomph::TElement< 2, NNODE_1D >::write_tecplot_zone_footer ( std::ostream &  outfile,
const unsigned &  nplot 
) const
inlinevirtual

Add tecplot zone "footer" to output stream (when plotting nplot points in each "coordinate direction). Empty by default – can be used, e.g., to add FE connectivity lists to elements that need it.

Reimplemented from oomph::FiniteElement.

1840  {
1841  // Output node lists for sub elements for Tecplot (node index
1842  // must start at 1)
1843  unsigned nod_count = 1;
1844  for (unsigned i = 0; i < nplot; i++)
1845  {
1846  for (unsigned j = 0; j < nplot - i; j++)
1847  {
1848  if (j < nplot - i - 1)
1849  {
1850  outfile << nod_count << " " << nod_count + 1 << " "
1851  << nod_count + nplot - i << std::endl;
1852  if (j < nplot - i - 2)
1853  {
1854  outfile << nod_count + 1 << " " << nod_count + nplot - i + 1
1855  << " " << nod_count + nplot - i << std::endl;
1856  }
1857  }
1858  ++nod_count;
1859  }
1860  }
1861  }

References i, and j.

Member Data Documentation

◆ Default_integration_scheme

template<unsigned NNODE_1D>
TGauss< 2, NNODE_1D > oomph::TElement< 2, NNODE_1D >::Default_integration_scheme
staticprivate

Default integration rule: Gaussian integration of same 'order' as the element

◆ Node_on_face

template<unsigned NNODE_1D>
const unsigned oomph::TElement< 2, NNODE_1D >::Node_on_face[3][NNODE_1D]
staticprivate

Nodal translation scheme for use when generating face elements.

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