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

#include <Telements.h>

+ Inheritance diagram for oomph::TElement< 1, 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
 
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 &i, 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
 
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 write_tecplot_zone_footer (std::ostream &outfile, const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (FILE *file_pt, const unsigned &nplot) const
 
virtual 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 unsigned get_bulk_node_number (const int &face_index, const unsigned &i) const
 
virtual int face_outer_unit_normal_sign (const int &face_index) const
 Get the sign of the outer unit normal on the face given by face_index. More...
 
virtual unsigned nnode_on_face () const
 
void face_node_number_error_check (const unsigned &i) const
 Range check for face node numbers. More...
 
virtual CoordinateMappingFctPt face_to_bulk_coordinate_fct_pt (const int &face_index) const
 
virtual BulkCoordinateDerivativesFctPt bulk_coordinate_derivatives_fct_pt (const int &face_index) const
 
- Public Member Functions inherited from oomph::GeneralisedElement
 GeneralisedElement ()
 Constructor: Initialise all pointers and all values to zero. More...
 
virtual ~GeneralisedElement ()
 Virtual destructor to clean up any memory allocated by the object. More...
 
 GeneralisedElement (const GeneralisedElement &)=delete
 Broken copy constructor. More...
 
void operator= (const GeneralisedElement &)=delete
 Broken assignment operator. More...
 
Data *& internal_data_pt (const unsigned &i)
 Return a pointer to i-th internal data object. More...
 
Data *const & internal_data_pt (const unsigned &i) const
 Return a pointer to i-th internal data object (const version) More...
 
Data *& external_data_pt (const unsigned &i)
 Return a pointer to i-th external data object. More...
 
Data *const & external_data_pt (const unsigned &i) const
 Return a pointer to i-th external data object (const version) More...
 
unsigned long eqn_number (const unsigned &ieqn_local) const
 
int local_eqn_number (const unsigned long &ieqn_global) const
 
unsigned add_external_data (Data *const &data_pt, const bool &fd=true)
 
bool external_data_fd (const unsigned &i) const
 
void exclude_external_data_fd (const unsigned &i)
 
void include_external_data_fd (const unsigned &i)
 
void flush_external_data ()
 Flush all external data. More...
 
void flush_external_data (Data *const &data_pt)
 Flush the object addressed by data_pt from the external data array. More...
 
unsigned ninternal_data () const
 Return the number of internal data objects. More...
 
unsigned nexternal_data () const
 Return the number of external data objects. More...
 
unsigned ndof () const
 Return the number of equations/dofs in the element. More...
 
void dof_vector (const unsigned &t, Vector< double > &dof)
 Return the vector of dof values at time level t. More...
 
void dof_pt_vector (Vector< double * > &dof_pt)
 Return the vector of pointers to dof values. More...
 
void set_internal_data_time_stepper (const unsigned &i, TimeStepper *const &time_stepper_pt, const bool &preserve_existing_data)
 
void assign_internal_eqn_numbers (unsigned long &global_number, Vector< double * > &Dof_pt)
 
void describe_dofs (std::ostream &out, const std::string &current_string) const
 
void add_internal_value_pt_to_map (std::map< unsigned, double * > &map_of_value_pt)
 
virtual void assign_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void complete_setup_of_dependencies ()
 
virtual void get_residuals (Vector< double > &residuals)
 
virtual void get_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
virtual void get_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &mass_matrix)
 
virtual void get_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
virtual void get_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
virtual void get_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
virtual void get_djacobian_and_dmass_matrix_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam)
 
virtual void get_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
virtual void get_inner_products (Vector< std::pair< unsigned, unsigned >> const &history_index, Vector< double > &inner_product)
 
virtual void get_inner_product_vectors (Vector< unsigned > const &history_index, Vector< Vector< double >> &inner_product_vector)
 
virtual void compute_norm (Vector< double > &norm)
 
virtual void compute_norm (double &norm)
 
virtual unsigned ndof_types () const
 
virtual void get_dof_numbers_for_unknowns (std::list< std::pair< unsigned long, unsigned >> &dof_lookup_list) const
 
- Public Member Functions inherited from oomph::GeomObject
 GeomObject ()
 Default constructor. More...
 
 GeomObject (const unsigned &ndim)
 
 GeomObject (const unsigned &nlagrangian, const unsigned &ndim)
 
 GeomObject (const unsigned &nlagrangian, const unsigned &ndim, TimeStepper *time_stepper_pt)
 
 GeomObject (const GeomObject &dummy)=delete
 Broken copy constructor. More...
 
void operator= (const GeomObject &)=delete
 Broken assignment operator. More...
 
virtual ~GeomObject ()
 (Empty) destructor More...
 
unsigned nlagrangian () const
 Access function to # of Lagrangian coordinates. More...
 
unsigned ndim () const
 Access function to # of Eulerian coordinates. More...
 
void set_nlagrangian_and_ndim (const unsigned &n_lagrangian, const unsigned &n_dim)
 Set # of Lagrangian and Eulerian coordinates. More...
 
TimeStepper *& time_stepper_pt ()
 
TimeStepper * time_stepper_pt () const
 
virtual void position (const double &t, const Vector< double > &zeta, Vector< double > &r) const
 
virtual void dposition (const Vector< double > &zeta, DenseMatrix< double > &drdzeta) const
 
virtual void d2position (const Vector< double > &zeta, RankThreeTensor< double > &ddrdzeta) const
 
virtual void d2position (const Vector< double > &zeta, Vector< double > &r, DenseMatrix< double > &drdzeta, RankThreeTensor< double > &ddrdzeta) const
 

Static Private Attributes

static TGauss< 1, NNODE_1D > Default_integration_scheme
 Assign the static integral. More...
 

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< 1, NNODE_1D >

General TElement class specialised to one spatial dimensions Ordering of nodes is 0 at local coordinate s[0] = 0, 1 at local coordinate s[0] = 1 and then filling in the intermediate values from s[0]=0 to 1.

Constructor & Destructor Documentation

◆ TElement() [1/2]

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

Constructor.

1233  {
1234  // Number of nodes
1235  switch (NNODE_1D)
1236  {
1237  case 2:
1238  case 3:
1239  case 4:
1240  break;
1241 
1242  default:
1243  std::string error_message =
1244  "One-dimensional TElements are currently only implemented for\n";
1245  error_message += "three and six nodes, i.e. NNODE_1D=2 , 3 , 4\n";
1246 
1247  throw OomphLibError(
1248  error_message, OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1249  }
1250 
1251  // Set the number of nodes
1252  unsigned n_node = NNODE_1D;
1253  this->set_n_node(n_node);
1254 
1255  // Set the elemental and nodal dimension
1256  set_dimension(1);
1257 
1258  // Assign default (full) spatial integration scheme
1259  set_integration_scheme(&Default_integration_scheme);
1260  }
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< 1, NNODE_1D >::Default_integration_scheme
static TGauss< 1, NNODE_1D > Default_integration_scheme
Assign the static integral.
Definition: Telements.h:1228
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< 1, NNODE_1D >::TElement ( const TElement< 1, NNODE_1D > &  )
delete

Broken copy constructor.

◆ ~TElement()

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

Broken assignment operator.

Destructor

1271 {}

Member Function Documentation

◆ build_face_element()

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

Build the lower-dimensional FaceElement (an element of type PointElement). The face index takes two values corresponding to the two possible faces: -1 (Left) s[0] = -1.0 +1 (Right) s[0] = 1.0

Function to setup geometrical information for lower-dimensional FaceElements (which are Point Elements).

Reimplemented from oomph::FiniteElement.

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

179  {
180  // Overload the nodal dimension
181  face_element_pt->set_nodal_dimension(nodal_dimension());
182 
183  // Set the pointer to the "bulk" element
184  face_element_pt->bulk_element_pt() = this;
185 
186 #ifdef OOMPH_HAS_MPI
187  // Pass on non-halo proc ID
188  face_element_pt->set_halo(Non_halo_proc_ID);
189 #endif
190 
191  // Resize the storage for the original number of values at the (one and
192  // only) node of the face element.
193  face_element_pt->nbulk_value_resize(1);
194 
195  // Resize the storage for the bulk node number corresponding to the (one
196  // and only) node of the face element
197  face_element_pt->bulk_node_number_resize(1);
198 
199  // Set the face index in the face element
200  face_element_pt->face_index() = face_index;
201 
202  // Now set up the node pointer
203  // The convention is that the "normal", should always point
204  // out of the element
205  switch (face_index)
206  {
207  // Bottom, normal sign is negative (coordinate points into element)
208  case (-1):
209  face_element_pt->node_pt(0) = node_pt(0);
210  face_element_pt->bulk_node_number(0) = 0;
211  face_element_pt->normal_sign() = -1;
212 
213  // Set the pointer to the function that determines the bulk coordinates
214  // in the face element
215  face_element_pt->face_to_bulk_coordinate_fct_pt() =
216  &TElement1FaceToBulkCoordinates::face0;
217 
218  // Set the pointer to the function that determines the mapping of
219  // derivatives
220  face_element_pt->bulk_coordinate_derivatives_fct_pt() =
221  &TElement1BulkCoordinateDerivatives::faces0;
222 
223  // Set the number of values stored when the node is part of the "bulk"
224  // element. The required_nvalue() must be used, rather than nvalue(),
225  // because otherwise nodes on boundaries will be resized multiple
226  // times. If you want any other behaviour, you MUST set nbulk_value()
227  // manually after construction of your specific face element.
228  face_element_pt->nbulk_value(0) = required_nvalue(0);
229  break;
230 
231  // Top, normal sign is positive (coordinate points out of element)
232  case (1):
233  face_element_pt->node_pt(0) = node_pt(NNODE_1D - 1);
234  face_element_pt->bulk_node_number(0) = NNODE_1D - 1;
235  face_element_pt->normal_sign() = +1;
236 
237  // Set the pointer to the function that determines the bulk coordinates
238  // in the face element
239  face_element_pt->face_to_bulk_coordinate_fct_pt() =
240  &TElement1FaceToBulkCoordinates::face1;
241 
242  // Set the pointer to the function that determines the mapping of
243  // derivatives
244  face_element_pt->bulk_coordinate_derivatives_fct_pt() =
245  &TElement1BulkCoordinateDerivatives::faces0;
246 
247 
248  // Set the number of values stored when the node is part of the "bulk"
249  // element.
250  face_element_pt->nbulk_value(0) = required_nvalue(NNODE_1D - 1);
251  break;
252 
253  // All other cases throw an error
254  default:
255  std::ostringstream error_message;
256  error_message << "Face_index should only take "
257  << "the values +/-1, not " << face_index << std::endl;
258 
259  throw OomphLibError(error_message.str(),
260  OOMPH_CURRENT_FUNCTION,
261  OOMPH_EXCEPTION_LOCATION);
262  }
263  }
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::TElement1BulkCoordinateDerivatives::faces0
void faces0(const Vector< double > &s, DenseMatrix< double > &dsbulk_dsface, unsigned &interior_direction)
Function for both faces – the bulk coordinate is fixed on both.
Definition: Telements.cc:160
oomph::TElement1FaceToBulkCoordinates::face0
void face0(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the face s0 = 0.0.
Definition: Telements.cc:139
oomph::TElement1FaceToBulkCoordinates::face1
void face1(const Vector< double > &s, Vector< double > &s_bulk)
The translation scheme for the face s0 = 1.0.
Definition: Telements.cc:145

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::TElement1FaceToBulkCoordinates::face0(), oomph::TElement1FaceToBulkCoordinates::face1(), oomph::FaceElement::face_index(), oomph::FaceElement::face_to_bulk_coordinate_fct_pt(), oomph::TElement1BulkCoordinateDerivatives::faces0(), oomph::FaceElement::nbulk_value(), oomph::FaceElement::nbulk_value_resize(), 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< 1, 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.

1338  {
1339  TElementShape<1, NNODE_1D>::d2shape_local(s, psi, dpsids, d2psids);
1340  }
s
RealScalar s
Definition: level1_cplx_impl.h:130

References s.

◆ dshape_local()

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

1325  {
1326  TElementShape<1, NNODE_1D>::dshape_local(s, psi, dpsids);
1327  }

References s.

◆ get_s_plot()

template<unsigned NNODE_1D>
void oomph::TElement< 1, NNODE_1D >::get_s_plot ( const unsigned &  i,
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.

1451  {
1452  if (nplot > 1)
1453  {
1454  s[0] = double(i) / double(nplot - 1);
1455 
1456  if (use_equally_spaced_interior_sample_points)
1457  {
1458  double range = 1.0;
1459  double dx_new = range / double(nplot);
1460  double range_new = double(nplot - 1) * dx_new;
1461  s[0] = 0.5 * dx_new + range_new * s[0] / range;
1462  }
1463  }
1464  else
1465  {
1466  s[0] = 0.5;
1467  }
1468  }
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9

References i, and s.

◆ invert_jacobian_mapping()

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

Reimplemented from oomph::FiniteElement.

1347  {
1348  return FiniteElement::invert_jacobian<1>(jacobian, inverse_jacobian);
1349  }

◆ local_coordinate_of_node()

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

Return local coordinates of node j.

Reimplemented from oomph::FiniteElement.

1366  {
1367  TElementShape<1, NNODE_1D>::local_coordinate_of_node(j, s);
1368  }
j
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References j, and s.

◆ nnode_1d()

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

Number of nodes along each element edge.

Reimplemented from oomph::FiniteElement.

1275  {
1276  return NNODE_1D;
1277  }

◆ nplot_points()

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

1482  {
1483  return nplot;
1484  }

◆ nplot_points_paraview()

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

1373  {
1374  return nplot;
1375  }

◆ nsub_elements_paraview()

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

1380  {
1381  return (nplot - 1);
1382  }

◆ nvertex_node()

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

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

Reimplemented from oomph::FiniteElement.

1283  {
1284  return 2;
1285  }

◆ output() [1/4]

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

C-style output.

C style output function for general 1D TElements.

Reimplemented from oomph::FiniteElement.

94  {
95  output(file_pt, NNODE_1D);
96  }
oomph::TElement< 1, NNODE_1D >::output
void output(std::ostream &output)
Output.
Definition: Telements.cc:53

References oomph::output().

◆ output() [2/4]

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

C_style output at n_plot points.

C style output function for n_plot points in each coordinate direction.

Reimplemented from oomph::FiniteElement.

103  {
104  // Local variables
105  Vector<double> s(1);
106 
107  // Tecplot header info
108  // outfile << "ZONE I=" << n_plot << std::endl;
109  fprintf(file_pt, "ZONE I=%i\n", n_plot);
110 
111  // Find the dimension of the first node
112  unsigned n_dim = this->nodal_dimension();
113 
114  // Loop over plot points
115  for (unsigned l = 0; l < n_plot; l++)
116  {
117  s[0] = l * 1.0 / (n_plot - 1);
118 
119  // Output the coordinates
120  for (unsigned i = 0; i < n_dim; i++)
121  {
122  // outfile << interpolated_x(s,i) << " " ;
123  fprintf(file_pt, "%g ", interpolated_x(s, i));
124  }
125  // outfile << std::endl;
126  fprintf(file_pt, "\n");
127  }
128  // outfile << std::endl;
129  fprintf(file_pt, "\n");
130  }
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

References i, and s.

◆ output() [3/4]

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

Output at specified number of plot points.

The output function for n_plot points in each coordinate direction.

Reimplemented from oomph::FiniteElement.

64  {
65  // Local variables
66  Vector<double> s(1);
67 
68  // Tecplot header info
69  outfile << "ZONE I=" << n_plot << std::endl;
70 
71  // Find the dimension of the nodes
72  unsigned n_dim = this->nodal_dimension();
73 
74  // Loop over plot points
75  for (unsigned l = 0; l < n_plot; l++)
76  {
77  s[0] = l * 1.0 / (n_plot - 1);
78  // Output the coordinates
79  for (unsigned i = 0; i < n_dim; i++)
80  {
81  outfile << interpolated_x(s, i) << " ";
82  }
83  outfile << std::endl;
84  }
85  outfile << std::endl;
86  }

References i, and s.

◆ output() [4/4]

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

Output.

The output function for general 1D TElements.

Reimplemented from oomph::FiniteElement.

54  {
55  output(outfile, NNODE_1D);
56  }

References oomph::output().

◆ s_max()

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

Max. value of local coordinate.

Reimplemented from oomph::FiniteElement.

1359  {
1360  return 1.0;
1361  }

◆ s_min()

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

Min. value of local coordinate.

Reimplemented from oomph::FiniteElement.

1353  {
1354  return 0.0;
1355  }

◆ shape()

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

Calculate the geometric shape functions at local coordinate s.

Implements oomph::FiniteElement.

1316  {
1317  TElementShape<1, NNODE_1D>::shape(s, psi);
1318  }
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< 1, 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.

1473  {
1474  std::ostringstream header;
1475  header << "ZONE I=" << nplot << "\n";
1476  return header.str();
1477  }

◆ vertex_node_pt()

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

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

Reimplemented from oomph::FiniteElement.

1289  {
1290  switch (j)
1291  {
1292  case 0:
1293 
1294  return node_pt(0);
1295  break;
1296 
1297  case 1:
1298 
1299  return node_pt(NNODE_1D - 1);
1300  break;
1301 
1302  default:
1303 
1304  std::ostringstream error_message;
1305  error_message
1306  << "Element only has two vertex nodes; called with node number "
1307  << j << std::endl;
1308  throw OomphLibError(error_message.str(),
1309  OOMPH_CURRENT_FUNCTION,
1310  OOMPH_EXCEPTION_LOCATION);
1311  }
1312  }

References j, OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

◆ write_paraview_offsets()

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

1421  {
1422  // Loop over all local elements and add its offset to the overall
1423  // offset_sum
1424  unsigned local_loop = nsub_elements_paraview(nplot);
1425  for (unsigned i = 0; i < local_loop; i++)
1426  {
1427  offset_sum += 2;
1428  file_out << offset_sum << std::endl;
1429  }
1430  }
oomph::TElement< 1, NNODE_1D >::nsub_elements_paraview
unsigned nsub_elements_paraview(const unsigned &nplot) const
Definition: Telements.h:1379

References i.

◆ write_paraview_output_offset_information()

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

1390  {
1391  // Number of local elements we want to plot over
1392  unsigned plot = nsub_elements_paraview(nplot);
1393 
1394  // loops over the i-th local element in parent element
1395  for (unsigned i = 0; i < plot; i++)
1396  {
1397  file_out << i + counter << " " << i + 1 + counter << std::endl;
1398  }
1399  counter += nplot_points_paraview(nplot);
1400  }
oomph::TElement< 1, NNODE_1D >::nplot_points_paraview
unsigned nplot_points_paraview(const unsigned &nplot) const
Definition: Telements.h:1372
PlanarWave::plot
void plot()
Plot.
Definition: sphere_scattering.cc:180

References i, and PlanarWave::plot().

◆ write_paraview_type()

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

1407  {
1408  unsigned local_loop = nsub_elements_paraview(nplot);
1409  for (unsigned i = 0; i < local_loop; i++)
1410  {
1411  file_out << "3" << std::endl;
1412  }
1413  }

References i.

Member Data Documentation

◆ Default_integration_scheme

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

Assign the static integral.

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

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