Public Member Functions | List of all members
oomph::RefineableQSpectralElement< 3 > Class Reference

#include <refineable_brick_spectral_element.h>

+ Inheritance diagram for oomph::RefineableQSpectralElement< 3 >:

Public Member Functions

 RefineableQSpectralElement ()
 Constructor. More...
 
 RefineableQSpectralElement (const RefineableQSpectralElement< 3 > &dummy)=delete
 Broken copy constructor. More...
 
virtual ~RefineableQSpectralElement ()
 Broken assignment operator. More...
 
void rebuild_from_sons (Mesh *&mesh_pt)
 The only thing to add is rebuild from sons. More...
 
virtual bool nodes_built ()
 Overload the nodes built function. More...
 
- Public Member Functions inherited from oomph::RefineableQElement< 3 >
 RefineableQElement ()
 Constructor: Pass refinement level (default 0 = root) More...
 
 RefineableQElement (const RefineableQElement< 3 > &dummy)=delete
 Broken copy constructor. More...
 
virtual ~RefineableQElement ()
 Broken assignment operator. More...
 
unsigned required_nsons () const
 A refineable brick element has eight sons. More...
 
virtual Node * node_created_by_neighbour (const Vector< double > &s_fraction, bool &is_periodic)
 
virtual Node * node_created_by_son_of_neighbour (const Vector< double > &s_fraction, bool &is_periodic)
 
virtual void build (Mesh *&mesh_pt, Vector< Node * > &new_node_pt, bool &was_already_built, std::ofstream &new_nodes_file)
 
void check_integrity (double &max_error)
 
void output_corners (std::ostream &outfile, const std::string &colour) const
 Print corner nodes, use colour. More...
 
OcTree * octree_pt ()
 Pointer to octree representation of this element. More...
 
OcTree * octree_pt () const
 Pointer to octree representation of this element. More...
 
void setup_hanging_nodes (Vector< std::ofstream * > &output_stream)
 
virtual void further_setup_hanging_nodes ()=0
 
- 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...
 
Tree * tree_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...
 
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)
 
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 ()
 
long number () const
 Element number (for debugging/plotting) More...
 
void set_number (const long &mynumber)
 Set element number (for debugging/plotting) More...
 
virtual unsigned ncont_interpolated_values () const =0
 
virtual void get_interpolated_values (const Vector< double > &s, Vector< double > &values)
 
virtual void get_interpolated_values (const unsigned &t, const Vector< double > &s, Vector< double > &values)=0
 
virtual Node * interpolating_node_pt (const unsigned &n, const int &value_id)
 
virtual double local_one_d_fraction_of_interpolating_node (const unsigned &n1d, const unsigned &i, const int &value_id)
 
virtual Node * get_interpolating_node_at_local_coordinate (const Vector< double > &s, const int &value_id)
 
virtual unsigned ninterpolating_node (const int &value_id)
 
virtual unsigned ninterpolating_node_1d (const int &value_id)
 
virtual void interpolating_basis (const Vector< double > &s, Shape &psi, const int &value_id) const
 
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 RefineableElement * root_element_pt ()
 
virtual RefineableElement * father_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 further_build ()
 Further build: e.g. deal with interpolation of internal values. More...
 
void get_dresidual_dnodal_coordinates (RankThreeTensor< double > &dresidual_dnodal_coordinates)
 
unsigned nshape_controlling_nodes ()
 
std::map< Node *, unsigned > shape_controlling_node_lookup ()
 
- Public Member Functions inherited from oomph::FiniteElement
void set_dimension (const unsigned &dim)
 
void set_nodal_dimension (const unsigned &nodal_dim)
 
void set_nnodal_position_type (const unsigned &nposition_type)
 Set the number of types required to interpolate the coordinate. More...
 
void set_n_node (const unsigned &n)
 
int nodal_local_eqn (const unsigned &n, const unsigned &i) const
 
double dJ_eulerian_at_knot (const unsigned &ipt, Shape &psi, DenseMatrix< double > &djacobian_dX) const
 
 FiniteElement ()
 Constructor. More...
 
virtual ~FiniteElement ()
 
 FiniteElement (const FiniteElement &)=delete
 Broken copy constructor. More...
 
void get_centre_of_gravity_and_max_radius_in_terms_of_zeta (Vector< double > &cog, double &max_radius) const
 
virtual void local_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)
 
MacroElement * macro_elem_pt ()
 Access function to pointer to macro element. More...
 
void get_x (const Vector< double > &s, Vector< double > &x) const
 
void get_x (const unsigned &t, const Vector< double > &s, Vector< double > &x)
 
virtual void 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
 
virtual void disable_ALE ()
 
virtual void enable_ALE ()
 
virtual unsigned required_nvalue (const unsigned &n) const
 
unsigned nnodal_position_type () const
 
bool has_hanging_nodes () const
 
unsigned nodal_dimension () const
 Return the required Eulerian dimension of the nodes in this element. More...
 
virtual Node * construct_node (const unsigned &n)
 
virtual Node * construct_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 
virtual Node * construct_boundary_node (const unsigned &n)
 
virtual Node * construct_boundary_node (const unsigned &n, TimeStepper *const &time_stepper_pt)
 
int get_node_number (Node *const &node_pt) const
 
virtual Node * get_node_at_local_coordinate (const Vector< double > &s) const
 
double raw_nodal_value (const unsigned &n, const unsigned &i) const
 
double raw_nodal_value (const unsigned &t, const unsigned &n, const unsigned &i) const
 
double nodal_value (const unsigned &n, const unsigned &i) const
 
double nodal_value (const unsigned &t, const unsigned &n, const unsigned &i) const
 
unsigned dim () const
 
virtual double interpolated_x (const Vector< double > &s, const unsigned &i) const
 Return FE interpolated coordinate x[i] at local coordinate s. More...
 
virtual double interpolated_x (const unsigned &t, const Vector< double > &s, const unsigned &i) const
 
virtual void interpolated_x (const Vector< double > &s, Vector< double > &x) const
 Return FE interpolated position x[] at local coordinate s as Vector. More...
 
virtual void interpolated_x (const unsigned &t, const Vector< double > &s, Vector< double > &x) const
 
virtual double interpolated_dxdt (const Vector< double > &s, const unsigned &i, const unsigned &t)
 
virtual void interpolated_dxdt (const Vector< double > &s, const unsigned &t, Vector< double > &dxdt)
 
unsigned ngeom_data () const
 
Data * geom_data_pt (const unsigned &j)
 
void position (const Vector< double > &zeta, Vector< double > &r) const
 
void position (const unsigned &t, const Vector< double > &zeta, Vector< double > &r) const
 
void dposition_dt (const Vector< double > &zeta, const unsigned &t, Vector< double > &drdt)
 
virtual double zeta_nodal (const unsigned &n, const unsigned &k, const unsigned &i) const
 
void interpolated_zeta (const Vector< double > &s, Vector< double > &zeta) const
 
void locate_zeta (const Vector< double > &zeta, GeomObject *&geom_object_pt, Vector< double > &s, const bool &use_coordinate_as_initial_guess=false)
 
virtual void node_update ()
 
virtual void identify_geometric_data (std::set< Data * > &geometric_data_pt)
 
virtual double s_min () const
 Min value of local coordinate. More...
 
virtual double s_max () const
 Max. value of local coordinate. More...
 
double size () const
 
virtual double compute_physical_size () const
 
virtual void point_output_data (const Vector< double > &s, Vector< double > &data)
 
void point_output (std::ostream &outfile, const Vector< double > &s)
 
virtual unsigned nplot_points_paraview (const unsigned &nplot) const
 
virtual unsigned nsub_elements_paraview (const unsigned &nplot) const
 
void output_paraview (std::ofstream &file_out, const unsigned &nplot) const
 
virtual void write_paraview_output_offset_information (std::ofstream &file_out, const unsigned &nplot, unsigned &counter) const
 
virtual void write_paraview_type (std::ofstream &file_out, const unsigned &nplot) const
 
virtual void write_paraview_offsets (std::ofstream &file_out, const unsigned &nplot, unsigned &offset_sum) const
 
virtual unsigned nscalar_paraview () const
 
virtual void scalar_value_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot) const
 
virtual void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt) const
 
virtual void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt) const
 
virtual std::string scalar_name_paraview (const unsigned &i) const
 
virtual void output (std::ostream &outfile)
 
virtual void output (std::ostream &outfile, const unsigned &n_plot)
 
virtual void output (const unsigned &t, std::ostream &outfile, const unsigned &n_plot) const
 
virtual void output (FILE *file_pt)
 
virtual void output (FILE *file_pt, const unsigned &n_plot)
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 Output an exact solution over the element. More...
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 Output a time-dependent exact solution over the element. More...
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, const SolutionFunctorBase &exact_soln) const
 Output a time-dependent exact solution over the element. More...
 
virtual void get_s_plot (const unsigned &i, const unsigned &nplot, Vector< double > &s, const bool &shifted_to_interior=false) const
 
virtual std::string tecplot_zone_string (const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (std::ostream &outfile, const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (FILE *file_pt, const unsigned &nplot) const
 
virtual unsigned nplot_points (const unsigned &nplot) const
 
virtual void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 Calculate the norm of the error and that of the exact solution. More...
 
virtual void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 Calculate the norm of the error and that of the exact solution. More...
 
virtual void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, Vector< double > &error, Vector< double > &norm)
 
virtual void compute_abs_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error)
 
void integrate_fct (FiniteElement::SteadyExactSolutionFctPt integrand_fct_pt, Vector< double > &integral)
 Integrate Vector-valued function over element. More...
 
void integrate_fct (FiniteElement::UnsteadyExactSolutionFctPt integrand_fct_pt, const double &time, Vector< double > &integral)
 Integrate Vector-valued time-dep function over element. More...
 
virtual void build_face_element (const int &face_index, FaceElement *face_element_pt)
 
virtual unsigned self_test ()
 
virtual unsigned get_bulk_node_number (const int &face_index, const unsigned &i) const
 
virtual int face_outer_unit_normal_sign (const int &face_index) const
 Get the sign of the outer unit normal on the face given by face_index. More...
 
void face_node_number_error_check (const unsigned &i) const
 Range check for face node numbers. More...
 
virtual CoordinateMappingFctPt face_to_bulk_coordinate_fct_pt (const int &face_index) const
 
virtual BulkCoordinateDerivativesFctPt bulk_coordinate_derivatives_fct_pt (const int &face_index) const
 
- Public Member Functions inherited from oomph::GeneralisedElement
 GeneralisedElement ()
 Constructor: Initialise all pointers and all values to zero. More...
 
virtual ~GeneralisedElement ()
 Virtual destructor to clean up any memory allocated by the object. More...
 
 GeneralisedElement (const GeneralisedElement &)=delete
 Broken copy constructor. More...
 
void operator= (const GeneralisedElement &)=delete
 Broken assignment operator. More...
 
Data *& internal_data_pt (const unsigned &i)
 Return a pointer to i-th internal data object. More...
 
Data *const & internal_data_pt (const unsigned &i) const
 Return a pointer to i-th internal data object (const version) More...
 
Data *& external_data_pt (const unsigned &i)
 Return a pointer to i-th external data object. More...
 
Data *const & external_data_pt (const unsigned &i) const
 Return a pointer to i-th external data object (const version) More...
 
unsigned long eqn_number (const unsigned &ieqn_local) const
 
int local_eqn_number (const unsigned long &ieqn_global) const
 
unsigned add_external_data (Data *const &data_pt, const bool &fd=true)
 
bool external_data_fd (const unsigned &i) const
 
void exclude_external_data_fd (const unsigned &i)
 
void include_external_data_fd (const unsigned &i)
 
void flush_external_data ()
 Flush all external data. More...
 
void flush_external_data (Data *const &data_pt)
 Flush the object addressed by data_pt from the external data array. More...
 
unsigned ninternal_data () const
 Return the number of internal data objects. More...
 
unsigned nexternal_data () const
 Return the number of external data objects. More...
 
unsigned ndof () const
 Return the number of equations/dofs in the element. More...
 
void dof_vector (const unsigned &t, Vector< double > &dof)
 Return the vector of dof values at time level t. More...
 
void dof_pt_vector (Vector< double * > &dof_pt)
 Return the vector of pointers to dof values. More...
 
void set_internal_data_time_stepper (const unsigned &i, TimeStepper *const &time_stepper_pt, const bool &preserve_existing_data)
 
void assign_internal_eqn_numbers (unsigned long &global_number, Vector< double * > &Dof_pt)
 
void describe_dofs (std::ostream &out, const std::string &current_string) const
 
void add_internal_value_pt_to_map (std::map< unsigned, double * > &map_of_value_pt)
 
virtual void assign_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void complete_setup_of_dependencies ()
 
virtual void get_residuals (Vector< double > &residuals)
 
virtual void get_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
virtual void get_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &mass_matrix)
 
virtual void get_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
virtual void get_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
virtual void get_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
virtual void get_djacobian_and_dmass_matrix_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam)
 
virtual void get_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
virtual void get_inner_products (Vector< std::pair< unsigned, unsigned >> const &history_index, Vector< double > &inner_product)
 
virtual void get_inner_product_vectors (Vector< unsigned > const &history_index, Vector< Vector< double >> &inner_product_vector)
 
virtual void compute_norm (Vector< double > &norm)
 
virtual void compute_norm (double &norm)
 
virtual unsigned ndof_types () const
 
virtual void get_dof_numbers_for_unknowns (std::list< std::pair< unsigned long, unsigned >> &dof_lookup_list) const
 
- Public Member Functions inherited from oomph::GeomObject
 GeomObject ()
 Default constructor. More...
 
 GeomObject (const unsigned &ndim)
 
 GeomObject (const unsigned &nlagrangian, const unsigned &ndim)
 
 GeomObject (const unsigned &nlagrangian, const unsigned &ndim, TimeStepper *time_stepper_pt)
 
 GeomObject (const GeomObject &dummy)=delete
 Broken copy constructor. More...
 
void operator= (const GeomObject &)=delete
 Broken assignment operator. More...
 
virtual ~GeomObject ()
 (Empty) destructor More...
 
unsigned nlagrangian () const
 Access function to # of Lagrangian coordinates. More...
 
unsigned ndim () const
 Access function to # of Eulerian coordinates. More...
 
void set_nlagrangian_and_ndim (const unsigned &n_lagrangian, const unsigned &n_dim)
 Set # of Lagrangian and Eulerian coordinates. More...
 
TimeStepper *& time_stepper_pt ()
 
TimeStepper * time_stepper_pt () const
 
virtual void position (const double &t, const Vector< double > &zeta, Vector< double > &r) const
 
virtual void dposition (const Vector< double > &zeta, DenseMatrix< double > &drdzeta) const
 
virtual void d2position (const Vector< double > &zeta, RankThreeTensor< double > &ddrdzeta) const
 
virtual void d2position (const Vector< double > &zeta, Vector< double > &r, DenseMatrix< double > &drdzeta, RankThreeTensor< double > &ddrdzeta) const
 
- Public Member Functions inherited from oomph::BrickElementBase
 BrickElementBase ()
 Constructor. Empty. More...
 
virtual unsigned nvertex_node () const =0
 Number of vertex nodes in the element. More...
 
virtual Node * vertex_node_pt (const unsigned &j) const =0
 Pointer to the j-th vertex node in the element. More...
 
- Public Member Functions inherited from oomph::QElementBase
 QElementBase ()
 Constructor: Initialise pointers to macro element reference coords. More...
 
 QElementBase (const QElementBase &)=delete
 Broken copy constructor. More...
 
virtual ~QElementBase ()
 Broken assignment operator. More...
 
bool local_coord_is_valid (const Vector< double > &s)
 Check whether the local coordinate are valid or not. More...
 
void move_local_coord_back_into_element (Vector< double > &s) const
 
virtual void set_macro_elem_pt (MacroElement *macro_elem_pt)
 
double & s_macro_ll (const unsigned &i)
 
double & s_macro_ur (const unsigned &i)
 
double s_macro_ll (const unsigned &i) const
 
double s_macro_ur (const unsigned &i) const
 
void get_x_from_macro_element (const Vector< double > &s, Vector< double > &x) const
 
void get_x_from_macro_element (const unsigned &t, const Vector< double > &s, Vector< double > &x)
 
unsigned nnode_on_face () const
 
ElementGeometry::ElementGeometry element_geometry () const
 It's a Q element! More...
 
- Public Member Functions inherited from oomph::QElementGeometricBase
 QElementGeometricBase ()
 Empty default constructor. More...
 
 QElementGeometricBase (const QElementGeometricBase &)=delete
 Broken copy constructor. More...
 

Additional Inherited Members

- Public Types inherited from oomph::RefineableQElement< 3 >
typedef void(RefineableQElement< 3 >::* VoidMemFctPt) ()
 
- 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::RefineableElement
static double & max_integrity_tolerance ()
 Max. allowed discrepancy in element integrity check. 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::RefineableQElement< 3 >
void setup_father_bounds ()
 
void get_face_bcs (const int &edge, Vector< int > &bound_cons) const
 
void get_boundaries (const int &edge, std::set< unsigned > &boundaries) const
 
void get_bcs (int bound, Vector< int > &bound_cons) const
 
void interpolated_zeta_on_face (const unsigned &boundary, const int &face, const Vector< double > &s, Vector< double > &zeta)
 
void setup_hang_for_value (const int &value_id)
 
virtual void oc_hang_helper (const int &value_id, const int &my_edge, std::ofstream &output_hangfile)
 
- 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)
 
- 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)
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 Zero-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 One-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 Two-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
- Protected Member Functions inherited from oomph::GeneralisedElement
unsigned add_internal_data (Data *const &data_pt, const bool &fd=true)
 
bool internal_data_fd (const unsigned &i) const
 
void exclude_internal_data_fd (const unsigned &i)
 
void include_internal_data_fd (const unsigned &i)
 
void clear_global_eqn_numbers ()
 
void add_global_eqn_numbers (std::deque< unsigned long > const &global_eqn_numbers, std::deque< double * > const &global_dof_pt)
 
virtual void assign_internal_and_external_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void assign_additional_local_eqn_numbers ()
 
int internal_local_eqn (const unsigned &i, const unsigned &j) const
 
int external_local_eqn (const unsigned &i, const unsigned &j)
 
virtual void fill_in_contribution_to_residuals (Vector< double > &residuals)
 
void fill_in_jacobian_from_internal_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_internal_by_fd (DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_external_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_external_by_fd (DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
virtual void update_before_internal_fd ()
 
virtual void reset_after_internal_fd ()
 
virtual void update_in_internal_fd (const unsigned &i)
 
virtual void reset_in_internal_fd (const unsigned &i)
 
virtual void update_before_external_fd ()
 
virtual void reset_after_external_fd ()
 
virtual void update_in_external_fd (const unsigned &i)
 
virtual void reset_in_external_fd (const unsigned &i)
 
virtual void fill_in_contribution_to_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &mass_matrix)
 
virtual void fill_in_contribution_to_jacobian_and_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &jacobian, DenseMatrix< double > &mass_matrix)
 
virtual void fill_in_contribution_to_dresiduals_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam)
 
virtual void fill_in_contribution_to_djacobian_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam)
 
virtual void fill_in_contribution_to_djacobian_and_dmass_matrix_dparameter (double *const &parameter_pt, Vector< double > &dres_dparam, DenseMatrix< double > &djac_dparam, DenseMatrix< double > &dmass_matrix_dparam)
 
virtual void fill_in_contribution_to_hessian_vector_products (Vector< double > const &Y, DenseMatrix< double > const &C, DenseMatrix< double > &product)
 
virtual void fill_in_contribution_to_inner_products (Vector< std::pair< unsigned, unsigned >> const &history_index, Vector< double > &inner_product)
 
virtual void fill_in_contribution_to_inner_product_vectors (Vector< unsigned > const &history_index, Vector< Vector< double >> &inner_product_vector)
 
- 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::RefineableElement
Tree * Tree_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...
 
- Protected Attributes inherited from oomph::FiniteElement
MacroElement * Macro_elem_pt
 Pointer to the element's macro element (NULL by default) More...
 
- Protected Attributes inherited from oomph::GeomObject
unsigned NLagrangian
 Number of Lagrangian (intrinsic) coordinates. More...
 
unsigned Ndim
 Number of Eulerian coordinates. More...
 
TimeStepper * Geom_object_time_stepper_pt
 
- Static Protected Attributes inherited from oomph::RefineableQElement< 3 >
static std::map< unsigned, DenseMatrix< int > > Father_bound
 
- Static Protected Attributes inherited from oomph::RefineableElement
static double Max_integrity_tolerance = 1.0e-8
 Max. allowed discrepancy in element integrity check. More...
 
- Static Protected Attributes inherited from oomph::FiniteElement
static const unsigned Default_Initial_Nvalue = 0
 Default value for the number of values at a node. More...
 
static const double Node_location_tolerance = 1.0e-14
 
static const unsigned N2deriv [] = {0, 1, 3, 6}
 
- Static Protected Attributes inherited from oomph::GeneralisedElement
static DenseMatrix< double > Dummy_matrix
 
static std::deque< double * > Dof_pt_deque
 

Detailed Description

Refineable version of QuadElements that add functionality for spectral Elements.

Constructor & Destructor Documentation

◆ RefineableQSpectralElement() [1/2]

oomph::RefineableQSpectralElement< 3 >::RefineableQSpectralElement ( )
inline

Constructor.

50  : RefineableElement()
51  {
52 #ifdef LEAK_CHECK
53  LeakCheckNames::RefineableQSpectralElement<3> _build += 1;
54 #endif
55  }
oomph::RefineableElement::RefineableElement
RefineableElement()
Definition: refineable_elements.h:188

◆ RefineableQSpectralElement() [2/2]

oomph::RefineableQSpectralElement< 3 >::RefineableQSpectralElement ( const RefineableQSpectralElement< 3 > &  dummy)
delete

Broken copy constructor.

◆ ~RefineableQSpectralElement()

virtual oomph::RefineableQSpectralElement< 3 >::~RefineableQSpectralElement ( )
inlinevirtual

Broken assignment operator.

Destructor

72  {
73 #ifdef LEAK_CHECK
74  LeakCheckNames::RefineableQSpectralElement<3> _build -= 1;
75 #endif
76  }

Member Function Documentation

◆ nodes_built()

virtual bool oomph::RefineableQSpectralElement< 3 >::nodes_built ( )
inlinevirtual

Overload the nodes built function.

Reimplemented from oomph::RefineableElement.

555  {
556  unsigned n_node = this->nnode();
557  for (unsigned n = 0; n < n_node; n++)
558  {
559  if (node_pt(n) == 0)
560  {
561  return false;
562  }
563  }
564  // If we get to here, OK
565  return true;
566  }
n
const unsigned n
Definition: CG3DPackingUnitTest.cpp:11
oomph::FiniteElement::node_pt
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175
oomph::FiniteElement::nnode
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210

References n.

◆ rebuild_from_sons()

void oomph::RefineableQSpectralElement< 3 >::rebuild_from_sons ( Mesh *&  mesh_pt)
inlinevirtual

The only thing to add is rebuild from sons.

Implements oomph::RefineableElement.

80  {
81  // The timestepper should be the same for all nodes and node 0 should
82  // never be deleted.
83  if (this->node_pt(0) == 0)
84  {
85  throw OomphLibError("The Corner node (0) does not exist",
86  OOMPH_CURRENT_FUNCTION,
87  OOMPH_EXCEPTION_LOCATION);
88  }
89 
90  TimeStepper* time_stepper_pt = this->node_pt(0)->time_stepper_pt();
91  unsigned ntstorage = time_stepper_pt->ntstorage();
92 
93  unsigned jnod = 0;
94  Vector<double> s_fraction(3), s(3);
95  // Loop over the nodes in the element
96  unsigned n_p = this->nnode_1d();
97  for (unsigned i0 = 0; i0 < n_p; i0++)
98  {
99  // Get the fractional position of the node
100  s_fraction[0] = this->local_one_d_fraction_of_node(i0, 0);
101  // Local coordinate
102  s[0] = -1.0 + 2.0 * s_fraction[0];
103 
104  for (unsigned i1 = 0; i1 < n_p; i1++)
105  {
106  // Get the fractional position of the node in the direction of s[1]
107  s_fraction[1] = this->local_one_d_fraction_of_node(i1, 1);
108  // Local coordinate in father element
109  s[1] = -1.0 + 2.0 * s_fraction[1];
110 
111  for (unsigned i2 = 0; i2 < n_p; i2++)
112  {
113  // Get the fractional position of the node in the direction of s[1]
114  s_fraction[2] = this->local_one_d_fraction_of_node(i2, 2);
115  // Local coordinate in father element
116  s[2] = -1.0 + 2.0 * s_fraction[2];
117 
118  // Set the local node number
119  jnod = i0 + n_p * i1 + n_p * n_p * i2;
120 
121  // If the node has not been built
122  if (this->node_pt(jnod) == 0)
123  {
124  // Has the node been created by one of its neighbours
125  bool is_periodic = false;
126  Node* created_node_pt =
127  this->node_created_by_neighbour(s_fraction, is_periodic);
128 
129  // If it has set the pointer
130  if (created_node_pt != 0)
131  {
132  // If the node is periodic
133  if (is_periodic)
134  {
135  throw OomphLibError("Cannot handle periodic nodes in "
136  "refineable spectral elements yet",
137  OOMPH_CURRENT_FUNCTION,
138  OOMPH_EXCEPTION_LOCATION);
139  }
140  // Non-periodic case, just set the pointer
141  else
142  {
143  this->node_pt(jnod) = created_node_pt;
144  }
145  }
146  // Otherwise, we need to build it
147  else
148  {
149  // First we need to find the pointer to the son that
150  // would contain the node
151 
152  // Find coordinates in the sons
153  Vector<double> s_in_son(3);
154  using namespace OcTreeNames;
155  int son = -10;
156  // If less than one-half on the left side
157  if (s_fraction[0] < 0.5)
158  {
159  // On the down side
160  if (s_fraction[1] < 0.5)
161  {
162  // On the back side
163  if (s_fraction[2] < 0.5)
164  {
165  // It's the left down back son
166  son = LDB;
167  s_in_son[0] = -1.0 + 4.0 * s_fraction[0];
168  s_in_son[1] = -1.0 + 4.0 * s_fraction[1];
169  s_in_son[2] = -1.0 + 4.0 * s_fraction[2];
170  }
171  // On the front side
172  else
173  {
174  // It's the left down front son
175  son = LDF;
176  s_in_son[0] = -1.0 + 4.0 * s_fraction[0];
177  s_in_son[1] = -1.0 + 4.0 * s_fraction[1];
178  s_in_son[2] = -1.0 + 4.0 * (s_fraction[2] - 0.5);
179  }
180  } // End of on down side
181  // else its on the top
182  else
183  {
184  // On the back
185  if (s_fraction[2] < 0.5)
186  {
187  // It's the left up back son
188  son = LUB;
189  s_in_son[0] = -1.0 + 4.0 * s_fraction[0];
190  s_in_son[1] = -1.0 + 4.0 * (s_fraction[1] - 0.5);
191  s_in_son[2] = -1.0 + 4.0 * s_fraction[2];
192  }
193  // On the front side
194  else
195  {
196  // It's the left up front son
197  son = LUF;
198  s_in_son[0] = -1.0 + 4.0 * s_fraction[0];
199  s_in_son[1] = -1.0 + 4.0 * (s_fraction[1] - 0.5);
200  s_in_son[2] = -1.0 + 4.0 * (s_fraction[2] - 0.5);
201  }
202  } // End of on top
203  } // End of on left
204  // Otherwise its on the right
205  else
206  {
207  // On the down side
208  if (s_fraction[1] < 0.5)
209  {
210  // On the back side
211  if (s_fraction[2] < 0.5)
212  {
213  // It's the right down back son
214  son = RDB;
215  s_in_son[0] = -1.0 + 4.0 * (s_fraction[0] - 0.5);
216  s_in_son[1] = -1.0 + 4.0 * s_fraction[1];
217  s_in_son[2] = -1.0 + 4.0 * s_fraction[2];
218  }
219  // On the front side
220  else
221  {
222  // It's the right down front son
223  son = RDF;
224  s_in_son[0] = -1.0 + 4.0 * (s_fraction[0] - 0.5);
225  s_in_son[1] = -1.0 + 4.0 * s_fraction[1];
226  s_in_son[2] = -1.0 + 4.0 * (s_fraction[2] - 0.5);
227  }
228  } // End of on down side
229  // else its on the top
230  else
231  {
232  // On the back
233  if (s_fraction[2] < 0.5)
234  {
235  // It's the right up back son
236  son = RUB;
237  s_in_son[0] = -1.0 + 4.0 * (s_fraction[0] - 0.5);
238  s_in_son[1] = -1.0 + 4.0 * (s_fraction[1] - 0.5);
239  s_in_son[2] = -1.0 + 4.0 * s_fraction[2];
240  }
241  // On the front side
242  else
243  {
244  // It's the right up front son
245  son = RUF;
246  s_in_son[0] = -1.0 + 4.0 * (s_fraction[0] - 0.5);
247  s_in_son[1] = -1.0 + 4.0 * (s_fraction[1] - 0.5);
248  s_in_son[2] = -1.0 + 4.0 * (s_fraction[2] - 0.5);
249  }
250  }
251  }
252 
253  // Get the pointer to the son element
254  RefineableQSpectralElement<3>* son_el_pt =
255  dynamic_cast<RefineableQSpectralElement<3>*>(
256  this->tree_pt()->son_pt(son)->object_pt());
257 
258  // If we are rebuilding, then worry about the boundary
259  // conditions Find the boundary of the node Initially none
260  int boundary = Tree::OMEGA;
261  // If we are on the left boundary
262  if (i0 == 0)
263  {
264  boundary = L;
265  }
266  // If we are on the right boundary
267  else if (i0 == n_p - 1)
268  {
269  boundary = R;
270  }
271 
272  // If we are on the lower boundary
273  if (i1 == 0)
274  {
275  // If we already have already set the boundary, we're on an
276  // edge
277  switch (boundary)
278  {
279  case L:
280  boundary = LD;
281  break;
282  case R:
283  boundary = RD;
284  break;
285  // Boundary not set
286  default:
287  boundary = D;
288  break;
289  }
290  }
291  // If we are the northern bounadry
292  else if (i1 == n_p - 1)
293  {
294  // If we already have a boundary
295  switch (boundary)
296  {
297  case L:
298  boundary = LU;
299  break;
300  case R:
301  boundary = RU;
302  break;
303  default:
304  boundary = U;
305  break;
306  }
307  }
308 
309  // If we are on the back face
310  if (i2 == 0)
311  {
312  // If we already have boundaries
313  switch (boundary)
314  {
315  case L:
316  boundary = LB;
317  break;
318  case R:
319  boundary = RB;
320  break;
321  case U:
322  boundary = UB;
323  break;
324  case D:
325  boundary = DB;
326  break;
327  case LD:
328  boundary = LDB;
329  break;
330  case RD:
331  boundary = RDB;
332  break;
333  case LU:
334  boundary = LUB;
335  break;
336  case RU:
337  boundary = RUB;
338  break;
339  default:
340  boundary = B;
341  break;
342  }
343  }
344  else if (i2 == n_p - 1)
345  {
346  // If we already have boundaries
347  switch (boundary)
348  {
349  case L:
350  boundary = LF;
351  break;
352  case R:
353  boundary = RF;
354  break;
355  case U:
356  boundary = UF;
357  break;
358  case D:
359  boundary = DF;
360  break;
361  case LD:
362  boundary = LDF;
363  break;
364  case RD:
365  boundary = RDF;
366  break;
367  case LU:
368  boundary = LUF;
369  break;
370  case RU:
371  boundary = RUF;
372  break;
373  default:
374  boundary = F;
375  break;
376  }
377  }
378 
379  // set of boundaries that this edge in the son lives on
380  std::set<unsigned> boundaries;
381  // If we are on a boundary of the Element, find the
382  // mesh boundaries on which we live
383  // The boundaries will be common to the son because there can be
384  // no rotations here
385  if (boundary != Tree::OMEGA)
386  {
387  son_el_pt->get_boundaries(boundary, boundaries);
388  }
389 
390  // If the node lives on a boundary:
391  // construct a boundary node,
392  // get boundary conditions and
393  // update both lookup schemes
394  if (boundaries.size() > 0)
395  {
396  // Construct the new node
397  this->node_pt(jnod) =
398  this->construct_boundary_node(jnod, time_stepper_pt);
399 
400  // Get the boundary conditions from the son
401  Vector<int> bound_cons(ncont_interpolated_values());
402  son_el_pt->get_bcs(boundary, bound_cons);
403 
404  // Loop over the values and pin if necessary
405  unsigned nval = this->node_pt(jnod)->nvalue();
406  for (unsigned k = 0; k < nval; k++)
407  {
408  if (bound_cons[k])
409  {
410  this->node_pt(jnod)->pin(k);
411  }
412  }
413 
414  // Solid node? If so, deal with the positional boundary
415  // conditions:
416  SolidNode* solid_node_pt =
417  dynamic_cast<SolidNode*>(this->node_pt(jnod));
418  if (solid_node_pt != 0)
419  {
420  // Get the positional boundary conditions from the father:
421  unsigned n_dim = this->node_pt(jnod)->ndim();
422  Vector<int> solid_bound_cons(n_dim);
423  RefineableSolidQElement<3>* son_solid_el_pt =
424  dynamic_cast<RefineableSolidQElement<3>*>(son_el_pt);
425 #ifdef PARANOID
426  if (son_solid_el_pt == 0)
427  {
428  std::string error_message = "We have a SolidNode outside "
429  "a refineable SolidElement\n";
430  error_message +=
431  "during mesh refinement -- this doesn't make sense\n";
432 
433  throw OomphLibError(error_message,
434  OOMPH_CURRENT_FUNCTION,
435  OOMPH_EXCEPTION_LOCATION);
436  }
437 #endif
438  son_solid_el_pt->get_solid_bcs(boundary, solid_bound_cons);
439 
440  // Loop over the positions and pin, if necessary
441  for (unsigned k = 0; k < n_dim; k++)
442  {
443  if (solid_bound_cons[k])
444  {
445  solid_node_pt->pin_position(k);
446  }
447  }
448  }
449 
450  // Next we update the boundary look-up schemes
451  // Loop over the boundaries stored in the set
452  for (std::set<unsigned>::iterator it = boundaries.begin();
453  it != boundaries.end();
454  ++it)
455  {
456  // Add the node to the boundary
457  mesh_pt->add_boundary_node(*it, this->node_pt(jnod));
458 
459  // If we have set an intrinsic coordinate on this
460  // mesh boundary then it must also be interpolated on
461  // the new node
462  // Now interpolate the intrinsic boundary coordinate
463  if (mesh_pt->boundary_coordinate_exists(*it) == true)
464  {
465  Vector<double> zeta(2);
466  son_el_pt->interpolated_zeta_on_face(
467  *it, boundary, s_in_son, zeta);
468 
469  this->node_pt(jnod)->set_coordinates_on_boundary(*it,
470  zeta);
471  }
472  }
473  }
474  // Otherwise we construct a normal "bulk" node
475  else
476  {
477  // Construct the new node
478  this->node_pt(jnod) =
479  this->construct_node(jnod, time_stepper_pt);
480  }
481 
482  // Now we set the position and values at the newly created node
483 
484  // In the first instance use macro element or FE representation
485  // to create past and present nodal positions.
486  // (THIS STEP SHOULD NOT BE SKIPPED FOR ALGEBRAIC
487  // ELEMENTS AS NOT ALL OF THEM NECESSARILY IMPLEMENT
488  // NONTRIVIAL NODE UPDATE FUNCTIONS. CALLING
489  // THE NODE UPDATE FOR SUCH ELEMENTS/NODES WILL LEAVE
490  // THEIR NODAL POSITIONS WHERE THEY WERE (THIS IS APPROPRIATE
491  // ONCE THEY HAVE BEEN GIVEN POSITIONS) BUT WILL
492  // NOT ASSIGN SENSIBLE INITIAL POSITONS!
493 
494  // Loop over # of history values
495  for (unsigned t = 0; t < ntstorage; t++)
496  {
497  // Get the position from the son
498  Vector<double> x_prev(3);
499 
500  // Now let's fill in the value
501  son_el_pt->get_x(t, s_in_son, x_prev);
502  for (unsigned i = 0; i < 3; i++)
503  {
504  this->node_pt(jnod)->x(t, i) = x_prev[i];
505  }
506  }
507 
508  // Now set the values
509  // Loop over all history values
510  for (unsigned t = 0; t < ntstorage; t++)
511  {
512  // Get values from father element
513  // Note: get_interpolated_values() sets Vector size itself.
514  Vector<double> prev_values;
515  son_el_pt->get_interpolated_values(t, s_in_son, prev_values);
516 
517  // Initialise the values at the new node
518  for (unsigned k = 0; k < this->node_pt(jnod)->nvalue(); k++)
519  {
520  this->node_pt(jnod)->set_value(t, k, prev_values[k]);
521  }
522  }
523 
524  // Add the node to the mesh
525  mesh_pt->add_node_pt(this->node_pt(jnod));
526 
527  } // Node has been constructed
528 
529  // Algebraic stuff here
530  // Check whether the element is an algebraic element
531  AlgebraicElementBase* alg_el_pt =
532  dynamic_cast<AlgebraicElementBase*>(this);
533 
534  // If we do have an algebraic element
535  if (alg_el_pt != 0)
536  {
537  std::string error_message =
538  "Have not implemented rebuilding from sons for";
539  error_message += "Algebraic Spectral elements yet\n";
540 
541  throw OomphLibError(
542  error_message,
543  "RefineableQSpectralElement::rebuild_from_sons()",
544  OOMPH_EXCEPTION_LOCATION);
545  }
546 
547  } // End of the case when the node was not built
548  }
549  }
550  }
551  }
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9
D
dominoes D
Definition: Domino.cpp:55
L
MatrixXd L
Definition: LLT_example.cpp:6
R
@ R
Definition: StatisticsVector.h:21
B
Matrix< SCALARB, Dynamic, Dynamic, opt_B > B
Definition: bench_gemm.cpp:48
oomph::Data::pin
void pin(const unsigned &i)
Pin the i-th stored variable.
Definition: nodes.h:385
oomph::Data::time_stepper_pt
TimeStepper *& time_stepper_pt()
Return the pointer to the timestepper.
Definition: nodes.h:238
oomph::Data::set_value
void set_value(const unsigned &i, const double &value_)
Definition: nodes.h:271
oomph::Data::nvalue
unsigned nvalue() const
Return number of values stored in data object (incl pinned ones).
Definition: nodes.h:483
oomph::FiniteElement::construct_node
virtual Node * construct_node(const unsigned &n)
Definition: elements.h:2509
oomph::FiniteElement::nnode_1d
virtual unsigned nnode_1d() const
Definition: elements.h:2218
oomph::FiniteElement::construct_boundary_node
virtual Node * construct_boundary_node(const unsigned &n)
Definition: elements.h:2538
oomph::FiniteElement::local_one_d_fraction_of_node
virtual double local_one_d_fraction_of_node(const unsigned &n1d, const unsigned &i)
Definition: elements.h:1858
oomph::GeomObject::time_stepper_pt
TimeStepper *& time_stepper_pt()
Definition: geom_objects.h:192
oomph::Node::x
double & x(const unsigned &i)
Return the i-th nodal coordinate.
Definition: nodes.h:1060
oomph::Node::ndim
unsigned ndim() const
Return (Eulerian) spatial dimension of the node.
Definition: nodes.h:1054
oomph::Node::set_coordinates_on_boundary
virtual void set_coordinates_on_boundary(const unsigned &b, const unsigned &k, const Vector< double > &boundary_zeta)
Definition: nodes.cc:2394
oomph::RefineableElement::tree_pt
Tree * tree_pt()
Access function: Pointer to quadtree representation of this element.
Definition: refineable_elements.h:211
oomph::RefineableElement::ncont_interpolated_values
virtual unsigned ncont_interpolated_values() const =0
oomph::RefineableQElement< 3 >::node_created_by_neighbour
virtual Node * node_created_by_neighbour(const Vector< double > &s_fraction, bool &is_periodic)
Definition: refineable_brick_element.cc:992
oomph::TimeStepper::ntstorage
unsigned ntstorage() const
Definition: timesteppers.h:601
oomph::Tree::object_pt
RefineableElement * object_pt() const
Definition: tree.h:88
oomph::Tree::OMEGA
static const int OMEGA
Default value for an unassigned neighbour.
Definition: tree.h:262
oomph::Tree::son_pt
Tree * son_pt(const int &son_index) const
Definition: tree.h:103
s
RealScalar s
Definition: level1_cplx_impl.h:130
k
char char char int int * k
Definition: level2_impl.h:374
Eigen::zeta
EIGEN_STRONG_INLINE const Eigen::CwiseBinaryOp< Eigen::internal::scalar_zeta_op< typename DerivedX::Scalar >, const DerivedX, const DerivedQ > zeta(const Eigen::ArrayBase< DerivedX > &x, const Eigen::ArrayBase< DerivedQ > &q)
Definition: SpecialFunctionsArrayAPI.h:152
RachelsAdvectionDiffusion::U
double U
Swimming speed.
Definition: two_d_variable_diff_adapt.cc:53
oomph::Global_string_for_annotation::string
std::string string(const unsigned &i)
Definition: oomph_definitions.cc:286
oomph::OcTreeNames::RDF
@ RDF
Definition: octree.h:54
oomph::OcTreeNames::DF
@ DF
Definition: octree.h:67
oomph::OcTreeNames::RD
@ RD
Definition: octree.h:62
oomph::OcTreeNames::RUB
@ RUB
Definition: octree.h:52
oomph::OcTreeNames::RB
@ RB
Definition: octree.h:58
oomph::OcTreeNames::UF
@ UF
Definition: octree.h:68
oomph::OcTreeNames::LUF
@ LUF
Definition: octree.h:55
oomph::OcTreeNames::RU
@ RU
Definition: octree.h:64
oomph::OcTreeNames::F
@ F
Definition: octree.h:74
oomph::OcTreeNames::LDF
@ LDF
Definition: octree.h:53
oomph::OcTreeNames::LF
@ LF
Definition: octree.h:65
oomph::OcTreeNames::RDB
@ RDB
Definition: octree.h:50
oomph::OcTreeNames::LD
@ LD
Definition: octree.h:61
oomph::OcTreeNames::LUB
@ LUB
Definition: octree.h:51
oomph::OcTreeNames::RUF
@ RUF
Definition: octree.h:56
oomph::OcTreeNames::RF
@ RF
Definition: octree.h:66
oomph::OcTreeNames::LDB
@ LDB
Definition: octree.h:49
oomph::OcTreeNames::LU
@ LU
Definition: octree.h:63
oomph::OcTreeNames::DB
@ DB
Definition: octree.h:59
oomph::OcTreeNames::UB
@ UB
Definition: octree.h:60
oomph::OcTreeNames::LB
@ LB
Definition: octree.h:57
plotPSD.t
t
Definition: plotPSD.py:36
OOMPH_EXCEPTION_LOCATION
#define OOMPH_EXCEPTION_LOCATION
Definition: oomph_definitions.h:61
OOMPH_CURRENT_FUNCTION
#define OOMPH_CURRENT_FUNCTION
Definition: oomph_definitions.h:86

References oomph::Mesh::add_boundary_node(), oomph::Mesh::add_node_pt(), oomph::Mesh::boundary_coordinate_exists(), D, oomph::OcTreeNames::DB, oomph::OcTreeNames::DF, oomph::OcTreeNames::F, oomph::RefineableQElement< 3 >::get_bcs(), oomph::RefineableQElement< 3 >::get_boundaries(), oomph::RefineableElement::get_interpolated_values(), oomph::RefineableSolidQElement< 3 >::get_solid_bcs(), oomph::FiniteElement::get_x(), i, oomph::RefineableQElement< 3 >::interpolated_zeta_on_face(), k, L, oomph::OcTreeNames::LB, oomph::OcTreeNames::LD, oomph::OcTreeNames::LDB, oomph::OcTreeNames::LDF, oomph::OcTreeNames::LF, oomph::OcTreeNames::LU, oomph::OcTreeNames::LUB, oomph::OcTreeNames::LUF, oomph::Node::ndim(), oomph::TimeStepper::ntstorage(), oomph::Tree::OMEGA, OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, oomph::SolidNode::pin_position(), R, oomph::OcTreeNames::RB, oomph::OcTreeNames::RD, oomph::OcTreeNames::RDB, oomph::OcTreeNames::RDF, oomph::OcTreeNames::RF, oomph::OcTreeNames::RU, oomph::OcTreeNames::RUB, oomph::OcTreeNames::RUF, s, oomph::Global_string_for_annotation::string(), plotPSD::t, RachelsAdvectionDiffusion::U, oomph::OcTreeNames::UB, oomph::OcTreeNames::UF, and Eigen::zeta().

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