Public Member Functions | List of all members
oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_ELEMENT > Class Template Reference

SpaceTimeUnsteadyHeat upgraded to become projectable. More...

#include <space_time_unsteady_heat_elements.h>

+ Inheritance diagram for oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_ELEMENT >:

Public Member Functions

 ProjectableUnsteadyHeatSpaceTimeElement ()
 
Vector< std::pair< Data *, unsigned > > data_values_of_field (const unsigned &fld)
 
unsigned nfields_for_projection ()
 Number of fields to be projected: Just one. More...
 
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: One per node. More...
 
int local_equation (const unsigned &fld, const unsigned &j)
 Return local equation number of value j in field fld. More...
 
void output (std::ostream &outfile, const unsigned &nplot)
 
 ProjectableUnsteadyHeatSpaceTimeElement ()
 
Vector< std::pair< Data *, unsigned > > data_values_of_field (const unsigned &fld)
 
unsigned nfields_for_projection ()
 Number of fields to be projected: Just one. More...
 
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: One per node. More...
 
int local_equation (const unsigned &fld, const unsigned &j)
 Return local equation number of value j in field fld. More...
 
void output (std::ostream &outfile, const unsigned &nplot)
 
- Public Member Functions inherited from oomph::ProjectableElement< UNSTEADY_HEAT_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 (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< UNSTEADY_HEAT_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 UNSTEADY_HEAT_ELEMENT>
class oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_ELEMENT >

SpaceTimeUnsteadyHeat upgraded to become projectable.

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

Constructor & Destructor Documentation

◆ ProjectableUnsteadyHeatSpaceTimeElement() [1/2]

template<class UNSTEADY_HEAT_ELEMENT >
oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_ELEMENT >::ProjectableUnsteadyHeatSpaceTimeElement ( )
inline

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

931 {}

◆ ProjectableUnsteadyHeatSpaceTimeElement() [2/2]

template<class UNSTEADY_HEAT_ELEMENT >
oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_ELEMENT >::ProjectableUnsteadyHeatSpaceTimeElement ( )
inline

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

1013 {}

Member Function Documentation

◆ data_values_of_field() [1/2]

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

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

Implements oomph::ProjectableElementBase.

938  {
939 #ifdef PARANOID
940  // If we're not dealing with the first field
941  if (fld != 0)
942  {
943  // Create a stringstream object to create an error message
944  std::stringstream error_stream;
945 
946  // Create the error string
947  error_stream << "SpaceTimeUnsteadyHeat elements only store a single "
948  << "field so fld must be 0 rather than " << fld
949  << std::endl;
950 
951  // Throw an error
952  throw OomphLibError(
953  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
954  }
955 #endif
956 
957  // The number of nodes in this element
958  unsigned nnod = this->nnode();
959 
960  // Storage for the pairs
961  Vector<std::pair<Data*, unsigned>> data_values(nnod);
962 
963  // Loop over all nodes
964  for (unsigned j = 0; j < nnod; j++)
965  {
966  // Add the data value and associated field: The node itself
967  data_values[j] = std::make_pair(this->node_pt(j), fld);
968  }
969 
970  // Return the vector
971  return data_values;
972  } // 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::nnode
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210
OOMPH_EXCEPTION_LOCATION
#define OOMPH_EXCEPTION_LOCATION
Definition: oomph_definitions.h:61
OOMPH_CURRENT_FUNCTION
#define OOMPH_CURRENT_FUNCTION
Definition: oomph_definitions.h:86
j
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

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

◆ data_values_of_field() [2/2]

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

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

Implements oomph::ProjectableElementBase.

1020  {
1021 #ifdef PARANOID
1022  // If we're not dealing with the first field
1023  if (fld != 0)
1024  {
1025  // Create a stringstream object to create an error message
1026  std::stringstream error_stream;
1027 
1028  // Create the error string
1029  error_stream << "SpaceTimeUnsteadyHeat elements only store a single "
1030  << "field so fld must be 0 rather than " << fld
1031  << std::endl;
1032 
1033  // Throw an error
1034  throw OomphLibError(
1035  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1036  }
1037 #endif
1038 
1039  // The number of nodes in this element
1040  unsigned nnod = this->nnode();
1041 
1042  // Storage for the pairs
1043  Vector<std::pair<Data*, unsigned>> data_values(nnod);
1044 
1045  // Loop over all nodes
1046  for (unsigned j = 0; j < nnod; j++)
1047  {
1048  // Add the data value and associated field: The node itself
1049  data_values[j] = std::make_pair(this->node_pt(j), fld);
1050  }
1051 
1052  // Return the vector
1053  return data_values;
1054  } // End of data_values_of_field

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

◆ get_field() [1/2]

template<class UNSTEADY_HEAT_ELEMENT >
double oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_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.

1074  {
1075 #ifdef PARANOID
1076  // If we're not dealing with the first field
1077  if (fld != 0)
1078  {
1079  // Create a stringstream object to create an error message
1080  std::stringstream error_stream;
1081 
1082  // Create the error string
1083  error_stream << "SpaceTimeUnsteadyHeat elements only store a single "
1084  << "field so fld must be 0 rather than " << fld
1085  << std::endl;
1086 
1087  // Throw an error
1088  throw OomphLibError(
1089  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1090  }
1091 #endif
1092 
1093  // Find the index at which the variable is stored
1094  unsigned u_nodal_index = this->u_index_ust_heat();
1095 
1096  // Get the number of nodes in the element
1097  unsigned n_node = this->nnode();
1098 
1099  // Local shape function
1100  Shape psi(n_node);
1101 
1102  // Find values of shape function
1103  this->shape(s, psi);
1104 
1105  // Initialise value of u
1106  double interpolated_u = 0.0;
1107 
1108  // Loop over the local nodes
1109  for (unsigned l = 0; l < n_node; l++)
1110  {
1111  // Update the interpolated solution value
1112  interpolated_u += this->nodal_value(t, l, u_nodal_index) * psi[l];
1113  }
1114 
1115  // Return the interpolated solution value
1116  return interpolated_u;
1117  } // End of get_field
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

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

◆ get_field() [2/2]

template<class UNSTEADY_HEAT_ELEMENT >
double oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_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.

1156  {
1157 #ifdef PARANOID
1158  // If we're not dealing with the first field
1159  if (fld != 0)
1160  {
1161  // Create a stringstream object to create an error message
1162  std::stringstream error_stream;
1163 
1164  // Create the error string
1165  error_stream << "SpaceTimeUnsteadyHeat elements only store a single "
1166  << "field so fld must be 0 rather than " << fld
1167  << std::endl;
1168 
1169  // Throw an error
1170  throw OomphLibError(
1171  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1172  }
1173 #endif
1174 
1175  // Find the index at which the variable is stored
1176  unsigned u_nodal_index = this->u_index_ust_heat();
1177 
1178  // Get the number of nodes in the element
1179  unsigned n_node = this->nnode();
1180 
1181  // Local shape function
1182  Shape psi(n_node);
1183 
1184  // Find values of shape function
1185  this->shape(s, psi);
1186 
1187  // Initialise value of u
1188  double interpolated_u = 0.0;
1189 
1190  // Loop over the local nodes
1191  for (unsigned l = 0; l < n_node; l++)
1192  {
1193  // Update the interpolated solution value
1194  interpolated_u += this->nodal_value(t, l, u_nodal_index) * psi[l];
1195  }
1196 
1197  // Return the interpolated solution value
1198  return interpolated_u;
1199  } // End of get_field

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

◆ jacobian_and_shape_of_field() [1/2]

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

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

Implements oomph::ProjectableElementBase.

1024  {
1025 #ifdef PARANOID
1026  // If we're not dealing with the first field
1027  if (fld != 0)
1028  {
1029  // Create a stringstream object to create an error message
1030  std::stringstream error_stream;
1031 
1032  // Create the error string
1033  error_stream << "SpaceTimeUnsteadyHeat elements only store a single "
1034  << "field so fld must be 0 rather than " << fld
1035  << std::endl;
1036 
1037  // Throw an error
1038  throw OomphLibError(
1039  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1040  }
1041 #endif
1042 
1043  // Get the number of dimensions in the element
1044  unsigned n_dim = this->dim();
1045 
1046  // Get the number of nodes in the element
1047  unsigned n_node = this->nnode();
1048 
1049  // Allocate space for the test functions
1050  Shape test(n_node);
1051 
1052  // Allocate space for the derivatives of the shape functions
1053  DShape dpsidx(n_node, n_dim);
1054 
1055  // Allocate space for the derivatives of the test functions
1056  DShape dtestdx(n_node, n_dim);
1057 
1058  // Calculate the shape functions and their derivatives at the local
1059  // coordinate s (and the same for the test functions). On top of this
1060  // calculate the determinant of the Jacobian
1061  double J =
1062  this->dshape_and_dtest_eulerian_ust_heat(s, psi, dpsidx, test, dtestdx);
1063 
1064  // Return the determinant of the Jacobian
1065  return J;
1066  } // End of jacobian_and_shape_of_field
J
JacobiRotation< float > J
Definition: Jacobi_makeJacobi.cpp:3
oomph::FiniteElement::dim
unsigned dim() const
Definition: elements.h:2611
s
RealScalar s
Definition: level1_cplx_impl.h:130
Eigen::test
squared absolute sa ArrayBase::abs2 DOXCOMMA MatrixBase::cwiseAbs2 sa Eigen::abs2 DOXCOMMA Eigen::pow DOXCOMMA ArrayBase::square nearest sa Eigen::floor DOXCOMMA Eigen::ceil DOXCOMMA ArrayBase::round nearest integer not less than the given sa Eigen::floor DOXCOMMA ArrayBase::ceil not a number test
Definition: GlobalFunctions.h:109
test
Definition: indexed_view.cpp:20

References oomph::FiniteElement::dim(), J, oomph::FiniteElement::nnode(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, s, and Eigen::test.

◆ jacobian_and_shape_of_field() [2/2]

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

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

Implements oomph::ProjectableElementBase.

1106  {
1107 #ifdef PARANOID
1108  // If we're not dealing with the first field
1109  if (fld != 0)
1110  {
1111  // Create a stringstream object to create an error message
1112  std::stringstream error_stream;
1113 
1114  // Create the error string
1115  error_stream << "SpaceTimeUnsteadyHeat elements only store a single "
1116  << "field so fld must be 0 rather than " << fld
1117  << std::endl;
1118 
1119  // Throw an error
1120  throw OomphLibError(
1121  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1122  }
1123 #endif
1124 
1125  // Get the number of dimensions in the element
1126  unsigned n_dim = this->dim();
1127 
1128  // Get the number of nodes in the element
1129  unsigned n_node = this->nnode();
1130 
1131  // Allocate space for the test functions
1132  Shape test(n_node);
1133 
1134  // Allocate space for the derivatives of the shape functions
1135  DShape dpsidx(n_node, n_dim);
1136 
1137  // Allocate space for the derivatives of the test functions
1138  DShape dtestdx(n_node, n_dim);
1139 
1140  // Calculate the shape functions and their derivatives at the local
1141  // coordinate s (and the same for the test functions). On top of this
1142  // calculate the determinant of the Jacobian
1143  double J =
1144  this->dshape_and_dtest_eulerian_ust_heat(s, psi, dpsidx, test, dtestdx);
1145 
1146  // Return the determinant of the Jacobian
1147  return J;
1148  } // End of jacobian_and_shape_of_field

References oomph::FiniteElement::dim(), J, oomph::FiniteElement::nnode(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, s, and Eigen::test.

◆ local_equation() [1/2]

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

Return local equation number of value j in field fld.

Implements oomph::ProjectableElementBase.

1148  {
1149 #ifdef PARANOID
1150  // If we're not dealing with the first field
1151  if (fld != 0)
1152  {
1153  // Create a stringstream object to create an error message
1154  std::stringstream error_stream;
1155 
1156  // Create the error string
1157  error_stream << "SpaceTimeUnsteadyHeat elements only store a single "
1158  << "field so fld must be 0 rather than " << fld
1159  << std::endl;
1160 
1161  // Throw an error
1162  throw OomphLibError(
1163  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1164  }
1165 #endif
1166 
1167  // Get the nodal index of the unknown
1168  const unsigned u_nodal_index = this->u_index_ust_heat();
1169 
1170  // Output the local equation number
1171  return this->nodal_local_eqn(j, u_nodal_index);
1172  } // End of local_equation
oomph::FiniteElement::nodal_local_eqn
int nodal_local_eqn(const unsigned &n, const unsigned &i) const
Definition: elements.h:1432

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

◆ local_equation() [2/2]

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

Return local equation number of value j in field fld.

Implements oomph::ProjectableElementBase.

1230  {
1231 #ifdef PARANOID
1232  // If we're not dealing with the first field
1233  if (fld != 0)
1234  {
1235  // Create a stringstream object to create an error message
1236  std::stringstream error_stream;
1237 
1238  // Create the error string
1239  error_stream << "SpaceTimeUnsteadyHeat elements only store a single "
1240  << "field so fld must be 0 rather than " << fld
1241  << std::endl;
1242 
1243  // Throw an error
1244  throw OomphLibError(
1245  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1246  }
1247 #endif
1248 
1249  // Get the nodal index of the unknown
1250  const unsigned u_nodal_index = this->u_index_ust_heat();
1251 
1252  // Output the local equation number
1253  return this->nodal_local_eqn(j, u_nodal_index);
1254  } // End of local_equation

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

◆ nfields_for_projection() [1/2]

template<class UNSTEADY_HEAT_ELEMENT >
unsigned oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_ELEMENT >::nfields_for_projection ( )
inlinevirtual

Number of fields to be projected: Just one.

Implements oomph::ProjectableElementBase.

977  {
978  // Return the appropriate value
979  return 1;
980  } // End of nfields_for_projection

◆ nfields_for_projection() [2/2]

template<class UNSTEADY_HEAT_ELEMENT >
unsigned oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_ELEMENT >::nfields_for_projection ( )
inlinevirtual

Number of fields to be projected: Just one.

Implements oomph::ProjectableElementBase.

1059  {
1060  // Return the appropriate value
1061  return 1;
1062  } // End of nfields_for_projection

◆ nhistory_values_for_coordinate_projection() [1/2]

template<class UNSTEADY_HEAT_ELEMENT >
unsigned oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_ELEMENT >::nhistory_values_for_coordinate_projection ( )
inlinevirtual

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

Implements oomph::ProjectableElementBase.

1013  {
1014  // Return the number of history values stored by the position timestepper
1015  return this->node_pt(0)->position_time_stepper_pt()->ntstorage();
1016  } // 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

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

◆ nhistory_values_for_coordinate_projection() [2/2]

template<class UNSTEADY_HEAT_ELEMENT >
unsigned oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_ELEMENT >::nhistory_values_for_coordinate_projection ( )
inlinevirtual

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

Implements oomph::ProjectableElementBase.

1095  {
1096  // Return the number of history values stored by the position timestepper
1097  return this->node_pt(0)->position_time_stepper_pt()->ntstorage();
1098  } // End of nhistory_values_for_coordinate_projection

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

◆ nhistory_values_for_projection() [1/2]

template<class UNSTEADY_HEAT_ELEMENT >
unsigned oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_ELEMENT >::nhistory_values_for_projection ( const unsigned &  fld)
inlinevirtual

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

Implements oomph::ProjectableElementBase.

986  {
987 #ifdef PARANOID
988  // If we're not dealing with the first field
989  if (fld != 0)
990  {
991  // Create a stringstream object to create an error message
992  std::stringstream error_stream;
993 
994  // Create the error string
995  error_stream << "SpaceTimeUnsteadyHeat elements only store a single "
996  << "field so fld must be 0 rather than " << fld
997  << std::endl;
998 
999  // Throw an error
1000  throw OomphLibError(
1001  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1002  }
1003 #endif
1004 
1005  // Return the number of stored values
1006  return this->node_pt(0)->ntstorage();
1007  } // End of nhistory_values_for_projection
oomph::Data::ntstorage
unsigned ntstorage() const
Definition: nodes.cc:879

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

◆ nhistory_values_for_projection() [2/2]

template<class UNSTEADY_HEAT_ELEMENT >
unsigned oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_ELEMENT >::nhistory_values_for_projection ( const unsigned &  fld)
inlinevirtual

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

Implements oomph::ProjectableElementBase.

1068  {
1069 #ifdef PARANOID
1070  // If we're not dealing with the first field
1071  if (fld != 0)
1072  {
1073  // Create a stringstream object to create an error message
1074  std::stringstream error_stream;
1075 
1076  // Create the error string
1077  error_stream << "SpaceTimeUnsteadyHeat elements only store a single "
1078  << "field so fld must be 0 rather than " << fld
1079  << std::endl;
1080 
1081  // Throw an error
1082  throw OomphLibError(
1083  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1084  }
1085 #endif
1086 
1087  // Return the number of stored values
1088  return this->node_pt(0)->ntstorage();
1089  } // End of nhistory_values_for_projection

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

◆ nvalue_of_field() [1/2]

template<class UNSTEADY_HEAT_ELEMENT >
unsigned oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_ELEMENT >::nvalue_of_field ( const unsigned &  fld)
inlinevirtual

Return number of values in field fld: One per node.

Implements oomph::ProjectableElementBase.

1122  {
1123 #ifdef PARANOID
1124  // If we're not dealing with the first field
1125  if (fld != 0)
1126  {
1127  // Create a stringstream object to create an error message
1128  std::stringstream error_stream;
1129 
1130  // Create the error string
1131  error_stream << "SpaceTimeUnsteadyHeat elements only store a single "
1132  << "field so fld must be 0 rather than " << fld
1133  << std::endl;
1134 
1135  // Throw an error
1136  throw OomphLibError(
1137  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1138  }
1139 #endif
1140 
1141  // Return the number of nodes in the element
1142  return this->nnode();
1143  } // End of nvalue_of_field

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

◆ nvalue_of_field() [2/2]

template<class UNSTEADY_HEAT_ELEMENT >
unsigned oomph::ProjectableUnsteadyHeatSpaceTimeElement< UNSTEADY_HEAT_ELEMENT >::nvalue_of_field ( const unsigned &  fld)
inlinevirtual

Return number of values in field fld: One per node.

Implements oomph::ProjectableElementBase.

1204  {
1205 #ifdef PARANOID
1206  // If we're not dealing with the first field
1207  if (fld != 0)
1208  {
1209  // Create a stringstream object to create an error message
1210  std::stringstream error_stream;
1211 
1212  // Create the error string
1213  error_stream << "SpaceTimeUnsteadyHeat elements only store a single "
1214  << "field so fld must be 0 rather than " << fld
1215  << std::endl;
1216 
1217  // Throw an error
1218  throw OomphLibError(
1219  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
1220  }
1221 #endif
1222 
1223  // Return the number of nodes in the element
1224  return this->nnode();
1225  } // End of nvalue_of_field

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

◆ output() [1/2]

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

Output FE representation of soln: x,t,u or x,y,t,u at n_plot^(SPATIAL_DIM+1) plot points

Reimplemented from oomph::FiniteElement.

1178  {
1179  // Get the dimension of the element
1180  unsigned el_dim = this->dim();
1181 
1182  // Vector of local coordinates
1183  Vector<double> s(el_dim, 0.0);
1184 
1185  // Tecplot header info
1186  outfile << this->tecplot_zone_string(nplot);
1187 
1188  // Get the number of plot points
1189  unsigned num_plot_points = this->nplot_points(nplot);
1190 
1191  // Loop over plot points
1192  for (unsigned iplot = 0; iplot < num_plot_points; iplot++)
1193  {
1194  // Get local coordinates of plot point
1195  this->get_s_plot(iplot, nplot, s);
1196 
1197  // Loop over the coordinate directions
1198  for (unsigned i = 0; i < el_dim; i++)
1199  {
1200  // Output the interpolated coordinates
1201  outfile << this->interpolated_x(s, i) << " ";
1202  }
1203 
1204  // Output the interpolated value of u(s)
1205  outfile << this->interpolated_u_ust_heat(s) << " ";
1206 
1207  // Output the interpolated value of du/dt(s)
1208  outfile << this->interpolated_du_dt_ust_heat(s) << " ";
1209 
1210  // History values of coordinates
1211  unsigned n_prev =
1212  this->node_pt(0)->position_time_stepper_pt()->ntstorage();
1213 
1214  // Loop over the previous timesteps
1215  for (unsigned t = 1; t < n_prev; t++)
1216  {
1217  // Loop over the coordinate directions
1218  for (unsigned i = 0; i < el_dim; i++)
1219  {
1220  // Output the coordinates
1221  outfile << this->interpolated_x(t, s, i) << " ";
1222  }
1223  } // for (unsigned t=1;t<n_prev;t++)
1224 
1225  // Number of history values of velocities
1226  n_prev = this->node_pt(0)->time_stepper_pt()->ntstorage();
1227 
1228  // Loop over the previous timesteps
1229  for (unsigned t = 1; t < n_prev; t++)
1230  {
1231  // Output the solution
1232  outfile << this->interpolated_u_ust_heat(t, s) << " ";
1233  }
1234 
1235  // Finish the line
1236  outfile << std::endl;
1237  } // for (unsigned iplot=0;iplot<num_plot_points;iplot++)
1238 
1239  // Write tecplot footer (e.g. FE connectivity lists)
1240  this->write_tecplot_zone_footer(outfile, nplot);
1241  } // End of output
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9
oomph::Data::time_stepper_pt
TimeStepper *& time_stepper_pt()
Return the pointer to the timestepper.
Definition: nodes.h:238
oomph::FiniteElement::tecplot_zone_string
virtual std::string tecplot_zone_string(const unsigned &nplot) const
Definition: elements.h:3161
oomph::FiniteElement::interpolated_x
virtual double interpolated_x(const Vector< double > &s, const unsigned &i) const
Return FE interpolated coordinate x[i] at local coordinate s.
Definition: elements.cc:3962
oomph::FiniteElement::get_s_plot
virtual void get_s_plot(const unsigned &i, const unsigned &nplot, Vector< double > &s, const bool &shifted_to_interior=false) const
Definition: elements.h:3148
oomph::FiniteElement::nplot_points
virtual unsigned nplot_points(const unsigned &nplot) const
Definition: elements.h:3186
oomph::FiniteElement::write_tecplot_zone_footer
virtual void write_tecplot_zone_footer(std::ostream &outfile, const unsigned &nplot) const
Definition: elements.h:3174
plotPSD.t
t
Definition: plotPSD.py:36
el_dim
unsigned el_dim
dimension
Definition: overloaded_cartesian_element_body.h:30

References oomph::FiniteElement::dim(), el_dim, oomph::FiniteElement::get_s_plot(), i, oomph::FiniteElement::interpolated_x(), oomph::FiniteElement::node_pt(), oomph::FiniteElement::nplot_points(), oomph::TimeStepper::ntstorage(), oomph::Node::position_time_stepper_pt(), s, plotPSD::t, oomph::FiniteElement::tecplot_zone_string(), oomph::Data::time_stepper_pt(), and oomph::FiniteElement::write_tecplot_zone_footer().

◆ output() [2/2]

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

Output FE representation of soln: x,t,u or x,y,t,u at n_plot^(SPATIAL_DIM+1) plot points

Reimplemented from oomph::FiniteElement.

1260  {
1261  // Get the dimension of the element
1262  unsigned el_dim = this->dim();
1263 
1264  // Vector of local coordinates
1265  Vector<double> s(el_dim, 0.0);
1266 
1267  // Tecplot header info
1268  outfile << this->tecplot_zone_string(nplot);
1269 
1270  // Get the number of plot points
1271  unsigned num_plot_points = this->nplot_points(nplot);
1272 
1273  // Loop over plot points
1274  for (unsigned iplot = 0; iplot < num_plot_points; iplot++)
1275  {
1276  // Get local coordinates of plot point
1277  this->get_s_plot(iplot, nplot, s);
1278 
1279  // Loop over the coordinate directions
1280  for (unsigned i = 0; i < el_dim; i++)
1281  {
1282  // Output the interpolated coordinates
1283  outfile << this->interpolated_x(s, i) << " ";
1284  }
1285 
1286  // Output the interpolated value of u(s)
1287  outfile << this->interpolated_u_ust_heat(s) << " ";
1288 
1289  // Output the interpolated value of du/dt(s)
1290  outfile << this->interpolated_du_dt_ust_heat(s) << " ";
1291 
1292  // History values of coordinates
1293  unsigned n_prev =
1294  this->node_pt(0)->position_time_stepper_pt()->ntstorage();
1295 
1296  // Loop over the previous timesteps
1297  for (unsigned t = 1; t < n_prev; t++)
1298  {
1299  // Loop over the coordinate directions
1300  for (unsigned i = 0; i < el_dim; i++)
1301  {
1302  // Output the coordinates
1303  outfile << this->interpolated_x(t, s, i) << " ";
1304  }
1305  } // for (unsigned t=1;t<n_prev;t++)
1306 
1307  // Number of history values of velocities
1308  n_prev = this->node_pt(0)->time_stepper_pt()->ntstorage();
1309 
1310  // Loop over the previous timesteps
1311  for (unsigned t = 1; t < n_prev; t++)
1312  {
1313  // Output the solution
1314  outfile << this->interpolated_u_ust_heat(t, s) << " ";
1315  }
1316 
1317  // Finish the line
1318  outfile << std::endl;
1319  } // for (unsigned iplot=0;iplot<num_plot_points;iplot++)
1320 
1321  // Write tecplot footer (e.g. FE connectivity lists)
1322  this->write_tecplot_zone_footer(outfile, nplot);
1323  } // End of output

References oomph::FiniteElement::dim(), el_dim, oomph::FiniteElement::get_s_plot(), i, oomph::FiniteElement::interpolated_x(), oomph::FiniteElement::node_pt(), oomph::FiniteElement::nplot_points(), oomph::TimeStepper::ntstorage(), oomph::Node::position_time_stepper_pt(), s, plotPSD::t, oomph::FiniteElement::tecplot_zone_string(), oomph::Data::time_stepper_pt(), and oomph::FiniteElement::write_tecplot_zone_footer().

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