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

#include <Telements.h>

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

Public Member Functions

 TElement ()
 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 [4][(NNODE_1D *(NNODE_1D+1))/2]
 Nodal translation scheme for use when generating face elements. More...
 
static TGauss< 3, 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< 3, NNODE_1D >

General TElement class specialised to three spatial dimensions (tet) Ordering of nodes inverted from Zienkiewizc sketches: When looking into the tet from vertex node 0. The vertex nodes on the opposite face are 1 - 2 - 3 in anticlockwise direction. Other nodes filled in edge by edge, then the face ones, then the internal ones.

Constructor & Destructor Documentation

◆ TElement() [1/2]

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

Constructor.

3057  {
3058  switch (NNODE_1D)
3059  {
3060  case 2:
3061  case 3:
3062  break;
3063 
3064  /* case 4: */
3065  /* n_node = 20; */
3066  /* break; */
3067 
3068  default:
3069  std::string error_message =
3070  "Tets are currently only implemented for\n";
3071  error_message += "four and ten nodes, i.e. NNODE_1D=2 , 3 \n";
3072 
3073  throw OomphLibError(
3074  error_message, OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
3075  }
3076 
3077 
3078  // Set the number of nodes
3079  unsigned n_node =
3080  (NNODE_1D * (NNODE_1D + 1)) / 2 + 1 + 3 * (NNODE_1D - 2);
3081  this->set_n_node(n_node);
3082 
3083  // Set the elemental and nodal dimensions
3084  set_dimension(3);
3085 
3086  // Assign default (full) spatial integration scheme
3087  set_integration_scheme(&Default_integration_scheme);
3088  }
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< 3, NNODE_1D >::Default_integration_scheme
static TGauss< 3, NNODE_1D > Default_integration_scheme
Definition: Telements.h:3052
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/2]

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

Broken copy constructor.

◆ ~TElement()

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

Broken assignment operator.

Destructor

3099 {}

Member Function Documentation

◆ build_face_element()

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

Build the lower-dimensional FaceElement (an element of type TElement<2,NNODE_1D>). The face index can take one of four values corresponding to the four possible faces: 0: (left) s[0] = 0.0 1: (bottom) s[1] = 0.0 2: (back) s[2] = 0.0 3: (sloping face) s[0] + s[1] + s[2] = 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< 3, NNODE_1D >.

735  {
736  // Set the nodal dimension
737  face_element_pt->set_nodal_dimension(nodal_dimension());
738 
739  // Set the pointer to the orginal "bulk" element
740  face_element_pt->bulk_element_pt() = this;
741 
742 #ifdef OOMPH_HAS_MPI
743  // Pass on non-halo proc ID
744  face_element_pt->set_halo(Non_halo_proc_ID);
745 #endif
746 
747  // Calculate the number of nodes in the face element
748  const unsigned n_face_nodes = (NNODE_1D * (NNODE_1D + 1)) / 2;
749 
750  // Resize storage for the number of values originally stored
751  // at the face element's nodes.
752  face_element_pt->nbulk_value_resize(n_face_nodes);
753 
754  // Resize storage for the bulk node numbers corresponding to
755  // the nodes of the face
756  face_element_pt->bulk_node_number_resize(n_face_nodes);
757 
758  // Set the face index in the face element
759  face_element_pt->face_index() = face_index;
760 
761  // So the faces are
762  // 0 : s_0 fixed
763  // 1 : s_1 fixed
764  // 2 : s_2 fixed
765  // 3 : sloping face
766 
767  // Copy nodes across to the face
768  for (unsigned i = 0; i < n_face_nodes; i++)
769  {
770  // The number is just offset by one
771  unsigned bulk_number = Node_on_face[face_index][i];
772  face_element_pt->node_pt(i) = node_pt(bulk_number);
773  face_element_pt->bulk_node_number(i) = bulk_number;
774  // set the number of values originally stored at this node
775  face_element_pt->nbulk_value(i) = required_nvalue(bulk_number);
776  }
777 
778  // Now set up the node pointers and the normal vectors
779  switch (face_index)
780  {
781  //
782  //-----The face s0 is constant
783  case 0:
784 
785  // Set the pointer to the function that determines the bulk
786  // coordinates in the face element
787  face_element_pt->face_to_bulk_coordinate_fct_pt() =
788  &TElement3FaceToBulkCoordinates::face0;
789 
790  // Outer unit normal is the negative cross product of two in plane
791  // tangent vectors
792  face_element_pt->normal_sign() = -1;
793 
794  break;
795 
796  //-----The face s1 is constant
797  case 1:
798 
799  // Set the pointer to the function that determines the bulk
800  // coordinates in the face element
801  face_element_pt->face_to_bulk_coordinate_fct_pt() =
802  &TElement3FaceToBulkCoordinates::face1;
803 
804  // Outer unit normal is the positive cross product of two in plane
805  // tangent vectors
806  face_element_pt->normal_sign() = +1;
807 
808  break;
809 
810  //
811  //-----The face s2 is constant
812  case 2:
813 
814  // Set the pointer to the function that determines the bulk
815  // coordinates in the face element
816  face_element_pt->face_to_bulk_coordinate_fct_pt() =
817  &TElement3FaceToBulkCoordinates::face2;
818 
819  // Outer unit normal is the negative cross product of two in plane
820  // tangent vectors
821  face_element_pt->normal_sign() = -1;
822 
823  break;
824 
825  //-----The sloping face of the tetrahedron
826  case 3:
827 
828  // Set the pointer to the function that determines the bulk
829  // coordinates in the face element
830  face_element_pt->face_to_bulk_coordinate_fct_pt() =
831  &TElement3FaceToBulkCoordinates::face3;
832 
833  // Outer unit normal is the positive cross product of two in plane
834  // tangent vectors
835  face_element_pt->normal_sign() = +1;
836 
837  break;
838 
839 
840  default:
841 
842  std::ostringstream error_message;
843  error_message << "Face_index should only take "
844  << "the values 0, 1, 2 or 3, not " << face_index
845  << std::endl;
846 
847  throw OomphLibError(error_message.str(),
848  OOMPH_CURRENT_FUNCTION,
849  OOMPH_EXCEPTION_LOCATION);
850  } // end switch
851  }
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< 3, NNODE_1D >::Node_on_face
static const unsigned Node_on_face[4][(NNODE_1D *(NNODE_1D+1))/2]
Nodal translation scheme for use when generating face elements.
Definition: Telements.h:3045
oomph::TElement3FaceToBulkCoordinates::face1
void face1(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the face s1 = 0.
Definition: Telements.cc:690
oomph::TElement3FaceToBulkCoordinates::face2
void face2(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the face s2 = 0.
Definition: Telements.cc:698
oomph::TElement3FaceToBulkCoordinates::face0
void face0(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the face s0 = 0.
Definition: Telements.cc:682
oomph::TElement3FaceToBulkCoordinates::face3
void face3(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the sloping face.
Definition: Telements.cc:706

References oomph::FaceElement::bulk_element_pt(), oomph::FaceElement::bulk_node_number(), oomph::FaceElement::bulk_node_number_resize(), oomph::TElement3FaceToBulkCoordinates::face0(), oomph::TElement3FaceToBulkCoordinates::face1(), oomph::TElement3FaceToBulkCoordinates::face2(), oomph::TElement3FaceToBulkCoordinates::face3(), 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< 3, 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_2^2 \) d2psids(i,3) = \( \partial^2 \psi_j / \partial s_0 \partial s_1 \) d2psids(i,4) = \( \partial^2 \psi_j / \partial s_0 \partial s_2 \) d2psids(i,5) = \( \partial^2 \psi_j / \partial s_1 \partial s_2 \)

Reimplemented from oomph::FiniteElement.

3168  {
3169  TElementShape<3, NNODE_1D>::d2shape_local(s, psi, dpsids, d2psids);
3170  }
s
RealScalar s
Definition: level1_cplx_impl.h:130

References s.

◆ dshape_local()

template<unsigned NNODE_1D>
void oomph::TElement< 3, 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.

3152  {
3153  TElementShape<3, NNODE_1D>::dshape_local(s, psi, dpsids);
3154  }

References s.

◆ get_bulk_node_number()

template<unsigned NNODE_1D>
unsigned oomph::TElement< 3, 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.

3118  {
3119  return Node_on_face[face_index][i];
3120  }

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

◆ get_s_plot()

template<unsigned NNODE_1D>
void oomph::TElement< 3, 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.

3359  {
3360  if (nplot > 1)
3361  {
3362  unsigned np = 0;
3363  for (unsigned i = 0; i < nplot; ++i)
3364  {
3365  for (unsigned j = 0; j < nplot - i; ++j)
3366  {
3367  for (unsigned k = 0; k < nplot - i - j; ++k)
3368  {
3369  if (np == iplot)
3370  {
3371  {
3372  s[0] = double(j) / double(nplot - 1);
3373  s[1] = double(i) / double(nplot - 1);
3374  s[2] = double(k) / double(nplot - 1);
3375  if (use_equally_spaced_interior_sample_points)
3376  {
3377  double range = 1.0;
3378  double dx_new = range / double(nplot + 1);
3379  double range_new = double(nplot - 1) * dx_new;
3380  s[0] = 0.5 * dx_new + range_new * s[0] / range;
3381  s[1] = 0.5 * dx_new + range_new * s[1] / range;
3382  s[2] = 0.5 * dx_new + range_new * s[2] / range;
3383  }
3384  return;
3385  }
3386  }
3387  np++;
3388  }
3389  }
3390  }
3391  }
3392  else
3393  {
3394  s[0] = 1.0 / 4.0;
3395  s[1] = 1.0 / 4.0;
3396  s[2] = 1.0 / 4.0;
3397  }
3398  }
k
char char char int int * k
Definition: level2_impl.h:374
j
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References i, j, k, and s.

◆ invert_jacobian_mapping()

template<unsigned NNODE_1D>
double oomph::TElement< 3, 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 three dimensional element, so use the 3D version.

Reimplemented from oomph::FiniteElement.

3177  {
3178  return FiniteElement::invert_jacobian<3>(jacobian, inverse_jacobian);
3179  }

◆ local_coordinate_of_node()

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

Return local coordinates of node j.

Reimplemented from oomph::FiniteElement.

3196  {
3197  TElementShape<3, NNODE_1D>::local_coordinate_of_node(j, s);
3198  }

References j, and s.

◆ nnode_1d()

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

Number of nodes along each element edge.

Reimplemented from oomph::FiniteElement.

3103  {
3104  return NNODE_1D;
3105  }

◆ nplot_points()

template<unsigned NNODE_1D>
unsigned oomph::TElement< 3, 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.

3530  {
3531  unsigned res = 0;
3532  if (nplot > 1)
3533  {
3534  res = 4;
3535  for (unsigned i = 3; i <= nplot; i++)
3536  {
3537  res += (i * (i + 1) / 2);
3538  }
3539  return res;
3540  }
3541  // Otherwise we return 1(?)
3542  return 1;
3543  }
res
cout<< "Here is the matrix m:"<< endl<< m<< endl;Matrix< ptrdiff_t, 3, 1 > res
Definition: PartialRedux_count.cpp:3

References i, and res.

◆ nplot_points_paraview()

template<unsigned NNODE_1D>
unsigned oomph::TElement< 3, 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.

3203  {
3204  unsigned node_sum = 0;
3205  for (unsigned j = 1; j <= nplot; j++)
3206  {
3207  for (unsigned i = 1; i <= j; i++)
3208  {
3209  node_sum += i;
3210  }
3211  }
3212  return node_sum;
3213  }

References i, and j.

◆ nsub_elements_paraview()

template<unsigned NNODE_1D>
unsigned oomph::TElement< 3, 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.

3218  {
3219  return (nplot - 1) * (nplot - 1) * (nplot - 1);
3220  }

◆ nvertex_node()

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

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

Reimplemented from oomph::FiniteElement.

3111  {
3112  return 4;
3113  }

◆ output() [1/4]

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

C-style output.

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

Reimplemented from oomph::FiniteElement.

635  {
636  output(file_pt, NNODE_1D);
637  }
oomph::TElement< 3, NNODE_1D >::output
void output(std::ostream &output)
Output.
Definition: Telements.cc:589

References oomph::output().

◆ output() [2/4]

template<unsigned NNODE_1D>
void oomph::TElement< 3, 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<3,NNODE_1D>

Reimplemented from oomph::FiniteElement.

645  {
646  // Vector of local coordinates
647  Vector<double> s(3);
648 
649  // Find the dimension of the nodes
650  unsigned n_dim = this->nodal_dimension();
651 
652  // Tecplot header info
653  fprintf(file_pt, "%s \n", tecplot_zone_string(nplot).c_str());
654 
655  // Loop over plot points
656  unsigned num_plot_points = nplot_points(nplot);
657  for (unsigned iplot = 0; iplot < num_plot_points; iplot++)
658  {
659  // Get local coordinates of plot point
660  get_s_plot(iplot, nplot, s);
661 
662  for (unsigned i = 0; i < n_dim; i++)
663  {
664  fprintf(file_pt, "%g ", interpolated_x(s, i));
665  // outfile << interpolated_x(s,i) << " ";
666  }
667  fprintf(file_pt, "\n");
668  // outfile << std::endl;
669  }
670 
671  // Write tecplot footer (e.g. FE connectivity lists)
672  write_tecplot_zone_footer(file_pt, nplot);
673  }
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< 3, NNODE_1D >::write_tecplot_zone_footer
void write_tecplot_zone_footer(std::ostream &outfile, const unsigned &nplot) const
Definition: Telements.h:3416
oomph::TElement< 3, NNODE_1D >::tecplot_zone_string
std::string tecplot_zone_string(const unsigned &nplot) const
Definition: Telements.h:3402
oomph::TElement< 3, NNODE_1D >::nplot_points
unsigned nplot_points(const unsigned &nplot) const
Definition: Telements.h:3529
oomph::TElement< 3, 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:3354

References i, and s.

◆ output() [3/4]

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

Output at specified number of plot points.

The output function for TElement<3,NNODE_1D>

Reimplemented from oomph::FiniteElement.

601  {
602  // Vector of local coordinates
603  Vector<double> s(3);
604 
605  // Find the dimension of the nodes
606  unsigned n_dim = this->nodal_dimension();
607 
608  // Tecplot header info
609  outfile << tecplot_zone_string(nplot);
610 
611  // Loop over plot points
612  unsigned num_plot_points = nplot_points(nplot);
613  for (unsigned iplot = 0; iplot < num_plot_points; iplot++)
614  {
615  // Get local coordinates of plot point
616  get_s_plot(iplot, nplot, s);
617 
618  for (unsigned i = 0; i < n_dim; i++)
619  {
620  outfile << interpolated_x(s, i) << " ";
621  }
622  outfile << "\n";
623  }
624 
625  // Write tecplot footer (e.g. FE connectivity lists)
626  write_tecplot_zone_footer(outfile, nplot);
627  }

References i, and s.

◆ output() [4/4]

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

Output.

The output function for TElement<3,NNODE_1D>

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

Reimplemented from oomph::FiniteElement.

590  {
591  output(outfile, NNODE_1D);
592  }

References oomph::output().

◆ s_max()

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

Max. value of local coordinate.

Reimplemented from oomph::FiniteElement.

3189  {
3190  return 1.0;
3191  }

◆ s_min()

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

Min. value of local coordinate.

Reimplemented from oomph::FiniteElement.

3183  {
3184  return 0.0;
3185  }

◆ shape()

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

Calculate the geometric shape functions at local coordinate s.

Implements oomph::FiniteElement.

3143  {
3144  TElementShape<3, NNODE_1D>::shape(s, psi);
3145  }
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< 3, 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.

3403  {
3404  std::ostringstream header;
3405  unsigned nel = 0;
3406  nel = (nplot - 1) * (nplot - 1) * (nplot - 1);
3407  header << "ZONE N=" << nplot_points(nplot) << ", E=" << nel
3408  << ", F=FEPOINT, ET=TETRAHEDRON\n";
3409  return header.str();
3410  }

◆ vertex_node_pt()

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

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

Reimplemented from oomph::FiniteElement.

3124  {
3125  // Vertex nodes come first:
3126 #ifdef PARANOID
3127  if (j > 3)
3128  {
3129  std::ostringstream error_message;
3130  error_message
3131  << "Element only has four vertex nodes; called with node number " << j
3132  << std::endl;
3133  throw OomphLibError(error_message.str(),
3134  OOMPH_CURRENT_FUNCTION,
3135  OOMPH_EXCEPTION_LOCATION);
3136  }
3137 #endif
3138  return node_pt(j);
3139  }

References j, OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

◆ write_paraview_offsets()

template<unsigned NNODE_1D>
void oomph::TElement< 3, 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.

3331  {
3332  unsigned local_loop = nsub_elements_paraview(nplot);
3333  for (unsigned i = 0; i < local_loop; i++)
3334  {
3335  offset_sum += 4;
3336  file_out << offset_sum << std::endl;
3337  }
3338  }
oomph::TElement< 3, NNODE_1D >::nsub_elements_paraview
unsigned nsub_elements_paraview(const unsigned &nplot) const
Definition: Telements.h:3217

References i.

◆ write_paraview_output_offset_information()

template<unsigned NNODE_1D>
void oomph::TElement< 3, 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.

3228  {
3229  // Output node lists for sub elements for Paraview (node index
3230  // must start at 0, fixed with magical counter-1)
3231 
3232  unsigned paraview_fix = counter - 1;
3233  unsigned nod_count = 1;
3234 
3235  for (unsigned i = 0; i < nplot; i++)
3236  {
3237  for (unsigned j = 0; j < nplot - i; j++)
3238  {
3239  for (unsigned k = 0; k < nplot - i - j; k++)
3240  {
3241  if (k < nplot - i - j - 1)
3242  {
3243  file_out << nod_count + paraview_fix << " "
3244  << nod_count + 1 + paraview_fix << " "
3245  << nod_count + nplot - i - j + paraview_fix << " "
3246  << nod_count + nplot - i - j +
3247  ((nplot - 1 - i) * (nplot - i) / 2) + paraview_fix
3248  << std::endl;
3249  if (k < nplot - i - j - 2)
3250  {
3251  file_out << nod_count + 1 + paraview_fix << " "
3252  << nod_count + nplot - i - j + paraview_fix << " "
3253  << nod_count + nplot - i - j +
3254  ((nplot - 1 - i) * (nplot - i) / 2) + paraview_fix
3255  << " "
3256  << nod_count + 2 * (nplot - i - j) - 1 +
3257  ((nplot - 1 - i) * (nplot - i) / 2) + paraview_fix
3258  << std::endl;
3259  file_out << nod_count + 1 + paraview_fix << " "
3260  << nod_count + nplot - i - j + paraview_fix << " "
3261  << nod_count + nplot - i - j + 1 + paraview_fix << " "
3262  << nod_count + 2 * (nplot - i - j) - 1 +
3263  ((nplot - 1 - i) * (nplot - i) / 2) + paraview_fix
3264  << std::endl;
3265  file_out << nod_count + 1 + paraview_fix << " "
3266  << nod_count + nplot - i - j +
3267  ((nplot - 1 - i) * (nplot - i) / 2) + paraview_fix
3268  << " "
3269  << nod_count + nplot - i - j +
3270  ((nplot - 1 - i) * (nplot - i) / 2) + 1 +
3271  paraview_fix
3272  << " "
3273  << nod_count + 2 * (nplot - i - j) - 1 +
3274  ((nplot - 1 - i) * (nplot - i) / 2) + paraview_fix
3275  << std::endl;
3276  file_out << nod_count + 1 + paraview_fix << " "
3277  << nod_count + nplot - i - j + 1 + paraview_fix << " "
3278  << nod_count + nplot - i - j +
3279  ((nplot - 1 - i) * (nplot - i) / 2) + 1 +
3280  paraview_fix
3281  << " "
3282  << nod_count + 2 * (nplot - i - j) - 1 +
3283  ((nplot - 1 - i) * (nplot - i) / 2) + paraview_fix
3284  << std::endl;
3285  }
3286  if (k > 1)
3287  {
3288  file_out << nod_count + nplot - i - j - 1 + paraview_fix << " "
3289  << nod_count + nplot - i - j +
3290  ((nplot - 1 - i) * (nplot - i) / 2) - 1 +
3291  paraview_fix
3292  << " "
3293  << nod_count + 2 * (nplot - i - j - 1) +
3294  ((nplot - 1 - i) * (nplot - i) / 2) - 1 +
3295  paraview_fix
3296  << " "
3297  << nod_count + 2 * (nplot - i - j - 1) +
3298  ((nplot - 1 - i) * (nplot - i) / 2) + paraview_fix
3299  << std::endl;
3300  }
3301  } // end if(k<nplot-i-j-1)
3302  ++nod_count;
3303  }
3304  }
3305  }
3306 
3307  // increment the counter to keep track of global connectivity
3308  counter += nplot_points_paraview(nplot);
3309 
3310  } // end of write Paraview_element...
oomph::TElement< 3, NNODE_1D >::nplot_points_paraview
unsigned nplot_points_paraview(const unsigned &nplot) const
Definition: Telements.h:3202

References i, j, and k.

◆ write_paraview_type()

template<unsigned NNODE_1D>
void oomph::TElement< 3, 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.

3317  {
3318  unsigned local_loop = nsub_elements_paraview(nplot);
3319  for (unsigned i = 0; i < local_loop; i++)
3320  {
3321  file_out << "10" << std::endl;
3322  }
3323  }

References i.

◆ write_tecplot_zone_footer() [1/2]

template<unsigned NNODE_1D>
void oomph::TElement< 3, 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.

3495  {
3496  // Output node lists for sub elements for Tecplot (node index
3497  // must start at 1)
3498  unsigned nod_count = 1;
3499  for (unsigned i = 0; i < nplot; i++)
3500  {
3501  for (unsigned j = 0; j < nplot - i; j++)
3502  {
3503  for (unsigned k = 0; k < nplot - i - j; k++)
3504  {
3505  if (j < nplot - i - 1)
3506  {
3507  fprintf(file_pt,
3508  "%i %i %i \n",
3509  nod_count,
3510  nod_count + 1,
3511  nod_count + nplot - i);
3512  if (j < nplot - i - 2)
3513  {
3514  fprintf(file_pt,
3515  "%i %i %i \n",
3516  nod_count + 1,
3517  nod_count + nplot - i + 1,
3518  nod_count + nplot - i);
3519  }
3520  }
3521  ++nod_count;
3522  }
3523  }
3524  }
3525  }

References i, j, and k.

◆ write_tecplot_zone_footer() [2/2]

template<unsigned NNODE_1D>
void oomph::TElement< 3, 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.

3418  {
3419  // Output node lists for sub elements for Tecplot (node index
3420  // must start at 1)
3421  unsigned nod_count = 1;
3422  for (unsigned i = 0; i < nplot; i++)
3423  {
3424  for (unsigned j = 0; j < nplot - i; j++)
3425  {
3426  for (unsigned k = 0; k < nplot - i - j; k++)
3427  {
3428  if (k < nplot - i - j - 1)
3429  {
3430  outfile << nod_count << " " << nod_count + 1 << " "
3431  << nod_count + nplot - i - j << " "
3432  << nod_count + nplot - i - j +
3433  ((nplot - 1 - i) * (nplot - i) / 2)
3434  << std::endl;
3435  if (k < nplot - i - j - 2)
3436  {
3437  outfile << nod_count + 1 << " " << nod_count + nplot - i - j
3438  << " "
3439  << nod_count + nplot - i - j +
3440  ((nplot - 1 - i) * (nplot - i) / 2)
3441  << " "
3442  << nod_count + 2 * (nplot - i - j) - 1 +
3443  ((nplot - 1 - i) * (nplot - i) / 2)
3444  << std::endl;
3445  outfile << nod_count + 1 << " " << nod_count + nplot - i - j
3446  << " " << nod_count + nplot - i - j + 1 << " "
3447  << nod_count + 2 * (nplot - i - j) - 1 +
3448  ((nplot - 1 - i) * (nplot - i) / 2)
3449  << std::endl;
3450  outfile << nod_count + 1 << " "
3451  << nod_count + nplot - i - j +
3452  ((nplot - 1 - i) * (nplot - i) / 2)
3453  << " "
3454  << nod_count + nplot - i - j +
3455  ((nplot - 1 - i) * (nplot - i) / 2) + 1
3456  << " "
3457  << nod_count + 2 * (nplot - i - j) - 1 +
3458  ((nplot - 1 - i) * (nplot - i) / 2)
3459  << std::endl;
3460  outfile << nod_count + 1 << " " << nod_count + nplot - i - j + 1
3461  << " "
3462  << nod_count + nplot - i - j +
3463  ((nplot - 1 - i) * (nplot - i) / 2) + 1
3464  << " "
3465  << nod_count + 2 * (nplot - i - j) - 1 +
3466  ((nplot - 1 - i) * (nplot - i) / 2)
3467  << std::endl;
3468  }
3469  if (k > 1)
3470  {
3471  outfile << nod_count + nplot - i - j - 1 << " "
3472  << nod_count + nplot - i - j +
3473  ((nplot - 1 - i) * (nplot - i) / 2) - 1
3474  << " "
3475  << nod_count + 2 * (nplot - i - j - 1) +
3476  ((nplot - 1 - i) * (nplot - i) / 2) - 1
3477  << " "
3478  << nod_count + 2 * (nplot - i - j - 1) +
3479  ((nplot - 1 - i) * (nplot - i) / 2)
3480  << std::endl;
3481  }
3482  } // end if(k<nplot-i-j-1)
3483  ++nod_count;
3484  }
3485  }
3486  }
3487  } // end of write tecplot...

References i, j, and k.

Member Data Documentation

◆ Default_integration_scheme

template<unsigned NNODE_1D>
TGauss< 3, NNODE_1D > oomph::TElement< 3, 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< 3, NNODE_1D >::Node_on_face[4][(NNODE_1D *(NNODE_1D+1))/2]
staticprivate

Nodal translation scheme for use when generating face elements.

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