oomph::RefineableQTaylorHoodSpaceTimeElement< DIM > Class Template Reference

#include <refineable_discontinuous_galerkin_equal_order_pressure_spacetime_navier_stokes_elements.h>

+ Inheritance diagram for oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >:

Public Member Functions

 RefineableQTaylorHoodSpaceTimeElement ()
 Constructor. More...
 
unsigned required_nvalue (const unsigned &n) const
 
unsigned ncont_interpolated_values () const
 Number of continuously interpolated values. More...
 
void rebuild_from_sons (Mesh *&mesh_pt)
 Rebuild from sons: empty. More...
 
unsigned nrecovery_order ()
 
unsigned nvertex_node () const
 Number of vertex nodes in the element. More...
 
Nodevertex_node_pt (const unsigned &j) const
 Pointer to the j-th vertex node in the element. More...
 
void get_interpolated_values (const Vector< double > &s, Vector< double > &values)
 
void get_interpolated_values (const unsigned &t, const Vector< double > &s, Vector< double > &values)
 
void further_setup_hanging_nodes ()
 
Nodepressure_node_pt (const unsigned &n_p)
 Pointer to n_p-th pressure node. More...
 
Nodeinterpolating_node_pt (const unsigned &n, const int &value_id)
 
double local_one_d_fraction_of_interpolating_node (const unsigned &n_1d, const unsigned &i, const int &value_id)
 
Nodeget_interpolating_node_at_local_coordinate (const Vector< double > &s, const int &value_id)
 
unsigned ninterpolating_node_1d (const int &value_id)
 
unsigned ninterpolating_node (const int &value_id)
 
void interpolating_basis (const Vector< double > &s, Shape &psi, const int &value_id) const
 
void build_fp_press_adv_diff_robin_bc_element (const unsigned &face_index)
 
void identify_load_data (std::set< std::pair< Data *, unsigned >> &paired_load_data)
 
 RefineableQTaylorHoodSpaceTimeElement ()
 Constructor. More...
 
unsigned required_nvalue (const unsigned &n) const
 
unsigned ncont_interpolated_values () const
 Number of continuously interpolated values. More...
 
void rebuild_from_sons (Mesh *&mesh_pt)
 Rebuild from sons: empty. More...
 
unsigned nrecovery_order ()
 
unsigned nvertex_node () const
 Number of vertex nodes in the element. More...
 
Nodevertex_node_pt (const unsigned &j) const
 Pointer to the j-th vertex node in the element. More...
 
void get_interpolated_values (const Vector< double > &s, Vector< double > &values)
 
void get_interpolated_values (const unsigned &t, const Vector< double > &s, Vector< double > &values)
 
void further_setup_hanging_nodes ()
 
Nodepressure_node_pt (const unsigned &n_p)
 Pointer to n_p-th pressure node. More...
 
Nodeinterpolating_node_pt (const unsigned &n, const int &value_id)
 
double local_one_d_fraction_of_interpolating_node (const unsigned &n_1d, const unsigned &i, const int &value_id)
 
Nodeget_interpolating_node_at_local_coordinate (const Vector< double > &s, const int &value_id)
 
unsigned ninterpolating_node_1d (const int &value_id)
 
unsigned ninterpolating_node (const int &value_id)
 
void interpolating_basis (const Vector< double > &s, Shape &psi, const int &value_id) const
 
void build_fp_press_adv_diff_robin_bc_element (const unsigned &face_index)
 
void identify_load_data (std::set< std::pair< Data *, unsigned >> &paired_load_data)
 
 RefineableQTaylorHoodSpaceTimeElement ()
 Constructor. More...
 
unsigned required_nvalue (const unsigned &n) const
 
unsigned ncont_interpolated_values () const
 Number of continuously interpolated values. More...
 
void rebuild_from_sons (Mesh *&mesh_pt)
 Rebuild from sons: empty. More...
 
unsigned nrecovery_order ()
 
unsigned nvertex_node () const
 Number of vertex nodes in the element. More...
 
Nodevertex_node_pt (const unsigned &j) const
 Pointer to the j-th vertex node in the element. More...
 
void get_interpolated_values (const Vector< double > &s, Vector< double > &values)
 
void get_interpolated_values (const unsigned &t, const Vector< double > &s, Vector< double > &values)
 
void further_setup_hanging_nodes ()
 
Nodepressure_node_pt (const unsigned &n_p)
 Pointer to n_p-th pressure node. More...
 
Nodeinterpolating_node_pt (const unsigned &n, const int &value_id)
 
double local_one_d_fraction_of_interpolating_node (const unsigned &n_1d, const unsigned &i, const int &value_id)
 
Nodeget_interpolating_node_at_local_coordinate (const Vector< double > &s, const int &value_id)
 
unsigned ninterpolating_node_1d (const int &value_id)
 
unsigned ninterpolating_node (const int &value_id)
 
void interpolating_basis (const Vector< double > &s, Shape &psi, const int &value_id) const
 
void build_fp_press_adv_diff_robin_bc_element (const unsigned &face_index)
 
void identify_load_data (std::set< std::pair< Data *, unsigned >> &paired_load_data)
 
- Public Member Functions inherited from oomph::QTaylorHoodSpaceTimeElement< DIM >
 QTaylorHoodSpaceTimeElement ()
 Constructor, no internal data points. More...
 
void pshape_nst (const Vector< double > &s, Shape &psi) const
 Pressure shape functions at local coordinate s. More...
 
void ptest_nst (const Vector< double > &s, Shape &psi) const
 Pressure test functions at local coordinate s. More...
 
void pshape_nst (const Vector< double > &s, Shape &psi, Shape &test) const
 Pressure shape and test functions at local coordinte s. More...
 
virtual int p_nodal_index_nst () const
 Set the value at which the pressure is stored in the nodes. More...
 
int p_local_eqn (const unsigned &n) const
 Return the local equation numbers for the pressure values. More...
 
double p_nst (const unsigned &n_p) const
 
double p_nst (const unsigned &t, const unsigned &n_p) const
 
unsigned npres_nst () const
 Return number of pressure values. More...
 
void fix_pressure (const unsigned &p_dof, const double &p_value)
 Pin p_dof-th pressure dof and set it to value specified by p_value. More...
 
void identify_pressure_data (std::set< std::pair< Data *, unsigned >> &paired_pressure_data)
 
void output (std::ostream &outfile)
 Redirect output to NavierStokesEquations output. More...
 
void output (std::ostream &outfile, const unsigned &nplot)
 Redirect output to NavierStokesEquations output. More...
 
void output (FILE *file_pt)
 Redirect output to NavierStokesEquations output. More...
 
void output (FILE *file_pt, const unsigned &nplot)
 Redirect output to NavierStokesEquations output. More...
 
unsigned ndof_types () const
 
void get_dof_numbers_for_unknowns (std::list< std::pair< unsigned long, unsigned >> &dof_lookup_list) const
 
 QTaylorHoodSpaceTimeElement ()
 Constructor, no internal data points. More...
 
 ~QTaylorHoodSpaceTimeElement ()
 
void pshape_nst (const Vector< double > &s, Shape &psi) const
 Pressure shape functions at local coordinate s. More...
 
void ptest_nst (const Vector< double > &s, Shape &psi) const
 Pressure test functions at local coordinate s. More...
 
void pshape_nst (const Vector< double > &s, Shape &psi, Shape &test) const
 Pressure shape and test functions at local coordinte s. More...
 
virtual int p_nodal_index_nst () const
 Set the value at which the pressure is stored in the nodes. More...
 
int p_local_eqn (const unsigned &n) const
 Return the local equation numbers for the pressure values. More...
 
double p_nst (const unsigned &n_p) const
 
double p_nst (const unsigned &t, const unsigned &n_p) const
 
unsigned npres_nst () const
 Return number of pressure values. More...
 
void fix_pressure (const unsigned &p_dof, const double &p_value)
 Pin p_dof-th pressure dof and set it to value specified by p_value. More...
 
void identify_pressure_data (std::set< std::pair< Data *, unsigned >> &paired_pressure_data)
 
void output (std::ostream &outfile)
 Redirect output to NavierStokesEquations output. More...
 
void output (std::ostream &outfile, const unsigned &nplot)
 Redirect output to NavierStokesEquations output. More...
 
void output (FILE *file_pt)
 Redirect output to NavierStokesEquations output. More...
 
void output (FILE *file_pt, const unsigned &nplot)
 Redirect output to NavierStokesEquations output. More...
 
 QTaylorHoodSpaceTimeElement ()
 Constructor, no internal data points. More...
 
void pshape_nst (const Vector< double > &s, Shape &psi) const
 Pressure shape functions at local coordinate s. More...
 
void pshape_nst (const Vector< double > &s, Shape &psi, Shape &test) const
 Pressure shape and test functions at local coordinte s. More...
 
virtual int p_nodal_index_nst () const
 Set the value at which the pressure is stored in the nodes. More...
 
int p_local_eqn (const unsigned &n) const
 Return the local equation numbers for the pressure values. More...
 
double p_nst (const unsigned &n_p) const
 
double p_nst (const unsigned &t, const unsigned &n_p) const
 
unsigned npres_nst () const
 Return number of pressure values. More...
 
void fix_pressure (const unsigned &p_dof, const double &p_value)
 Pin p_dof-th pressure dof and set it to value specified by p_value. More...
 
void build_fp_press_adv_diff_robin_bc_element (const unsigned &face_index)
 
void identify_load_data (std::set< std::pair< Data *, unsigned >> &paired_load_data)
 
void identify_pressure_data (std::set< std::pair< Data *, unsigned >> &paired_pressure_data)
 
void output (std::ostream &outfile)
 Redirect output to NavierStokesEquations output. More...
 
void output (std::ostream &outfile, const unsigned &nplot)
 Redirect output to NavierStokesEquations output. More...
 
void output (FILE *file_pt)
 Redirect output to NavierStokesEquations output. More...
 
void output (FILE *file_pt, const unsigned &nplot)
 Redirect output to NavierStokesEquations output. More...
 
void pshape_nst (const Vector< double > &s, Shape &psi) const
 2D (in space): Pressure shape functions More...
 
void ptest_nst (const Vector< double > &s, Shape &test) const
 2D (in space): Pressure shape functions More...
 
void pshape_nst (const Vector< double > &s, Shape &psi) const
 2D (in space): Pressure shape functions More...
 
void ptest_nst (const Vector< double > &s, Shape &test) const
 2D (in space): Pressure shape functions More...
 
void pshape_nst (const Vector< double > &s, Shape &psi) const
 
- Public Member Functions inherited from oomph::SpaceTimeNavierStokesEquations< DIM >
 SpaceTimeNavierStokesEquations ()
 
void store_strouhal_as_external_data (Data *strouhal_data_pt)
 Function that tells us whether the period is stored as external data. More...
 
const doublere () const
 Reynolds number. More...
 
double *& re_pt ()
 Pointer to Reynolds number. More...
 
bool is_strouhal_stored_as_external_data () const
 Are we storing the Strouhal number as external data? More...
 
const doublest () const
 Strouhal parameter (const. version) More...
 
doublest_pt () const
 Pointer to Strouhal parameter (const. version) More...
 
double *& st_pt ()
 Pointer to Strouhal number (can only assign to private member data) More...
 
double re_st () const
 Product of Reynolds and Strouhal number (=Womersley number) More...
 
const doubleviscosity_ratio () const
 
double *& viscosity_ratio_pt ()
 Pointer to Viscosity Ratio. More...
 
const doubledensity_ratio () const
 
double *& density_ratio_pt ()
 Pointer to Density ratio. More...
 
const doublere_invfr () const
 Global inverse Froude number. More...
 
double *& re_invfr_pt ()
 Pointer to global inverse Froude number. More...
 
const Vector< double > & g () const
 Vector of gravitational components. More...
 
Vector< double > *& g_pt ()
 Pointer to Vector of gravitational components. More...
 
NavierStokesBodyForceFctPtbody_force_fct_pt ()
 Access function for the body-force pointer. More...
 
NavierStokesBodyForceFctPt body_force_fct_pt () const
 Access function for the body-force pointer. Const version. More...
 
NavierStokesSourceFctPtsource_fct_pt ()
 Access function for the source-function pointer. More...
 
NavierStokesSourceFctPt source_fct_pt () const
 Access function for the source-function pointer. Const version. More...
 
NavierStokesPressureAdvDiffSourceFctPtsource_fct_for_pressure_adv_diff ()
 
NavierStokesPressureAdvDiffSourceFctPt source_fct_for_pressure_adv_diff () const
 
intpinned_fp_pressure_eqn ()
 
double u_nst (const unsigned &n, const unsigned &i) const
 
double u_nst (const unsigned &t, const unsigned &n, const unsigned &i) const
 
virtual unsigned u_index_nst (const unsigned &i) const
 
unsigned n_u_nst () const
 
double get_du_dt (const unsigned &n, const unsigned &i) const
 
double du_dt_nst (const unsigned &n, const unsigned &i) const
 
double interpolated_du_dt_nst (const Vector< double > &s, const unsigned &i) const
 
void disable_ALE ()
 
void enable_ALE ()
 
double pressure_integral () const
 Integral of pressure over element. More...
 
double dissipation () const
 Return integral of dissipation over element. More...
 
double dissipation (const Vector< double > &s) const
 Return dissipation at local coordinate s. More...
 
void get_vorticity (const Vector< double > &s, Vector< double > &vorticity) const
 Compute the vorticity vector at local coordinate s. More...
 
void get_vorticity (const Vector< double > &s, double &vorticity) const
 Compute the vorticity vector at local coordinate s. More...
 
double kin_energy () const
 Get integral of kinetic energy over element. More...
 
double d_kin_energy_dt () const
 Get integral of time derivative of kinetic energy over element. More...
 
void strain_rate (const Vector< double > &s, DenseMatrix< double > &strain_rate) const
 Strain-rate tensor: 1/2 (du_i/dx_j+du_j/dx_i) More...
 
void get_traction (const Vector< double > &s, const Vector< double > &N, Vector< double > &traction)
 
void get_traction (const Vector< double > &s, const Vector< double > &N, Vector< double > &traction_p, Vector< double > &traction_visc_n, Vector< double > &traction_visc_t)
 
void get_load (const Vector< double > &s, const Vector< double > &N, Vector< double > &load)
 
unsigned nscalar_paraview () const
 
void scalar_value_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot) const
 
void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt) const
 
std::string scalar_name_paraview (const unsigned &i) const
 
void full_output (std::ostream &outfile)
 
void full_output (std::ostream &outfile, const unsigned &n_plot)
 
void output_veloc (std::ostream &outfile, const unsigned &nplot, const unsigned &t)
 
void output_vorticity (std::ostream &outfile, const unsigned &nplot)
 
void output_fct (std::ostream &outfile, const unsigned &nplot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 
void output_fct (std::ostream &outfile, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 
void compute_norm (Vector< double > &norm)
 Compute the vector norm of the FEM solution. More...
 
void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 
void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 
void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 
void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Compute the element's residual Vector. More...
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_contribution_to_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
void fill_in_contribution_to_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
void fill_in_contribution_to_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
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)
 
void fill_in_pressure_advection_diffusion_residuals (Vector< double > &residuals)
 
void fill_in_pressure_advection_diffusion_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void pin_all_non_pressure_dofs (std::map< Data *, std::vector< int >> &eqn_number_backup)
 Pin all non-pressure dofs and backup eqn numbers. More...
 
void output_pressure_advection_diffusion_robin_elements (std::ostream &outfile)
 
void delete_pressure_advection_diffusion_robin_elements ()
 
void interpolated_u_nst (const Vector< double > &s, Vector< double > &velocity) const
 Compute vector of FE interpolated velocity u at local coordinate s. More...
 
double interpolated_u_nst (const Vector< double > &s, const unsigned &i) const
 Return FE interpolated velocity u[i] at local coordinate s. More...
 
double interpolated_u_nst (const unsigned &t, const Vector< double > &s, const unsigned &i) const
 
virtual double interpolated_p_nst (const Vector< double > &s) const
 Return FE interpolated pressure at local coordinate s. More...
 
double interpolated_p_nst (const unsigned &t, const Vector< double > &s) const
 Return FE interpolated pressure at local coordinate s at time level t. More...
 
void point_output_data (const Vector< double > &s, Vector< double > &data)
 
 SpaceTimeNavierStokesEquations ()
 
void store_strouhal_as_external_data (Data *strouhal_data_pt)
 Function that tells us whether the period is stored as external data. More...
 
const doublere () const
 Reynolds number. More...
 
double *& re_pt ()
 Pointer to Reynolds number. More...
 
bool is_strouhal_stored_as_external_data () const
 Are we storing the Strouhal number as external data? More...
 
const doublest () const
 Strouhal parameter (const. version) More...
 
doublest_pt () const
 Pointer to Strouhal parameter (const. version) More...
 
double *& st_pt ()
 Pointer to Strouhal number (can only assign to private member data) More...
 
double re_st () const
 Product of Reynolds and Strouhal number (=Womersley number) More...
 
const doubleviscosity_ratio () const
 
double *& viscosity_ratio_pt ()
 Pointer to Viscosity Ratio. More...
 
const doubledensity_ratio () const
 
double *& density_ratio_pt ()
 Pointer to Density ratio. More...
 
const doublere_invfr () const
 Global inverse Froude number. More...
 
double *& re_invfr_pt ()
 Pointer to global inverse Froude number. More...
 
const Vector< double > & g () const
 Vector of gravitational components. More...
 
Vector< double > *& g_pt ()
 Pointer to Vector of gravitational components. More...
 
NavierStokesBodyForceFctPtbody_force_fct_pt ()
 Access function for the body-force pointer. More...
 
NavierStokesBodyForceFctPt body_force_fct_pt () const
 Access function for the body-force pointer. Const version. More...
 
NavierStokesSourceFctPtsource_fct_pt ()
 Access function for the source-function pointer. More...
 
NavierStokesSourceFctPt source_fct_pt () const
 Access function for the source-function pointer. Const version. More...
 
NavierStokesPressureAdvDiffSourceFctPtsource_fct_for_pressure_adv_diff ()
 
NavierStokesPressureAdvDiffSourceFctPt source_fct_for_pressure_adv_diff () const
 
intpinned_fp_pressure_eqn ()
 
double u_nst (const unsigned &n, const unsigned &i) const
 
double u_nst (const unsigned &t, const unsigned &n, const unsigned &i) const
 
virtual unsigned u_index_nst (const unsigned &i) const
 
unsigned n_u_nst () const
 
double get_du_dt (const unsigned &n, const unsigned &i) const
 
double du_dt_nst (const unsigned &n, const unsigned &i) const
 
double interpolated_du_dt_nst (const Vector< double > &s, const unsigned &i) const
 
void disable_ALE ()
 
void enable_ALE ()
 
double pressure_integral () const
 Integral of pressure over element. More...
 
double dissipation () const
 Return integral of dissipation over element. More...
 
double dissipation (const Vector< double > &s) const
 Return dissipation at local coordinate s. More...
 
void get_vorticity (const Vector< double > &s, Vector< double > &vorticity) const
 Compute the vorticity vector at local coordinate s. More...
 
void get_vorticity (const Vector< double > &s, double &vorticity) const
 Compute the vorticity vector at local coordinate s. More...
 
double kin_energy () const
 Get integral of kinetic energy over element. More...
 
double d_kin_energy_dt () const
 Get integral of time derivative of kinetic energy over element. More...
 
void strain_rate (const Vector< double > &s, DenseMatrix< double > &strain_rate) const
 Strain-rate tensor: 1/2 (du_i/dx_j+du_j/dx_i) More...
 
void get_traction (const Vector< double > &s, const Vector< double > &N, Vector< double > &traction)
 
void get_traction (const Vector< double > &s, const Vector< double > &N, Vector< double > &traction_p, Vector< double > &traction_visc_n, Vector< double > &traction_visc_t)
 
void get_load (const Vector< double > &s, const Vector< double > &N, Vector< double > &load)
 
unsigned nscalar_paraview () const
 
void scalar_value_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot) const
 
void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt) const
 
std::string scalar_name_paraview (const unsigned &i) const
 
void full_output (std::ostream &outfile)
 
void full_output (std::ostream &outfile, const unsigned &n_plot)
 
void output_veloc (std::ostream &outfile, const unsigned &nplot, const unsigned &t)
 
void output_vorticity (std::ostream &outfile, const unsigned &nplot)
 
void output_fct (std::ostream &outfile, const unsigned &nplot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 
void output_fct (std::ostream &outfile, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 
void compute_norm (Vector< double > &norm)
 Compute the vector norm of the FEM solution. More...
 
void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 
void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 
void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 
void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Compute the element's residual Vector. More...
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_contribution_to_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
void fill_in_contribution_to_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
void fill_in_contribution_to_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
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)
 
void fill_in_pressure_advection_diffusion_residuals (Vector< double > &residuals)
 
void fill_in_pressure_advection_diffusion_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void pin_all_non_pressure_dofs (std::map< Data *, std::vector< int >> &eqn_number_backup)
 Pin all non-pressure dofs and backup eqn numbers. More...
 
void output_pressure_advection_diffusion_robin_elements (std::ostream &outfile)
 
void delete_pressure_advection_diffusion_robin_elements ()
 
void interpolated_u_nst (const Vector< double > &s, Vector< double > &velocity) const
 Compute vector of FE interpolated velocity u at local coordinate s. More...
 
double interpolated_u_nst (const Vector< double > &s, const unsigned &i) const
 Return FE interpolated velocity u[i] at local coordinate s. More...
 
double interpolated_u_nst (const unsigned &t, const Vector< double > &s, const unsigned &i) const
 
virtual double interpolated_p_nst (const Vector< double > &s) const
 Return FE interpolated pressure at local coordinate s. More...
 
double interpolated_p_nst (const unsigned &t, const Vector< double > &s) const
 Return FE interpolated pressure at local coordinate s at time level t. More...
 
void point_output_data (const Vector< double > &s, Vector< double > &data)
 
 SpaceTimeNavierStokesEquations ()
 
void store_strouhal_as_external_data (Data *strouhal_data_pt)
 Function that tells us whether the period is stored as external data. More...
 
bool is_reynolds_strouhal_stored_as_external_data () const
 Are we storing the Strouhal number as external data? More...
 
const doublere () const
 Reynolds number. More...
 
double *& re_pt ()
 Pointer to Reynolds number. More...
 
const doublere_st () const
 ReSt parameter (const. version) More...
 
doublere_st_pt () const
 Pointer to Strouhal parameter (const. version) More...
 
double *& re_st_pt ()
 Pointer to ReSt number (can only assign to private member data) More...
 
const doubleviscosity_ratio () const
 
double *& viscosity_ratio_pt ()
 Pointer to Viscosity Ratio. More...
 
const doubledensity_ratio () const
 
double *& density_ratio_pt ()
 Pointer to Density ratio. More...
 
const doublere_invfr () const
 Global inverse Froude number. More...
 
double *& re_invfr_pt ()
 Pointer to global inverse Froude number. More...
 
const Vector< double > & g () const
 Vector of gravitational components. More...
 
Vector< double > *& g_pt ()
 Pointer to Vector of gravitational components. More...
 
NavierStokesBodyForceFctPtbody_force_fct_pt ()
 Access function for the body-force pointer. More...
 
NavierStokesBodyForceFctPt body_force_fct_pt () const
 Access function for the body-force pointer. Const version. More...
 
NavierStokesSourceFctPtsource_fct_pt ()
 Access function for the source-function pointer. More...
 
NavierStokesSourceFctPt source_fct_pt () const
 Access function for the source-function pointer. Const version. More...
 
NavierStokesPressureAdvDiffSourceFctPtsource_fct_for_pressure_adv_diff ()
 
NavierStokesPressureAdvDiffSourceFctPt source_fct_for_pressure_adv_diff () const
 
intpinned_fp_pressure_eqn ()
 
double u_nst (const unsigned &n, const unsigned &i) const
 
double u_nst (const unsigned &t, const unsigned &n, const unsigned &i) const
 
virtual unsigned u_index_nst (const unsigned &i) const
 
unsigned n_u_nst () const
 
double get_du_dt (const unsigned &n, const unsigned &i) const
 
double du_dt_nst (const unsigned &n, const unsigned &i) const
 
double interpolated_du_dt_nst (const Vector< double > &s, const unsigned &i) const
 
void disable_ALE ()
 
void enable_ALE ()
 
double pressure_integral () const
 Integral of pressure over element. More...
 
double dissipation () const
 Return integral of dissipation over element. More...
 
double dissipation (const Vector< double > &s) const
 Return dissipation at local coordinate s. More...
 
void get_vorticity (const Vector< double > &s, Vector< double > &vorticity) const
 Compute the vorticity vector at local coordinate s. More...
 
void get_vorticity (const Vector< double > &s, double &vorticity) const
 Compute the vorticity vector at local coordinate s. More...
 
double kin_energy () const
 Get integral of kinetic energy over element. More...
 
double d_kin_energy_dt () const
 Get integral of time derivative of kinetic energy over element. More...
 
void strain_rate (const Vector< double > &s, DenseMatrix< double > &strain_rate) const
 Strain-rate tensor: 1/2 (du_i/dx_j+du_j/dx_i) More...
 
void get_traction (const Vector< double > &s, const Vector< double > &N, Vector< double > &traction)
 
void get_traction (const Vector< double > &s, const Vector< double > &N, Vector< double > &traction_p, Vector< double > &traction_visc_n, Vector< double > &traction_visc_t)
 
void get_load (const Vector< double > &s, const Vector< double > &N, Vector< double > &load)
 
void get_pressure_and_velocity_mass_matrix_diagonal (Vector< double > &press_mass_diag, Vector< double > &veloc_mass_diag, const unsigned &which_one=0)
 
unsigned nscalar_paraview () const
 
void scalar_value_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot) const
 
void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt) const
 
std::string scalar_name_paraview (const unsigned &i) const
 
void full_output (std::ostream &outfile)
 
void full_output (std::ostream &outfile, const unsigned &n_plot)
 
void output_veloc (std::ostream &outfile, const unsigned &nplot, const unsigned &t)
 
void output_vorticity (std::ostream &outfile, const unsigned &nplot)
 
void output_fct (std::ostream &outfile, const unsigned &nplot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 
void output_fct (std::ostream &outfile, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 
void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 
void compute_norm (Vector< double > &norm)
 Compute the vector norm of the FEM solution. More...
 
void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 
void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 
void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Compute the element's residual Vector. More...
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void fill_in_contribution_to_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
void fill_in_contribution_to_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
void fill_in_contribution_to_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
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)
 
void fill_in_pressure_advection_diffusion_residuals (Vector< double > &residuals)
 
void fill_in_pressure_advection_diffusion_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
void pin_all_non_pressure_dofs (std::map< Data *, std::vector< int >> &eqn_number_backup)
 Pin all non-pressure dofs and backup eqn numbers. More...
 
void output_pressure_advection_diffusion_robin_elements (std::ostream &outfile)
 
void delete_pressure_advection_diffusion_robin_elements ()
 
void interpolated_u_nst (const Vector< double > &s, Vector< double > &velocity) const
 Compute vector of FE interpolated velocity u at local coordinate s. More...
 
double interpolated_u_nst (const Vector< double > &s, const unsigned &i) const
 Return FE interpolated velocity u[i] at local coordinate s. More...
 
double interpolated_u_nst (const unsigned &t, const Vector< double > &s, const unsigned &i) const
 
virtual double interpolated_p_nst (const Vector< double > &s) const
 Return FE interpolated pressure at local coordinate s. More...
 
double interpolated_p_nst (const unsigned &t, const Vector< double > &s) const
 Return FE interpolated pressure at local coordinate s at time level t. More...
 
void point_output_data (const Vector< double > &s, Vector< double > &data)
 
void get_vorticity (const Vector< double > &s, Vector< double > &vorticity) const
 
void get_vorticity (const Vector< double > &s, double &vorticity) const
 
void get_vorticity (const Vector< double > &s, Vector< double > &vorticity) const
 
void get_vorticity (const Vector< double > &s, double &vorticity) const
 
void get_vorticity (const Vector< double > &s, Vector< double > &vorticity) const
 
void get_vorticity (const Vector< double > &s, double &vorticity) const
 
- Public Member Functions inherited from oomph::FSIFluidElement
 FSIFluidElement ()
 Constructor. More...
 
 FSIFluidElement (const FSIFluidElement &)=delete
 Broken copy constructor. More...
 
void operator= (const FSIFluidElement &)=delete
 Broken assignment operator. More...
 
- Public Member Functions inherited from oomph::FiniteElement
void set_dimension (const unsigned &dim)
 
void set_nodal_dimension (const unsigned &nodal_dim)
 
void set_nnodal_position_type (const unsigned &nposition_type)
 Set the number of types required to interpolate the coordinate. More...
 
void set_n_node (const unsigned &n)
 
int nodal_local_eqn (const unsigned &n, const unsigned &i) const
 
double dJ_eulerian_at_knot (const unsigned &ipt, Shape &psi, DenseMatrix< double > &djacobian_dX) const
 
 FiniteElement ()
 Constructor. More...
 
virtual ~FiniteElement ()
 
 FiniteElement (const FiniteElement &)=delete
 Broken copy constructor. More...
 
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)
 
MacroElementmacro_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 describe_local_dofs (std::ostream &out, const std::string &current_string) const
 
virtual void describe_nodal_local_dofs (std::ostream &out, const std::string &current_string) const
 
virtual void assign_all_generic_local_eqn_numbers (const bool &store_local_dof_pt)
 
Node *& node_pt (const unsigned &n)
 Return a pointer to the local node n. More...
 
Node *const & node_pt (const unsigned &n) const
 Return a pointer to the local node n (const version) More...
 
unsigned nnode () const
 Return the number of nodes. More...
 
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
 
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 Nodeconstruct_node (const unsigned &n)
 
virtual Nodeconstruct_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 
virtual Nodeconstruct_boundary_node (const unsigned &n)
 
virtual Nodeconstruct_boundary_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 
int get_node_number (Node *const &node_pt) const
 
virtual Nodeget_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
 
Datageom_data_pt (const unsigned &j)
 
void position (const Vector< double > &zeta, Vector< double > &r) const
 
void position (const unsigned &t, const Vector< double > &zeta, Vector< double > &r) const
 
void dposition_dt (const Vector< double > &zeta, const unsigned &t, Vector< double > &drdt)
 
virtual double zeta_nodal (const unsigned &n, const unsigned &k, const unsigned &i) const
 
void interpolated_zeta (const Vector< double > &s, Vector< double > &zeta) const
 
void locate_zeta (const Vector< double > &zeta, GeomObject *&geom_object_pt, Vector< double > &s, const bool &use_coordinate_as_initial_guess=false)
 
virtual void node_update ()
 
virtual void identify_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
 
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 void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt) const
 
virtual void output (const unsigned &t, std::ostream &outfile, const unsigned &n_plot) const
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, 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, 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, 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 (double &norm)
 
- 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 ()
 
TimeSteppertime_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
 
- Public Member Functions inherited from oomph::TemplateFreeSpaceTimeNavierStokesEquationsBase
 TemplateFreeSpaceTimeNavierStokesEquationsBase ()
 Constructor (empty) More...
 
virtual ~TemplateFreeSpaceTimeNavierStokesEquationsBase ()
 Virtual destructor (empty) More...
 
 TemplateFreeSpaceTimeNavierStokesEquationsBase ()
 Constructor (empty) More...
 
virtual ~TemplateFreeSpaceTimeNavierStokesEquationsBase ()
 Virtual destructor (empty) More...
 
 TemplateFreeSpaceTimeNavierStokesEquationsBase ()
 Constructor (empty) More...
 
virtual ~TemplateFreeSpaceTimeNavierStokesEquationsBase ()
 Virtual destructor (empty) More...
 
- Public Member Functions inherited from oomph::NavierStokesElementWithDiagonalMassMatrices
 NavierStokesElementWithDiagonalMassMatrices ()
 Empty constructor. More...
 
virtual ~NavierStokesElementWithDiagonalMassMatrices ()
 Virtual destructor. More...
 
 NavierStokesElementWithDiagonalMassMatrices (const NavierStokesElementWithDiagonalMassMatrices &)=delete
 Broken copy constructor. More...
 
void operator= (const NavierStokesElementWithDiagonalMassMatrices &)=delete
 Broken assignment operator. More...
 
- Public Member Functions inherited from oomph::RefineableSpaceTimeNavierStokesEquations< DIM >
 RefineableSpaceTimeNavierStokesEquations ()
 Constructor. More...
 
void get_pressure_and_velocity_mass_matrix_diagonal (Vector< double > &press_mass_diag, Vector< double > &veloc_mass_diag, const unsigned &which_one=0)
 
unsigned num_Z2_flux_terms ()
 Number of 'flux' terms for Z2 error estimation. More...
 
void get_Z2_flux (const Vector< double > &s, Vector< double > &flux)
 
void further_build ()
 Further build, pass the pointers down to the sons. More...
 
void dinterpolated_u_nst_ddata (const Vector< double > &s, const unsigned &i, Vector< double > &du_ddata, Vector< unsigned > &global_eqn_number)
 
 RefineableSpaceTimeNavierStokesEquations ()
 Constructor. More...
 
void get_pressure_and_velocity_mass_matrix_diagonal (Vector< double > &press_mass_diag, Vector< double > &veloc_mass_diag, const unsigned &which_one=0)
 
unsigned num_Z2_flux_terms ()
 Number of 'flux' terms for Z2 error estimation. More...
 
void get_Z2_flux (const Vector< double > &s, Vector< double > &flux)
 
void further_build ()
 Further build, pass the pointers down to the sons. More...
 
void dinterpolated_u_nst_ddata (const Vector< double > &s, const unsigned &i, Vector< double > &du_ddata, Vector< unsigned > &global_eqn_number)
 
 RefineableSpaceTimeNavierStokesEquations ()
 Constructor. More...
 
void get_pressure_and_velocity_mass_matrix_diagonal (Vector< double > &press_mass_diag, Vector< double > &veloc_mass_diag, const unsigned &which_one=0)
 
unsigned num_Z2_flux_terms ()
 Number of 'flux' terms for Z2 error estimation. More...
 
void get_Z2_flux (const Vector< double > &s, Vector< double > &flux)
 
void further_build ()
 Further build, pass the pointers down to the sons. More...
 
void dinterpolated_u_nst_ddata (const Vector< double > &s, const unsigned &i, Vector< double > &du_ddata, Vector< unsigned > &global_eqn_number)
 
- Public Member Functions inherited from oomph::RefineableElement
 RefineableElement ()
 
virtual ~RefineableElement ()
 Destructor, delete the allocated storage for the hanging equations. More...
 
 RefineableElement (const RefineableElement &)=delete
 Broken copy constructor. More...
 
void operator= (const RefineableElement &)=delete
 Broken assignment operator. More...
 
Treetree_pt ()
 Access function: Pointer to quadtree representation of this element. More...
 
void set_tree_pt (Tree *my_tree_pt)
 Set pointer to quadtree representation of this element. More...
 
virtual unsigned required_nsons () const
 
bool refinement_is_enabled ()
 Flag to indicate suppression of any refinement. More...
 
void disable_refinement ()
 Suppress of any refinement for this element. More...
 
void enable_refinement ()
 Emnable refinement for this element. More...
 
template<class ELEMENT >
void split (Vector< ELEMENT * > &son_pt) const
 
int local_hang_eqn (Node *const &node_pt, const unsigned &i)
 
virtual void build (Mesh *&mesh_pt, Vector< Node * > &new_node_pt, bool &was_already_built, std::ofstream &new_nodes_file)=0
 
void set_refinement_level (const int &refine_level)
 Set the refinement level. More...
 
unsigned refinement_level () const
 Return the Refinement level. More...
 
void select_for_refinement ()
 Select the element for refinement. More...
 
void deselect_for_refinement ()
 Deselect the element for refinement. More...
 
void select_sons_for_unrefinement ()
 Unrefinement will be performed by merging the four sons of this element. More...
 
void deselect_sons_for_unrefinement ()
 
bool to_be_refined ()
 Has the element been selected for refinement? More...
 
bool sons_to_be_unrefined ()
 Has the element been selected for unrefinement? More...
 
virtual void unbuild ()
 
virtual void deactivate_element ()
 
virtual bool nodes_built ()
 Return true if all the nodes have been built, false if not. More...
 
long number () const
 Element number (for debugging/plotting) More...
 
void set_number (const long &mynumber)
 Set element number (for debugging/plotting) More...
 
virtual void check_integrity (double &max_error)=0
 
void identify_field_data_for_interactions (std::set< std::pair< Data *, unsigned >> &paired_field_data)
 
void assign_nodal_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual RefineableElementroot_element_pt ()
 
virtual RefineableElementfather_element_pt () const
 Return a pointer to the father element. More...
 
void get_father_at_refinement_level (unsigned &refinement_level, RefineableElement *&father_at_reflevel_pt)
 
virtual void initial_setup (Tree *const &adopted_father_pt=0, const unsigned &initial_p_order=0)
 
virtual void pre_build (Mesh *&mesh_pt, Vector< Node * > &new_node_pt)
 Pre-build the element. More...
 
virtual void setup_hanging_nodes (Vector< std::ofstream * > &output_stream)
 
unsigned nshape_controlling_nodes ()
 
std::map< Node *, unsignedshape_controlling_node_lookup ()
 
- Public Member Functions inherited from oomph::ElementWithZ2ErrorEstimator
 ElementWithZ2ErrorEstimator ()
 Default empty constructor. More...
 
 ElementWithZ2ErrorEstimator (const ElementWithZ2ErrorEstimator &)=delete
 Broken copy constructor. More...
 
void operator= (const ElementWithZ2ErrorEstimator &)=delete
 Broken assignment operator. More...
 
virtual unsigned ncompound_fluxes ()
 
virtual void compute_exact_Z2_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_flux_pt, double &error, double &norm)
 
virtual void get_Z2_compound_flux_indices (Vector< unsigned > &flux_index)
 
virtual double geometric_jacobian (const Vector< double > &x)
 
- Public Member Functions inherited from oomph::RefineableQElement< DIM+1 >
 RefineableQElement ()
 Empty constuctor. More...
 

Private Member Functions

void unpin_elemental_pressure_dofs ()
 Unpin all pressure dofs. More...
 
void pin_elemental_redundant_nodal_pressure_dofs ()
 Pin all nodal pressure dofs that are not required. More...
 
void unpin_elemental_pressure_dofs ()
 Unpin all pressure dofs. More...
 
void pin_elemental_redundant_nodal_pressure_dofs ()
 Pin all nodal pressure dofs that are not required. More...
 
void unpin_elemental_pressure_dofs ()
 Unpin all pressure dofs. More...
 
void pin_elemental_redundant_nodal_pressure_dofs ()
 Pin all nodal pressure dofs that are not required. More...
 

Additional Inherited Members

- Public Types inherited from oomph::SpaceTimeNavierStokesEquations< DIM >
typedef void(* NavierStokesBodyForceFctPt) (const double &time, const Vector< double > &x, Vector< double > &body_force)
 
typedef double(* NavierStokesSourceFctPt) (const double &time, const Vector< double > &x)
 Function pointer to source function fct(t,x) (x is a Vector!) More...
 
typedef double(* NavierStokesPressureAdvDiffSourceFctPt) (const Vector< double > &x)
 
typedef void(* NavierStokesBodyForceFctPt) (const double &time, const Vector< double > &x, Vector< double > &body_force)
 
typedef double(* NavierStokesSourceFctPt) (const double &time, const Vector< double > &x)
 Function pointer to source function fct(t,x) (x is a Vector!) More...
 
typedef double(* NavierStokesPressureAdvDiffSourceFctPt) (const Vector< double > &x)
 
typedef void(* NavierStokesBodyForceFctPt) (const double &time, const Vector< double > &x, Vector< double > &body_force)
 
typedef double(* NavierStokesSourceFctPt) (const double &time, const Vector< double > &x)
 Function pointer to source function fct(t,x) (x is a Vector!) More...
 
typedef double(* NavierStokesPressureAdvDiffSourceFctPt) (const Vector< double > &x)
 
- 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 Member Functions inherited from oomph::RefineableSpaceTimeNavierStokesEquations< DIM >
static void pin_redundant_nodal_pressures (const Vector< GeneralisedElement * > &element_pt)
 
static void unpin_all_pressure_dofs (const Vector< GeneralisedElement * > &element_pt)
 Unpin all pressure dofs in elements listed in vector. More...
 
static void pin_redundant_nodal_pressures (const Vector< GeneralisedElement * > &element_pt)
 
static void unpin_all_pressure_dofs (const Vector< GeneralisedElement * > &element_pt)
 Unpin all pressure dofs in elements listed in vector. More...
 
static void pin_redundant_nodal_pressures (const Vector< GeneralisedElement * > &element_pt)
 
static void unpin_all_pressure_dofs (const Vector< GeneralisedElement * > &element_pt)
 Unpin all pressure dofs in elements listed in vector. More...
 
- Static Public Member Functions inherited from oomph::RefineableElement
static doublemax_integrity_tolerance ()
 Max. allowed discrepancy in element integrity check. More...
 
- Static Public Attributes inherited from oomph::SpaceTimeNavierStokesEquations< DIM >
static Vector< doubleGamma
 Navier-Stokes equations static data. More...
 
- Static Public Attributes inherited from oomph::FiniteElement
static double Tolerance_for_singular_jacobian = 1.0e-16
 Tolerance below which the jacobian is considered singular. More...
 
static bool Accept_negative_jacobian = false
 
static bool Suppress_output_while_checking_for_inverted_elements
 
- Static Public Attributes inherited from oomph::GeneralisedElement
static bool Suppress_warning_about_repeated_internal_data
 
static bool Suppress_warning_about_repeated_external_data = true
 
static double Default_fd_jacobian_step = 1.0e-8
 
- Protected Member Functions inherited from oomph::QTaylorHoodSpaceTimeElement< DIM >
double dshape_and_dtest_eulerian_nst (const Vector< double > &s, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_at_knot_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_at_knot_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, RankFourTensor< double > &d_dpsidx_dX, Shape &test, DShape &dtestdx, RankFourTensor< double > &d_dtestdx_dX, DenseMatrix< double > &djacobian_dX) const
 
double dpshape_eulerian (const Vector< double > &s, Shape &ppsi, DShape &dppsidx) const
 
double dptest_eulerian (const Vector< double > &s, Shape &ptest, DShape &dptestdx) const
 
double dpshape_and_dptest_eulerian_nst (const Vector< double > &s, Shape &ppsi, DShape &dppsidx, Shape &ptest, DShape &dptestdx) const
 
double dshape_and_dtest_eulerian_nst (const Vector< double > &s, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_at_knot_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_at_knot_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, RankFourTensor< double > &d_dpsidx_dX, Shape &test, DShape &dtestdx, RankFourTensor< double > &d_dtestdx_dX, DenseMatrix< double > &djacobian_dX) const
 
double dpshape_eulerian (const Vector< double > &s, Shape &ppsi, DShape &dppsidx) const
 
double dptest_eulerian (const Vector< double > &s, Shape &ptest, DShape &dptestdx) const
 
double dpshape_and_dptest_eulerian_nst (const Vector< double > &s, Shape &ppsi, DShape &dppsidx, Shape &ptest, DShape &dptestdx) const
 
double dshape_and_dtest_eulerian_nst (const Vector< double > &s, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_at_knot_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const
 
double dshape_and_dtest_eulerian_at_knot_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, RankFourTensor< double > &d_dpsidx_dX, Shape &test, DShape &dtestdx, RankFourTensor< double > &d_dtestdx_dX, DenseMatrix< double > &djacobian_dX) const
 
double dpshape_and_dptest_eulerian_nst (const Vector< double > &s, Shape &ppsi, DShape &dppsidx, Shape &ptest, DShape &dptestdx) const
 
const unsigned Pconv [12]
 
double dpshape_eulerian (const Vector< double > &s, Shape &ppsi, DShape &dppsidx) const
 
double dptest_eulerian (const Vector< double > &s, Shape &ptest, DShape &dptestdx) const
 
double dpshape_and_dptest_eulerian_nst (const Vector< double > &s, Shape &ppsi, DShape &dppsidx, Shape &ptest, DShape &dptestdx) const
 
double dshape_and_dtest_eulerian_at_knot_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, RankFourTensor< double > &d_dpsidx_dX, Shape &test, DShape &dtestdx, RankFourTensor< double > &d_dtestdx_dX, DenseMatrix< double > &djacobian_dX) const
 
const unsigned Pconv [8]
 
double dpshape_eulerian (const Vector< double > &s, Shape &ppsi, DShape &dppsidx) const
 
double dptest_eulerian (const Vector< double > &s, Shape &ptest, DShape &dptestdx) const
 
double dpshape_and_dptest_eulerian_nst (const Vector< double > &s, Shape &ppsi, DShape &dppsidx, Shape &ptest, DShape &dptestdx) const
 
double dshape_and_dtest_eulerian_at_knot_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, RankFourTensor< double > &d_dpsidx_dX, Shape &test, DShape &dtestdx, RankFourTensor< double > &d_dtestdx_dX, DenseMatrix< double > &djacobian_dX) const
 
const unsigned Pconv [8]
 
double dpshape_and_dptest_eulerian_nst (const Vector< double > &s, Shape &ppsi, DShape &dppsidx, Shape &ptest, DShape &dptestdx) const
 
double dshape_and_dtest_eulerian_at_knot_nst (const unsigned &ipt, Shape &psi, DShape &dpsidx, RankFourTensor< double > &d_dpsidx_dX, Shape &test, DShape &dtestdx, RankFourTensor< double > &d_dtestdx_dX, DenseMatrix< double > &djacobian_dX) const
 
- Protected Member Functions inherited from oomph::SpaceTimeNavierStokesEquations< DIM >
virtual void get_body_force_nst (const double &time, const unsigned &ipt, const Vector< double > &s, const Vector< double > &x, Vector< double > &result)
 
virtual void get_body_force_gradient_nst (const double &time, const unsigned &ipt, const Vector< double > &s, const Vector< double > &x, DenseMatrix< double > &d_body_force_dx)
 
virtual double get_source_nst (const double &time, const unsigned &ipt, const Vector< double > &x)
 Calculate the source fct at a given time and Eulerian position. More...
 
virtual void get_source_gradient_nst (const double &time, const unsigned &ipt, const Vector< double > &x, Vector< double > &gradient)
 
virtual void fill_in_generic_dresidual_contribution_nst (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam, const unsigned &flag)
 
void fill_in_contribution_to_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
virtual void get_body_force_nst (const double &time, const unsigned &ipt, const Vector< double > &s, const Vector< double > &x, Vector< double > &result)
 
virtual void get_body_force_gradient_nst (const double &time, const unsigned &ipt, const Vector< double > &s, const Vector< double > &x, DenseMatrix< double > &d_body_force_dx)
 
virtual double get_source_nst (const double &time, const unsigned &ipt, const Vector< double > &x)
 Calculate the source fct at a given time and Eulerian position. More...
 
virtual void get_source_gradient_nst (const double &time, const unsigned &ipt, const Vector< double > &x, Vector< double > &gradient)
 
virtual void fill_in_generic_dresidual_contribution_nst (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam, const unsigned &flag)
 
void fill_in_contribution_to_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
virtual void get_body_force_nst (const double &time, const unsigned &ipt, const Vector< double > &s, const Vector< double > &x, Vector< double > &result)
 
virtual void get_body_force_gradient_nst (const double &time, const unsigned &ipt, const Vector< double > &s, const Vector< double > &x, DenseMatrix< double > &d_body_force_dx)
 
virtual double get_source_nst (const double &time, const unsigned &ipt, const Vector< double > &x)
 Calculate the source fct at a given time and Eulerian position. More...
 
virtual void get_source_gradient_nst (const double &time, const unsigned &ipt, const Vector< double > &x, Vector< double > &gradient)
 
virtual void fill_in_generic_dresidual_contribution_nst (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam, const unsigned &flag)
 
void fill_in_contribution_to_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
- Protected Member Functions inherited from oomph::FiniteElement
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 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
 
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_internal_and_external_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void assign_additional_local_eqn_numbers ()
 
int internal_local_eqn (const unsigned &i, const unsigned &j) const
 
int external_local_eqn (const unsigned &i, const unsigned &j)
 
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_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 Member Functions inherited from oomph::RefineableSpaceTimeNavierStokesEquations< DIM >
void fill_in_generic_residual_contribution_nst (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix, const unsigned &compute_jacobian_flag)
 
void fill_in_generic_pressure_advection_diffusion_contribution_nst (Vector< double > &residuals, DenseMatrix< double > &jacobian, const unsigned &compute_jacobian_flag)
 
virtual void get_dresidual_dnodal_coordinates (RankThreeTensor< double > &dresidual_dnodal_coordinates)
 
void fill_in_generic_residual_contribution_nst (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix, const unsigned &compute_jacobian_flag)
 
void fill_in_generic_pressure_advection_diffusion_contribution_nst (Vector< double > &residuals, DenseMatrix< double > &jacobian, const unsigned &compute_jacobian_flag)
 
virtual void get_dresidual_dnodal_coordinates (RankThreeTensor< double > &dresidual_dnodal_coordinates)
 
void fill_in_generic_residual_contribution_nst (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix, const unsigned &compute_jacobian_flag)
 
void fill_in_generic_pressure_advection_diffusion_contribution_nst (Vector< double > &residuals, DenseMatrix< double > &jacobian, const unsigned &compute_jacobian_flag)
 
virtual void get_dresidual_dnodal_coordinates (RankThreeTensor< double > &dresidual_dnodal_coordinates)
 
- Protected Member Functions inherited from oomph::RefineableElement
void assemble_local_to_eulerian_jacobian (const DShape &dpsids, DenseMatrix< double > &jacobian) const
 
void assemble_local_to_eulerian_jacobian2 (const DShape &d2psids, DenseMatrix< double > &jacobian2) const
 
void assemble_eulerian_base_vectors (const DShape &dpsids, DenseMatrix< double > &interpolated_G) const
 
double local_to_eulerian_mapping_diagonal (const DShape &dpsids, DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
void assign_hanging_local_eqn_numbers (const bool &store_local_dof_pt)
 Assign the local equation numbers for hanging node variables. More...
 
virtual void fill_in_jacobian_from_nodal_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
- Static Protected Member Functions inherited from oomph::RefineableElement
static void check_value_id (const int &n_continuously_interpolated_values, const int &value_id)
 
- Protected Attributes inherited from oomph::SpaceTimeNavierStokesEquations< DIM >
doubleViscosity_Ratio_pt
 
doubleDensity_Ratio_pt
 
doubleRe_pt
 Pointer to global Reynolds number. More...
 
doubleSt_pt
 Pointer to global Reynolds number x Strouhal number (=Womersley) More...
 
bool Strouhal_is_stored_as_external_data
 
doubleReInvFr_pt
 
Vector< double > * G_pt
 Pointer to global gravity Vector. More...
 
NavierStokesBodyForceFctPt Body_force_fct_pt
 Pointer to body force function. More...
 
NavierStokesSourceFctPt Source_fct_pt
 Pointer to volumetric source function. More...
 
NavierStokesPressureAdvDiffSourceFctPt Press_adv_diff_source_fct_pt
 
bool ALE_is_disabled
 
Vector< FpPressureAdvDiffRobinBCSpaceTimeElementBase * > Pressure_advection_diffusion_robin_element_pt
 
int Pinned_fp_pressure_eqn
 
doubleReSt_pt
 Pointer to global Reynolds number x Strouhal number (= Womersley) More...
 
bool ReynoldsStrouhal_is_stored_as_external_data
 
- Protected Attributes inherited from oomph::FiniteElement
MacroElementMacro_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...
 
TimeStepperGeom_object_time_stepper_pt
 
- Protected Attributes inherited from oomph::RefineableElement
TreeTree_pt
 A pointer to a general tree object. More...
 
unsigned Refine_level
 Refinement level. More...
 
bool To_be_refined
 Flag for refinement. More...
 
bool Refinement_is_enabled
 Flag to indicate suppression of any refinement. More...
 
bool Sons_to_be_unrefined
 Flag for unrefinement. More...
 
long Number
 Global element number – for plotting/validation purposes. More...
 
- Static Protected Attributes inherited from oomph::QTaylorHoodSpaceTimeElement< DIM >
static const unsigned Pconv []
 
- 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< doubleDummy_matrix
 
static std::deque< double * > Dof_pt_deque
 
- Static Protected Attributes inherited from oomph::RefineableElement
static double Max_integrity_tolerance = 1.0e-8
 Max. allowed discrepancy in element integrity check. More...
 

Detailed Description

template<unsigned DIM>
class oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >

Refineable version of Taylor Hood elements. These classes can be written in total generality.

Constructor & Destructor Documentation

◆ RefineableQTaylorHoodSpaceTimeElement() [1/3]

Constructor.

766  : RefineableElement(),
767  RefineableSpaceTimeNavierStokesEquations<DIM>(),
768  RefineableQElement<DIM + 1>(),
769  QTaylorHoodSpaceTimeElement<DIM>()
770  {
771  }
RefineableElement()
Definition: refineable_elements.h:188

◆ RefineableQTaylorHoodSpaceTimeElement() [2/3]

Constructor.

766  : RefineableElement(),
767  RefineableSpaceTimeNavierStokesEquations<DIM>(),
768  RefineableQElement<DIM + 1>(),
769  QTaylorHoodSpaceTimeElement<DIM>()
770  {
771  }

◆ RefineableQTaylorHoodSpaceTimeElement() [3/3]

Constructor.

766  : RefineableElement(),
767  RefineableSpaceTimeNavierStokesEquations<DIM>(),
768  RefineableQElement<DIM + 1>(),
769  QTaylorHoodSpaceTimeElement<DIM>()
770  {
771  }

Member Function Documentation

◆ build_fp_press_adv_diff_robin_bc_element() [1/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::build_fp_press_adv_diff_robin_bc_element ( const unsigned face_index)
inlinevirtual

Build FaceElements that apply the Robin boundary condition to the pressure advection diffusion problem required by Fp preconditioner

Reimplemented from oomph::QTaylorHoodSpaceTimeElement< DIM >.

1074  {
1076  new RefineableFpPressureAdvDiffRobinBCSpaceTimeElement<
1077  RefineableQTaylorHoodSpaceTimeElement<DIM>>(this, face_index));
1078  } // End of build_fp_press_adv_diff_robin_bc_element
Vector< FpPressureAdvDiffRobinBCSpaceTimeElementBase * > Pressure_advection_diffusion_robin_element_pt
Definition: discontinuous_galerkin_equal_order_pressure_spacetime_navier_stokes_elements.h:548

References oomph::SpaceTimeNavierStokesEquations< DIM >::Pressure_advection_diffusion_robin_element_pt.

◆ build_fp_press_adv_diff_robin_bc_element() [2/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::build_fp_press_adv_diff_robin_bc_element ( const unsigned face_index)
inlinevirtual

Build FaceElements that apply the Robin boundary condition to the pressure advection diffusion problem required by Fp preconditioner

Reimplemented from oomph::QTaylorHoodSpaceTimeElement< DIM >.

1074  {
1076  new RefineableFpPressureAdvDiffRobinBCSpaceTimeElement<
1077  RefineableQTaylorHoodSpaceTimeElement<DIM>>(this, face_index));
1078  } // End of build_fp_press_adv_diff_robin_bc_element

References oomph::SpaceTimeNavierStokesEquations< DIM >::Pressure_advection_diffusion_robin_element_pt.

◆ build_fp_press_adv_diff_robin_bc_element() [3/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::build_fp_press_adv_diff_robin_bc_element ( const unsigned face_index)
inlinevirtual

Build FaceElements that apply the Robin boundary condition to the pressure advection diffusion problem required by Fp preconditioner

Reimplemented from oomph::QTaylorHoodSpaceTimeElement< DIM >.

1074  {
1076  new RefineableFpPressureAdvDiffRobinBCSpaceTimeElement<
1077  RefineableQTaylorHoodSpaceTimeElement<DIM>>(this, face_index));
1078  } // End of build_fp_press_adv_diff_robin_bc_element

References oomph::SpaceTimeNavierStokesEquations< DIM >::Pressure_advection_diffusion_robin_element_pt.

◆ further_setup_hanging_nodes() [1/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::further_setup_hanging_nodes ( )
inlinevirtual

Perform additional hanging node procedures for variables that are not interpolated by all nodes. The pressures are stored at the p_nodal_index_nst-th location in each node

Reimplemented from oomph::RefineableElement.

878  {
879  // Set up the hang info for the pressure nodes
880  this->setup_hang_for_value(this->p_nodal_index_nst());
881  } // End of further_setup_hanging_nodes
virtual int p_nodal_index_nst() const
Set the value at which the pressure is stored in the nodes.
Definition: discontinuous_galerkin_equal_order_pressure_spacetime_navier_stokes_elements.h:2221

References oomph::QTaylorHoodSpaceTimeElement< DIM >::p_nodal_index_nst().

◆ further_setup_hanging_nodes() [2/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::further_setup_hanging_nodes ( )
inlinevirtual

Perform additional hanging node procedures for variables that are not interpolated by all nodes. The pressures are stored at the p_nodal_index_nst-th location in each node

Reimplemented from oomph::RefineableElement.

878  {
879  // Set up the hang info for the pressure nodes
880  this->setup_hang_for_value(this->p_nodal_index_nst());
881  } // End of further_setup_hanging_nodes

References oomph::QTaylorHoodSpaceTimeElement< DIM >::p_nodal_index_nst().

◆ further_setup_hanging_nodes() [3/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::further_setup_hanging_nodes ( )
inlinevirtual

Perform additional hanging node procedures for variables that are not interpolated by all nodes. The pressures are stored at the p_nodal_index_nst-th location in each node

Reimplemented from oomph::RefineableElement.

878  {
879  // Set up the hang info for the pressure nodes
880  this->setup_hang_for_value(this->p_nodal_index_nst());
881  } // End of further_setup_hanging_nodes

References oomph::QTaylorHoodSpaceTimeElement< DIM >::p_nodal_index_nst().

◆ get_interpolated_values() [1/6]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::get_interpolated_values ( const unsigned t,
const Vector< double > &  s,
Vector< double > &  values 
)
inlinevirtual

Get the function value u in Vector. Note: Given the generality of the interface (this function is usually called from black-box documentation or interpolation routines), the values Vector sets its own size in here.

Implements oomph::RefineableElement.

850  {
851  // The only time this makes sense (if you can even say that...)
852  if (t == 0)
853  {
854  // Call the other function
855  get_interpolated_values(s, values);
856  }
857  else
858  {
859  // Create an output stream
860  std::ostringstream error_message_stream;
861 
862  // Create an error message
863  error_message_stream << "History values don't make sense in "
864  << "space-time elements!" << std::endl;
865 
866  // Throw an error message
867  throw OomphLibError(error_message_stream.str(),
870  }
871  } // End of get_interpolated_values
void get_interpolated_values(const Vector< double > &s, Vector< double > &values)
Definition: refineable_discontinuous_galerkin_equal_order_pressure_spacetime_navier_stokes_elements.h:825
RealScalar s
Definition: level1_cplx_impl.h:130
t
Definition: plotPSD.py:36
#define OOMPH_EXCEPTION_LOCATION
Definition: oomph_definitions.h:61
#define OOMPH_CURRENT_FUNCTION
Definition: oomph_definitions.h:86

References oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::get_interpolated_values(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, s, and plotPSD::t.

◆ get_interpolated_values() [2/6]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::get_interpolated_values ( const unsigned t,
const Vector< double > &  s,
Vector< double > &  values 
)
inlinevirtual

Get the function value u in Vector. Note: Given the generality of the interface (this function is usually called from black-box documentation or interpolation routines), the values Vector sets its own size in here.

Implements oomph::RefineableElement.

850  {
851  // The only time this makes sense (if you can even say that...)
852  if (t == 0)
853  {
854  // Call the other function
855  get_interpolated_values(s, values);
856  }
857  else
858  {
859  // Create an output stream
860  std::ostringstream error_message_stream;
861 
862  // Create an error message
863  error_message_stream << "History values don't make sense in "
864  << "space-time elements!" << std::endl;
865 
866  // Throw an error message
867  throw OomphLibError(error_message_stream.str(),
870  }
871  } // End of get_interpolated_values

References oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::get_interpolated_values(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, s, and plotPSD::t.

◆ get_interpolated_values() [3/6]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::get_interpolated_values ( const unsigned t,
const Vector< double > &  s,
Vector< double > &  values 
)
inlinevirtual

Get the function value u in Vector. Note: Given the generality of the interface (this function is usually called from black-box documentation or interpolation routines), the values Vector sets its own size in here.

Implements oomph::RefineableElement.

850  {
851  // The only time this makes sense (if you can even say that...)
852  if (t == 0)
853  {
854  // Call the other function
855  get_interpolated_values(s, values);
856  }
857  else
858  {
859  // Create an output stream
860  std::ostringstream error_message_stream;
861 
862  // Create an error message
863  error_message_stream << "History values don't make sense in "
864  << "space-time elements!" << std::endl;
865 
866  // Throw an error message
867  throw OomphLibError(error_message_stream.str(),
870  }
871  } // End of get_interpolated_values

References oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::get_interpolated_values(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, s, and plotPSD::t.

◆ get_interpolated_values() [4/6]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::get_interpolated_values ( const Vector< double > &  s,
Vector< double > &  values 
)
inlinevirtual

Get the function value u in Vector. Note: Given the generality of the interface (this function is usually called from black-box documentation or interpolation routines), the values Vector sets its own size in here.

Reimplemented from oomph::RefineableElement.

827  {
828  // Set size of Vector: u,v,p and initialise to zero
829  values.resize(DIM + 1, 0.0);
830 
831  // Calculate velocities: values[0],...
832  for (unsigned i = 0; i < DIM; i++)
833  {
834  // Calculate the i-th velocity component at local coordinates s
835  values[i] = this->interpolated_u_nst(s, i);
836  }
837 
838  // Calculate the pressure field at local coordinates s
839  values[DIM] = this->interpolated_p_nst(s);
840  } // End of get_interpolated_values
int i
Definition: BiCGSTAB_step_by_step.cpp:9
virtual double interpolated_p_nst(const Vector< double > &s) const
Return FE interpolated pressure at local coordinate s.
Definition: discontinuous_galerkin_equal_order_pressure_spacetime_navier_stokes_elements.h:2037
void interpolated_u_nst(const Vector< double > &s, Vector< double > &velocity) const
Compute vector of FE interpolated velocity u at local coordinate s.
Definition: discontinuous_galerkin_equal_order_pressure_spacetime_navier_stokes_elements.h:1892
#define DIM
Definition: linearised_navier_stokes_elements.h:44

References DIM, i, oomph::SpaceTimeNavierStokesEquations< DIM >::interpolated_p_nst(), and oomph::SpaceTimeNavierStokesEquations< DIM >::interpolated_u_nst().

Referenced by oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::get_interpolated_values().

◆ get_interpolated_values() [5/6]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::get_interpolated_values ( const Vector< double > &  s,
Vector< double > &  values 
)
inlinevirtual

Get the function value u in Vector. Note: Given the generality of the interface (this function is usually called from black-box documentation or interpolation routines), the values Vector sets its own size in here.

Reimplemented from oomph::RefineableElement.

827  {
828  // Set size of Vector: u,v,p and initialise to zero
829  values.resize(DIM + 1, 0.0);
830 
831  // Calculate velocities: values[0],...
832  for (unsigned i = 0; i < DIM; i++)
833  {
834  // Calculate the i-th velocity component at local coordinates s
835  values[i] = this->interpolated_u_nst(s, i);
836  }
837 
838  // Calculate the pressure field at local coordinates s
839  values[DIM] = this->interpolated_p_nst(s);
840  } // End of get_interpolated_values

References DIM, i, oomph::SpaceTimeNavierStokesEquations< DIM >::interpolated_p_nst(), and oomph::SpaceTimeNavierStokesEquations< DIM >::interpolated_u_nst().

◆ get_interpolated_values() [6/6]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::get_interpolated_values ( const Vector< double > &  s,
Vector< double > &  values 
)
inlinevirtual

Get the function value u in Vector. Note: Given the generality of the interface (this function is usually called from black-box documentation or interpolation routines), the values Vector sets its own size in here.

Reimplemented from oomph::RefineableElement.

827  {
828  // Set size of Vector: u,v,p and initialise to zero
829  values.resize(DIM + 1, 0.0);
830 
831  // Calculate velocities: values[0],...
832  for (unsigned i = 0; i < DIM; i++)
833  {
834  // Calculate the i-th velocity component at local coordinates s
835  values[i] = this->interpolated_u_nst(s, i);
836  }
837 
838  // Calculate the pressure field at local coordinates s
839  values[DIM] = this->interpolated_p_nst(s);
840  } // End of get_interpolated_values

References DIM, i, oomph::SpaceTimeNavierStokesEquations< DIM >::interpolated_p_nst(), and oomph::SpaceTimeNavierStokesEquations< DIM >::interpolated_u_nst().

◆ get_interpolating_node_at_local_coordinate() [1/3]

template<unsigned DIM>
Node* oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::get_interpolating_node_at_local_coordinate ( const Vector< double > &  s,
const int value_id 
)
inlinevirtual

The velocity nodes are the same as the geometric nodes. The pressure nodes must be calculated by using the same methods as the geometric nodes, but by recalling that there are only two pressure nodes per edge.

Reimplemented from oomph::RefineableElement.

941  {
942  // If we are calculating pressure nodes
943  if (value_id == DIM)
944  {
945  // Storage for the index of the pressure node
946  unsigned total_index = 0;
947 
948  // The number of nodes along each 1D edge is 2.
949  unsigned nnode_1d = 2;
950 
951  // Storage for the index along each boundary
952  Vector<int> index(DIM + 1);
953 
954  // Loop over the coordinate directions
955  for (unsigned i = 0; i < DIM + 1; i++)
956  {
957  // If we are at the lower limit, the index is zero
958  if (s[i] == -1.0)
959  {
960  // We're at the first node
961  index[i] = 0;
962  }
963  // If we are at the upper limit, the index is the number of nodes
964  // minus 1
965  else if (s[i] == 1.0)
966  {
967  // We're on the last node
968  index[i] = nnode_1d - 1;
969  }
970  // Otherwise, we have to calculate the index in general
971  else
972  {
973  // For uniformly spaced nodes this should produce an integer when
974  // s[i] is associated with a nodal location
975  double float_index = 0.5 * (1.0 + s[i]) * (nnode_1d - 1);
976 
977  // Take the integer part of the float_index value
978  index[i] = int(float_index);
979 
980  // What is the excess. This should be safe because the taking the
981  // integer part rounds down
982  double excess = float_index - index[i];
983 
984  // If the excess is bigger than our tolerance there is no node
986  ((1.0 - excess) > FiniteElement::Node_location_tolerance))
987  {
988  // As there is no node, return null
989  return 0;
990  }
991  } // if (s[i]==-1.0)
992 
993  // Construct the general pressure index from the components.
994  total_index +=
995  index[i] * static_cast<unsigned>(pow(static_cast<float>(nnode_1d),
996  static_cast<int>(i)));
997  } // for (unsigned i=0;i<DIM;i++)
998 
999  // If we've got here we have a node, so let's return a pointer to it
1000  return this->pressure_node_pt(total_index);
1001  }
1002  // Otherwise velocity nodes are the same as pressure nodes
1003  else
1004  {
1005  // Call the regular helper function to find the right node
1006  return this->get_node_at_local_coordinate(s);
1007  } // if (value_id==DIM)
1008  } // End of get_interpolating_node_at_local_coordinate
static const double Node_location_tolerance
Definition: elements.h:1374
virtual Node * get_node_at_local_coordinate(const Vector< double > &s) const
Definition: elements.cc:3882
virtual unsigned nnode_1d() const
Definition: elements.h:2218
Node * pressure_node_pt(const unsigned &n_p)
Pointer to n_p-th pressure node.
Definition: refineable_discontinuous_galerkin_equal_order_pressure_spacetime_navier_stokes_elements.h:885
return int(ret)+1
EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC bfloat16 pow(const bfloat16 &a, const bfloat16 &b)
Definition: BFloat16.h:625

References DIM, oomph::FiniteElement::get_node_at_local_coordinate(), i, int(), oomph::FiniteElement::nnode_1d(), oomph::FiniteElement::Node_location_tolerance, Eigen::bfloat16_impl::pow(), oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt(), and s.

◆ get_interpolating_node_at_local_coordinate() [2/3]

template<unsigned DIM>
Node* oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::get_interpolating_node_at_local_coordinate ( const Vector< double > &  s,
const int value_id 
)
inlinevirtual

The velocity nodes are the same as the geometric nodes. The pressure nodes must be calculated by using the same methods as the geometric nodes, but by recalling that there are only two pressure nodes per edge.

Reimplemented from oomph::RefineableElement.

941  {
942  // If we are calculating pressure nodes
943  if (value_id == DIM)
944  {
945  // Storage for the index of the pressure node
946  unsigned total_index = 0;
947 
948  // The number of nodes along each 1D edge is 2.
949  unsigned nnode_1d = 2;
950 
951  // Storage for the index along each boundary
952  Vector<int> index(DIM + 1);
953 
954  // Loop over the coordinate directions
955  for (unsigned i = 0; i < DIM + 1; i++)
956  {
957  // If we are at the lower limit, the index is zero
958  if (s[i] == -1.0)
959  {
960  // We're at the first node
961  index[i] = 0;
962  }
963  // If we are at the upper limit, the index is the number of nodes
964  // minus 1
965  else if (s[i] == 1.0)
966  {
967  // We're on the last node
968  index[i] = nnode_1d - 1;
969  }
970  // Otherwise, we have to calculate the index in general
971  else
972  {
973  // For uniformly spaced nodes this should produce an integer when
974  // s[i] is associated with a nodal location
975  double float_index = 0.5 * (1.0 + s[i]) * (nnode_1d - 1);
976 
977  // Take the integer part of the float_index value
978  index[i] = int(float_index);
979 
980  // What is the excess. This should be safe because the taking the
981  // integer part rounds down
982  double excess = float_index - index[i];
983 
984  // If the excess is bigger than our tolerance there is no node
986  ((1.0 - excess) > FiniteElement::Node_location_tolerance))
987  {
988  // As there is no node, return null
989  return 0;
990  }
991  } // if (s[i]==-1.0)
992 
993  // Construct the general pressure index from the components.
994  total_index +=
995  index[i] * static_cast<unsigned>(pow(static_cast<float>(nnode_1d),
996  static_cast<int>(i)));
997  } // for (unsigned i=0;i<DIM;i++)
998 
999  // If we've got here we have a node, so let's return a pointer to it
1000  return this->pressure_node_pt(total_index);
1001  }
1002  // Otherwise velocity nodes are the same as pressure nodes
1003  else
1004  {
1005  // Call the regular helper function to find the right node
1006  return this->get_node_at_local_coordinate(s);
1007  } // if (value_id==DIM)
1008  } // End of get_interpolating_node_at_local_coordinate

References DIM, oomph::FiniteElement::get_node_at_local_coordinate(), i, int(), oomph::FiniteElement::nnode_1d(), oomph::FiniteElement::Node_location_tolerance, Eigen::bfloat16_impl::pow(), oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt(), and s.

◆ get_interpolating_node_at_local_coordinate() [3/3]

template<unsigned DIM>
Node* oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::get_interpolating_node_at_local_coordinate ( const Vector< double > &  s,
const int value_id 
)
inlinevirtual

The velocity nodes are the same as the geometric nodes. The pressure nodes must be calculated by using the same methods as the geometric nodes, but by recalling that there are only two pressure nodes per edge.

Reimplemented from oomph::RefineableElement.

941  {
942  // If we are calculating pressure nodes
943  if (value_id == DIM)
944  {
945  // Storage for the index of the pressure node
946  unsigned total_index = 0;
947 
948  // The number of nodes along each 1D edge is 2.
949  unsigned nnode_1d = 2;
950 
951  // Storage for the index along each boundary
952  Vector<int> index(DIM + 1);
953 
954  // Loop over the coordinate directions
955  for (unsigned i = 0; i < DIM + 1; i++)
956  {
957  // If we are at the lower limit, the index is zero
958  if (s[i] == -1.0)
959  {
960  // We're at the first node
961  index[i] = 0;
962  }
963  // If we are at the upper limit, the index is the number of nodes
964  // minus 1
965  else if (s[i] == 1.0)
966  {
967  // We're on the last node
968  index[i] = nnode_1d - 1;
969  }
970  // Otherwise, we have to calculate the index in general
971  else
972  {
973  // For uniformly spaced nodes this should produce an integer when
974  // s[i] is associated with a nodal location
975  double float_index = 0.5 * (1.0 + s[i]) * (nnode_1d - 1);
976 
977  // Take the integer part of the float_index value
978  index[i] = int(float_index);
979 
980  // What is the excess. This should be safe because the taking the
981  // integer part rounds down
982  double excess = float_index - index[i];
983 
984  // If the excess is bigger than our tolerance there is no node
986  ((1.0 - excess) > FiniteElement::Node_location_tolerance))
987  {
988  // As there is no node, return null
989  return 0;
990  }
991  } // if (s[i]==-1.0)
992 
993  // Construct the general pressure index from the components.
994  total_index +=
995  index[i] * static_cast<unsigned>(pow(static_cast<float>(nnode_1d),
996  static_cast<int>(i)));
997  } // for (unsigned i=0;i<DIM;i++)
998 
999  // If we've got here we have a node, so let's return a pointer to it
1000  return this->pressure_node_pt(total_index);
1001  }
1002  // Otherwise velocity nodes are the same as pressure nodes
1003  else
1004  {
1005  // Call the regular helper function to find the right node
1006  return this->get_node_at_local_coordinate(s);
1007  } // if (value_id==DIM)
1008  } // End of get_interpolating_node_at_local_coordinate

References DIM, oomph::FiniteElement::get_node_at_local_coordinate(), i, int(), oomph::FiniteElement::nnode_1d(), oomph::FiniteElement::Node_location_tolerance, Eigen::bfloat16_impl::pow(), oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt(), and s.

◆ identify_load_data() [1/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::identify_load_data ( std::set< std::pair< Data *, unsigned >> &  paired_load_data)
inlinevirtual

Add to the set paired_load_data pairs containing

  • the pointer to a Data object and
  • the index of the value in that Data object

for all values (pressures, velocities) that affect the load computed in the get_load(...) function. (Overloads non-refineable version and takes hanging nodes into account)

Reimplemented from oomph::QTaylorHoodSpaceTimeElement< DIM >.

1092  {
1093  // Allocate space for the velocity component indices
1094  unsigned u_index[DIM];
1095 
1096  // Loop over the velocity components
1097  for (unsigned i = 0; i < DIM; i++)
1098  {
1099  // Get the nodal indices at which the velocities are stored
1100  u_index[i] = this->u_index_nst(i);
1101  }
1102 
1103  // Get the number of nodes in the element
1104  unsigned n_node = this->nnode();
1105 
1106  // Loop over the nodes
1107  for (unsigned n = 0; n < n_node; n++)
1108  {
1109  // Pointer to current node
1110  Node* nod_pt = this->node_pt(n);
1111 
1112  // Check if it's hanging:
1113  if (nod_pt->is_hanging())
1114  {
1115  // It's hanging -- get number of master nodes
1116  unsigned nmaster = nod_pt->hanging_pt()->nmaster();
1117 
1118  // Loop over master nodes
1119  for (unsigned j = 0; j < nmaster; j++)
1120  {
1121  // Create a pointer to the j-th master node
1122  Node* master_nod_pt = nod_pt->hanging_pt()->master_node_pt(j);
1123 
1124  // Loop over the velocity components and add a pointer to their data
1125  // and indices to the vectors
1126  for (unsigned i = 0; i < DIM; i++)
1127  {
1128  // Create a pair and add it to the storage
1129  paired_load_data.insert(
1130  std::make_pair(master_nod_pt, u_index[i]));
1131  }
1132  } // for (unsigned j=0;j<nmaster;j++)
1133  }
1134  // If the node is not hanging
1135  else
1136  {
1137  // Loop over the velocity components and add pointer to their data
1138  // and indices to the vectors
1139  for (unsigned i = 0; i < DIM; i++)
1140  {
1141  // Create a pair and add it to the storage
1142  paired_load_data.insert(
1143  std::make_pair(this->node_pt(n), u_index[i]));
1144  }
1145  } // if (nod_pt->is_hanging())
1146  } // for (unsigned n=0;n<n_node;n++)
1147 
1148  // Get the nodal index at which the pressure is stored
1149  int p_index = this->p_nodal_index_nst();
1150 
1151  // Get the number of pressure degrees of freedom
1152  unsigned n_pres = this->npres_nst();
1153 
1154  // Loop over the pressure data
1155  for (unsigned l = 0; l < n_pres; l++)
1156  {
1157  // Get the pointer to the l-th pressure node
1158  Node* pres_node_pt = this->pressure_node_pt(l);
1159 
1160  // Check if the pressure dof is hanging
1161  if (pres_node_pt->is_hanging(p_index))
1162  {
1163  // Get the pointer to the hang info object (pressure is stored
1164  // as p_index-th nodal dof).
1165  HangInfo* hang_info_pt = pres_node_pt->hanging_pt(p_index);
1166 
1167  // Get number of pressure master nodes (pressure is stored
1168  unsigned nmaster = hang_info_pt->nmaster();
1169 
1170  // Loop over pressure master nodes
1171  for (unsigned m = 0; m < nmaster; m++)
1172  {
1173  // The p_index-th entry in each nodal data is the pressure, which
1174  // affects the traction
1175  paired_load_data.insert(
1176  std::make_pair(hang_info_pt->master_node_pt(m), p_index));
1177  }
1178  }
1179  // If the pressure dof is not hanging
1180  else
1181  {
1182  // The p_index-th entry in each nodal data is the pressure, which
1183  // affects the traction
1184  paired_load_data.insert(std::make_pair(pres_node_pt, p_index));
1185  } // if (pres_node_pt->is_hanging(p_index))
1186  } // for (unsigned l=0;l<n_pres;l++)
1187  } // End of identify_load_data
const unsigned n
Definition: CG3DPackingUnitTest.cpp:11
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210
unsigned npres_nst() const
Return number of pressure values.
Definition: discontinuous_galerkin_equal_order_pressure_spacetime_navier_stokes_elements.h:2255
virtual unsigned u_index_nst(const unsigned &i) const
Definition: discontinuous_galerkin_equal_order_pressure_spacetime_navier_stokes_elements.h:1102
int * m
Definition: level2_cplx_impl.h:294
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References DIM, oomph::Node::hanging_pt(), i, oomph::Node::is_hanging(), j, m, oomph::HangInfo::master_node_pt(), n, oomph::HangInfo::nmaster(), oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), oomph::QTaylorHoodSpaceTimeElement< DIM >::npres_nst(), oomph::QTaylorHoodSpaceTimeElement< DIM >::p_nodal_index_nst(), oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt(), and oomph::SpaceTimeNavierStokesEquations< DIM >::u_index_nst().

◆ identify_load_data() [2/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::identify_load_data ( std::set< std::pair< Data *, unsigned >> &  paired_load_data)
inlinevirtual

Add to the set paired_load_data pairs containing

  • the pointer to a Data object and
  • the index of the value in that Data object

for all values (pressures, velocities) that affect the load computed in the get_load(...) function. (Overloads non-refineable version and takes hanging nodes into account)

Reimplemented from oomph::QTaylorHoodSpaceTimeElement< DIM >.

1092  {
1093  // Allocate space for the velocity component indices
1094  unsigned u_index[DIM];
1095 
1096  // Loop over the velocity components
1097  for (unsigned i = 0; i < DIM; i++)
1098  {
1099  // Get the nodal indices at which the velocities are stored
1100  u_index[i] = this->u_index_nst(i);
1101  }
1102 
1103  // Get the number of nodes in the element
1104  unsigned n_node = this->nnode();
1105 
1106  // Loop over the nodes
1107  for (unsigned n = 0; n < n_node; n++)
1108  {
1109  // Pointer to current node
1110  Node* nod_pt = this->node_pt(n);
1111 
1112  // Check if it's hanging:
1113  if (nod_pt->is_hanging())
1114  {
1115  // It's hanging -- get number of master nodes
1116  unsigned nmaster = nod_pt->hanging_pt()->nmaster();
1117 
1118  // Loop over master nodes
1119  for (unsigned j = 0; j < nmaster; j++)
1120  {
1121  // Create a pointer to the j-th master node
1122  Node* master_nod_pt = nod_pt->hanging_pt()->master_node_pt(j);
1123 
1124  // Loop over the velocity components and add a pointer to their data
1125  // and indices to the vectors
1126  for (unsigned i = 0; i < DIM; i++)
1127  {
1128  // Create a pair and add it to the storage
1129  paired_load_data.insert(
1130  std::make_pair(master_nod_pt, u_index[i]));
1131  }
1132  } // for (unsigned j=0;j<nmaster;j++)
1133  }
1134  // If the node is not hanging
1135  else
1136  {
1137  // Loop over the velocity components and add pointer to their data
1138  // and indices to the vectors
1139  for (unsigned i = 0; i < DIM; i++)
1140  {
1141  // Create a pair and add it to the storage
1142  paired_load_data.insert(
1143  std::make_pair(this->node_pt(n), u_index[i]));
1144  }
1145  } // if (nod_pt->is_hanging())
1146  } // for (unsigned n=0;n<n_node;n++)
1147 
1148  // Get the nodal index at which the pressure is stored
1149  int p_index = this->p_nodal_index_nst();
1150 
1151  // Get the number of pressure degrees of freedom
1152  unsigned n_pres = this->npres_nst();
1153 
1154  // Loop over the pressure data
1155  for (unsigned l = 0; l < n_pres; l++)
1156  {
1157  // Get the pointer to the l-th pressure node
1158  Node* pres_node_pt = this->pressure_node_pt(l);
1159 
1160  // Check if the pressure dof is hanging
1161  if (pres_node_pt->is_hanging(p_index))
1162  {
1163  // Get the pointer to the hang info object (pressure is stored
1164  // as p_index-th nodal dof).
1165  HangInfo* hang_info_pt = pres_node_pt->hanging_pt(p_index);
1166 
1167  // Get number of pressure master nodes (pressure is stored
1168  unsigned nmaster = hang_info_pt->nmaster();
1169 
1170  // Loop over pressure master nodes
1171  for (unsigned m = 0; m < nmaster; m++)
1172  {
1173  // The p_index-th entry in each nodal data is the pressure, which
1174  // affects the traction
1175  paired_load_data.insert(
1176  std::make_pair(hang_info_pt->master_node_pt(m), p_index));
1177  }
1178  }
1179  // If the pressure dof is not hanging
1180  else
1181  {
1182  // The p_index-th entry in each nodal data is the pressure, which
1183  // affects the traction
1184  paired_load_data.insert(std::make_pair(pres_node_pt, p_index));
1185  } // if (pres_node_pt->is_hanging(p_index))
1186  } // for (unsigned l=0;l<n_pres;l++)
1187  } // End of identify_load_data

References DIM, oomph::Node::hanging_pt(), i, oomph::Node::is_hanging(), j, m, oomph::HangInfo::master_node_pt(), n, oomph::HangInfo::nmaster(), oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), oomph::QTaylorHoodSpaceTimeElement< DIM >::npres_nst(), oomph::QTaylorHoodSpaceTimeElement< DIM >::p_nodal_index_nst(), oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt(), and oomph::SpaceTimeNavierStokesEquations< DIM >::u_index_nst().

◆ identify_load_data() [3/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::identify_load_data ( std::set< std::pair< Data *, unsigned >> &  paired_load_data)
inlinevirtual

Add to the set paired_load_data pairs containing

  • the pointer to a Data object and
  • the index of the value in that Data object

for all values (pressures, velocities) that affect the load computed in the get_load(...) function. (Overloads non-refineable version and takes hanging nodes into account)

Reimplemented from oomph::QTaylorHoodSpaceTimeElement< DIM >.

1092  {
1093  // Allocate space for the velocity component indices
1094  unsigned u_index[DIM];
1095 
1096  // Loop over the velocity components
1097  for (unsigned i = 0; i < DIM; i++)
1098  {
1099  // Get the nodal indices at which the velocities are stored
1100  u_index[i] = this->u_index_nst(i);
1101  }
1102 
1103  // Get the number of nodes in the element
1104  unsigned n_node = this->nnode();
1105 
1106  // Loop over the nodes
1107  for (unsigned n = 0; n < n_node; n++)
1108  {
1109  // Pointer to current node
1110  Node* nod_pt = this->node_pt(n);
1111 
1112  // Check if it's hanging:
1113  if (nod_pt->is_hanging())
1114  {
1115  // It's hanging -- get number of master nodes
1116  unsigned nmaster = nod_pt->hanging_pt()->nmaster();
1117 
1118  // Loop over master nodes
1119  for (unsigned j = 0; j < nmaster; j++)
1120  {
1121  // Create a pointer to the j-th master node
1122  Node* master_nod_pt = nod_pt->hanging_pt()->master_node_pt(j);
1123 
1124  // Loop over the velocity components and add a pointer to their data
1125  // and indices to the vectors
1126  for (unsigned i = 0; i < DIM; i++)
1127  {
1128  // Create a pair and add it to the storage
1129  paired_load_data.insert(
1130  std::make_pair(master_nod_pt, u_index[i]));
1131  }
1132  } // for (unsigned j=0;j<nmaster;j++)
1133  }
1134  // If the node is not hanging
1135  else
1136  {
1137  // Loop over the velocity components and add pointer to their data
1138  // and indices to the vectors
1139  for (unsigned i = 0; i < DIM; i++)
1140  {
1141  // Create a pair and add it to the storage
1142  paired_load_data.insert(
1143  std::make_pair(this->node_pt(n), u_index[i]));
1144  }
1145  } // if (nod_pt->is_hanging())
1146  } // for (unsigned n=0;n<n_node;n++)
1147 
1148  // Get the nodal index at which the pressure is stored
1149  int p_index = this->p_nodal_index_nst();
1150 
1151  // Get the number of pressure degrees of freedom
1152  unsigned n_pres = this->npres_nst();
1153 
1154  // Loop over the pressure data
1155  for (unsigned l = 0; l < n_pres; l++)
1156  {
1157  // Get the pointer to the l-th pressure node
1158  Node* pres_node_pt = this->pressure_node_pt(l);
1159 
1160  // Check if the pressure dof is hanging
1161  if (pres_node_pt->is_hanging(p_index))
1162  {
1163  // Get the pointer to the hang info object (pressure is stored
1164  // as p_index-th nodal dof).
1165  HangInfo* hang_info_pt = pres_node_pt->hanging_pt(p_index);
1166 
1167  // Get number of pressure master nodes (pressure is stored
1168  unsigned nmaster = hang_info_pt->nmaster();
1169 
1170  // Loop over pressure master nodes
1171  for (unsigned m = 0; m < nmaster; m++)
1172  {
1173  // The p_index-th entry in each nodal data is the pressure, which
1174  // affects the traction
1175  paired_load_data.insert(
1176  std::make_pair(hang_info_pt->master_node_pt(m), p_index));
1177  }
1178  }
1179  // If the pressure dof is not hanging
1180  else
1181  {
1182  // The p_index-th entry in each nodal data is the pressure, which
1183  // affects the traction
1184  paired_load_data.insert(std::make_pair(pres_node_pt, p_index));
1185  } // if (pres_node_pt->is_hanging(p_index))
1186  } // for (unsigned l=0;l<n_pres;l++)
1187  } // End of identify_load_data

References DIM, oomph::Node::hanging_pt(), i, oomph::Node::is_hanging(), j, m, oomph::HangInfo::master_node_pt(), n, oomph::HangInfo::nmaster(), oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), oomph::QTaylorHoodSpaceTimeElement< DIM >::npres_nst(), oomph::QTaylorHoodSpaceTimeElement< DIM >::p_nodal_index_nst(), oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt(), and oomph::SpaceTimeNavierStokesEquations< DIM >::u_index_nst().

◆ interpolating_basis() [1/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::interpolating_basis ( const Vector< double > &  s,
Shape psi,
const int value_id 
) const
inlinevirtual

The basis interpolating the pressure is given by pshape(). / The basis interpolating the velocity is shape().

Reimplemented from oomph::RefineableElement.

1054  {
1055  // If we want the pressure basis functions
1056  if (value_id == DIM)
1057  {
1058  // Call the pressure shape function
1059  return this->pshape_nst(s, psi);
1060  }
1061  // If we want the velocity basis functions
1062  else
1063  {
1064  // Call the velocity shape function
1065  return this->shape(s, psi);
1066  }
1067  } // End of interpolating_basis
virtual void shape(const Vector< double > &s, Shape &psi) const =0
void pshape_nst(const Vector< double > &s, Shape &psi) const
Pressure shape functions at local coordinate s.

References DIM, oomph::QTaylorHoodSpaceTimeElement< DIM >::pshape_nst(), and oomph::FiniteElement::shape().

◆ interpolating_basis() [2/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::interpolating_basis ( const Vector< double > &  s,
Shape psi,
const int value_id 
) const
inlinevirtual

The basis interpolating the pressure is given by pshape(). / The basis interpolating the velocity is shape().

Reimplemented from oomph::RefineableElement.

1054  {
1055  // If we want the pressure basis functions
1056  if (value_id == DIM)
1057  {
1058  // Call the pressure shape function
1059  return this->pshape_nst(s, psi);
1060  }
1061  // If we want the velocity basis functions
1062  else
1063  {
1064  // Call the velocity shape function
1065  return this->shape(s, psi);
1066  }
1067  } // End of interpolating_basis

References DIM, oomph::QTaylorHoodSpaceTimeElement< DIM >::pshape_nst(), and oomph::FiniteElement::shape().

◆ interpolating_basis() [3/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::interpolating_basis ( const Vector< double > &  s,
Shape psi,
const int value_id 
) const
inlinevirtual

The basis interpolating the pressure is given by pshape(). / The basis interpolating the velocity is shape().

Reimplemented from oomph::RefineableElement.

1054  {
1055  // If we want the pressure basis functions
1056  if (value_id == DIM)
1057  {
1058  // Call the pressure shape function
1059  return this->pshape_nst(s, psi);
1060  }
1061  // If we want the velocity basis functions
1062  else
1063  {
1064  // Call the velocity shape function
1065  return this->shape(s, psi);
1066  }
1067  } // End of interpolating_basis

References DIM, oomph::QTaylorHoodSpaceTimeElement< DIM >::pshape_nst(), and oomph::FiniteElement::shape().

◆ interpolating_node_pt() [1/3]

template<unsigned DIM>
Node* oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::interpolating_node_pt ( const unsigned n,
const int value_id 
)
inlinevirtual

The velocities are isoparametric and so the "nodes" interpolating the velocities are the geometric nodes. The pressure "nodes" are a subset of the nodes, so when value_id==DIM, the n-th pressure node is returned.

Reimplemented from oomph::RefineableElement.

898  {
899  // The only different nodes are the pressure nodes
900  if (value_id == DIM)
901  {
902  // Return a pointer to the n-th pressure node
903  return this->pressure_node_pt(n);
904  }
905  // The other variables are interpolated via the usual nodes
906  else
907  {
908  // Return a pointer to the n-th regular node
909  return this->node_pt(n);
910  }
911  } // End of interpolating_node_pt

References DIM, oomph::FiniteElement::node_pt(), and oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt().

◆ interpolating_node_pt() [2/3]

template<unsigned DIM>
Node* oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::interpolating_node_pt ( const unsigned n,
const int value_id 
)
inlinevirtual

The velocities are isoparametric and so the "nodes" interpolating the velocities are the geometric nodes. The pressure "nodes" are a subset of the nodes, so when value_id==DIM, the n-th pressure node is returned.

Reimplemented from oomph::RefineableElement.

898  {
899  // The only different nodes are the pressure nodes
900  if (value_id == DIM)
901  {
902  // Return a pointer to the n-th pressure node
903  return this->pressure_node_pt(n);
904  }
905  // The other variables are interpolated via the usual nodes
906  else
907  {
908  // Return a pointer to the n-th regular node
909  return this->node_pt(n);
910  }
911  } // End of interpolating_node_pt

References DIM, oomph::FiniteElement::node_pt(), and oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt().

◆ interpolating_node_pt() [3/3]

template<unsigned DIM>
Node* oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::interpolating_node_pt ( const unsigned n,
const int value_id 
)
inlinevirtual

The velocities are isoparametric and so the "nodes" interpolating the velocities are the geometric nodes. The pressure "nodes" are a subset of the nodes, so when value_id==DIM, the n-th pressure node is returned.

Reimplemented from oomph::RefineableElement.

898  {
899  // The only different nodes are the pressure nodes
900  if (value_id == DIM)
901  {
902  // Return a pointer to the n-th pressure node
903  return this->pressure_node_pt(n);
904  }
905  // The other variables are interpolated via the usual nodes
906  else
907  {
908  // Return a pointer to the n-th regular node
909  return this->node_pt(n);
910  }
911  } // End of interpolating_node_pt

References DIM, oomph::FiniteElement::node_pt(), and oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt().

◆ local_one_d_fraction_of_interpolating_node() [1/3]

template<unsigned DIM>
double oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::local_one_d_fraction_of_interpolating_node ( const unsigned n_1d,
const unsigned i,
const int value_id 
)
inlinevirtual

The pressure nodes are the corner nodes, so when n_value==DIM, the fraction is the same as the 1D node number, 0 or 1.

Reimplemented from oomph::RefineableElement.

919  {
920  // If we're dealing with a pressure node
921  if (value_id == DIM)
922  {
923  // The pressure nodes are just located on the boundaries at 0 or 1
924  return double(n_1d);
925  }
926  // Otherwise we're dealing with a velocity node
927  else
928  {
929  // The velocity nodes are the same as the geometric ones
930  return this->local_one_d_fraction_of_node(n_1d, i);
931  }
932  } // End of local_one_d_fraction_of_interpolating_node
virtual double local_one_d_fraction_of_node(const unsigned &n1d, const unsigned &i)
Definition: elements.h:1858

References DIM, i, and oomph::FiniteElement::local_one_d_fraction_of_node().

◆ local_one_d_fraction_of_interpolating_node() [2/3]

template<unsigned DIM>
double oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::local_one_d_fraction_of_interpolating_node ( const unsigned n_1d,
const unsigned i,
const int value_id 
)
inlinevirtual

The pressure nodes are the corner nodes, so when n_value==DIM, the fraction is the same as the 1D node number, 0 or 1.

Reimplemented from oomph::RefineableElement.

919  {
920  // If we're dealing with a pressure node
921  if (value_id == DIM)
922  {
923  // The pressure nodes are just located on the boundaries at 0 or 1
924  return double(n_1d);
925  }
926  // Otherwise we're dealing with a velocity node
927  else
928  {
929  // The velocity nodes are the same as the geometric ones
930  return this->local_one_d_fraction_of_node(n_1d, i);
931  }
932  } // End of local_one_d_fraction_of_interpolating_node

References DIM, i, and oomph::FiniteElement::local_one_d_fraction_of_node().

◆ local_one_d_fraction_of_interpolating_node() [3/3]

template<unsigned DIM>
double oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::local_one_d_fraction_of_interpolating_node ( const unsigned n_1d,
const unsigned i,
const int value_id 
)
inlinevirtual

The pressure nodes are the corner nodes, so when n_value==DIM, the fraction is the same as the 1D node number, 0 or 1.

Reimplemented from oomph::RefineableElement.

919  {
920  // If we're dealing with a pressure node
921  if (value_id == DIM)
922  {
923  // The pressure nodes are just located on the boundaries at 0 or 1
924  return double(n_1d);
925  }
926  // Otherwise we're dealing with a velocity node
927  else
928  {
929  // The velocity nodes are the same as the geometric ones
930  return this->local_one_d_fraction_of_node(n_1d, i);
931  }
932  } // End of local_one_d_fraction_of_interpolating_node

References DIM, i, and oomph::FiniteElement::local_one_d_fraction_of_node().

◆ ncont_interpolated_values() [1/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::ncont_interpolated_values ( ) const
inlinevirtual

Number of continuously interpolated values.

Implements oomph::RefineableElement.

786  {
787  // There are DIM velocities and 1 pressure component
788  return DIM + 1;
789  } // End of ncont_interpolated_values

References DIM.

◆ ncont_interpolated_values() [2/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::ncont_interpolated_values ( ) const
inlinevirtual

Number of continuously interpolated values.

Implements oomph::RefineableElement.

786  {
787  // There are DIM velocities and 1 pressure component
788  return DIM + 1;
789  } // End of ncont_interpolated_values

References DIM.

◆ ncont_interpolated_values() [3/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::ncont_interpolated_values ( ) const
inlinevirtual

Number of continuously interpolated values.

Implements oomph::RefineableElement.

786  {
787  // There are DIM velocities and 1 pressure component
788  return DIM + 1;
789  } // End of ncont_interpolated_values

References DIM.

◆ ninterpolating_node() [1/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::ninterpolating_node ( const int value_id)
inlinevirtual

The number of pressure nodes is 2^DIM. The number of velocity nodes is the same as the number of geometric nodes.

Reimplemented from oomph::RefineableElement.

1033  {
1034  // If we want the pressure basis functions
1035  if (value_id == DIM)
1036  {
1037  // There are 2^{DIM+1} pressure dofs (also interpolating in time)
1038  return static_cast<unsigned>(pow(2.0, static_cast<int>(DIM + 1)));
1039  }
1040  // If we want the velocity basis functions
1041  else
1042  {
1043  // Every node is a velocity interpolating node
1044  return this->nnode();
1045  }
1046  } // End if ninterpolating_node

References DIM, oomph::FiniteElement::nnode(), and Eigen::bfloat16_impl::pow().

◆ ninterpolating_node() [2/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::ninterpolating_node ( const int value_id)
inlinevirtual

The number of pressure nodes is 2^DIM. The number of velocity nodes is the same as the number of geometric nodes.

Reimplemented from oomph::RefineableElement.

1033  {
1034  // If we want the pressure basis functions
1035  if (value_id == DIM)
1036  {
1037  // There are 2^{DIM+1} pressure dofs (also interpolating in time)
1038  return static_cast<unsigned>(pow(2.0, static_cast<int>(DIM + 1)));
1039  }
1040  // If we want the velocity basis functions
1041  else
1042  {
1043  // Every node is a velocity interpolating node
1044  return this->nnode();
1045  }
1046  } // End if ninterpolating_node

References DIM, oomph::FiniteElement::nnode(), and Eigen::bfloat16_impl::pow().

◆ ninterpolating_node() [3/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::ninterpolating_node ( const int value_id)
inlinevirtual

The number of pressure nodes is 2^DIM. The number of velocity nodes is the same as the number of geometric nodes.

Reimplemented from oomph::RefineableElement.

1033  {
1034  // If we want the pressure basis functions
1035  if (value_id == DIM)
1036  {
1037  // There are 2^{DIM+1} pressure dofs (also interpolating in time)
1038  return static_cast<unsigned>(pow(2.0, static_cast<int>(DIM + 1)));
1039  }
1040  // If we want the velocity basis functions
1041  else
1042  {
1043  // Every node is a velocity interpolating node
1044  return this->nnode();
1045  }
1046  } // End if ninterpolating_node

References DIM, oomph::FiniteElement::nnode(), and Eigen::bfloat16_impl::pow().

◆ ninterpolating_node_1d() [1/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::ninterpolating_node_1d ( const int value_id)
inlinevirtual

The number of 1D pressure nodes is 2, the number of 1D velocity nodes is the same as the number of 1D geometric nodes.

Reimplemented from oomph::RefineableElement.

1014  {
1015  // If we're dealing with the pressure nodes
1016  if (value_id == DIM)
1017  {
1018  // Using linear interpolation for the pressure
1019  return 2;
1020  }
1021  // If we're dealing with the velocity nodes
1022  else
1023  {
1024  // Every node is a velocity interpolating node
1025  return this->nnode_1d();
1026  }
1027  } // End of ninterpolating_node_1d

References DIM, and oomph::FiniteElement::nnode_1d().

◆ ninterpolating_node_1d() [2/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::ninterpolating_node_1d ( const int value_id)
inlinevirtual

The number of 1D pressure nodes is 2, the number of 1D velocity nodes is the same as the number of 1D geometric nodes.

Reimplemented from oomph::RefineableElement.

1014  {
1015  // If we're dealing with the pressure nodes
1016  if (value_id == DIM)
1017  {
1018  // Using linear interpolation for the pressure
1019  return 2;
1020  }
1021  // If we're dealing with the velocity nodes
1022  else
1023  {
1024  // Every node is a velocity interpolating node
1025  return this->nnode_1d();
1026  }
1027  } // End of ninterpolating_node_1d

References DIM, and oomph::FiniteElement::nnode_1d().

◆ ninterpolating_node_1d() [3/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::ninterpolating_node_1d ( const int value_id)
inlinevirtual

The number of 1D pressure nodes is 2, the number of 1D velocity nodes is the same as the number of 1D geometric nodes.

Reimplemented from oomph::RefineableElement.

1014  {
1015  // If we're dealing with the pressure nodes
1016  if (value_id == DIM)
1017  {
1018  // Using linear interpolation for the pressure
1019  return 2;
1020  }
1021  // If we're dealing with the velocity nodes
1022  else
1023  {
1024  // Every node is a velocity interpolating node
1025  return this->nnode_1d();
1026  }
1027  } // End of ninterpolating_node_1d

References DIM, and oomph::FiniteElement::nnode_1d().

◆ nrecovery_order() [1/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::nrecovery_order ( )
inlinevirtual

Order of recovery shape functions for Z2 error estimation: Same order as shape functions.

Implements oomph::ElementWithZ2ErrorEstimator.

799  {
800  // Using quadratic interpolation
801  return 2;
802  } // End of nrecovery_order

◆ nrecovery_order() [2/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::nrecovery_order ( )
inlinevirtual

Order of recovery shape functions for Z2 error estimation: Same order as shape functions.

Implements oomph::ElementWithZ2ErrorEstimator.

799  {
800  // Using quadratic interpolation
801  return 2;
802  } // End of nrecovery_order

◆ nrecovery_order() [3/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::nrecovery_order ( )
inlinevirtual

Order of recovery shape functions for Z2 error estimation: Same order as shape functions.

Implements oomph::ElementWithZ2ErrorEstimator.

799  {
800  // Using quadratic interpolation
801  return 2;
802  } // End of nrecovery_order

◆ nvertex_node() [1/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::nvertex_node ( ) const
inlinevirtual

Number of vertex nodes in the element.

Implements oomph::ElementWithZ2ErrorEstimator.

807  {
808  // Call the base class implementation of the function
810  } // End of nvertex_node
virtual unsigned nvertex_node() const
Definition: elements.h:2491

References oomph::FiniteElement::nvertex_node().

◆ nvertex_node() [2/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::nvertex_node ( ) const
inlinevirtual

Number of vertex nodes in the element.

Implements oomph::ElementWithZ2ErrorEstimator.

807  {
808  // Call the base class implementation of the function
810  } // End of nvertex_node

References oomph::FiniteElement::nvertex_node().

◆ nvertex_node() [3/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::nvertex_node ( ) const
inlinevirtual

Number of vertex nodes in the element.

Implements oomph::ElementWithZ2ErrorEstimator.

807  {
808  // Call the base class implementation of the function
810  } // End of nvertex_node

References oomph::FiniteElement::nvertex_node().

◆ pin_elemental_redundant_nodal_pressure_dofs() [1/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pin_elemental_redundant_nodal_pressure_dofs ( )
inlineprivatevirtual

Pin all nodal pressure dofs that are not required.

Reimplemented from oomph::RefineableSpaceTimeNavierStokesEquations< DIM >.

1210  {
1211  // Find the pressure index
1212  int p_index = this->p_nodal_index_nst();
1213 
1214  // Loop over all nodes
1215  unsigned n_node = this->nnode();
1216 
1217  // Loop over all nodes and pin all the nodal pressures
1218  for (unsigned n = 0; n < n_node; n++)
1219  {
1220  // Pin the p_index-th dof at node n
1221  this->node_pt(n)->pin(p_index);
1222  }
1223 
1224  // Loop over all actual pressure nodes and unpin if they're not hanging
1225  unsigned n_pres = this->npres_nst();
1226 
1227  // Loop over all pressure nodes
1228  for (unsigned l = 0; l < n_pres; l++)
1229  {
1230  // Get a pointer to the l-th pressure node
1231  Node* nod_pt = this->pressure_node_pt(l);
1232 
1233  // If the node isn't hanging
1234  if (!nod_pt->is_hanging(p_index))
1235  {
1236  // Unpin the p_index-th dof at this node
1237  nod_pt->unpin(p_index);
1238  }
1239  } // for (unsigned l=0;l<n_pres;l++)
1240  } // End of pin_elemental_redundant_nodal_pressure_dofs
void pin(const unsigned &i)
Pin the i-th stored variable.
Definition: nodes.h:385

References oomph::Node::is_hanging(), n, oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), oomph::QTaylorHoodSpaceTimeElement< DIM >::npres_nst(), oomph::QTaylorHoodSpaceTimeElement< DIM >::p_nodal_index_nst(), oomph::Data::pin(), oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt(), and oomph::Data::unpin().

◆ pin_elemental_redundant_nodal_pressure_dofs() [2/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pin_elemental_redundant_nodal_pressure_dofs ( )
inlineprivatevirtual

Pin all nodal pressure dofs that are not required.

Reimplemented from oomph::RefineableSpaceTimeNavierStokesEquations< DIM >.

1210  {
1211  // Find the pressure index
1212  int p_index = this->p_nodal_index_nst();
1213 
1214  // Loop over all nodes
1215  unsigned n_node = this->nnode();
1216 
1217  // Loop over all nodes and pin all the nodal pressures
1218  for (unsigned n = 0; n < n_node; n++)
1219  {
1220  // Pin the p_index-th dof at node n
1221  this->node_pt(n)->pin(p_index);
1222  }
1223 
1224  // Loop over all actual pressure nodes and unpin if they're not hanging
1225  unsigned n_pres = this->npres_nst();
1226 
1227  // Loop over all pressure nodes
1228  for (unsigned l = 0; l < n_pres; l++)
1229  {
1230  // Get a pointer to the l-th pressure node
1231  Node* nod_pt = this->pressure_node_pt(l);
1232 
1233  // If the node isn't hanging
1234  if (!nod_pt->is_hanging(p_index))
1235  {
1236  // Unpin the p_index-th dof at this node
1237  nod_pt->unpin(p_index);
1238  }
1239  } // for (unsigned l=0;l<n_pres;l++)
1240  } // End of pin_elemental_redundant_nodal_pressure_dofs

References oomph::Node::is_hanging(), n, oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), oomph::QTaylorHoodSpaceTimeElement< DIM >::npres_nst(), oomph::QTaylorHoodSpaceTimeElement< DIM >::p_nodal_index_nst(), oomph::Data::pin(), oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt(), and oomph::Data::unpin().

◆ pin_elemental_redundant_nodal_pressure_dofs() [3/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pin_elemental_redundant_nodal_pressure_dofs ( )
inlineprivatevirtual

Pin all nodal pressure dofs that are not required.

Reimplemented from oomph::RefineableSpaceTimeNavierStokesEquations< DIM >.

1210  {
1211  // Find the pressure index
1212  int p_index = this->p_nodal_index_nst();
1213 
1214  // Loop over all nodes
1215  unsigned n_node = this->nnode();
1216 
1217  // Loop over all nodes and pin all the nodal pressures
1218  for (unsigned n = 0; n < n_node; n++)
1219  {
1220  // Pin the p_index-th dof at node n
1221  this->node_pt(n)->pin(p_index);
1222  }
1223 
1224  // Loop over all actual pressure nodes and unpin if they're not hanging
1225  unsigned n_pres = this->npres_nst();
1226 
1227  // Loop over all pressure nodes
1228  for (unsigned l = 0; l < n_pres; l++)
1229  {
1230  // Get a pointer to the l-th pressure node
1231  Node* nod_pt = this->pressure_node_pt(l);
1232 
1233  // If the node isn't hanging
1234  if (!nod_pt->is_hanging(p_index))
1235  {
1236  // Unpin the p_index-th dof at this node
1237  nod_pt->unpin(p_index);
1238  }
1239  } // for (unsigned l=0;l<n_pres;l++)
1240  } // End of pin_elemental_redundant_nodal_pressure_dofs

References oomph::Node::is_hanging(), n, oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), oomph::QTaylorHoodSpaceTimeElement< DIM >::npres_nst(), oomph::QTaylorHoodSpaceTimeElement< DIM >::p_nodal_index_nst(), oomph::Data::pin(), oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt(), and oomph::Data::unpin().

◆ pressure_node_pt() [1/3]

template<unsigned DIM>
Node* oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt ( const unsigned n_p)
inlinevirtual

Pointer to n_p-th pressure node.

Reimplemented from oomph::RefineableSpaceTimeNavierStokesEquations< DIM >.

886  {
887  // Return a pointer to the n_p-th pressure node
888  return this->node_pt(this->Pconv[n_p]);
889  } // End of pressure node_pt
static const unsigned Pconv[]
Definition: discontinuous_galerkin_equal_order_pressure_spacetime_navier_stokes_elements.h:2128

References oomph::FiniteElement::node_pt(), and oomph::QTaylorHoodSpaceTimeElement< DIM >::Pconv.

Referenced by oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::get_interpolating_node_at_local_coordinate(), oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::identify_load_data(), oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::interpolating_node_pt(), and oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pin_elemental_redundant_nodal_pressure_dofs().

◆ pressure_node_pt() [2/3]

template<unsigned DIM>
Node* oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt ( const unsigned n_p)
inlinevirtual

Pointer to n_p-th pressure node.

Reimplemented from oomph::RefineableSpaceTimeNavierStokesEquations< DIM >.

886  {
887  // Return a pointer to the n_p-th pressure node
888  return this->node_pt(this->Pconv[n_p]);
889  } // End of pressure node_pt

References oomph::FiniteElement::node_pt(), and oomph::QTaylorHoodSpaceTimeElement< DIM >::Pconv.

◆ pressure_node_pt() [3/3]

template<unsigned DIM>
Node* oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::pressure_node_pt ( const unsigned n_p)
inlinevirtual

Pointer to n_p-th pressure node.

Reimplemented from oomph::RefineableSpaceTimeNavierStokesEquations< DIM >.

886  {
887  // Return a pointer to the n_p-th pressure node
888  return this->node_pt(this->Pconv[n_p]);
889  } // End of pressure node_pt

References oomph::FiniteElement::node_pt(), and oomph::QTaylorHoodSpaceTimeElement< DIM >::Pconv.

◆ rebuild_from_sons() [1/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::rebuild_from_sons ( Mesh *&  mesh_pt)
inlinevirtual

Rebuild from sons: empty.

Implements oomph::RefineableElement.

793 {}

◆ rebuild_from_sons() [2/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::rebuild_from_sons ( Mesh *&  mesh_pt)
inlinevirtual

Rebuild from sons: empty.

Implements oomph::RefineableElement.

793 {}

◆ rebuild_from_sons() [3/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::rebuild_from_sons ( Mesh *&  mesh_pt)
inlinevirtual

Rebuild from sons: empty.

Implements oomph::RefineableElement.

793 {}

◆ required_nvalue() [1/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::required_nvalue ( const unsigned n) const
inlinevirtual

Number of values required at local node n. In order to simplify matters, we allocate storage for pressure variables at all the nodes and then pin those that are not used.

Reimplemented from oomph::QTaylorHoodSpaceTimeElement< DIM >.

778  {
779  // There are DIM velocity components and 1 pressure component
780  return DIM + 1;
781  } // End of required_nvalue

References DIM.

◆ required_nvalue() [2/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::required_nvalue ( const unsigned n) const
inlinevirtual

Number of values required at local node n. In order to simplify matters, we allocate storage for pressure variables at all the nodes and then pin those that are not used.

Reimplemented from oomph::QTaylorHoodSpaceTimeElement< DIM >.

778  {
779  // There are DIM velocity components and 1 pressure component
780  return DIM + 1;
781  } // End of required_nvalue

References DIM.

◆ required_nvalue() [3/3]

template<unsigned DIM>
unsigned oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::required_nvalue ( const unsigned n) const
inlinevirtual

Number of values required at local node n. In order to simplify matters, we allocate storage for pressure variables at all the nodes and then pin those that are not used.

Reimplemented from oomph::QTaylorHoodSpaceTimeElement< DIM >.

778  {
779  // There are DIM velocity components and 1 pressure component
780  return DIM + 1;
781  } // End of required_nvalue

References DIM.

◆ unpin_elemental_pressure_dofs() [1/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::unpin_elemental_pressure_dofs ( )
inlineprivatevirtual

Unpin all pressure dofs.

Implements oomph::RefineableSpaceTimeNavierStokesEquations< DIM >.

1192  {
1193  // Find the index at which the pressure is stored
1194  int p_index = this->p_nodal_index_nst();
1195 
1196  // Get the number of nodes in the element
1197  unsigned n_node = this->nnode();
1198 
1199  // Loop over nodes
1200  for (unsigned n = 0; n < n_node; n++)
1201  {
1202  // Unpin the p_index-th dof at node n
1203  this->node_pt(n)->unpin(p_index);
1204  }
1205  } // End of unpin_elemental_pressure_dofs
void unpin(const unsigned &i)
Unpin the i-th stored variable.
Definition: nodes.h:391

References n, oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), oomph::QTaylorHoodSpaceTimeElement< DIM >::p_nodal_index_nst(), and oomph::Data::unpin().

◆ unpin_elemental_pressure_dofs() [2/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::unpin_elemental_pressure_dofs ( )
inlineprivatevirtual

Unpin all pressure dofs.

Implements oomph::RefineableSpaceTimeNavierStokesEquations< DIM >.

1192  {
1193  // Find the index at which the pressure is stored
1194  int p_index = this->p_nodal_index_nst();
1195 
1196  // Get the number of nodes in the element
1197  unsigned n_node = this->nnode();
1198 
1199  // Loop over nodes
1200  for (unsigned n = 0; n < n_node; n++)
1201  {
1202  // Unpin the p_index-th dof at node n
1203  this->node_pt(n)->unpin(p_index);
1204  }
1205  } // End of unpin_elemental_pressure_dofs

References n, oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), oomph::QTaylorHoodSpaceTimeElement< DIM >::p_nodal_index_nst(), and oomph::Data::unpin().

◆ unpin_elemental_pressure_dofs() [3/3]

template<unsigned DIM>
void oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::unpin_elemental_pressure_dofs ( )
inlineprivatevirtual

Unpin all pressure dofs.

Implements oomph::RefineableSpaceTimeNavierStokesEquations< DIM >.

1192  {
1193  // Find the index at which the pressure is stored
1194  int p_index = this->p_nodal_index_nst();
1195 
1196  // Get the number of nodes in the element
1197  unsigned n_node = this->nnode();
1198 
1199  // Loop over nodes
1200  for (unsigned n = 0; n < n_node; n++)
1201  {
1202  // Unpin the p_index-th dof at node n
1203  this->node_pt(n)->unpin(p_index);
1204  }
1205  } // End of unpin_elemental_pressure_dofs

References n, oomph::FiniteElement::nnode(), oomph::FiniteElement::node_pt(), oomph::QTaylorHoodSpaceTimeElement< DIM >::p_nodal_index_nst(), and oomph::Data::unpin().

◆ vertex_node_pt() [1/3]

template<unsigned DIM>
Node* oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::vertex_node_pt ( const unsigned j) const
inlinevirtual

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

Implements oomph::ElementWithZ2ErrorEstimator.

815  {
816  // Call the base class implementation of the function
818  } // End of vertex_node_pt
virtual Node * vertex_node_pt(const unsigned &j) const
Definition: elements.h:2500

References j, and oomph::FiniteElement::vertex_node_pt().

◆ vertex_node_pt() [2/3]

template<unsigned DIM>
Node* oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::vertex_node_pt ( const unsigned j) const
inlinevirtual

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

Implements oomph::ElementWithZ2ErrorEstimator.

815  {
816  // Call the base class implementation of the function
818  } // End of vertex_node_pt

References j, and oomph::FiniteElement::vertex_node_pt().

◆ vertex_node_pt() [3/3]

template<unsigned DIM>
Node* oomph::RefineableQTaylorHoodSpaceTimeElement< DIM >::vertex_node_pt ( const unsigned j) const
inlinevirtual

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

Implements oomph::ElementWithZ2ErrorEstimator.

815  {
816  // Call the base class implementation of the function
818  } // End of vertex_node_pt

References j, and oomph::FiniteElement::vertex_node_pt().


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