Public Member Functions | Private Member Functions | List of all members
oomph::RefineablePMLHelmholtzEquations< DIM > Class Template Reference

#include <refineable_pml_helmholtz_elements.h>

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

Public Member Functions

 RefineablePMLHelmholtzEquations ()
 Constructor, simply call other constructors. More...
 
 RefineablePMLHelmholtzEquations (const RefineablePMLHelmholtzEquations< DIM > &dummy)=delete
 Broken copy constructor. More...
 
unsigned num_Z2_flux_terms ()
 Broken assignment operator. More...
 
void get_Z2_flux (const Vector< double > &s, Vector< double > &flux)
 
void get_interpolated_values (const Vector< double > &s, Vector< double > &values)
 
void get_interpolated_values (const unsigned &t, const Vector< double > &s, Vector< double > &values)
 
void further_build ()
 Further build: Copy source function pointer from father element. More...
 
- Public Member Functions inherited from oomph::PMLHelmholtzEquations< DIM >
 PMLHelmholtzEquations ()
 Constructor. More...
 
 PMLHelmholtzEquations (const PMLHelmholtzEquations &dummy)=delete
 Broken copy constructor. More...
 
virtual std::complex< unsigned > u_index_helmholtz () const
 Broken assignment operator. More...
 
double *& k_squared_pt ()
 Get pointer to k_squared. More...
 
double k_squared ()
 Get the square of wavenumber. More...
 
const double & alpha () const
 Alpha, wavenumber complex shift. More...
 
double *& alpha_pt ()
 Pointer to Alpha, wavenumber complex shift. More...
 
unsigned nscalar_paraview () const
 
void scalar_value_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot) const
 
void scalar_value_fct_paraview (std::ofstream &file_out, const unsigned &i, const unsigned &nplot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt) const
 
std::string scalar_name_paraview (const unsigned &i) const
 
void output (std::ostream &outfile)
 Output with default number of plot points. More...
 
void output (std::ostream &outfile, const unsigned &n_plot)
 
void output_total_real (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt incoming_wave_fct_pt, const double &phi, const unsigned &nplot)
 
void output_real (std::ostream &outfile, const double &phi, const unsigned &n_plot)
 
void output_imag (std::ostream &outfile, const double &phi, const unsigned &n_plot)
 
void output (FILE *file_pt)
 C_style output with default number of plot points. More...
 
void output (FILE *file_pt, const unsigned &n_plot)
 
void output_fct (std::ostream &outfile, const unsigned &n_plot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt)
 
void output_real_fct (std::ostream &outfile, const double &phi, const unsigned &n_plot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 
void output_imag_fct (std::ostream &outfile, const double &phi, const unsigned &n_plot, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt)
 
void compute_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, double &error, double &norm)
 Get error against and norm of exact solution. More...
 
void compute_error (std::ostream &outfile, FiniteElement::UnsteadyExactSolutionFctPt exact_soln_pt, const double &time, double &error, double &norm)
 Dummy, time dependent error checker. More...
 
void compute_norm (double &norm)
 Compute norm of fe solution. More...
 
PMLHelmholtzSourceFctPt & source_fct_pt ()
 Access function: Pointer to source function. More...
 
PMLHelmholtzSourceFctPt source_fct_pt () const
 Access function: Pointer to source function. Const version. More...
 
virtual void get_source_helmholtz (const unsigned &ipt, const Vector< double > &x, std::complex< double > &source) const
 
void values_to_be_pinned_on_outer_pml_boundary (Vector< unsigned > &values_to_pin)
 
void get_flux (const Vector< double > &s, Vector< std::complex< double >> &flux) const
 Get flux: flux[i] = du/dx_i for real and imag part. More...
 
void fill_in_contribution_to_residuals (Vector< double > &residuals)
 Add the element's contribution to its residual vector (wrapper) More...
 
void fill_in_contribution_to_jacobian (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
std::complex< double > interpolated_u_pml_helmholtz (const Vector< double > &s) const
 
unsigned self_test ()
 Self-test: Return 0 for OK. More...
 
void compute_pml_coefficients (const unsigned &ipt, const Vector< double > &x, Vector< std::complex< double >> &pml_laplace_factor, std::complex< double > &pml_k_squared_factor)
 
PMLMapping *& pml_mapping_pt ()
 Return a pointer to the PML Mapping object. More...
 
PMLMapping *const & pml_mapping_pt () const
 Return a pointer to the PML Mapping object (const version) More...
 
unsigned ndof_types () const
 
void get_dof_numbers_for_unknowns (std::list< std::pair< unsigned long, unsigned >> &dof_lookup_list) const
 
- Public Member Functions inherited from oomph::PMLElementBase< DIM >
 PMLElementBase ()
 Constructor. More...
 
virtual ~PMLElementBase ()
 Virtual destructor. More...
 
void disable_pml ()
 
void enable_pml (const int &direction, const double &interface_border_value, const double &outer_domain_border_value)
 
- 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 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 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)
 
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 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 void output (const unsigned &t, std::ostream &outfile, const unsigned &n_plot) const
 
virtual void output_fct (std::ostream &outfile, const unsigned &n_plot, const double &time, const SolutionFunctorBase &exact_soln) const
 Output a time-dependent exact solution over the element. More...
 
virtual void get_s_plot (const unsigned &i, const unsigned &nplot, Vector< double > &s, const bool &shifted_to_interior=false) const
 
virtual std::string tecplot_zone_string (const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (std::ostream &outfile, const unsigned &nplot) const
 
virtual void write_tecplot_zone_footer (FILE *file_pt, const unsigned &nplot) const
 
virtual unsigned nplot_points (const unsigned &nplot) const
 
virtual void compute_error (FiniteElement::SteadyExactSolutionFctPt exact_soln_pt, 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, 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 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)
 
- 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::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...
 
virtual unsigned required_nsons () const
 
bool refinement_is_enabled ()
 Flag to indicate suppression of any refinement. More...
 
void disable_refinement ()
 Suppress of any refinement for this element. More...
 
void enable_refinement ()
 Emnable refinement for this element. More...
 
template<class ELEMENT >
void split (Vector< ELEMENT * > &son_pt) const
 
int local_hang_eqn (Node *const &node_pt, const unsigned &i)
 
virtual void build (Mesh *&mesh_pt, Vector< Node * > &new_node_pt, bool &was_already_built, std::ofstream &new_nodes_file)=0
 
void set_refinement_level (const int &refine_level)
 Set the refinement level. More...
 
unsigned refinement_level () const
 Return the Refinement level. More...
 
void select_for_refinement ()
 Select the element for refinement. More...
 
void deselect_for_refinement ()
 Deselect the element for refinement. More...
 
void select_sons_for_unrefinement ()
 Unrefinement will be performed by merging the four sons of this element. More...
 
void deselect_sons_for_unrefinement ()
 
bool to_be_refined ()
 Has the element been selected for refinement? More...
 
bool sons_to_be_unrefined ()
 Has the element been selected for unrefinement? More...
 
virtual void rebuild_from_sons (Mesh *&mesh_pt)=0
 
virtual void unbuild ()
 
virtual void deactivate_element ()
 
virtual bool nodes_built ()
 Return true if all the nodes have been built, false if not. More...
 
long number () const
 Element number (for debugging/plotting) More...
 
void set_number (const long &mynumber)
 Set element number (for debugging/plotting) More...
 
virtual unsigned ncont_interpolated_values () const =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
 
virtual void check_integrity (double &max_error)=0
 
void identify_field_data_for_interactions (std::set< std::pair< Data *, unsigned >> &paired_field_data)
 
void assign_nodal_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual 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 setup_hanging_nodes (Vector< std::ofstream * > &output_stream)
 
virtual void further_setup_hanging_nodes ()
 
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::ElementWithZ2ErrorEstimator
 ElementWithZ2ErrorEstimator ()
 Default empty constructor. More...
 
 ElementWithZ2ErrorEstimator (const ElementWithZ2ErrorEstimator &)=delete
 Broken copy constructor. More...
 
void operator= (const ElementWithZ2ErrorEstimator &)=delete
 Broken assignment operator. More...
 
virtual unsigned ncompound_fluxes ()
 
virtual void compute_exact_Z2_error (std::ostream &outfile, FiniteElement::SteadyExactSolutionFctPt exact_flux_pt, double &error, double &norm)
 
virtual void get_Z2_compound_flux_indices (Vector< unsigned > &flux_index)
 
virtual unsigned nvertex_node () const =0
 Number of vertex nodes in the element. More...
 
virtual Node * vertex_node_pt (const unsigned &j) const =0
 
virtual unsigned nrecovery_order ()=0
 Order of recovery shape functions. More...
 
virtual double geometric_jacobian (const Vector< double > &x)
 

Private Member Functions

void fill_in_generic_residual_contribution_helmholtz (Vector< double > &residuals, DenseMatrix< double > &jacobian, const unsigned &flag)
 

Additional Inherited Members

- Public Types inherited from oomph::PMLHelmholtzEquations< DIM >
typedef void(* PMLHelmholtzSourceFctPt) (const Vector< double > &x, std::complex< double > &f)
 
- 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::PMLHelmholtzEquations< DIM >
static BermudezPMLMapping Default_pml_mapping
 
- 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::PMLHelmholtzEquations< DIM >
virtual double dshape_and_dtest_eulerian_helmholtz (const Vector< double > &s, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const =0
 
virtual double dshape_and_dtest_eulerian_at_knot_helmholtz (const unsigned &ipt, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const =0
 
- Protected Member Functions inherited from oomph::FiniteElement
template<unsigned DIM>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
virtual double invert_jacobian_mapping (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
virtual double local_to_eulerian_mapping (const DShape &dpsids, DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
double local_to_eulerian_mapping (const DShape &dpsids, DenseMatrix< double > &inverse_jacobian) const
 
virtual void dJ_eulerian_dnodal_coordinates (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<unsigned DIM>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
virtual void d_dshape_eulerian_dnodal_coordinates (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<unsigned DIM>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
virtual void transform_derivatives (const DenseMatrix< double > &inverse_jacobian, DShape &dbasis) const
 
void transform_derivatives_diagonal (const DenseMatrix< double > &inverse_jacobian, DShape &dbasis) const
 
virtual void transform_second_derivatives (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<unsigned DIM>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<unsigned DIM>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
void fill_in_jacobian_from_nodal_by_fd (DenseMatrix< double > &jacobian)
 
virtual void update_before_nodal_fd ()
 
virtual void reset_after_nodal_fd ()
 
virtual void update_in_nodal_fd (const unsigned &i)
 
virtual void reset_in_nodal_fd (const unsigned &i)
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 Zero-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 One-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 Two-d specialisation of function to calculate inverse of jacobian mapping. More...
 
template<>
double invert_jacobian (const DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void dJ_eulerian_dnodal_coordinates_templated_helper (const DenseMatrix< double > &jacobian, const DShape &dpsids, DenseMatrix< double > &djacobian_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void d_dshape_eulerian_dnodal_coordinates_templated_helper (const double &det_jacobian, const DenseMatrix< double > &jacobian, const DenseMatrix< double > &djacobian_dX, const DenseMatrix< double > &inverse_jacobian, const DShape &dpsids, RankFourTensor< double > &d_dpsidx_dX) const
 
template<>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_template (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
template<>
void transform_second_derivatives_diagonal (const DenseMatrix< double > &jacobian, const DenseMatrix< double > &inverse_jacobian, const DenseMatrix< double > &jacobian2, DShape &dbasis, DShape &d2basis) const
 
- Protected Member Functions inherited from oomph::GeneralisedElement
unsigned add_internal_data (Data *const &data_pt, const bool &fd=true)
 
bool internal_data_fd (const unsigned &i) const
 
void exclude_internal_data_fd (const unsigned &i)
 
void include_internal_data_fd (const unsigned &i)
 
void clear_global_eqn_numbers ()
 
void add_global_eqn_numbers (std::deque< unsigned long > const &global_eqn_numbers, std::deque< double * > const &global_dof_pt)
 
virtual void assign_internal_and_external_local_eqn_numbers (const bool &store_local_dof_pt)
 
virtual void assign_additional_local_eqn_numbers ()
 
int internal_local_eqn (const unsigned &i, const unsigned &j) const
 
int external_local_eqn (const unsigned &i, const unsigned &j)
 
void fill_in_jacobian_from_internal_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_internal_by_fd (DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_external_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
void fill_in_jacobian_from_external_by_fd (DenseMatrix< double > &jacobian, const bool &fd_all_data=false)
 
virtual void update_before_internal_fd ()
 
virtual void reset_after_internal_fd ()
 
virtual void update_in_internal_fd (const unsigned &i)
 
virtual void reset_in_internal_fd (const unsigned &i)
 
virtual void update_before_external_fd ()
 
virtual void reset_after_external_fd ()
 
virtual void update_in_external_fd (const unsigned &i)
 
virtual void reset_in_external_fd (const unsigned &i)
 
virtual void fill_in_contribution_to_mass_matrix (Vector< double > &residuals, DenseMatrix< double > &mass_matrix)
 
virtual void fill_in_contribution_to_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 Member Functions inherited from oomph::RefineableElement
void assemble_local_to_eulerian_jacobian (const DShape &dpsids, DenseMatrix< double > &jacobian) const
 
void assemble_local_to_eulerian_jacobian2 (const DShape &d2psids, DenseMatrix< double > &jacobian2) const
 
void assemble_eulerian_base_vectors (const DShape &dpsids, DenseMatrix< double > &interpolated_G) const
 
double local_to_eulerian_mapping_diagonal (const DShape &dpsids, DenseMatrix< double > &jacobian, DenseMatrix< double > &inverse_jacobian) const
 
void assign_hanging_local_eqn_numbers (const bool &store_local_dof_pt)
 Assign the local equation numbers for hanging node variables. More...
 
virtual void fill_in_jacobian_from_nodal_by_fd (Vector< double > &residuals, DenseMatrix< double > &jacobian)
 
- Static Protected Member Functions inherited from oomph::RefineableElement
static void check_value_id (const int &n_continuously_interpolated_values, const int &value_id)
 
- Protected Attributes inherited from oomph::PMLHelmholtzEquations< DIM >
double * Alpha_pt
 Pointer to wavenumber complex shift. More...
 
PMLHelmholtzSourceFctPt Source_fct_pt
 Pointer to source function: More...
 
double * K_squared_pt
 Pointer to wave number (must be set!) More...
 
PMLMapping * Pml_mapping_pt
 
- Protected Attributes inherited from oomph::PMLElementBase< DIM >
bool Pml_is_enabled
 Boolean indicating if element is used in pml mode. More...
 
std::vector< bool > Pml_direction_active
 
Vector< double > Pml_inner_boundary
 
Vector< double > Pml_outer_boundary
 
- 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
 
- 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...
 
- Static Protected Attributes inherited from oomph::PMLHelmholtzEquations< DIM >
static double Default_Physical_Constant_Value
 Static default value for the physical constants (initialised to zero) 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
 
- Static Protected Attributes inherited from oomph::RefineableElement
static double Max_integrity_tolerance = 1.0e-8
 Max. allowed discrepancy in element integrity check. More...
 

Detailed Description

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

//////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////// Refineable version of PMLHelmholtz equations

Constructor & Destructor Documentation

◆ RefineablePMLHelmholtzEquations() [1/2]

template<unsigned DIM>
oomph::RefineablePMLHelmholtzEquations< DIM >::RefineablePMLHelmholtzEquations ( )
inline

Constructor, simply call other constructors.

64  : PMLHelmholtzEquations<DIM>(),
65  RefineableElement(),
66  ElementWithZ2ErrorEstimator()
67  {
68  }
oomph::ElementWithZ2ErrorEstimator::ElementWithZ2ErrorEstimator
ElementWithZ2ErrorEstimator()
Default empty constructor.
Definition: error_estimator.h:82
oomph::RefineableElement::RefineableElement
RefineableElement()
Definition: refineable_elements.h:188

◆ RefineablePMLHelmholtzEquations() [2/2]

template<unsigned DIM>
oomph::RefineablePMLHelmholtzEquations< DIM >::RefineablePMLHelmholtzEquations ( const RefineablePMLHelmholtzEquations< DIM > &  dummy)
delete

Broken copy constructor.

Member Function Documentation

◆ fill_in_generic_residual_contribution_helmholtz()

template<unsigned DIM>
void oomph::RefineablePMLHelmholtzEquations< DIM >::fill_in_generic_residual_contribution_helmholtz ( Vector< double > &  residuals,
DenseMatrix< double > &  jacobian,
const unsigned &  flag 
)
privatevirtual

Add element's contribution to elemental residual vector and/or Jacobian matrix flag=1: compute both flag=0: compute only residual vector

Compute element residual Vector and/or element Jacobian matrix

flag=1: compute both flag=0: compute only residual Vector

Pure version without hanging nodes

Reimplemented from oomph::PMLHelmholtzEquations< DIM >.

45  {
46  // Find out how many nodes there are
47  const unsigned n_node = nnode();
48 
49  // Set up memory for the shape and test functions
50  Shape psi(n_node), test(n_node);
51  DShape dpsidx(n_node, DIM), dtestdx(n_node, DIM);
52 
53  // Local storage for pointers to hang_info objects
54  HangInfo *hang_info_pt = 0, *hang_info2_pt = 0;
55 
56  // Set the value of n_intpt
57  const unsigned n_intpt = integral_pt()->nweight();
58 
59  // Integers to store the local equation and unknown numbers
60  int local_eqn_real = 0, local_unknown_real = 0;
61  int local_eqn_imag = 0, local_unknown_imag = 0;
62 
63  // Loop over the integration points
64  for (unsigned ipt = 0; ipt < n_intpt; ipt++)
65  {
66  // Get the integral weight
67  double w = integral_pt()->weight(ipt);
68 
69  // Call the derivatives of the shape and test functions
70  double J = this->dshape_and_dtest_eulerian_at_knot_helmholtz(
71  ipt, psi, dpsidx, test, dtestdx);
72 
73  // Premultiply the weights and the Jacobian
74  double W = w * J;
75 
76  // Calculate local values of unknown
77  // Allocate and initialise to zero
78  std::complex<double> interpolated_u(0.0, 0.0);
79  Vector<double> interpolated_x(DIM, 0.0);
80  Vector<std::complex<double>> interpolated_dudx(DIM);
81 
82  // Calculate function value and derivatives:
83  //-----------------------------------------
84  // Loop over nodes
85  for (unsigned l = 0; l < n_node; l++)
86  {
87  // Loop over directions
88  for (unsigned j = 0; j < DIM; j++)
89  {
90  interpolated_x[j] += nodal_position(l, j) * psi(l);
91  }
92 
93  // Get the nodal value of the helmholtz unknown
94  const std::complex<double> u_value(
95  this->nodal_value(l, this->u_index_helmholtz().real()),
96  this->nodal_value(l, this->u_index_helmholtz().imag()));
97 
98  // Add to the interpolated value
99  interpolated_u += u_value * psi(l);
100 
101  // Loop over directions
102  for (unsigned j = 0; j < DIM; j++)
103  {
104  interpolated_dudx[j] += u_value * dpsidx(l, j);
105  }
106  }
107 
108  // Get source function
109  //-------------------
110  std::complex<double> source(0.0, 0.0);
111  this->get_source_helmholtz(ipt, interpolated_x, source);
112 
113 
114  // Declare a vector of complex numbers for pml weights on the Laplace bit
115  Vector<std::complex<double>> pml_laplace_factor(DIM);
116  // Declare a complex number for pml weights on the mass matrix bit
117  std::complex<double> pml_k_squared_factor =
118  std::complex<double>(1.0, 0.0);
119 
120  // All the PML weights that participate in the assemby process
121  // are computed here. pml_laplace_factor will contain the entries
122  // for the Laplace bit, while pml_k_squared_factor contains the
123  // contributions to the Helmholtz bit. Both default to 1.0, should the PML
124  // not be enabled via enable_pml.
125  this->compute_pml_coefficients(
126  ipt, interpolated_x, pml_laplace_factor, pml_k_squared_factor);
127 
128  // Alpha adjusts the pml factors, the imaginary part produces cross terms
129  std::complex<double> alpha_pml_k_squared_factor =
130  std::complex<double>(pml_k_squared_factor.real() -
131  this->alpha() * pml_k_squared_factor.imag(),
132  this->alpha() * pml_k_squared_factor.real() +
133  pml_k_squared_factor.imag());
134 
135 
136  // std::complex<double> alpha_pml_k_squared_factor
137  // if(alpha_pt() == 0)
138  // {
139  // std::complex<double> alpha_pml_k_squared_factor =
140  // std::complex<double>(
141  // pml_k_squared_factor.real() - alpha() *
142  // pml_k_squared_factor.imag(), alpha() * pml_k_squared_factor.real() +
143  // pml_k_squared_factor.imag()
144  // );
145  // }
146  // Assemble residuals and Jacobian
147  //--------------------------------
148  // Loop over the test functions
149  for (unsigned l = 0; l < n_node; l++)
150  {
151  // Local variables used to store the number of master nodes and the
152  // weight associated with the shape function if the node is hanging
153  unsigned n_master = 1;
154  double hang_weight = 1.0;
155 
156  // Local bool (is the node hanging)
157  bool is_node_hanging = this->node_pt(l)->is_hanging();
158 
159  // If the node is hanging, get the number of master nodes
160  if (is_node_hanging)
161  {
162  hang_info_pt = this->node_pt(l)->hanging_pt();
163  n_master = hang_info_pt->nmaster();
164  }
165  // Otherwise there is just one master node, the node itself
166  else
167  {
168  n_master = 1;
169  }
170 
171  // Loop over the master nodes
172  for (unsigned m = 0; m < n_master; m++)
173  {
174  // Get the local equation number and hang_weight
175  // If the node is hanging
176  if (is_node_hanging)
177  {
178  // Read out the local equation number from the m-th master node
179  local_eqn_real =
180  this->local_hang_eqn(hang_info_pt->master_node_pt(m),
181  this->u_index_helmholtz().real());
182 
183  local_eqn_imag =
184  this->local_hang_eqn(hang_info_pt->master_node_pt(m),
185  this->u_index_helmholtz().imag());
186 
187  // Read out the weight from the master node
188  hang_weight = hang_info_pt->master_weight(m);
189  }
190  // If the node is not hanging
191  else
192  {
193  // The local equation number comes from the node itself
194  local_eqn_real =
195  this->nodal_local_eqn(l, this->u_index_helmholtz().real());
196  local_eqn_imag =
197  this->nodal_local_eqn(l, this->u_index_helmholtz().imag());
198 
199  // The hang weight is one
200  hang_weight = 1.0;
201  }
202 
203  // first, compute the real part contribution
204  //-------------------------------------------
205 
206  /*IF it's not a boundary condition*/
207  if (local_eqn_real >= 0)
208  {
209  // Add body force/source term and Helmholtz bit
210  residuals[local_eqn_real] +=
211  (source.real() - (alpha_pml_k_squared_factor.real() *
212  this->k_squared() * interpolated_u.real() -
213  alpha_pml_k_squared_factor.imag() *
214  this->k_squared() * interpolated_u.imag())) *
215  test(l) * W * hang_weight;
216 
217  // The Laplace bit
218  for (unsigned k = 0; k < DIM; k++)
219  {
220  residuals[local_eqn_real] +=
221  (pml_laplace_factor[k].real() * interpolated_dudx[k].real() -
222  pml_laplace_factor[k].imag() * interpolated_dudx[k].imag()) *
223  dtestdx(l, k) * W * hang_weight;
224  }
225 
226  // Calculate the jacobian
227  //-----------------------
228  if (flag)
229  {
230  // Local variables to store the number of master nodes
231  // and the weights associated with each hanging node
232  unsigned n_master2 = 1;
233  double hang_weight2 = 1.0;
234 
235  // Loop over the nodes for the variables
236  for (unsigned l2 = 0; l2 < n_node; l2++)
237  {
238  // Local bool (is the node hanging)
239  bool is_node2_hanging = this->node_pt(l2)->is_hanging();
240 
241  // If the node is hanging, get the number of master nodes
242  if (is_node2_hanging)
243  {
244  hang_info2_pt = this->node_pt(l2)->hanging_pt();
245  n_master2 = hang_info2_pt->nmaster();
246  }
247  // Otherwise there is one master node, the node itself
248  else
249  {
250  n_master2 = 1;
251  }
252 
253  // Loop over the master nodes
254  for (unsigned m2 = 0; m2 < n_master2; m2++)
255  {
256  // Get the local unknown and weight
257  // If the node is hanging
258  if (is_node2_hanging)
259  {
260  // Read out the local unknown from the master node
261  local_unknown_real =
262  this->local_hang_eqn(hang_info2_pt->master_node_pt(m2),
263  this->u_index_helmholtz().real());
264  local_unknown_imag =
265  this->local_hang_eqn(hang_info2_pt->master_node_pt(m2),
266  this->u_index_helmholtz().imag());
267 
268  // Read out the hanging weight from the master node
269  hang_weight2 = hang_info2_pt->master_weight(m2);
270  }
271  // If the node is not hanging
272  else
273  {
274  // The local unknown number comes from the node
275  local_unknown_real = this->nodal_local_eqn(
276  l2, this->u_index_helmholtz().real());
277 
278  local_unknown_imag = this->nodal_local_eqn(
279  l2, this->u_index_helmholtz().imag());
280 
281  // The hang weight is one
282  hang_weight2 = 1.0;
283  }
284 
285 
286  // If at a non-zero degree of freedom add in the entry
287  if (local_unknown_real >= 0)
288  {
289  // Add contribution to Elemental Matrix
290  for (unsigned i = 0; i < DIM; i++)
291  {
292  jacobian(local_eqn_real, local_unknown_real) +=
293  pml_laplace_factor[i].real() * dpsidx(l2, i) *
294  dtestdx(l, i) * W * hang_weight * hang_weight2;
295  }
296  // Add the helmholtz contribution
297  jacobian(local_eqn_real, local_unknown_real) +=
298  -alpha_pml_k_squared_factor.real() * this->k_squared() *
299  psi(l2) * test(l) * W * hang_weight * hang_weight2;
300  }
301  // If at a non-zero degree of freedom add in the entry
302  if (local_unknown_imag >= 0)
303  {
304  // Add contribution to Elemental Matrix
305  for (unsigned i = 0; i < DIM; i++)
306  {
307  jacobian(local_eqn_real, local_unknown_imag) -=
308  pml_laplace_factor[i].imag() * dpsidx(l2, i) *
309  dtestdx(l, i) * W * hang_weight * hang_weight2;
310  }
311  // Add the helmholtz contribution
312  jacobian(local_eqn_real, local_unknown_imag) +=
313  alpha_pml_k_squared_factor.imag() * this->k_squared() *
314  psi(l2) * test(l) * W * hang_weight * hang_weight2;
315  }
316  }
317  }
318  }
319  }
320 
321  // Second, compute the imaginary part contribution
322  //------------------------------------------------
323 
324  /*IF it's not a boundary condition*/
325  if (local_eqn_imag >= 0)
326  {
327  // Add body force/source term and Helmholtz bit
328  residuals[local_eqn_imag] +=
329  (source.imag() - (alpha_pml_k_squared_factor.imag() *
330  this->k_squared() * interpolated_u.real() +
331  alpha_pml_k_squared_factor.real() *
332  this->k_squared() * interpolated_u.imag())) *
333  test(l) * W * hang_weight;
334 
335  // The Laplace bit
336  for (unsigned k = 0; k < DIM; k++)
337  {
338  residuals[local_eqn_imag] +=
339  (pml_laplace_factor[k].imag() * interpolated_dudx[k].real() +
340  pml_laplace_factor[k].real() * interpolated_dudx[k].imag()) *
341  dtestdx(l, k) * W * hang_weight;
342  }
343 
344  // Calculate the jacobian
345  //-----------------------
346  if (flag)
347  {
348  // Local variables to store the number of master nodes
349  // and the weights associated with each hanging node
350  unsigned n_master2 = 1;
351  double hang_weight2 = 1.0;
352 
353  // Loop over the nodes for the variables
354  for (unsigned l2 = 0; l2 < n_node; l2++)
355  {
356  // Local bool (is the node hanging)
357  bool is_node2_hanging = this->node_pt(l2)->is_hanging();
358 
359  // If the node is hanging, get the number of master nodes
360  if (is_node2_hanging)
361  {
362  hang_info2_pt = this->node_pt(l2)->hanging_pt();
363  n_master2 = hang_info2_pt->nmaster();
364  }
365  // Otherwise there is one master node, the node itself
366  else
367  {
368  n_master2 = 1;
369  }
370 
371  // Loop over the master nodes
372  for (unsigned m2 = 0; m2 < n_master2; m2++)
373  {
374  // Get the local unknown and weight
375  // If the node is hanging
376  if (is_node2_hanging)
377  {
378  // Read out the local unknown from the master node
379  local_unknown_real =
380  this->local_hang_eqn(hang_info2_pt->master_node_pt(m2),
381  this->u_index_helmholtz().real());
382  local_unknown_imag =
383  this->local_hang_eqn(hang_info2_pt->master_node_pt(m2),
384  this->u_index_helmholtz().imag());
385 
386  // Read out the hanging weight from the master node
387  hang_weight2 = hang_info2_pt->master_weight(m2);
388  }
389  // If the node is not hanging
390  else
391  {
392  // The local unknown number comes from the node
393  local_unknown_real = this->nodal_local_eqn(
394  l2, this->u_index_helmholtz().real());
395 
396  local_unknown_imag = this->nodal_local_eqn(
397  l2, this->u_index_helmholtz().imag());
398 
399  // The hang weight is one
400  hang_weight2 = 1.0;
401  }
402 
403  // If at a non-zero degree of freedom add in the entry
404  if (local_unknown_imag >= 0)
405  {
406  // Add contribution to Elemental Matrix
407  for (unsigned i = 0; i < DIM; i++)
408  {
409  jacobian(local_eqn_imag, local_unknown_imag) +=
410  pml_laplace_factor[i].real() * dpsidx(l2, i) *
411  dtestdx(l, i) * W * hang_weight * hang_weight2;
412  }
413  // Add the helmholtz contribution
414  jacobian(local_eqn_imag, local_unknown_imag) +=
415  -alpha_pml_k_squared_factor.real() * this->k_squared() *
416  psi(l2) * test(l) * W * hang_weight * hang_weight2;
417  }
418  if (local_unknown_real >= 0)
419  {
420  // Add contribution to Elemental Matrix
421  for (unsigned i = 0; i < DIM; i++)
422  {
423  jacobian(local_eqn_imag, local_unknown_real) +=
424  pml_laplace_factor[i].imag() * dpsidx(l2, i) *
425  dtestdx(l, i) * W * hang_weight * hang_weight2;
426  }
427  // Add the helmholtz contribution
428  jacobian(local_eqn_imag, local_unknown_real) +=
429  -alpha_pml_k_squared_factor.imag() * this->k_squared() *
430  psi(l2) * test(l) * W * hang_weight * hang_weight2;
431  }
432  }
433  }
434  }
435  }
436  }
437  }
438 
439  } // End of loop over integration points
440  }
imag
AnnoyingScalar imag(const AnnoyingScalar &)
Definition: AnnoyingScalar.h:132
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9
J
JacobiRotation< float > J
Definition: Jacobi_makeJacobi.cpp:3
w
RowVector3d w
Definition: Matrix_resize_int.cpp:3
m2
MatrixType m2(n_dims)
oomph::FiniteElement::node_pt
Node *& node_pt(const unsigned &n)
Return a pointer to the local node n.
Definition: elements.h:2175
oomph::FiniteElement::nodal_value
double nodal_value(const unsigned &n, const unsigned &i) const
Definition: elements.h:2593
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::nodal_local_eqn
int nodal_local_eqn(const unsigned &n, const unsigned &i) const
Definition: elements.h:1432
oomph::FiniteElement::nnode
unsigned nnode() const
Return the number of nodes.
Definition: elements.h:2210
oomph::FiniteElement::integral_pt
Integral *const & integral_pt() const
Return the pointer to the integration scheme (const version)
Definition: elements.h:1963
oomph::FiniteElement::nodal_position
double nodal_position(const unsigned &n, const unsigned &i) const
Definition: elements.h:2317
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::Node::hanging_pt
HangInfo *const & hanging_pt() const
Definition: nodes.h:1228
oomph::Node::is_hanging
bool is_hanging() const
Test whether the node is geometrically hanging.
Definition: nodes.h:1285
oomph::PMLHelmholtzEquations::dshape_and_dtest_eulerian_at_knot_helmholtz
virtual double dshape_and_dtest_eulerian_at_knot_helmholtz(const unsigned &ipt, Shape &psi, DShape &dpsidx, Shape &test, DShape &dtestdx) const =0
oomph::PMLHelmholtzEquations::get_source_helmholtz
virtual void get_source_helmholtz(const unsigned &ipt, const Vector< double > &x, std::complex< double > &source) const
Definition: pml_helmholtz_elements.h:416
oomph::PMLHelmholtzEquations::u_index_helmholtz
virtual std::complex< unsigned > u_index_helmholtz() const
Broken assignment operator.
Definition: pml_helmholtz_elements.h:91
oomph::PMLHelmholtzEquations::k_squared
double k_squared()
Get the square of wavenumber.
Definition: pml_helmholtz_elements.h:104
oomph::PMLHelmholtzEquations::compute_pml_coefficients
void compute_pml_coefficients(const unsigned &ipt, const Vector< double > &x, Vector< std::complex< double >> &pml_laplace_factor, std::complex< double > &pml_k_squared_factor)
Definition: pml_helmholtz_elements.h:546
oomph::RefineableElement::local_hang_eqn
int local_hang_eqn(Node *const &node_pt, const unsigned &i)
Definition: refineable_elements.h:278
real
float real
Definition: datatypes.h:10
m
int * m
Definition: level2_cplx_impl.h:294
k
char char char int int * k
Definition: level2_impl.h:374
DIM
#define DIM
Definition: linearised_navier_stokes_elements.h:44
Eigen::test
squared absolute sa ArrayBase::abs2 DOXCOMMA MatrixBase::cwiseAbs2 sa Eigen::abs2 DOXCOMMA Eigen::pow DOXCOMMA ArrayBase::square nearest sa Eigen::floor DOXCOMMA Eigen::ceil DOXCOMMA ArrayBase::round nearest integer not less than the given sa Eigen::floor DOXCOMMA ArrayBase::ceil not a number test
Definition: GlobalFunctions.h:109
TestProblem::source
void source(const Vector< double > &x, Vector< double > &f)
Source function.
Definition: unstructured_two_d_circle.cc:46
oomph::QuadTreeNames::W
@ W
Definition: quadtree.h:63
test
Definition: indexed_view.cpp:20
j
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2

References DIM, i, imag(), J, j, k, m, m2(), oomph::HangInfo::master_node_pt(), oomph::HangInfo::master_weight(), oomph::HangInfo::nmaster(), TestProblem::source(), Eigen::test, w, and oomph::QuadTreeNames::W.

◆ further_build()

template<unsigned DIM>
void oomph::RefineablePMLHelmholtzEquations< DIM >::further_build ( )
inlinevirtual

Further build: Copy source function pointer from father element.

Reimplemented from oomph::RefineableElement.

170  {
171  this->Source_fct_pt = dynamic_cast<RefineablePMLHelmholtzEquations<DIM>*>(
172  this->father_element_pt())
173  ->source_fct_pt();
174  }
oomph::PMLHelmholtzEquations::Source_fct_pt
PMLHelmholtzSourceFctPt Source_fct_pt
Pointer to source function:
Definition: pml_helmholtz_elements.h:727
oomph::PMLHelmholtzEquations::source_fct_pt
PMLHelmholtzSourceFctPt & source_fct_pt()
Access function: Pointer to source function.
Definition: pml_helmholtz_elements.h:400
oomph::RefineableElement::father_element_pt
virtual RefineableElement * father_element_pt() const
Return a pointer to the father element.
Definition: refineable_elements.h:539

References oomph::RefineableElement::father_element_pt(), oomph::PMLHelmholtzEquations< DIM >::source_fct_pt(), and oomph::PMLHelmholtzEquations< DIM >::Source_fct_pt.

◆ get_interpolated_values() [1/2]

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

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

Implements oomph::RefineableElement.

147  {
148  if (t != 0)
149  {
150  std::string error_message =
151  "Time-dependent version of get_interpolated_values() ";
152  error_message += "not implemented for this element \n";
153  throw OomphLibError(
154  error_message,
155  "RefineablePMLHelmholtzEquations::get_interpolated_values()",
156  OOMPH_EXCEPTION_LOCATION);
157  }
158  else
159  {
160  // Make sure that we call this particular object's steady
161  // get_interpolated_values (it could get overloaded lower down)
162  RefineablePMLHelmholtzEquations<DIM>::get_interpolated_values(s,
163  values);
164  }
165  }
oomph::RefineablePMLHelmholtzEquations::get_interpolated_values
void get_interpolated_values(const Vector< double > &s, Vector< double > &values)
Definition: refineable_pml_helmholtz_elements.h:109
s
RealScalar s
Definition: level1_cplx_impl.h:130
oomph::Global_string_for_annotation::string
std::string string(const unsigned &i)
Definition: oomph_definitions.cc:286
plotPSD.t
t
Definition: plotPSD.py:36
OOMPH_EXCEPTION_LOCATION
#define OOMPH_EXCEPTION_LOCATION
Definition: oomph_definitions.h:61

References oomph::RefineablePMLHelmholtzEquations< DIM >::get_interpolated_values(), OOMPH_EXCEPTION_LOCATION, s, oomph::Global_string_for_annotation::string(), and plotPSD::t.

◆ get_interpolated_values() [2/2]

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

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

Reimplemented from oomph::RefineableElement.

111  {
112  // Set size of Vector: u
113  values.resize(2);
114 
115  // Find number of nodes
116  unsigned n_node = nnode();
117 
118  // Local shape function
119  Shape psi(n_node);
120 
121  // Find values of shape function
122  shape(s, psi);
123 
124  // Initialise value of u
125  values[0] = 0.0;
126  values[1] = 0.0;
127 
128  // Find the index at which the pml_helmholtz unknown is stored
129  std::complex<unsigned> u_nodal_index = this->u_index_helmholtz();
130 
131  // Loop over the local nodes and sum up the values
132  for (unsigned l = 0; l < n_node; l++)
133  {
134  values[0] += this->nodal_value(l, u_nodal_index.real()) * psi[l];
135  values[1] += this->nodal_value(l, u_nodal_index.imag()) * psi[l];
136  }
137  }
oomph::FiniteElement::shape
virtual void shape(const Vector< double > &s, Shape &psi) const =0

References oomph::FiniteElement::nnode(), oomph::FiniteElement::nodal_value(), s, oomph::FiniteElement::shape(), and oomph::PMLHelmholtzEquations< DIM >::u_index_helmholtz().

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

◆ get_Z2_flux()

template<unsigned DIM>
void oomph::RefineablePMLHelmholtzEquations< DIM >::get_Z2_flux ( const Vector< double > &  s,
Vector< double > &  flux 
)
inlinevirtual

Get 'flux' for Z2 error recovery: Complex flux from PMLHelmholtz equations, strung together

Implements oomph::ElementWithZ2ErrorEstimator.

91  {
92  Vector<std::complex<double>> complex_flux(DIM);
93  this->get_flux(s, complex_flux);
94  unsigned count = 0;
95  for (unsigned i = 0; i < DIM; i++)
96  {
97  flux[count] = complex_flux[i].real();
98  count++;
99  flux[count] = complex_flux[i].imag();
100  count++;
101  }
102  }
oomph::PMLHelmholtzEquations::get_flux
void get_flux(const Vector< double > &s, Vector< std::complex< double >> &flux) const
Get flux: flux[i] = du/dx_i for real and imag part.
Definition: pml_helmholtz_elements.h:446
ProblemParameters::flux
void flux(const double &time, const Vector< double > &x, double &flux)
Get flux applied along boundary x=0.
Definition: pretend_melt.cc:59

References DIM, ProblemParameters::flux(), oomph::PMLHelmholtzEquations< DIM >::get_flux(), and i.

◆ num_Z2_flux_terms()

template<unsigned DIM>
unsigned oomph::RefineablePMLHelmholtzEquations< DIM >::num_Z2_flux_terms ( )
inlinevirtual

Broken assignment operator.

Number of 'flux' terms for Z2 error estimation

Implements oomph::ElementWithZ2ErrorEstimator.

84  {
85  return 2 * DIM;
86  }

References DIM.

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