Public Member Functions | List of all members
oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT > Class Template Reference

Taylor Hood upgraded to become projectable. More...

#include <discontinuous_galerkin_equal_order_pressure_spacetime_navier_stokes_elements.h>

+ Inheritance diagram for oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >:

Public Member Functions

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

Additional Inherited Members

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

Detailed Description

template<class TAYLOR_HOOD_ELEMENT>
class oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >

Taylor Hood upgraded to become projectable.

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

Constructor & Destructor Documentation

◆ ProjectableTaylorHoodSpaceTimeElement() [1/3]

template<class TAYLOR_HOOD_ELEMENT >
oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::ProjectableTaylorHoodSpaceTimeElement ( )
inline

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

2934 {}

◆ ProjectableTaylorHoodSpaceTimeElement() [2/3]

template<class TAYLOR_HOOD_ELEMENT >
oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::ProjectableTaylorHoodSpaceTimeElement ( )
inline

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

2951 {}

◆ ProjectableTaylorHoodSpaceTimeElement() [3/3]

template<class TAYLOR_HOOD_ELEMENT >
oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::ProjectableTaylorHoodSpaceTimeElement ( )
inline

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

2634 {}

Member Function Documentation

◆ data_values_of_field() [1/3]

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

Specify the values associated with field fld. The information is returned in a vector of pairs which comprise the Data object and the value within it, that correspond to field fld. In the underlying Taylor Hood elements the fld-th velocities are stored at the fld-th value of the nodes; the pressures (the dim-th field) are the dim-th values at the vertex nodes etc.

Implements oomph::ProjectableElementBase.

2944  {
2945  // Create the vector
2946  Vector<std::pair<Data*, unsigned>> data_values;
2947 
2948  // If we're dealing with the velocity dofs
2949  if (fld < this->dim() - 1)
2950  {
2951  // How many nodes in the element?
2952  unsigned nnod = this->nnode();
2953 
2954  // Loop over all nodes
2955  for (unsigned j = 0; j < nnod; j++)
2956  {
2957  // Add the data value associated with the velocity components
2958  data_values.push_back(std::make_pair(this->node_pt(j), fld));
2959  }
2960  }
2961  // If we're dealing with the pressure dof
2962  else
2963  {
2964  // How many pressure dofs are there?
2965  // DRAIG: Shouldn't there be more?
2966  unsigned Pconv_size = this->dim();
2967 
2968  // Loop over all vertex nodes
2969  for (unsigned j = 0; j < Pconv_size; j++)
2970  {
2971  // Get the vertex index associated with the j-th pressure dof
2972  unsigned vertex_index = this->Pconv[j];
2973 
2974  // Add the data value associated with the pressure components
2975  data_values.push_back(
2976  std::make_pair(this->node_pt(vertex_index), fld));
2977  }
2978  } // if (fld<this->dim())
2979 
2980  // Return the vector
2981  return data_values;
2982  } // End of data_values_of_field
oomph::FiniteElement::node_pt
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175
oomph::FiniteElement::dim
unsigned dim() const
Definition: elements.h:2611
oomph::FiniteElement::nnode
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210
j
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References j.

◆ data_values_of_field() [2/3]

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

Specify the values associated with field fld. The information is returned in a vector of pairs which comprise the Data object and the value within it, that correspond to field fld. In the underlying Taylor Hood elements the fld-th velocities are stored at the fld-th value of the nodes; the pressures (the dim-th field) are the dim-th values at the vertex nodes etc.

Implements oomph::ProjectableElementBase.

2961  {
2962  // Create the vector
2963  Vector<std::pair<Data*, unsigned>> data_values;
2964 
2965  // If we're dealing with the velocity dofs
2966  if (fld < this->dim() - 1)
2967  {
2968  // How many nodes in the element?
2969  unsigned nnod = this->nnode();
2970 
2971  // Loop over all nodes
2972  for (unsigned j = 0; j < nnod; j++)
2973  {
2974  // Add the data value associated with the velocity components
2975  data_values.push_back(std::make_pair(this->node_pt(j), fld));
2976  }
2977  }
2978  // If we're dealing with the pressure dof
2979  else
2980  {
2981  // How many pressure dofs are there?
2982  // DRAIG: Shouldn't there be more?
2983  unsigned Pconv_size = this->dim();
2984 
2985  // Loop over all vertex nodes
2986  for (unsigned j = 0; j < Pconv_size; j++)
2987  {
2988  // Get the vertex index associated with the j-th pressure dof
2989  unsigned vertex_index = this->Pconv[j];
2990 
2991  // Add the data value associated with the pressure components
2992  data_values.push_back(
2993  std::make_pair(this->node_pt(vertex_index), fld));
2994  }
2995  } // if (fld<this->dim())
2996 
2997  // Return the vector
2998  return data_values;
2999  } // End of data_values_of_field

References j.

◆ data_values_of_field() [3/3]

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

Specify the values associated with field fld. The information is returned in a vector of pairs which comprise the Data object and the value within it, that correspond to field fld. In the underlying Taylor Hood elements the fld-th velocities are stored at the fld-th value of the nodes; the pressures (the dim-th field) are the dim-th values at the vertex nodes etc.

Implements oomph::ProjectableElementBase.

2643  {
2644  // Create the vector
2645  Vector<std::pair<Data*, unsigned>> data_values;
2646 
2647  // If we're dealing with the velocity dofs
2648  if (fld < this->dim() - 1)
2649  {
2650  // How many nodes in the element?
2651  unsigned nnod = this->nnode();
2652 
2653  // Loop over all nodes
2654  for (unsigned j = 0; j < nnod; j++)
2655  {
2656  // Add the data value associated with the velocity components
2657  data_values.push_back(std::make_pair(this->node_pt(j), fld));
2658  }
2659  }
2660  // If we're dealing with the pressure dof
2661  else
2662  {
2663  // How many pressure dofs are there?
2664  // DRAIG: Shouldn't there be more?
2665  unsigned Pconv_size = this->dim();
2666 
2667  // Loop over all vertex nodes
2668  for (unsigned j = 0; j < Pconv_size; j++)
2669  {
2670  // Get the vertex index associated with the j-th pressure dof
2671  unsigned vertex_index = this->Pconv[j];
2672 
2673  // Add the data value associated with the pressure components
2674  data_values.push_back(
2675  std::make_pair(this->node_pt(vertex_index), fld));
2676  }
2677  } // if (fld<this->dim())
2678 
2679  // Return the vector
2680  return data_values;
2681  } // End of data_values_of_field

References j.

◆ get_field() [1/3]

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

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

Implements oomph::ProjectableElementBase.

3087  {
3088  unsigned n_dim = this->dim();
3089  unsigned n_node = this->nnode();
3090 
3091  // If fld=n_dim, we deal with the pressure
3092  if (fld == n_dim)
3093  {
3094  return this->interpolated_p_nst(t, s);
3095  }
3096  // Velocity
3097  else
3098  {
3099  // Find the index at which the variable is stored
3100  unsigned u_nodal_index = this->u_index_nst(fld);
3101 
3102  // Local shape function
3103  Shape psi(n_node);
3104 
3105  // Find values of shape function
3106  this->shape(s, psi);
3107 
3108  // Initialise value of u
3109  double interpolated_u = 0.0;
3110 
3111  // Sum over the local nodes at that time
3112  for (unsigned l = 0; l < n_node; l++)
3113  {
3114  interpolated_u += this->nodal_value(t, l, u_nodal_index) * psi[l];
3115  }
3116  return interpolated_u;
3117  }
3118  }
oomph::FiniteElement::nodal_value
double nodal_value(const unsigned &n, const unsigned &i) const
Definition: elements.h:2593
oomph::FiniteElement::shape
virtual void shape(const Vector< double > &s, Shape &psi) const =0
s
RealScalar s
Definition: level1_cplx_impl.h:130
plotPSD.t
t
Definition: plotPSD.py:36

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

◆ get_field() [2/3]

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

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

Implements oomph::ProjectableElementBase.

3104  {
3105  unsigned n_dim = this->dim();
3106  unsigned n_node = this->nnode();
3107 
3108  // If fld=n_dim, we deal with the pressure
3109  if (fld == n_dim)
3110  {
3111  return this->interpolated_p_nst(t, s);
3112  }
3113  // Velocity
3114  else
3115  {
3116  // Find the index at which the variable is stored
3117  unsigned u_nodal_index = this->u_index_nst(fld);
3118 
3119  // Local shape function
3120  Shape psi(n_node);
3121 
3122  // Find values of shape function
3123  this->shape(s, psi);
3124 
3125  // Initialise value of u
3126  double interpolated_u = 0.0;
3127 
3128  // Sum over the local nodes at that time
3129  for (unsigned l = 0; l < n_node; l++)
3130  {
3131  interpolated_u += this->nodal_value(t, l, u_nodal_index) * psi[l];
3132  }
3133  return interpolated_u;
3134  }
3135  }

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

◆ get_field() [3/3]

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

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

Implements oomph::ProjectableElementBase.

2786  {
2787  unsigned n_dim = this->dim();
2788  unsigned n_node = this->nnode();
2789 
2790  // If fld=n_dim, we deal with the pressure
2791  if (fld == n_dim)
2792  {
2793  return this->interpolated_p_nst(t, s);
2794  }
2795  // Velocity
2796  else
2797  {
2798  // Find the index at which the variable is stored
2799  unsigned u_nodal_index = this->u_index_nst(fld);
2800 
2801  // Local shape function
2802  Shape psi(n_node);
2803 
2804  // Find values of shape function
2805  this->shape(s, psi);
2806 
2807  // Initialise value of u
2808  double interpolated_u = 0.0;
2809 
2810  // Sum over the local nodes at that time
2811  for (unsigned l = 0; l < n_node; l++)
2812  {
2813  interpolated_u += this->nodal_value(t, l, u_nodal_index) * psi[l];
2814  }
2815  return interpolated_u;
2816  }
2817  }

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

◆ jacobian_and_shape_of_field() [1/3]

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

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

Implements oomph::ProjectableElementBase.

3028  {
3029  // How many dimensions in this element?
3030  unsigned n_dim = this->dim();
3031 
3032  // How many nodes in this element?
3033  unsigned n_node = this->nnode();
3034 
3035  // If we're on the pressure dof
3036  if (fld == n_dim)
3037  {
3038  // Call the pressure interpolation function
3039  this->pshape_nst(s, psi);
3040 
3041  // Allocate space for the pressure shape function
3042  Shape psif(n_node);
3043 
3044  // Allocate space for the pressure test function
3045  Shape testf(n_node);
3046 
3047  // Allocate space for the pressure shape function derivatives
3048  DShape dpsifdx(n_node, n_dim);
3049 
3050  // Allocate space for the pressure test function derivatives
3051  DShape dtestfdx(n_node, n_dim);
3052 
3053  // Calculate the Jacobian of the mapping
3054  double J = this->dshape_and_dtest_eulerian_nst(
3055  s, psif, dpsifdx, testf, dtestfdx);
3056 
3057  // Return the Jacobian
3058  return J;
3059  }
3060  // If we're on the velocity components
3061  else
3062  {
3063  // Allocate space for the test functions
3064  Shape testf(n_node);
3065 
3066  // Allocate space for the shape function derivatives
3067  DShape dpsifdx(n_node, n_dim);
3068 
3069  // Allocate space for the test function derivatives
3070  DShape dtestfdx(n_node, n_dim);
3071 
3072  // Calculate the Jacobian of the mapping
3073  double J =
3074  this->dshape_and_dtest_eulerian_nst(s, psi, dpsifdx, testf, dtestfdx);
3075 
3076  // Return the Jacobian
3077  return J;
3078  }
3079  } // End of jacobian_and_shape_of_field
J
JacobiRotation< float > J
Definition: Jacobi_makeJacobi.cpp:3

References J, and s.

◆ jacobian_and_shape_of_field() [2/3]

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

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

Implements oomph::ProjectableElementBase.

3045  {
3046  // How many dimensions in this element?
3047  unsigned n_dim = this->dim();
3048 
3049  // How many nodes in this element?
3050  unsigned n_node = this->nnode();
3051 
3052  // If we're on the pressure dof
3053  if (fld == n_dim)
3054  {
3055  // Call the pressure interpolation function
3056  this->pshape_nst(s, psi);
3057 
3058  // Allocate space for the pressure shape function
3059  Shape psif(n_node);
3060 
3061  // Allocate space for the pressure test function
3062  Shape testf(n_node);
3063 
3064  // Allocate space for the pressure shape function derivatives
3065  DShape dpsifdx(n_node, n_dim);
3066 
3067  // Allocate space for the pressure test function derivatives
3068  DShape dtestfdx(n_node, n_dim);
3069 
3070  // Calculate the Jacobian of the mapping
3071  double J = this->dshape_and_dtest_eulerian_nst(
3072  s, psif, dpsifdx, testf, dtestfdx);
3073 
3074  // Return the Jacobian
3075  return J;
3076  }
3077  // If we're on the velocity components
3078  else
3079  {
3080  // Allocate space for the test functions
3081  Shape testf(n_node);
3082 
3083  // Allocate space for the shape function derivatives
3084  DShape dpsifdx(n_node, n_dim);
3085 
3086  // Allocate space for the test function derivatives
3087  DShape dtestfdx(n_node, n_dim);
3088 
3089  // Calculate the Jacobian of the mapping
3090  double J =
3091  this->dshape_and_dtest_eulerian_nst(s, psi, dpsifdx, testf, dtestfdx);
3092 
3093  // Return the Jacobian
3094  return J;
3095  }
3096  } // End of jacobian_and_shape_of_field

References J, and s.

◆ jacobian_and_shape_of_field() [3/3]

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

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

Implements oomph::ProjectableElementBase.

2727  {
2728  // How many dimensions in this element?
2729  unsigned n_dim = this->dim();
2730 
2731  // How many nodes in this element?
2732  unsigned n_node = this->nnode();
2733 
2734  // If we're on the pressure dof
2735  if (fld == n_dim)
2736  {
2737  // Call the pressure interpolation function
2738  this->pshape_nst(s, psi);
2739 
2740  // Allocate space for the pressure shape function
2741  Shape psif(n_node);
2742 
2743  // Allocate space for the pressure test function
2744  Shape testf(n_node);
2745 
2746  // Allocate space for the pressure shape function derivatives
2747  DShape dpsifdx(n_node, n_dim);
2748 
2749  // Allocate space for the pressure test function derivatives
2750  DShape dtestfdx(n_node, n_dim);
2751 
2752  // Calculate the Jacobian of the mapping
2753  double J = this->dshape_and_dtest_eulerian_nst(
2754  s, psif, dpsifdx, testf, dtestfdx);
2755 
2756  // Return the Jacobian
2757  return J;
2758  }
2759  // If we're on the velocity components
2760  else
2761  {
2762  // Allocate space for the test functions
2763  Shape testf(n_node);
2764 
2765  // Allocate space for the shape function derivatives
2766  DShape dpsifdx(n_node, n_dim);
2767 
2768  // Allocate space for the test function derivatives
2769  DShape dtestfdx(n_node, n_dim);
2770 
2771  // Calculate the Jacobian of the mapping
2772  double J =
2773  this->dshape_and_dtest_eulerian_nst(s, psi, dpsifdx, testf, dtestfdx);
2774 
2775  // Return the Jacobian
2776  return J;
2777  }
2778  } // End of jacobian_and_shape_of_field

References J, and s.

◆ local_equation() [1/3]

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

Return local equation number of value j in field fld.

Implements oomph::ProjectableElementBase.

3137  {
3138  if (fld == this->dim())
3139  {
3140  return this->p_local_eqn(j);
3141  }
3142  else
3143  {
3144  const unsigned u_nodal_index = this->u_index_nst(fld);
3145  return this->nodal_local_eqn(j, u_nodal_index);
3146  }
3147  }
oomph::FiniteElement::nodal_local_eqn
int nodal_local_eqn(const unsigned &n, const unsigned &i) const
Definition: elements.h:1432

References j.

◆ local_equation() [2/3]

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

Return local equation number of value j in field fld.

Implements oomph::ProjectableElementBase.

3154  {
3155  if (fld == this->dim())
3156  {
3157  return this->p_local_eqn(j);
3158  }
3159  else
3160  {
3161  const unsigned u_nodal_index = this->u_index_nst(fld);
3162  return this->nodal_local_eqn(j, u_nodal_index);
3163  }
3164  }

References j.

◆ local_equation() [3/3]

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

Return local equation number of value j in field fld.

Implements oomph::ProjectableElementBase.

2836  {
2837  if (fld == this->dim())
2838  {
2839  return this->p_local_eqn(j);
2840  }
2841  else
2842  {
2843  const unsigned u_nodal_index = this->u_index_nst(fld);
2844  return this->nodal_local_eqn(j, u_nodal_index);
2845  }
2846  }

References j.

◆ nfields_for_projection() [1/3]

template<class TAYLOR_HOOD_ELEMENT >
unsigned oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::nfields_for_projection ( )
inlinevirtual

Number of fields to be projected: dim+1, corresponding to velocity components and pressure

Implements oomph::ProjectableElementBase.

2988  {
2989  // There are dim velocity dofs and 1 pressure dof
2990  return this->dim();
2991  } // End of nfields_for_projection

◆ nfields_for_projection() [2/3]

template<class TAYLOR_HOOD_ELEMENT >
unsigned oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::nfields_for_projection ( )
inlinevirtual

Number of fields to be projected: dim+1, corresponding to velocity components and pressure

Implements oomph::ProjectableElementBase.

3005  {
3006  // There are dim velocity dofs and 1 pressure dof
3007  return this->dim();
3008  } // End of nfields_for_projection

◆ nfields_for_projection() [3/3]

template<class TAYLOR_HOOD_ELEMENT >
unsigned oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::nfields_for_projection ( )
inlinevirtual

Number of fields to be projected: dim+1, corresponding to velocity components and pressure

Implements oomph::ProjectableElementBase.

2687  {
2688  // There are dim velocity dofs and 1 pressure dof
2689  return this->dim();
2690  } // End of nfields_for_projection

◆ nhistory_values_for_coordinate_projection() [1/3]

template<class TAYLOR_HOOD_ELEMENT >
unsigned oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::nhistory_values_for_coordinate_projection ( )
inlinevirtual

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

Implements oomph::ProjectableElementBase.

3017  {
3018  // Return the number of positional history values
3019  return this->node_pt(0)->position_time_stepper_pt()->ntstorage();
3020  } // End of nhistory_values_for_coordinate_projection
oomph::Node::position_time_stepper_pt
TimeStepper *& position_time_stepper_pt()
Return a pointer to the position timestepper.
Definition: nodes.h:1022
oomph::TimeStepper::ntstorage
unsigned ntstorage() const
Definition: timesteppers.h:601

◆ nhistory_values_for_coordinate_projection() [2/3]

template<class TAYLOR_HOOD_ELEMENT >
unsigned oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::nhistory_values_for_coordinate_projection ( )
inlinevirtual

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

Implements oomph::ProjectableElementBase.

3034  {
3035  // Return the number of positional history values
3036  return this->node_pt(0)->position_time_stepper_pt()->ntstorage();
3037  } // End of nhistory_values_for_coordinate_projection

◆ nhistory_values_for_coordinate_projection() [3/3]

template<class TAYLOR_HOOD_ELEMENT >
unsigned oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::nhistory_values_for_coordinate_projection ( )
inlinevirtual

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

Implements oomph::ProjectableElementBase.

2716  {
2717  // Return the number of positional history values
2718  return this->node_pt(0)->position_time_stepper_pt()->ntstorage();
2719  } // End of nhistory_values_for_coordinate_projection

◆ nhistory_values_for_projection() [1/3]

template<class TAYLOR_HOOD_ELEMENT >
unsigned oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::nhistory_values_for_projection ( const unsigned &  fld)
inlinevirtual

Number of history values to be stored for fld-th field. Whatever the timestepper has set up for the velocity components and none for the pressure field. NOTE: The count includes the current value!

Implements oomph::ProjectableElementBase.

2998  {
2999  // If we're dealing with the pressure dof
3000  if (fld == this->dim())
3001  {
3002  // Pressure doesn't have history values
3003  return this->node_pt(0)->ntstorage();
3004  }
3005  // If we're dealing with the velocity dofs
3006  else
3007  {
3008  // The velocity dofs have ntstorage() history values
3009  return this->node_pt(0)->ntstorage();
3010  }
3011  } // End of nhistory_values_for_projection
oomph::Data::ntstorage
unsigned ntstorage() const
Definition: nodes.cc:879

◆ nhistory_values_for_projection() [2/3]

template<class TAYLOR_HOOD_ELEMENT >
unsigned oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::nhistory_values_for_projection ( const unsigned &  fld)
inlinevirtual

Number of history values to be stored for fld-th field. Whatever the timestepper has set up for the velocity components and none for the pressure field. NOTE: The count includes the current value!

Implements oomph::ProjectableElementBase.

3015  {
3016  // If we're dealing with the pressure dof
3017  if (fld == this->dim())
3018  {
3019  // Pressure doesn't have history values
3020  return this->node_pt(0)->ntstorage();
3021  }
3022  // If we're dealing with the velocity dofs
3023  else
3024  {
3025  // The velocity dofs have ntstorage() history values
3026  return this->node_pt(0)->ntstorage();
3027  }
3028  } // End of nhistory_values_for_projection

◆ nhistory_values_for_projection() [3/3]

template<class TAYLOR_HOOD_ELEMENT >
unsigned oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::nhistory_values_for_projection ( const unsigned &  fld)
inlinevirtual

Number of history values to be stored for fld-th field. Whatever the timestepper has set up for the velocity components and none for the pressure field. NOTE: The count includes the current value!

Implements oomph::ProjectableElementBase.

2697  {
2698  // If we're dealing with the pressure dof
2699  if (fld == this->dim())
2700  {
2701  // Pressure doesn't have history values
2702  return this->node_pt(0)->ntstorage();
2703  }
2704  // If we're dealing with the velocity dofs
2705  else
2706  {
2707  // The velocity dofs have ntstorage() history values
2708  return this->node_pt(0)->ntstorage();
2709  }
2710  } // End of nhistory_values_for_projection

◆ nvalue_of_field() [1/3]

template<class TAYLOR_HOOD_ELEMENT >
unsigned oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::nvalue_of_field ( const unsigned &  fld)
inlinevirtual

Return number of values in field fld.

Implements oomph::ProjectableElementBase.

3123  {
3124  if (fld == this->dim())
3125  {
3126  return this->npres_nst();
3127  }
3128  else
3129  {
3130  return this->nnode();
3131  }
3132  }

◆ nvalue_of_field() [2/3]

template<class TAYLOR_HOOD_ELEMENT >
unsigned oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::nvalue_of_field ( const unsigned &  fld)
inlinevirtual

Return number of values in field fld.

Implements oomph::ProjectableElementBase.

3140  {
3141  if (fld == this->dim())
3142  {
3143  return this->npres_nst();
3144  }
3145  else
3146  {
3147  return this->nnode();
3148  }
3149  }

◆ nvalue_of_field() [3/3]

template<class TAYLOR_HOOD_ELEMENT >
unsigned oomph::ProjectableTaylorHoodSpaceTimeElement< TAYLOR_HOOD_ELEMENT >::nvalue_of_field ( const unsigned &  fld)
inlinevirtual

Return number of values in field fld.

Implements oomph::ProjectableElementBase.

2822  {
2823  if (fld == this->dim())
2824  {
2825  return this->npres_nst();
2826  }
2827  else
2828  {
2829  return this->nnode();
2830  }
2831  }
The documentation for this class was generated from the following files: