Public Member Functions | List of all members
oomph::ProjectableDarcyElement< DARCY_ELEMENT > Class Template Reference

Darcy upgraded to become projectable. More...

#include <darcy_elements.h>

+ Inheritance diagram for oomph::ProjectableDarcyElement< DARCY_ELEMENT >:

Public Member Functions

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

Additional Inherited Members

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

Detailed Description

template<class DARCY_ELEMENT>
class oomph::ProjectableDarcyElement< DARCY_ELEMENT >

Darcy upgraded to become projectable.

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

Constructor & Destructor Documentation

◆ ProjectableDarcyElement()

template<class DARCY_ELEMENT >
oomph::ProjectableDarcyElement< DARCY_ELEMENT >::ProjectableDarcyElement ( )
inline

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

524 {}

Member Function Documentation

◆ data_values_of_field()

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

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

Implements oomph::ProjectableElementBase.

530  {
531 #ifdef PARANOID
532  if (fld > 1)
533  {
534  std::stringstream error_stream;
535  error_stream << "Darcy elements only store 2 fields so fld = " << fld
536  << " is illegal \n";
537  throw OomphLibError(
538  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
539  }
540 #endif
541 
542  // Create the vector
543  Vector<std::pair<Data*, unsigned>> data_values;
544 
545  // Pressure
546  if (fld == 0)
547  {
548  Data* dat_pt = this->p_data_pt();
549  unsigned nvalue = dat_pt->nvalue();
550  for (unsigned i = 0; i < nvalue; i++)
551  {
552  data_values.push_back(std::make_pair(dat_pt, i));
553  }
554  }
555  else
556  {
557  Vector<Data*> edge_dat_pt = this->q_edge_data_pt();
558  unsigned n = edge_dat_pt.size();
559  for (unsigned j = 0; j < n; j++)
560  {
561  unsigned nvalue = edge_dat_pt[j]->nvalue();
562  for (unsigned i = 0; i < nvalue; i++)
563  {
564  data_values.push_back(std::make_pair(edge_dat_pt[j], i));
565  }
566  }
567  if (this->nq_basis_internal() > 0)
568  {
569  Data* int_dat_pt = this->q_internal_data_pt();
570  unsigned nvalue = int_dat_pt->nvalue();
571  for (unsigned i = 0; i < nvalue; i++)
572  {
573  data_values.push_back(std::make_pair(int_dat_pt, i));
574  }
575  }
576  }
577 
578  // Return the vector
579  return data_values;
580  }
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9
n
const unsigned n
Definition: CG3DPackingUnitTest.cpp:11
OOMPH_EXCEPTION_LOCATION
#define OOMPH_EXCEPTION_LOCATION
Definition: oomph_definitions.h:61
OOMPH_CURRENT_FUNCTION
#define OOMPH_CURRENT_FUNCTION
Definition: oomph_definitions.h:86
j
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References i, j, n, oomph::Data::nvalue(), OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

◆ get_field()

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

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

Implements oomph::ProjectableElementBase.

682  {
683 #ifdef PARANOID
684  if (fld > 1)
685  {
686  std::stringstream error_stream;
687  error_stream << "Darcy elements only store two fields so fld = " << fld
688  << " is illegal\n";
689  throw OomphLibError(
690  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
691  }
692 #endif
693 
694  double return_value = 0.0;
695 
696  // Pressure
697  if (fld == 0)
698  {
699  this->interpolated_p(s, return_value);
700  }
701  else
702  {
703  throw OomphLibError("Don't call this function for Darcy!",
704  OOMPH_CURRENT_FUNCTION,
705  OOMPH_EXCEPTION_LOCATION);
706  }
707 
708  return return_value;
709  }
s
RealScalar s
Definition: level1_cplx_impl.h:130

References OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, and s.

Referenced by oomph::ProjectableDarcyElement< DARCY_ELEMENT >::residual_for_projection().

◆ jacobian_and_shape_of_field()

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

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

Implements oomph::ProjectableElementBase.

617  {
618 #ifdef PARANOID
619  if (fld > 1)
620  {
621  std::stringstream error_stream;
622  error_stream << "Darcy elements only store two fields so fld = " << fld
623  << " is illegal.\n";
624  throw OomphLibError(
625  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
626  }
627 #endif
628 
629 
630  // Get the number of geometric nodes, total number of basis functions,
631  // and number of edges basis functions
632  const unsigned n_dim = this->dim();
633  const unsigned n_node = this->nnode();
634  const unsigned n_q_basis = this->nq_basis();
635  const unsigned n_p_basis = this->np_basis();
636 
637  // Storage for the geometric and computational bases and the test
638  // functions
639  Shape psi_geom(n_node), q_basis(n_q_basis, n_dim),
640  q_test(n_q_basis, n_dim), p_basis(n_p_basis), p_test(n_p_basis),
641  div_q_basis_ds(n_q_basis), div_q_test_ds(n_q_basis);
642  double J = this->shape_basis_test_local(s,
643  psi_geom,
644  q_basis,
645  q_test,
646  p_basis,
647  p_test,
648  div_q_basis_ds,
649  div_q_test_ds);
650  // Pressure basis functions
651  if (fld == 0)
652  {
653  unsigned n = p_basis.nindex1();
654  for (unsigned i = 0; i < n; i++)
655  {
656  psi[i] = p_basis[i];
657  }
658  }
659  // Flux basis functions
660  else
661  {
662  unsigned n = q_basis.nindex1();
663  unsigned m = q_basis.nindex2();
664  for (unsigned i = 0; i < n; i++)
665  {
666  for (unsigned j = 0; j < m; j++)
667  {
668  psi(i, j) = q_basis(i, j);
669  }
670  }
671  }
672 
673  return J;
674  }
J
JacobiRotation< float > J
Definition: Jacobi_makeJacobi.cpp:3
oomph::FiniteElement::dim
unsigned dim() const
Definition: elements.h:2611
oomph::FiniteElement::nnode
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210
m
int * m
Definition: level2_cplx_impl.h:294

References oomph::FiniteElement::dim(), i, J, j, m, n, oomph::FiniteElement::nnode(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, and s.

Referenced by oomph::ProjectableDarcyElement< DARCY_ELEMENT >::residual_for_projection().

◆ local_equation()

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

Return local equation number of value j in field fld.

Implements oomph::ProjectableElementBase.

742  {
743 #ifdef PARANOID
744  if (fld > 1)
745  {
746  std::stringstream error_stream;
747  error_stream << "Darcy elements only store two fields so fld = " << fld
748  << " is illegal\n";
749  throw OomphLibError(
750  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
751  }
752 #endif
753 
754  int return_value = 0;
755  // Pressure
756  if (fld == 0)
757  {
758  return_value = this->p_local_eqn(j);
759  }
760  else
761  {
762  unsigned nedge = this->nq_basis_edge();
763  if (j < nedge)
764  {
765  return_value = this->q_edge_local_eqn(j);
766  }
767  else
768  {
769  return_value = this->q_internal_local_eqn(j - nedge);
770  }
771  }
772 
773  return return_value;
774  }

References j, OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

Referenced by oomph::ProjectableDarcyElement< DARCY_ELEMENT >::residual_for_projection().

◆ nfields_for_projection()

template<class DARCY_ELEMENT >
unsigned oomph::ProjectableDarcyElement< DARCY_ELEMENT >::nfields_for_projection ( )
inlinevirtual

Number of fields to be projected: 2 (pressure and flux)

Implements oomph::ProjectableElementBase.

584  {
585  return 2;
586  }

◆ nhistory_values_for_coordinate_projection()

template<class DARCY_ELEMENT >
unsigned oomph::ProjectableDarcyElement< DARCY_ELEMENT >::nhistory_values_for_coordinate_projection ( )
inlinevirtual

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

Implements oomph::ProjectableElementBase.

608  {
609  return this->node_pt(0)->position_time_stepper_pt()->ntstorage();
610  }
oomph::FiniteElement::node_pt
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175
oomph::Node::position_time_stepper_pt
TimeStepper *& position_time_stepper_pt()
Return a pointer to the position timestepper.
Definition: nodes.h:1022
oomph::TimeStepper::ntstorage
unsigned ntstorage() const
Definition: timesteppers.h:601

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

◆ nhistory_values_for_projection()

template<class DARCY_ELEMENT >
unsigned oomph::ProjectableDarcyElement< DARCY_ELEMENT >::nhistory_values_for_projection ( const unsigned &  fld)
inlinevirtual

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

Implements oomph::ProjectableElementBase.

591  {
592 #ifdef PARANOID
593  if (fld > 1)
594  {
595  std::stringstream error_stream;
596  error_stream << "Darcy elements only store two fields so fld = " << fld
597  << " is illegal\n";
598  throw OomphLibError(
599  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
600  }
601 #endif
602  return this->node_pt(0)->ntstorage();
603  }
oomph::Data::ntstorage
unsigned ntstorage() const
Definition: nodes.cc:879

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

◆ nvalue_of_field()

template<class DARCY_ELEMENT >
unsigned oomph::ProjectableDarcyElement< DARCY_ELEMENT >::nvalue_of_field ( const unsigned &  fld)
inlinevirtual

Return number of values in field fld.

Implements oomph::ProjectableElementBase.

714  {
715 #ifdef PARANOID
716  if (fld > 1)
717  {
718  std::stringstream error_stream;
719  error_stream << "Darcy elements only store two fields so fld = " << fld
720  << " is illegal\n";
721  throw OomphLibError(
722  error_stream.str(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION);
723  }
724 #endif
725 
726  unsigned return_value = 0;
727  if (fld == 0)
728  {
729  return_value = this->np_basis();
730  }
731  else
732  {
733  return_value = this->nq_basis();
734  }
735 
736  return return_value;
737  }

References OOMPH_CURRENT_FUNCTION, and OOMPH_EXCEPTION_LOCATION.

Referenced by oomph::ProjectableDarcyElement< DARCY_ELEMENT >::residual_for_projection().

◆ output()

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

Output FE representation of soln as in underlying element.

Reimplemented from oomph::FiniteElement.

778  {
779  DARCY_ELEMENT::output(outfile, nplot);
780  }
output
void output(std::ostream &outfile, const unsigned &nplot)
Overload output function.
Definition: overloaded_element_body.h:490

References output().

◆ pin_superfluous_darcy_dofs()

template<class DARCY_ELEMENT >
void oomph::ProjectableDarcyElement< DARCY_ELEMENT >::pin_superfluous_darcy_dofs ( )
inline

Helper function to pin superfluous dofs; required because we introduce at least one dof per node to allow projection during unstructured refinement)

793  {
794  // Pin dofs at vertex nodes (because they're only used for projection)
795  for (unsigned j = 0; j < 3; j++)
796  {
797  this->node_pt(j)->pin(0);
798  }
799  }
oomph::Data::pin
void pin(const unsigned &i)
Pin the i-th stored variable.
Definition: nodes.h:385

References j, oomph::FiniteElement::node_pt(), and oomph::Data::pin().

◆ required_nvalue()

template<class DARCY_ELEMENT >
unsigned oomph::ProjectableDarcyElement< DARCY_ELEMENT >::required_nvalue ( const unsigned &  n) const
inlinevirtual

Overload required_nvalue to return at least one value.

Reimplemented from oomph::FiniteElement.

784  {
785  return std::max(this->Initial_Nvalue[n], unsigned(1));
786  }
max
#define max(a, b)
Definition: datatypes.h:23

References max, and n.

◆ residual_for_projection()

template<class DARCY_ELEMENT >
void oomph::ProjectableDarcyElement< DARCY_ELEMENT >::residual_for_projection ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian,
const unsigned &  flag 
)
inlinevirtual

Residual for the projection step. Flag indicates if we want the Jacobian (1) or not (0). Virtual so it can be overloaded if necessary

Reimplemented from oomph::ProjectableElement< DARCY_ELEMENT >.

807  {
808  // Am I a solid element
809  SolidFiniteElement* solid_el_pt = dynamic_cast<SolidFiniteElement*>(this);
810 
811  unsigned n_dim = this->dim();
812 
813  // Allocate storage for local coordinates
814  Vector<double> s(n_dim);
815 
816  // Current field
817  unsigned fld = this->Projected_field;
818 
819  // Number of nodes
820  const unsigned n_node = this->nnode();
821  // Number of positional dofs
822  const unsigned n_position_type = this->nnodal_position_type();
823 
824  // Number of dof for current field
825  const unsigned n_value = nvalue_of_field(fld);
826 
827  // Set the value of n_intpt
828  const unsigned n_intpt = this->integral_pt()->nweight();
829 
830  // Loop over the integration points
831  for (unsigned ipt = 0; ipt < n_intpt; ipt++)
832  {
833  // Get the local coordinates of Gauss point
834  for (unsigned i = 0; i < n_dim; i++)
835  s[i] = this->integral_pt()->knot(ipt, i);
836 
837  // Get the integral weight
838  double w = this->integral_pt()->weight(ipt);
839 
840  // Find same point in base mesh using external storage
841  FiniteElement* other_el_pt = 0;
842  other_el_pt = this->external_element_pt(0, ipt);
843  Vector<double> other_s(n_dim);
844  other_s = this->external_element_local_coord(0, ipt);
845 
846  ProjectableElement<DARCY_ELEMENT>* cast_el_pt =
847  dynamic_cast<ProjectableElement<DARCY_ELEMENT>*>(other_el_pt);
848 
849  // Storage for the local equation and local unknown
850  int local_eqn = 0, local_unknown = 0;
851 
852  // Now set up the three different projection problems
853  switch (Projection_type)
854  {
855  case Lagrangian:
856  {
857  // If we don't have a solid element there is a problem
858  if (solid_el_pt == 0)
859  {
860  throw OomphLibError("Trying to project Lagrangian coordinates in "
861  "non-SolidElement\n",
862  OOMPH_CURRENT_FUNCTION,
863  OOMPH_EXCEPTION_LOCATION);
864  }
865 
866  // Find the position shape functions
867  Shape psi(n_node, n_position_type);
868  // Get the position shape functions
869  this->shape(s, psi);
870  // Get the jacobian
871  double J = this->J_eulerian(s);
872 
873  // Premultiply the weights and the Jacobian
874  double W = w * J;
875 
876  // Get the value of the current position of the 0th coordinate
877  // in the current element
878  // at the current time level (which is the unkown)
879  double interpolated_xi_proj = this->interpolated_x(s, 0);
880 
881  // Get the Lagrangian position in the other element
882  double interpolated_xi_bar =
883  dynamic_cast<SolidFiniteElement*>(cast_el_pt)
884  ->interpolated_xi(other_s, Projected_lagrangian);
885 
886  // Now loop over the nodes and position dofs
887  for (unsigned l = 0; l < n_node; ++l)
888  {
889  // Loop over position unknowns
890  for (unsigned k = 0; k < n_position_type; ++k)
891  {
892  // The local equation is going to be the positional local
893  // equation
894  local_eqn = solid_el_pt->position_local_eqn(l, k, 0);
895 
896  // Now assemble residuals
897  if (local_eqn >= 0)
898  {
899  // calculate residuals
900  residuals[local_eqn] +=
901  (interpolated_xi_proj - interpolated_xi_bar) * psi(l, k) *
902  W;
903 
904  // Calculate the jacobian
905  if (flag == 1)
906  {
907  for (unsigned l2 = 0; l2 < n_node; l2++)
908  {
909  // Loop over position dofs
910  for (unsigned k2 = 0; k2 < n_position_type; k2++)
911  {
912  local_unknown =
913  solid_el_pt->position_local_eqn(l2, k2, 0);
914  if (local_unknown >= 0)
915  {
916  jacobian(local_eqn, local_unknown) +=
917  psi(l2, k2) * psi(l, k) * W;
918  }
919  }
920  }
921  } // end of jacobian
922  }
923  }
924  }
925  } // End of Lagrangian coordinate case
926 
927  break;
928 
929  // Now the coordinate history case
930  case Coordinate:
931  {
932  // Find the position shape functions
933  Shape psi(n_node, n_position_type);
934  // Get the position shape functions
935  this->shape(s, psi);
936  // Get the jacobian
937  double J = this->J_eulerian(s);
938 
939  // Premultiply the weights and the Jacobian
940  double W = w * J;
941 
942  // Get the value of the current position in the current element
943  // at the current time level (which is the unkown)
944  double interpolated_x_proj = 0.0;
945  // If we are a solid element read it out directly from the data
946  if (solid_el_pt != 0)
947  {
948  interpolated_x_proj =
949  this->interpolated_x(s, Projected_coordinate);
950  }
951  // Otherwise it's the field value at the current time
952  else
953  {
954  interpolated_x_proj = this->get_field(0, fld, s);
955  }
956 
957  // Get the position in the other element
958  double interpolated_x_bar = cast_el_pt->interpolated_x(
959  Time_level_for_projection, other_s, Projected_coordinate);
960 
961  // Now loop over the nodes and position dofs
962  for (unsigned l = 0; l < n_node; ++l)
963  {
964  // Loop over position unknowns
965  for (unsigned k = 0; k < n_position_type; ++k)
966  {
967  // If I'm a solid element
968  if (solid_el_pt != 0)
969  {
970  // The local equation is going to be the positional local
971  // equation
972  local_eqn =
973  solid_el_pt->position_local_eqn(l, k, Projected_coordinate);
974  }
975  // Otherwise just pick the local unknown of a field to
976  // project into
977  else
978  {
979  // Complain if using generalised position types
980  // but this is not a solid element, because the storage
981  // is then not clear
982  if (n_position_type != 1)
983  {
984  throw OomphLibError("Trying to project generalised "
985  "positions not in SolidElement\n",
986  OOMPH_CURRENT_FUNCTION,
987  OOMPH_EXCEPTION_LOCATION);
988  }
989  local_eqn = local_equation(fld, l);
990  }
991 
992  // Now assemble residuals
993  if (local_eqn >= 0)
994  {
995  // calculate residuals
996  residuals[local_eqn] +=
997  (interpolated_x_proj - interpolated_x_bar) * psi(l, k) * W;
998 
999  // Calculate the jacobian
1000  if (flag == 1)
1001  {
1002  for (unsigned l2 = 0; l2 < n_node; l2++)
1003  {
1004  // Loop over position dofs
1005  for (unsigned k2 = 0; k2 < n_position_type; k2++)
1006  {
1007  // If I'm a solid element
1008  if (solid_el_pt != 0)
1009  {
1010  local_unknown = solid_el_pt->position_local_eqn(
1011  l2, k2, Projected_coordinate);
1012  }
1013  else
1014  {
1015  local_unknown = local_equation(fld, l2);
1016  }
1017 
1018  if (local_unknown >= 0)
1019  {
1020  jacobian(local_eqn, local_unknown) +=
1021  psi(l2, k2) * psi(l, k) * W;
1022  }
1023  }
1024  }
1025  } // end of jacobian
1026  }
1027  }
1028  }
1029  } // End of coordinate case
1030  break;
1031 
1032  // Now the value case
1033  case Value:
1034  {
1035  // Pressure -- "normal" procedure
1036  if (fld == 0)
1037  {
1038  // Field shape function
1039  Shape psi(n_value);
1040 
1041  // Calculate jacobian and shape functions for that field
1042  double J = jacobian_and_shape_of_field(fld, s, psi);
1043 
1044  // Premultiply the weights and the Jacobian
1045  double W = w * J;
1046 
1047  // Value of field in current element at current time level
1048  //(the unknown)
1049  unsigned now = 0;
1050  double interpolated_value_proj = this->get_field(now, fld, s);
1051 
1052  // Value of the interpolation of element located in base mesh
1053  double interpolated_value_bar =
1054  cast_el_pt->get_field(Time_level_for_projection, fld, other_s);
1055 
1056  // Loop over dofs of field fld
1057  for (unsigned l = 0; l < n_value; l++)
1058  {
1059  local_eqn = local_equation(fld, l);
1060  if (local_eqn >= 0)
1061  {
1062  // calculate residuals
1063  residuals[local_eqn] +=
1064  (interpolated_value_proj - interpolated_value_bar) *
1065  psi[l] * W;
1066 
1067  // Calculate the jacobian
1068  if (flag == 1)
1069  {
1070  for (unsigned l2 = 0; l2 < n_value; l2++)
1071  {
1072  local_unknown = local_equation(fld, l2);
1073  if (local_unknown >= 0)
1074  {
1075  jacobian(local_eqn, local_unknown) +=
1076  psi[l2] * psi[l] * W;
1077  }
1078  }
1079  } // end of jacobian
1080  }
1081  }
1082  }
1083  // Flux -- need inner product
1084  else if (fld == 1)
1085  {
1086  // Field shape function
1087  Shape psi(n_value, n_dim);
1088 
1089  // Calculate jacobian and shape functions for that field
1090  double J = jacobian_and_shape_of_field(fld, s, psi);
1091 
1092  // Premultiply the weights and the Jacobian
1093  double W = w * J;
1094 
1095  // Value of flux in current element at current time level
1096  //(the unknown)
1097  Vector<double> q_proj(n_dim);
1098  this->interpolated_q(s, q_proj);
1099 
1100  // Value of the interpolation of element located in base mesh
1101  Vector<double> q_bar(n_dim);
1102  cast_el_pt->interpolated_q(other_s, q_bar);
1103 
1104 #ifdef PARANOID
1105  if (Time_level_for_projection != 0)
1106  {
1107  std::stringstream error_stream;
1108  error_stream << "Darcy elements are steady!\n";
1109  throw OomphLibError(error_stream.str(),
1110  OOMPH_CURRENT_FUNCTION,
1111  OOMPH_EXCEPTION_LOCATION);
1112  }
1113 #endif
1114 
1115  // Loop over dofs of field fld
1116  for (unsigned l = 0; l < n_value; l++)
1117  {
1118  local_eqn = local_equation(fld, l);
1119  if (local_eqn >= 0)
1120  {
1121  // Loop over spatial dimension for dot product
1122  for (unsigned i = 0; i < n_dim; i++)
1123  {
1124  // Add to residuals
1125  residuals[local_eqn] +=
1126  (q_proj[i] - q_bar[i]) * psi(l, i) * W;
1127 
1128  // Calculate the jacobian
1129  if (flag == 1)
1130  {
1131  for (unsigned l2 = 0; l2 < n_value; l2++)
1132  {
1133  local_unknown = local_equation(fld, l2);
1134  if (local_unknown >= 0)
1135  {
1136  jacobian(local_eqn, local_unknown) +=
1137  psi(l2, i) * psi(l, i) * W;
1138  }
1139  }
1140  } // end of jacobian
1141  }
1142  }
1143  }
1144  }
1145  else
1146  {
1147  throw OomphLibError(
1148  "Wrong field specified. This should never happen\n",
1149  OOMPH_CURRENT_FUNCTION,
1150  OOMPH_EXCEPTION_LOCATION);
1151  }
1152 
1153 
1154  break;
1155 
1156  default:
1157  throw OomphLibError("Wrong type specificied in Projection_type. "
1158  "This should never happen\n",
1159  OOMPH_CURRENT_FUNCTION,
1160  OOMPH_EXCEPTION_LOCATION);
1161  }
1162  } // End of the switch statement
1163 
1164  } // End of loop over ipt
1165 
1166  } // End of residual_for_projection function
w
RowVector3d w
Definition: Matrix_resize_int.cpp:3
oomph::ElementWithExternalElement::external_element_local_coord
Vector< double > & external_element_local_coord(const unsigned &interaction_index, const unsigned &ipt)
Definition: element_with_external_element.h:136
oomph::ElementWithExternalElement::external_element_pt
FiniteElement *& external_element_pt(const unsigned &interaction_index, const unsigned &ipt)
Definition: element_with_external_element.h:107
oomph::FiniteElement::nnodal_position_type
unsigned nnodal_position_type() const
Definition: elements.h:2463
oomph::FiniteElement::interpolated_x
virtual double interpolated_x(const Vector< double > &s, const unsigned &i) const
Return FE interpolated coordinate x[i] at local coordinate s.
Definition: elements.cc:3962
oomph::FiniteElement::shape
virtual void shape(const Vector< double > &s, Shape &psi) const =0
oomph::FiniteElement::integral_pt
Integral *const & integral_pt() const
Return the pointer to the integration scheme (const version)
Definition: elements.h:1963
oomph::FiniteElement::FiniteElement
FiniteElement()
Constructor.
Definition: elements.h:1782
oomph::FiniteElement::J_eulerian
virtual double J_eulerian(const Vector< double > &s) const
Definition: elements.cc:4103
oomph::Integral::knot
virtual double knot(const unsigned &i, const unsigned &j) const =0
Return local coordinate s[j] of i-th integration point.
oomph::Integral::nweight
virtual unsigned nweight() const =0
Return the number of integration points of the scheme.
oomph::Integral::weight
virtual double weight(const unsigned &i) const =0
Return weight of i-th integration point.
oomph::ProjectableDarcyElement::jacobian_and_shape_of_field
double jacobian_and_shape_of_field(const unsigned &fld, const Vector< double > &s, Shape &psi)
Definition: darcy_elements.h:614
oomph::ProjectableDarcyElement::get_field
double get_field(const unsigned &t, const unsigned &fld, const Vector< double > &s)
Definition: darcy_elements.h:679
oomph::ProjectableDarcyElement::nvalue_of_field
unsigned nvalue_of_field(const unsigned &fld)
Return number of values in field fld.
Definition: darcy_elements.h:713
oomph::ProjectableDarcyElement::local_equation
int local_equation(const unsigned &fld, const unsigned &j)
Return local equation number of value j in field fld.
Definition: darcy_elements.h:741
oomph::ProjectableElementBase::Projected_lagrangian
unsigned Projected_lagrangian
Definition: projection.h:78
oomph::ProjectableElementBase::Projection_type
Projection_Type Projection_type
Definition: projection.h:83
oomph::ProjectableElementBase::Time_level_for_projection
unsigned Time_level_for_projection
Time level we are projecting (0=current values; >0: history values)
Definition: projection.h:70
oomph::ProjectableElementBase::Value
@ Value
Definition: projection.h:63
oomph::ProjectableElementBase::Lagrangian
@ Lagrangian
Definition: projection.h:62
oomph::ProjectableElementBase::Coordinate
@ Coordinate
Definition: projection.h:61
oomph::ProjectableElementBase::Projected_field
unsigned Projected_field
Field that is currently being projected.
Definition: projection.h:67
oomph::ProjectableElementBase::Projected_coordinate
unsigned Projected_coordinate
Definition: projection.h:74
k
char char char int int * k
Definition: level2_impl.h:374
oomph::QuadTreeNames::W
@ W
Definition: quadtree.h:63

References oomph::ProjectableElementBase::Coordinate, oomph::FiniteElement::dim(), oomph::ElementWithExternalElement::external_element_local_coord(), oomph::ElementWithExternalElement::external_element_pt(), oomph::ProjectableDarcyElement< DARCY_ELEMENT >::get_field(), oomph::ProjectableElementBase::get_field(), i, oomph::FiniteElement::integral_pt(), oomph::FiniteElement::interpolated_x(), J, oomph::FiniteElement::J_eulerian(), oomph::ProjectableDarcyElement< DARCY_ELEMENT >::jacobian_and_shape_of_field(), k, oomph::Integral::knot(), oomph::ProjectableElementBase::Lagrangian, oomph::ProjectableDarcyElement< DARCY_ELEMENT >::local_equation(), oomph::FiniteElement::nnodal_position_type(), oomph::FiniteElement::nnode(), oomph::ProjectableDarcyElement< DARCY_ELEMENT >::nvalue_of_field(), oomph::Integral::nweight(), OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, oomph::SolidFiniteElement::position_local_eqn(), oomph::ProjectableElementBase::Projected_coordinate, oomph::ProjectableElementBase::Projected_field, oomph::ProjectableElementBase::Projected_lagrangian, oomph::ProjectableElementBase::Projection_type, s, oomph::FiniteElement::shape(), oomph::ProjectableElementBase::Time_level_for_projection, oomph::ProjectableElementBase::Value, w, oomph::QuadTreeNames::W, and oomph::Integral::weight().

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