Class definition. More...

Inheritance diagram for UnstructuredFvKProblem< ELEMENT >:

Public Member Functions
	UnstructuredFvKProblem (double element_area=0.2)
	Constructor. More...

	~UnstructuredFvKProblem ()
	Destructor. More...

void	actions_after_newton_solve ()
	Update after solve (empty) More...

void	actions_before_newton_solve ()
	Update the problem specs before solve (empty) More...

void	doc_solution ()
	Doc the solution. More...

	UnstructuredFvKProblem (double element_area=0.2)
	Constructor. More...

	~UnstructuredFvKProblem ()
	Destructor. More...

void	actions_after_newton_solve ()
	Update after solve (empty) More...

void	actions_before_newton_solve ()
	Update the problem specs before solve (empty) More...

void	doc_solution ()
	Doc the solution. More...

	UnstructuredFvKProblem (double element_area=0.2)
	Constructor. More...

	~UnstructuredFvKProblem ()
	Destructor. More...

void	actions_before_adapt ()

void	actions_after_adapt ()

void	actions_after_newton_solve ()
	Update after solve (empty) More...

void	actions_before_newton_solve ()
	Update the problem specs before solve: Re-apply boundary conditons. More...

void	doc_solution (const std::string &comment="")
	Doc the solution. More...

Public Member Functions inherited from oomph::Problem
virtual void	debug_hook_fct (const unsigned &i)

void	set_analytic_dparameter (double *const &parameter_pt)

void	unset_analytic_dparameter (double *const &parameter_pt)

bool	is_dparameter_calculated_analytically (double *const &parameter_pt)

void	set_analytic_hessian_products ()

void	unset_analytic_hessian_products ()

bool	are_hessian_products_calculated_analytically ()

void	set_pinned_values_to_zero ()

bool	distributed () const

OomphCommunicator *	communicator_pt ()
	access function to the oomph-lib communicator More...

const OomphCommunicator *	communicator_pt () const
	access function to the oomph-lib communicator, const version More...

	Problem ()

	Problem (const Problem &dummy)=delete
	Broken copy constructor. More...

void	operator= (const Problem &)=delete
	Broken assignment operator. More...

virtual	~Problem ()
	Virtual destructor to clean up memory. More...

Mesh *&	mesh_pt ()
	Return a pointer to the global mesh. More...

Mesh *const &	mesh_pt () const
	Return a pointer to the global mesh (const version) More...

Mesh *&	mesh_pt (const unsigned &imesh)

Mesh *const &	mesh_pt (const unsigned &imesh) const
	Return a pointer to the i-th submesh (const version) More...

unsigned	nsub_mesh () const
	Return number of submeshes. More...

unsigned	add_sub_mesh (Mesh *const &mesh_pt)

void	flush_sub_meshes ()

void	build_global_mesh ()

void	rebuild_global_mesh ()

LinearSolver *&	linear_solver_pt ()
	Return a pointer to the linear solver object. More...

LinearSolver *const &	linear_solver_pt () const
	Return a pointer to the linear solver object (const version) More...

LinearSolver *&	mass_matrix_solver_for_explicit_timestepper_pt ()

LinearSolver *	mass_matrix_solver_for_explicit_timestepper_pt () const

EigenSolver *&	eigen_solver_pt ()
	Return a pointer to the eigen solver object. More...

EigenSolver *const &	eigen_solver_pt () const
	Return a pointer to the eigen solver object (const version) More...

Time *&	time_pt ()
	Return a pointer to the global time object. More...

Time *	time_pt () const
	Return a pointer to the global time object (const version). More...

double &	time ()
	Return the current value of continuous time. More...

double	time () const
	Return the current value of continuous time (const version) More...

TimeStepper *&	time_stepper_pt ()

const TimeStepper *	time_stepper_pt () const

TimeStepper *&	time_stepper_pt (const unsigned &i)
	Return a pointer to the i-th timestepper. More...

ExplicitTimeStepper *&	explicit_time_stepper_pt ()
	Return a pointer to the explicit timestepper. More...

unsigned long	set_timestepper_for_all_data (TimeStepper *const &time_stepper_pt, const bool &preserve_existing_data=false)

virtual void	shift_time_values ()
	Shift all values along to prepare for next timestep. More...

AssemblyHandler *&	assembly_handler_pt ()
	Return a pointer to the assembly handler object. More...

AssemblyHandler *const &	assembly_handler_pt () const
	Return a pointer to the assembly handler object (const version) More...

double &	minimum_dt ()
	Access function to min timestep in adaptive timestepping. More...

double &	maximum_dt ()
	Access function to max timestep in adaptive timestepping. More...

unsigned &	max_newton_iterations ()
	Access function to max Newton iterations before giving up. More...

void	problem_is_nonlinear (const bool &prob_lin)
	Access function to Problem_is_nonlinear. More...

double &	max_residuals ()

bool &	time_adaptive_newton_crash_on_solve_fail ()
	Access function for Time_adaptive_newton_crash_on_solve_fail. More...

double &	newton_solver_tolerance ()

void	add_time_stepper_pt (TimeStepper *const &time_stepper_pt)

void	set_explicit_time_stepper_pt (ExplicitTimeStepper *const &explicit_time_stepper_pt)

void	initialise_dt (const double &dt)

void	initialise_dt (const Vector< double > &dt)

Data *&	global_data_pt (const unsigned &i)
	Return a pointer to the the i-th global data object. More...

void	add_global_data (Data *const &global_data_pt)

void	flush_global_data ()

LinearAlgebraDistribution *const &	dof_distribution_pt () const
	Return the pointer to the dof distribution (read-only) More...

unsigned long	ndof () const
	Return the number of dofs. More...

unsigned	ntime_stepper () const
	Return the number of time steppers. More...

unsigned	nglobal_data () const
	Return the number of global data values. More...

unsigned	self_test ()
	Self-test: Check meshes and global data. Return 0 for OK. More...

void	enable_store_local_dof_pt_in_elements ()

void	disable_store_local_dof_pt_in_elements ()

unsigned long	assign_eqn_numbers (const bool &assign_local_eqn_numbers=true)

void	describe_dofs (std::ostream &out= *(oomph_info.stream_pt())) const

void	enable_discontinuous_formulation ()

void	disable_discontinuous_formulation ()

void	get_dofs (DoubleVector &dofs) const

void	get_dofs (const unsigned &t, DoubleVector &dofs) const
	Return vector of the t'th history value of all dofs. More...

void	set_dofs (const DoubleVector &dofs)
	Set the values of the dofs. More...

void	set_dofs (const unsigned &t, DoubleVector &dofs)
	Set the history values of the dofs. More...

void	set_dofs (const unsigned &t, Vector< double * > &dof_pt)

void	add_to_dofs (const double &lambda, const DoubleVector &increment_dofs)
	Add lambda x incremenet_dofs[l] to the l-th dof. More...

double *	global_dof_pt (const unsigned &i)

double &	dof (const unsigned &i)
	i-th dof in the problem More...

double	dof (const unsigned &i) const
	i-th dof in the problem (const version) More...

double *&	dof_pt (const unsigned &i)
	Pointer to i-th dof in the problem. More...

double *	dof_pt (const unsigned &i) const
	Pointer to i-th dof in the problem (const version) More...

virtual void	get_inverse_mass_matrix_times_residuals (DoubleVector &Mres)

virtual void	get_dvaluesdt (DoubleVector &f)

virtual void	get_residuals (DoubleVector &residuals)
	Get the total residuals Vector for the problem. More...

virtual void	get_jacobian (DoubleVector &residuals, DenseDoubleMatrix &jacobian)

virtual void	get_jacobian (DoubleVector &residuals, CRDoubleMatrix &jacobian)

virtual void	get_jacobian (DoubleVector &residuals, CCDoubleMatrix &jacobian)

virtual void	get_jacobian (DoubleVector &residuals, SumOfMatrices &jacobian)

void	get_fd_jacobian (DoubleVector &residuals, DenseMatrix< double > &jacobian)
	Get the full Jacobian by finite differencing. More...

void	get_derivative_wrt_global_parameter (double *const &parameter_pt, DoubleVector &result)

void	get_hessian_vector_products (DoubleVectorWithHaloEntries const &Y, Vector< DoubleVectorWithHaloEntries > const &C, Vector< DoubleVectorWithHaloEntries > &product)

void	solve_eigenproblem (const unsigned &n_eval, Vector< std::complex< double >> &eigenvalue, Vector< DoubleVector > &eigenvector, const bool &steady=true)
	Solve the eigenproblem. More...

void	solve_eigenproblem (const unsigned &n_eval, Vector< std::complex< double >> &eigenvalue, const bool &steady=true)

virtual void	get_eigenproblem_matrices (CRDoubleMatrix &mass_matrix, CRDoubleMatrix &main_matrix, const double &shift=0.0)

void	assign_eigenvector_to_dofs (DoubleVector &eigenvector)
	Assign the eigenvector passed to the function to the dofs. More...

void	add_eigenvector_to_dofs (const double &epsilon, const DoubleVector &eigenvector)

void	store_current_dof_values ()
	Store the current values of the degrees of freedom. More...

void	restore_dof_values ()
	Restore the stored values of the degrees of freedom. More...

void	enable_jacobian_reuse ()

void	disable_jacobian_reuse ()
	Disable recycling of Jacobian in Newton iteration. More...

bool	jacobian_reuse_is_enabled ()
	Is recycling of Jacobian in Newton iteration enabled? More...

bool &	use_predictor_values_as_initial_guess ()

void	newton_solve ()
	Use Newton method to solve the problem. More...

void	enable_globally_convergent_newton_method ()
	enable globally convergent Newton method More...

void	disable_globally_convergent_newton_method ()
	disable globally convergent Newton method More...

void	newton_solve (unsigned const &max_adapt)

void	steady_newton_solve (unsigned const &max_adapt=0)

void	copy (Problem *orig_problem_pt)

virtual Problem *	make_copy ()

virtual void	read (std::ifstream &restart_file, bool &unsteady_restart)

virtual void	read (std::ifstream &restart_file)

virtual void	dump (std::ofstream &dump_file) const

void	dump (const std::string &dump_file_name) const

void	delete_all_external_storage ()

virtual void	symmetrise_eigenfunction_for_adaptive_pitchfork_tracking ()

double *	bifurcation_parameter_pt () const

void	get_bifurcation_eigenfunction (Vector< DoubleVector > &eigenfunction)

void	activate_fold_tracking (double *const &parameter_pt, const bool &block_solve=true)

void	activate_bifurcation_tracking (double *const &parameter_pt, const DoubleVector &eigenvector, const bool &block_solve=true)

void	activate_bifurcation_tracking (double *const &parameter_pt, const DoubleVector &eigenvector, const DoubleVector &normalisation, const bool &block_solve=true)

void	activate_pitchfork_tracking (double *const &parameter_pt, const DoubleVector &symmetry_vector, const bool &block_solve=true)

void	activate_hopf_tracking (double *const &parameter_pt, const bool &block_solve=true)

void	activate_hopf_tracking (double *const &parameter_pt, const double &omega, const DoubleVector &null_real, const DoubleVector &null_imag, const bool &block_solve=true)

void	deactivate_bifurcation_tracking ()

void	reset_assembly_handler_to_default ()
	Reset the system to the standard non-augemented state. More...

double	arc_length_step_solve (double *const &parameter_pt, const double &ds, const unsigned &max_adapt=0)

double	arc_length_step_solve (Data *const &data_pt, const unsigned &data_index, const double &ds, const unsigned &max_adapt=0)

void	reset_arc_length_parameters ()

int &	sign_of_jacobian ()

void	explicit_timestep (const double &dt, const bool &shift_values=true)
	Take an explicit timestep of size dt. More...

void	unsteady_newton_solve (const double &dt)

void	unsteady_newton_solve (const double &dt, const bool &shift_values)

void	unsteady_newton_solve (const double &dt, const unsigned &max_adapt, const bool &first, const bool &shift=true)

double	doubly_adaptive_unsteady_newton_solve (const double &dt, const double &epsilon, const unsigned &max_adapt, const bool &first, const bool &shift=true)

double	doubly_adaptive_unsteady_newton_solve (const double &dt, const double &epsilon, const unsigned &max_adapt, const unsigned &suppress_resolve_after_spatial_adapt_flag, const bool &first, const bool &shift=true)

double	adaptive_unsteady_newton_solve (const double &dt_desired, const double &epsilon)

double	adaptive_unsteady_newton_solve (const double &dt_desired, const double &epsilon, const bool &shift_values)

void	assign_initial_values_impulsive ()

void	assign_initial_values_impulsive (const double &dt)

void	calculate_predictions ()
	Calculate predictions. More...

void	enable_mass_matrix_reuse ()

void	disable_mass_matrix_reuse ()

bool	mass_matrix_reuse_is_enabled ()
	Return whether the mass matrix is being reused. More...

void	refine_uniformly (const Vector< unsigned > &nrefine_for_mesh)

void	refine_uniformly (const Vector< unsigned > &nrefine_for_mesh, DocInfo &doc_info)

void	refine_uniformly_and_prune (const Vector< unsigned > &nrefine_for_mesh)

void	refine_uniformly_and_prune (const Vector< unsigned > &nrefine_for_mesh, DocInfo &doc_info)

void	refine_uniformly (DocInfo &doc_info)

void	refine_uniformly_and_prune (DocInfo &doc_info)

void	refine_uniformly ()

void	refine_uniformly (const unsigned &i_mesh, DocInfo &doc_info)
	Do uniform refinement for submesh i_mesh with documentation. More...

void	refine_uniformly (const unsigned &i_mesh)
	Do uniform refinement for submesh i_mesh without documentation. More...

void	p_refine_uniformly (const Vector< unsigned > &nrefine_for_mesh)

void	p_refine_uniformly (const Vector< unsigned > &nrefine_for_mesh, DocInfo &doc_info)

void	p_refine_uniformly_and_prune (const Vector< unsigned > &nrefine_for_mesh)

void	p_refine_uniformly_and_prune (const Vector< unsigned > &nrefine_for_mesh, DocInfo &doc_info)

void	p_refine_uniformly (DocInfo &doc_info)

void	p_refine_uniformly_and_prune (DocInfo &doc_info)

void	p_refine_uniformly ()

void	p_refine_uniformly (const unsigned &i_mesh, DocInfo &doc_info)
	Do uniform p-refinement for submesh i_mesh with documentation. More...

void	p_refine_uniformly (const unsigned &i_mesh)
	Do uniform p-refinement for submesh i_mesh without documentation. More...

void	refine_selected_elements (const Vector< unsigned > &elements_to_be_refined)

void	refine_selected_elements (const Vector< RefineableElement * > &elements_to_be_refined_pt)

void	refine_selected_elements (const unsigned &i_mesh, const Vector< unsigned > &elements_to_be_refined)

void	refine_selected_elements (const unsigned &i_mesh, const Vector< RefineableElement * > &elements_to_be_refined_pt)

void	refine_selected_elements (const Vector< Vector< unsigned >> &elements_to_be_refined)

void	refine_selected_elements (const Vector< Vector< RefineableElement * >> &elements_to_be_refined_pt)

void	p_refine_selected_elements (const Vector< unsigned > &elements_to_be_refined)

void	p_refine_selected_elements (const Vector< PRefineableElement * > &elements_to_be_refined_pt)

void	p_refine_selected_elements (const unsigned &i_mesh, const Vector< unsigned > &elements_to_be_refined)

void	p_refine_selected_elements (const unsigned &i_mesh, const Vector< PRefineableElement * > &elements_to_be_refined_pt)

void	p_refine_selected_elements (const Vector< Vector< unsigned >> &elements_to_be_refined)

void	p_refine_selected_elements (const Vector< Vector< PRefineableElement * >> &elements_to_be_refined_pt)

unsigned	unrefine_uniformly ()

unsigned	unrefine_uniformly (const unsigned &i_mesh)

void	p_unrefine_uniformly (DocInfo &doc_info)

void	p_unrefine_uniformly (const unsigned &i_mesh, DocInfo &doc_info)
	Do uniform p-unrefinement for submesh i_mesh without documentation. More...

void	adapt (unsigned &n_refined, unsigned &n_unrefined)

void	adapt ()

void	p_adapt (unsigned &n_refined, unsigned &n_unrefined)

void	p_adapt ()

void	adapt_based_on_error_estimates (unsigned &n_refined, unsigned &n_unrefined, Vector< Vector< double >> &elemental_error)

void	adapt_based_on_error_estimates (Vector< Vector< double >> &elemental_error)

void	get_all_error_estimates (Vector< Vector< double >> &elemental_error)

void	doc_errors (DocInfo &doc_info)
	Get max and min error for all elements in submeshes. More...

void	doc_errors ()
	Get max and min error for all elements in submeshes. More...

void	enable_info_in_newton_solve ()

void	disable_info_in_newton_solve ()
	Disable the output of information when in the newton solver. More...

Public Member Functions inherited from oomph::ExplicitTimeSteppableObject
	ExplicitTimeSteppableObject ()
	Empty constructor. More...

	ExplicitTimeSteppableObject (const ExplicitTimeSteppableObject &)=delete
	Broken copy constructor. More...

void	operator= (const ExplicitTimeSteppableObject &)=delete
	Broken assignment operator. More...

virtual	~ExplicitTimeSteppableObject ()
	Empty destructor. More...

virtual void	actions_before_explicit_stage ()

virtual void	actions_after_explicit_stage ()

Private Types
enum	{ Outer_boundary0 = 0 , Outer_boundary1 = 1 , Inner_boundary0 = 2 , Inner_boundary1 = 3 , Inner_boundary2 = 4 , Inner_boundary3 = 5 , Inner_boundary4 = 6 , Inner_boundary5 = 7 }

enum	{ Outer_boundary0 = 0 , Outer_boundary1 = 1 , Inner_boundary0 = 2 , Inner_boundary1 = 3 , Inner_boundary2 = 4 , Inner_boundary3 = 5 , Inner_boundary4 = 6 , Inner_boundary5 = 7 }

enum	{ Outer_boundary0 = 0 , Outer_boundary1 = 1 , Inner_boundary0 = 2 , Inner_boundary1 = 3 , Inner_boundary2 = 4 , Inner_boundary3 = 5 , Inner_boundary4 = 6 , Inner_boundary5 = 7 }

Private Member Functions
void	complete_problem_setup ()

void	complete_problem_setup ()

void	apply_boundary_conditions ()
	Helper function to apply boundary conditions. More...

void	complete_problem_setup ()

Private Attributes
DocInfo	Doc_info
	Doc info object for labeling output. More...

TriangleMesh< ELEMENT > *	My_mesh_pt
	Pointers to specific mesh. More...

ofstream	Trace_file
	Trace file to document norm of solution. More...

double	Element_area
	Element area. More...

MeshAsGeomObject *	Mesh_as_geom_object_pt
	Geom object made of mesh. More...

Vector< std::pair< GeomObject *, Vector< double > > >	Radial_sample_point_pt
	Plot points along radial line are located in this element/local coord. More...

RefineableTriangleMesh< ELEMENT > *	My_mesh_pt
	Pointers to specific mesh. More...

Mesh *	Volume_constraint_mesh_pt
	Mesh for the volume constraint element. More...

FoepplvonKarmanVolumeConstraintElement< ELEMENT, RefineableTriangleMesh > *	Volume_constraint_element_pt
	Single volume constraint element instance. More...

double *	Temp_pressure_pt

Vector< unsigned >	Bubble_regions

Additional Inherited Members
Public Types inherited from oomph::Problem
typedef void(*	SpatialErrorEstimatorFctPt) (Mesh *&mesh_pt, Vector< double > &elemental_error)
	Function pointer for spatial error estimator. More...

typedef void(*	SpatialErrorEstimatorWithDocFctPt) (Mesh *&mesh_pt, Vector< double > &elemental_error, DocInfo &doc_info)
	Function pointer for spatial error estimator with doc. More...

Public Attributes inherited from oomph::Problem
bool	Shut_up_in_newton_solve

Static Public Attributes inherited from oomph::Problem
static bool	Suppress_warning_about_actions_before_read_unstructured_meshes

Protected Types inherited from oomph::Problem
enum	Assembly_method { Perform_assembly_using_vectors_of_pairs , Perform_assembly_using_two_vectors , Perform_assembly_using_maps , Perform_assembly_using_lists , Perform_assembly_using_two_arrays }
	Enumerated flags to determine which sparse assembly method is used. More...

Protected Member Functions inherited from oomph::Problem
unsigned	setup_element_count_per_dof ()

virtual void	sparse_assemble_row_or_column_compressed (Vector< int * > &column_or_row_index, Vector< int * > &row_or_column_start, Vector< double * > &value, Vector< unsigned > &nnz, Vector< double * > &residual, bool compressed_row_flag)

virtual void	actions_before_newton_convergence_check ()

virtual void	actions_before_newton_step ()

virtual void	actions_after_newton_step ()

virtual void	actions_before_implicit_timestep ()

virtual void	actions_after_implicit_timestep ()

virtual void	actions_after_implicit_timestep_and_error_estimation ()

virtual void	actions_before_explicit_timestep ()
	Actions that should be performed before each explicit time step. More...

virtual void	actions_after_explicit_timestep ()
	Actions that should be performed after each explicit time step. More...

virtual void	actions_before_read_unstructured_meshes ()

virtual void	actions_after_read_unstructured_meshes ()

virtual void	actions_after_change_in_global_parameter (double *const &parameter_pt)

virtual void	actions_after_change_in_bifurcation_parameter ()

virtual void	actions_after_parameter_increase (double *const &parameter_pt)

double &	dof_derivative (const unsigned &i)

double &	dof_current (const unsigned &i)

virtual void	set_initial_condition ()

virtual double	global_temporal_error_norm ()

unsigned	newton_solve_continuation (double *const &parameter_pt)

unsigned	newton_solve_continuation (double *const &parameter_pt, DoubleVector &z)

void	calculate_continuation_derivatives (double *const &parameter_pt)

void	calculate_continuation_derivatives (const DoubleVector &z)

void	calculate_continuation_derivatives_fd (double *const &parameter_pt)

bool	does_pointer_correspond_to_problem_data (double *const &parameter_pt)

void	set_consistent_pinned_values_for_continuation ()

Protected Attributes inherited from oomph::Problem
Vector< Problem * >	Copy_of_problem_pt

std::map< double *, bool >	Calculate_dparameter_analytic

bool	Calculate_hessian_products_analytic

LinearAlgebraDistribution *	Dof_distribution_pt

Vector< double * >	Dof_pt
	Vector of pointers to dofs. More...

DoubleVectorWithHaloEntries	Element_count_per_dof

double	Relaxation_factor

double	Newton_solver_tolerance

unsigned	Max_newton_iterations
	Maximum number of Newton iterations. More...

unsigned	Nnewton_iter_taken

Vector< double >	Max_res
	Maximum residuals at start and after each newton iteration. More...

double	Max_residuals

bool	Time_adaptive_newton_crash_on_solve_fail

bool	Jacobian_reuse_is_enabled
	Is re-use of Jacobian in Newton iteration enabled? Default: false. More...

bool	Jacobian_has_been_computed

bool	Problem_is_nonlinear

bool	Pause_at_end_of_sparse_assembly

bool	Doc_time_in_distribute

unsigned	Sparse_assembly_method

unsigned	Sparse_assemble_with_arrays_initial_allocation

unsigned	Sparse_assemble_with_arrays_allocation_increment

Vector< Vector< unsigned > >	Sparse_assemble_with_arrays_previous_allocation

double	Numerical_zero_for_sparse_assembly

double	FD_step_used_in_get_hessian_vector_products

bool	Mass_matrix_reuse_is_enabled

bool	Mass_matrix_has_been_computed

bool	Discontinuous_element_formulation

double	Minimum_dt
	Minimum desired dt: if dt falls below this value, exit. More...

double	Maximum_dt
	Maximum desired dt. More...

double	DTSF_max_increase

double	DTSF_min_decrease

double	Minimum_dt_but_still_proceed

bool	Scale_arc_length
	Boolean to control whether arc-length should be scaled. More...

double	Desired_proportion_of_arc_length
	Proportion of the arc-length to taken by the parameter. More...

double	Theta_squared

int	Sign_of_jacobian
	Storage for the sign of the global Jacobian. More...

double	Continuation_direction

double	Parameter_derivative
	Storage for the derivative of the global parameter wrt arc-length. More...

double	Parameter_current
	Storage for the present value of the global parameter. More...

bool	Use_continuation_timestepper
	Boolean to control original or new storage of dof stuff. More...

unsigned	Dof_derivative_offset

unsigned	Dof_current_offset

Vector< double >	Dof_derivative
	Storage for the derivative of the problem variables wrt arc-length. More...

Vector< double >	Dof_current
	Storage for the present values of the variables. More...

double	Ds_current
	Storage for the current step value. More...

unsigned	Desired_newton_iterations_ds

double	Minimum_ds
	Minimum desired value of arc-length. More...

bool	Bifurcation_detection
	Boolean to control bifurcation detection via determinant of Jacobian. More...

bool	Bisect_to_find_bifurcation
	Boolean to control wheter bisection is used to located bifurcation. More...

bool	First_jacobian_sign_change
	Boolean to indicate whether a sign change has occured in the Jacobian. More...

bool	Arc_length_step_taken
	Boolean to indicate whether an arc-length step has been taken. More...

bool	Use_finite_differences_for_continuation_derivatives

OomphCommunicator *	Communicator_pt
	The communicator for this problem. More...

bool	Always_take_one_newton_step

double	Timestep_reduction_factor_after_nonconvergence

bool	Keep_temporal_error_below_tolerance

Static Protected Attributes inherited from oomph::Problem
static ContinuationStorageScheme	Continuation_time_stepper
	Storage for the single static continuation timestorage object. More...

Detailed Description

template<class ELEMENT>
class UnstructuredFvKProblem< ELEMENT >

Class definition.

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

Member Enumeration Documentation

◆ anonymous enum

template<class ELEMENT >

anonymous enum

private

Enumerator
Outer_boundary0
Outer_boundary1
Inner_boundary0
Inner_boundary1
Inner_boundary2
Inner_boundary3
Inner_boundary4
Inner_boundary5

   {
    Outer_boundary0 = 0,
    Outer_boundary1 = 1,
    Inner_boundary0 = 2,
    Inner_boundary1 = 3,
    Inner_boundary2 = 4,
    Inner_boundary3 = 5,
    Inner_boundary4 = 6,
    Inner_boundary5 = 7
   };

◆ anonymous enum

template<class ELEMENT >

anonymous enum

private

Enumerator
Outer_boundary0
Outer_boundary1
Inner_boundary0
Inner_boundary1
Inner_boundary2
Inner_boundary3
Inner_boundary4
Inner_boundary5

   {
    Outer_boundary0 = 0,
    Outer_boundary1 = 1,
    Inner_boundary0 = 2,
    Inner_boundary1 = 3,
    Inner_boundary2 = 4,
    Inner_boundary3 = 5,
    Inner_boundary4 = 6,
    Inner_boundary5 = 7
   };

◆ anonymous enum

template<class ELEMENT >

anonymous enum

private

Enumerator
Outer_boundary0
Outer_boundary1
Inner_boundary0
Inner_boundary1
Inner_boundary2
Inner_boundary3
Inner_boundary4
Inner_boundary5

   {
    Outer_boundary0 = 0,
    Outer_boundary1 = 1,
    Inner_boundary0 = 2,
    Inner_boundary1 = 3,
    Inner_boundary2 = 4,
    Inner_boundary3 = 5,
    Inner_boundary4 = 6,
    Inner_boundary5 = 7
   };

Constructor & Destructor Documentation

◆ UnstructuredFvKProblem() [1/3]

template<class ELEMENT >

UnstructuredFvKProblem< ELEMENT >::UnstructuredFvKProblem ( double element_area = 0.2 )

Constructor.

Constructor: Pass in element area.

Create the initial volume constraint element

Set the problem pressure to be that of the volume constraint element

Set initial pressure guess

Create the volume constraint mesh and add the element to the mesh

Add the sub meshes and build the global mesh

  : Element_area(element_area)
 {
  Vector<double> zeta(1);
  Vector<double> posn(2);
  
  //Outer boundary
  //--------------
  Ellipse* outer_boundary_ellipse_pt = new Ellipse(1.0,1.0);
  
  TriangleMeshClosedCurve* outer_boundary_pt = 0;
  
  Vector<TriangleMeshCurveSection*> outer_curvilinear_boundary_pt(2);
  
  //First bit
  double zeta_start = 0.0;
  double zeta_end = MathematicalConstants::Pi;
  unsigned nsegment = (int)(MathematicalConstants::Pi/sqrt(element_area));
  outer_curvilinear_boundary_pt[0] =
   new TriangleMeshCurviLine(outer_boundary_ellipse_pt, zeta_start,
     zeta_end, nsegment, Outer_boundary0);
  
  //Second bit
  zeta_start = MathematicalConstants::Pi;
  zeta_end = 2.0*MathematicalConstants::Pi;
  nsegment = (int)(MathematicalConstants::Pi/sqrt(element_area));
  outer_curvilinear_boundary_pt[1] =
   new TriangleMeshCurviLine(outer_boundary_ellipse_pt, zeta_start,
     zeta_end, nsegment, Outer_boundary1);
  
  outer_boundary_pt =
   new TriangleMeshClosedCurve(outer_curvilinear_boundary_pt);
  
  
  //Inner boundary
  //--------------
   Ellipse* inner_boundary_ellipse_pt = new Ellipse(GlobalParameters::R_b,
                                                   GlobalParameters::R_b);
  
  Vector<TriangleMeshClosedCurve*> inner_boundary_pt(1);
  Vector<TriangleMeshCurveSection*> inner_curvilinear_boundary_pt(2);
  
  //First part
  zeta_start = 0.0;
  zeta_end = MathematicalConstants::Pi;
  nsegment = (int)(GlobalParameters::R_b*
                   MathematicalConstants::Pi/sqrt(element_area));
  inner_curvilinear_boundary_pt[0] =
   new TriangleMeshCurviLine(inner_boundary_ellipse_pt, zeta_start,
     zeta_end, nsegment, Inner_boundary0);
  
  //Second part
  zeta_start = MathematicalConstants::Pi;
  zeta_end = 2.0*MathematicalConstants::Pi;
  nsegment = (int)(GlobalParameters::R_b*
                   MathematicalConstants::Pi/sqrt(element_area));
  inner_curvilinear_boundary_pt[1] =
   new TriangleMeshCurviLine(inner_boundary_ellipse_pt, zeta_start,
     zeta_end, nsegment, Inner_boundary1);
  
  //Combine to internal curvilinear boundary
  Vector<double> left_region_coords(2);
  Vector<double> right_region_coords(2);
  left_region_coords[0] = -0.5*GlobalParameters::R_b;
  left_region_coords[1] = 0.0;
  
  right_region_coords[0] = 0.5*GlobalParameters::R_b;
  right_region_coords[1] = 0.0;
  
  inner_boundary_pt[0] =
   new TriangleMeshClosedCurve(inner_curvilinear_boundary_pt);
  
  //Diameter boundary line
  Vector<TriangleMeshOpenCurve*> inner_open_boundary_pt(4);
  
  Vector<Vector<double> > vertex_coord(2, Vector<double>(2));
  
  //Section 1 of line
  vertex_coord[0][0] = 0.0;
  vertex_coord[0][1] = -1.0;
  vertex_coord[1][0] = 0.0;
  vertex_coord[1][1] = -GlobalParameters::R_b;
  
  TriangleMeshPolyLine *inner_open_polyline1_pt =
   new TriangleMeshPolyLine(vertex_coord, Inner_boundary2);
  
  inner_open_polyline1_pt->connect_initial_vertex_to_curviline(
    dynamic_cast<TriangleMeshCurviLine *>(outer_curvilinear_boundary_pt[1]),
    1.5*MathematicalConstants::Pi);
  
  inner_open_polyline1_pt->connect_final_vertex_to_curviline(
    dynamic_cast<TriangleMeshCurviLine *>(inner_curvilinear_boundary_pt[1]),
    1.5*MathematicalConstants::Pi);
  
  Vector<TriangleMeshCurveSection *> inner_boundary_line1_pt(1);
  inner_boundary_line1_pt[0] = inner_open_polyline1_pt;
  
  inner_open_boundary_pt[0] =
   new TriangleMeshOpenCurve(inner_boundary_line1_pt);
  
  //Section 2 of line
  vertex_coord[0][0] = 0.0;
  vertex_coord[0][1] = -GlobalParameters::R_b;
  vertex_coord[1][0] = 0.0;
  vertex_coord[1][1] = 0.0;
  
  TriangleMeshPolyLine *inner_open_polyline2_pt =
   new TriangleMeshPolyLine(vertex_coord, Inner_boundary3);
  
  inner_open_polyline2_pt->connect_initial_vertex_to_curviline(
    dynamic_cast<TriangleMeshCurviLine *>(inner_curvilinear_boundary_pt[1]),
    1.5*MathematicalConstants::Pi);
  
  Vector<TriangleMeshCurveSection *> inner_boundary_line2_pt(1);
  inner_boundary_line2_pt[0] = inner_open_polyline2_pt;
  
  inner_open_boundary_pt[1] =
   new TriangleMeshOpenCurve(inner_boundary_line2_pt);
  
  //Section 3 of line
  vertex_coord[0][0] = 0.0;
  vertex_coord[0][1] = 0.0;
  vertex_coord[1][0] = 0.0;
  vertex_coord[1][1] = GlobalParameters::R_b;
  
  TriangleMeshPolyLine *inner_open_polyline3_pt =
   new TriangleMeshPolyLine(vertex_coord, Inner_boundary4);
  
  inner_open_polyline3_pt->connect_initial_vertex_to_polyline(
    inner_open_polyline2_pt,1);
  
  inner_open_polyline3_pt->connect_final_vertex_to_curviline(
    dynamic_cast<TriangleMeshCurviLine *>(inner_curvilinear_boundary_pt[0]),
    0.5*MathematicalConstants::Pi);
  
  Vector<TriangleMeshCurveSection *> inner_boundary_line3_pt(1);
  inner_boundary_line3_pt[0] = inner_open_polyline3_pt;
  
  inner_open_boundary_pt[2] =
   new TriangleMeshOpenCurve(inner_boundary_line3_pt);
  
  //Section 4 of line
  vertex_coord[0][0] = 0.0;
  vertex_coord[0][1] = GlobalParameters::R_b;
  vertex_coord[1][0] = 0.0;
  vertex_coord[1][1] = 1.0;
  
  TriangleMeshPolyLine *inner_open_polyline4_pt =
   new TriangleMeshPolyLine(vertex_coord, Inner_boundary5);
  
  inner_open_polyline4_pt->connect_initial_vertex_to_curviline(
    dynamic_cast<TriangleMeshCurviLine *>(inner_curvilinear_boundary_pt[0]),
    0.5*MathematicalConstants::Pi);
  
  inner_open_polyline4_pt->connect_final_vertex_to_curviline(
    dynamic_cast<TriangleMeshCurviLine *>(outer_curvilinear_boundary_pt[0]),
    0.5*MathematicalConstants::Pi);
  
  Vector<TriangleMeshCurveSection *> inner_boundary_line4_pt(1);
  inner_boundary_line4_pt[0] = inner_open_polyline4_pt;
  
  inner_open_boundary_pt[3] =
   new TriangleMeshOpenCurve(inner_boundary_line4_pt);
  
  //Create the mesh
  //---------------
  
  //Create mesh parameters object
  TriangleMeshParameters mesh_parameters(outer_boundary_pt);
  mesh_parameters.internal_closed_curve_pt() = inner_boundary_pt;
  mesh_parameters.internal_open_curves_pt() = inner_open_boundary_pt;
  mesh_parameters.add_region_coordinates(1,left_region_coords);
  mesh_parameters.add_region_coordinates(2,right_region_coords);
  mesh_parameters.element_area() = element_area;
  
  // Build the bloody thing
  My_mesh_pt = new TriangleMesh<ELEMENT>(mesh_parameters);
  
  
  complete_problem_setup();
  add_sub_mesh(My_mesh_pt);
  build_global_mesh();
  
  char filename[100];
  sprintf(filename, "RESLT/trace.dat");
  Trace_file.open(filename);
  
  oomph_info << "Number of equations: "
             << this->assign_eqn_numbers() << '\n';
 }

References oomph::TriangleMeshParameters::add_region_coordinates(), oomph::Problem::add_sub_mesh(), oomph::Problem::assign_eqn_numbers(), oomph::Problem::build_global_mesh(), UnstructuredFvKProblem< ELEMENT >::complete_problem_setup(), oomph::TriangleMeshCurveSection::connect_final_vertex_to_curviline(), oomph::TriangleMeshCurveSection::connect_initial_vertex_to_curviline(), oomph::TriangleMeshCurveSection::connect_initial_vertex_to_polyline(), oomph::TriangleMeshParameters::element_area(), MergeRestartFiles::filename, UnstructuredFvKProblem< ELEMENT >::Inner_boundary0, UnstructuredFvKProblem< ELEMENT >::Inner_boundary1, UnstructuredFvKProblem< ELEMENT >::Inner_boundary2, UnstructuredFvKProblem< ELEMENT >::Inner_boundary3, UnstructuredFvKProblem< ELEMENT >::Inner_boundary4, UnstructuredFvKProblem< ELEMENT >::Inner_boundary5, int(), oomph::TriangleMeshParameters::internal_closed_curve_pt(), oomph::TriangleMeshParameters::internal_open_curves_pt(), UnstructuredFvKProblem< ELEMENT >::My_mesh_pt, oomph::oomph_info, UnstructuredFvKProblem< ELEMENT >::Outer_boundary0, UnstructuredFvKProblem< ELEMENT >::Outer_boundary1, BiharmonicTestFunctions2::Pi, GlobalParameters::R_b, sqrt(), UnstructuredFvKProblem< ELEMENT >::Trace_file, and Eigen::zeta().

◆ ~UnstructuredFvKProblem() [1/3]

template<class ELEMENT >

UnstructuredFvKProblem< ELEMENT >::~UnstructuredFvKProblem ( )

inline

Destructor.

   {
    delete My_mesh_pt;
   }

◆ UnstructuredFvKProblem() [2/3]

template<class ELEMENT >

UnstructuredFvKProblem< ELEMENT >::UnstructuredFvKProblem ( double element_area = 0.2 )

Constructor.

◆ ~UnstructuredFvKProblem() [2/3]

template<class ELEMENT >

UnstructuredFvKProblem< ELEMENT >::~UnstructuredFvKProblem ( )

inline

Destructor.

   {
    delete My_mesh_pt;
   }

◆ UnstructuredFvKProblem() [3/3]

template<class ELEMENT >

UnstructuredFvKProblem< ELEMENT >::UnstructuredFvKProblem ( double element_area = 0.2 )

Constructor.

◆ ~UnstructuredFvKProblem() [3/3]

template<class ELEMENT >

UnstructuredFvKProblem< ELEMENT >::~UnstructuredFvKProblem ( )

inline

Destructor.

   {
    delete My_mesh_pt;
   };

Member Function Documentation

◆ actions_after_adapt()

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::actions_after_adapt ( )

inlinevirtual

Actions after adapt: Setup the problem again – remember that the mesh has been completely rebuilt and its element's don't have any pointers to source fcts etc. yet Also create new volume constraint element with previous pressure value

Reimplemented from oomph::Problem.

   {
    Volume_constraint_element_pt =
     new FoepplvonKarmanVolumeConstraintElement
     <ELEMENT, RefineableTriangleMesh> (My_mesh_pt,
                                        Bubble_regions,
                                        *Temp_pressure_pt);
    Volume_constraint_element_pt->set_prescribed_volume(
      &TestSoln::prescribed_volume);
    Volume_constraint_mesh_pt->add_element_pt(Volume_constraint_element_pt);
    TestSoln::p_b_pt = Volume_constraint_element_pt->pressure_data_pt();
    rebuild_global_mesh();
    complete_problem_setup();
   }

References oomph::Mesh::add_element_pt(), TestSoln::p_b_pt, and TestSoln::prescribed_volume.

◆ actions_after_newton_solve() [1/3]

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::actions_after_newton_solve ( )

inlinevirtual

Update after solve (empty)

Reimplemented from oomph::Problem.

101 {}

◆ actions_after_newton_solve() [2/3]

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::actions_after_newton_solve ( )

inlinevirtual

Update after solve (empty)

Reimplemented from oomph::Problem.

101 {}

◆ actions_after_newton_solve() [3/3]

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::actions_after_newton_solve ( )

inlinevirtual

Update after solve (empty)

Reimplemented from oomph::Problem.

139 {

140 }

◆ actions_before_adapt()

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::actions_before_adapt ( )

inlinevirtual

Actions before adapt. Delete old volume constraint element and update global mesh

Reimplemented from oomph::Problem.

   {
    *Temp_pressure_pt =
     Volume_constraint_element_pt->pressure_data_pt()->value(0);
    Volume_constraint_mesh_pt->flush_element_and_node_storage();
    delete Volume_constraint_element_pt;
    rebuild_global_mesh();
   }

◆ actions_before_newton_solve() [1/3]

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::actions_before_newton_solve ( )

inlinevirtual

Update the problem specs before solve (empty)

Reimplemented from oomph::Problem.

104 {}

◆ actions_before_newton_solve() [2/3]

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::actions_before_newton_solve ( )

inlinevirtual

Update the problem specs before solve (empty)

Reimplemented from oomph::Problem.

104 {}

◆ actions_before_newton_solve() [3/3]

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::actions_before_newton_solve ( )

inlinevirtual

Update the problem specs before solve: Re-apply boundary conditons.

Reimplemented from oomph::Problem.

   {
    apply_boundary_conditions();
   }

◆ apply_boundary_conditions()

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::apply_boundary_conditions

private

Helper function to apply boundary conditions.

 {
  
  // Loop over all boundary nodes
  //unsigned nbound=this->My_mesh_pt->nboundary();
  //Just loop over outer boundary since inner boundary doesn't have boundary
  //conditions
  unsigned nbound = Outer_boundary1 + 1;
  
  for(unsigned ibound=0;ibound<nbound;ibound++)
   {
    unsigned num_nod=this->My_mesh_pt->nboundary_node(ibound);
    for (unsigned inod=0;inod<num_nod;inod++)
     {
      // Get node
      Node* nod_pt=this->My_mesh_pt->boundary_node_pt(ibound,inod);
      
      // Extract nodal coordinates from node:
      Vector<double> x(2);
      x[0]=nod_pt->x(0);
      x[1]=nod_pt->x(1);
     }
   } 
  
 } // end set bc

References plotDoE::x, and oomph::Node::x().

◆ complete_problem_setup() [1/3]

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::complete_problem_setup

private

Helper function to (re-)set boundary condition and complete the build of all elements

Set boundary condition exactly, and complete the build of all elements

Number of sample points

 {   
  
  // Set the boundary conditions for problem: All nodes are
  // free by default -- just pin the ones that have Dirichlet 
  // conditions here. 
  unsigned nbound = Outer_boundary1 + 1;
  for(unsigned ibound=0;ibound<nbound;ibound++)
   {
    unsigned num_nod=My_mesh_pt->nboundary_node(ibound);
    for (unsigned inod=0;inod<num_nod;inod++)
     {
      // Get node
      Node* nod_pt=My_mesh_pt->boundary_node_pt(ibound,inod);
      
      // Pin unknown values
      nod_pt->pin(0);  // pin vertical displacement; Laplacian of w unpinned
                       // imposes clamping
      nod_pt->pin(2);  // pin Airy stress fct 
     }   
   } // end loop over boundaries
  
  
  
  // Complete the build of all elements so they are fully functional
  unsigned n_element = My_mesh_pt->nelement();
  for(unsigned e=0;e<n_element;e++)
   {
    // Upcast from GeneralisedElement to the present element
    ELEMENT* el_pt = dynamic_cast<ELEMENT*>(My_mesh_pt->element_pt(e));
    
    //Set the pressure function pointers and the physical constants
    el_pt->eta_pt() = &GlobalParameters::Eta;
    el_pt->pressure_fct_pt() = &GlobalParameters::get_pressure;
  
   }
  
   // Create the mesh as Geom Object
   Mesh_as_geom_object_pt=new MeshAsGeomObject(My_mesh_pt);
  
  
   
   unsigned n_sample=50;
   const double dr = 1.0/double(n_sample);
  
   // Get 'em
   Vector<double> x(2,0.0);
   Radial_sample_point_pt.resize(n_sample);
   for (unsigned j=0;j<n_sample;j++)
    {
     Radial_sample_point_pt[j].second.resize(2);
     x[0]=double(j)*dr;
     
     // Get the element and its local coordinates
     Mesh_as_geom_object_pt->locate_zeta(x,
                                         Radial_sample_point_pt[j].first,
                                         Radial_sample_point_pt[j].second);
    }
 }

References e(), GlobalParameters::Eta, GlobalParameters::get_pressure(), j, oomph::Data::pin(), and plotDoE::x.

Referenced by UnstructuredFvKProblem< ELEMENT >::UnstructuredFvKProblem().

◆ complete_problem_setup() [2/3]

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::complete_problem_setup ( )

private

Helper function to (re-)set boundary condition and complete the build of all elements

◆ complete_problem_setup() [3/3]

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::complete_problem_setup ( )

private

Helper function to (re-)set boundary condition and complete the build of all elements

◆ doc_solution() [1/3]

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::doc_solution

Doc the solution.

 { 
  ofstream some_file;
  char filename[100];
  
  // Number of plot points
  unsigned npts = 5;
  
  sprintf(filename,"RESLT/soln%i.dat",Doc_info.number());
  some_file.open(filename);
  this->My_mesh_pt->output(some_file,npts); 
  some_file.close();
  
  // Output boundaries
  //------------------
  sprintf(filename,"RESLT/boundaries%i.dat",Doc_info.number());
  some_file.open(filename);
  My_mesh_pt->output_boundaries(some_file);
  some_file.close();
  
  // Find solution at r = 0
  // ----------------------
  unsigned n_boundary3_node = My_mesh_pt->nboundary_node(Inner_boundary3);
  unsigned n_boundary4_node = My_mesh_pt->nboundary_node(Inner_boundary4);
  double w_0 = 0;
  for(unsigned inode=0; inode < n_boundary3_node; inode++)
   {
    for(unsigned jnode=0; jnode < n_boundary4_node; jnode++)
     {
      if(My_mesh_pt->boundary_node_pt(Inner_boundary3,inode)
        == My_mesh_pt->boundary_node_pt(Inner_boundary4,jnode))
       {
        w_0 = My_mesh_pt->boundary_node_pt(Inner_boundary3,inode)->value(0);
       }
     }
   }
  
  // Plot solution along radial line
  sprintf(filename,"RESLT/soln_along_radial_line%i.dat",Doc_info.number());
  some_file.open(filename);
  Vector<double> x(2);
  unsigned nplot=Radial_sample_point_pt.size();
  for (unsigned j=0;j<nplot;j++)
   {
    ELEMENT* el_pt=dynamic_cast<ELEMENT*>(Radial_sample_point_pt[j].first);
    Vector<double> s=Radial_sample_point_pt[j].second;
    el_pt->interpolated_x(s,x);
    some_file << x[0] << " " 
              << x[1] << " " 
              << el_pt->interpolated_w_fvk(s) << " " // w 
              << std::endl;
   }
  some_file.close();
  
  
  // Doc error and return of the square of the L2 error
  //---------------------------------------------------
  double dummy_error,zero_norm;
  sprintf(filename,"RESLT/norm%i.dat",Doc_info.number());
  some_file.open(filename);
  My_mesh_pt->compute_error(some_file,GlobalParameters::zero,
                            dummy_error,zero_norm);
  some_file.close();
  
  // Doc L2 error and norm of solution
  oomph_info << "Norm of computed solution: " 
             << sqrt(dummy_error) << std::endl;
  
  Trace_file << GlobalParameters::Pressure << " "
             << w_0 << '\n';
  
  // Increment the doc_info number
  Doc_info.number()++;
  
 } // end of doc

References GlobalParameters::Doc_info, MergeRestartFiles::filename, j, oomph::DocInfo::number(), oomph::oomph_info, GlobalParameters::Pressure, s, sqrt(), oomph::Problem_Parameter::Trace_file, plotDoE::x, and GlobalParameters::zero().

◆ doc_solution() [2/3]

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::doc_solution ( )

Doc the solution.

◆ doc_solution() [3/3]

template<class ELEMENT >

void UnstructuredFvKProblem< ELEMENT >::doc_solution ( const std::string & comment = "" )

Doc the solution.

 { 
  ofstream some_file;
  char filename[100];
  
  // Number of plot points
  unsigned npts = 5;
  
  sprintf(filename,"RESLT/soln%i-%f.dat",Doc_info.number(),Element_area);
  some_file.open(filename);
  this->My_mesh_pt->output(some_file,npts); 
  some_file << "TEXT X = 22, Y = 92, CS=FRAME T = \"" 
            << comment << "\"\n";
  some_file.close();
  
  // Output boundaries
  //------------------
  sprintf(filename,"RESLT/boundaries%i-%f.dat",Doc_info.number(),Element_area);
  some_file.open(filename);
  My_mesh_pt->output_boundaries(some_file);
  some_file.close();
  
  
  // Output along radial boundary
  // ----------------------------
  sprintf(filename,"RESLT/radial_soln%i-%f.dat",Doc_info.number(),Element_area);
  some_file.open(filename);
  for(unsigned b = Inner_boundary2; b <= Inner_boundary5; b++)
   {
    unsigned n_boundary_el = My_mesh_pt->nboundary_element(b);
    for(unsigned e = 0; e < n_boundary_el; e++)
     {
      My_mesh_pt->boundary_element_pt(b,e)->output(some_file,npts);
     }
   }
  some_file.close();
  
  // Find solution at r = 0
  // ----------------------
  unsigned n_boundary3_node = My_mesh_pt->nboundary_node(Inner_boundary3);
  unsigned n_boundary4_node = My_mesh_pt->nboundary_node(Inner_boundary4);
  double w_0 = 0;
  for(unsigned inode=0; inode < n_boundary3_node; inode++)
   {
    for(unsigned jnode=0; jnode < n_boundary4_node; jnode++)
     {
      if(My_mesh_pt->boundary_node_pt(Inner_boundary3,inode)
        == My_mesh_pt->boundary_node_pt(Inner_boundary4,jnode))
       {
        w_0 = My_mesh_pt->boundary_node_pt(Inner_boundary3,inode)->value(0);
       }
     }
   }
  
  
  // Output regions
  unsigned n_region = My_mesh_pt->nregion();
  if (n_region > 1)
   {
    for (unsigned r = 0; r < n_region; r++)
     {
      //Attempt to output elements in different regions
      sprintf(filename,"RESLT/region%i%i-%f.dat",r,Doc_info.number(),
        Element_area);
      some_file.open(filename);
      unsigned nel = My_mesh_pt->nregion_element(r);
      for (unsigned e = 0; e < nel; e++)
       {
        My_mesh_pt->region_element_pt(r,e)->output(some_file,npts);
       }
      some_file.close();
     }
   }
  
  // Calculate the integral (i.e. volume) over the inner region
  double bubble_volume = 0;
  if (n_region > 1)
   {
    sprintf(filename,"RESLT/bubble_volume%i.dat",Doc_info.number());
    some_file.open(filename);
    unsigned n_inner_el;
    unsigned n_bubble_regions = Bubble_regions.size();
  
    for(unsigned r = 0;r < n_bubble_regions; r++)
     {
      n_inner_el = My_mesh_pt->nregion_element(Bubble_regions[r]);
      for(unsigned e = 0; e < n_inner_el; e++)
       {
        ELEMENT* el_pt = dynamic_cast<ELEMENT*>(
          My_mesh_pt->region_element_pt(Bubble_regions[r],e));
        if(el_pt != 0)
         {
          bubble_volume += el_pt->get_bounded_volume();
         }
       }
     }
    some_file << bubble_volume << '\n';
    some_file.close();
   }
  
  
  // Doc error and return of the square of the L2 error
  //---------------------------------------------------
  //double error,norm,dummy_error,zero_norm;
  double dummy_error,zero_norm;
  sprintf(filename,"RESLT/error%i-%f.dat",Doc_info.number(),Element_area);
  some_file.open(filename);
  
  My_mesh_pt->compute_error(some_file,TestSoln::zero,
                            dummy_error,zero_norm);
  some_file.close();
  
  // Doc L2 error and norm of solution
  oomph_info << "Norm of computed solution: " << sqrt(dummy_error) << std::endl;
  
  Trace_file << TestSoln::p_b_pt->value(0) << " "
             << TestSoln::p_0_pt->value(0) << " "
             << w_0 << " " << bubble_volume << '\n';
  
  // Increment the doc_info number
  Doc_info.number()++;
  
 } // end of doc

private

Trace file to document norm of solution.

Referenced by UnstructuredFvKProblem< ELEMENT >::UnstructuredFvKProblem().

◆ Volume_constraint_element_pt

template<class ELEMENT >

FoepplvonKarmanVolumeConstraintElement<ELEMENT, RefineableTriangleMesh>* UnstructuredFvKProblem< ELEMENT >::Volume_constraint_element_pt

private

Single volume constraint element instance.

◆ Volume_constraint_mesh_pt

template<class ELEMENT >

Mesh* UnstructuredFvKProblem< ELEMENT >::Volume_constraint_mesh_pt

private

Mesh for the volume constraint element.

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

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

template<class ELEMENT> class UnstructuredFvKProblem< ELEMENT >

Member Enumeration Documentation

◆ anonymous enum

◆ anonymous enum

◆ anonymous enum

Constructor & Destructor Documentation

◆ UnstructuredFvKProblem() [1/3]

◆ ~UnstructuredFvKProblem() [1/3]

◆ UnstructuredFvKProblem() [2/3]

◆ ~UnstructuredFvKProblem() [2/3]

◆ UnstructuredFvKProblem() [3/3]

◆ ~UnstructuredFvKProblem() [3/3]

Member Function Documentation

◆ actions_after_adapt()

◆ actions_after_newton_solve() [1/3]

◆ actions_after_newton_solve() [2/3]

◆ actions_after_newton_solve() [3/3]

◆ actions_before_adapt()

◆ actions_before_newton_solve() [1/3]

◆ actions_before_newton_solve() [2/3]

◆ actions_before_newton_solve() [3/3]

◆ apply_boundary_conditions()

◆ complete_problem_setup() [1/3]

◆ complete_problem_setup() [2/3]

◆ complete_problem_setup() [3/3]

◆ doc_solution() [1/3]

◆ doc_solution() [2/3]

◆ doc_solution() [3/3]

Member Data Documentation

◆ Bubble_regions

◆ Doc_info

◆ Element_area

◆ Mesh_as_geom_object_pt

◆ My_mesh_pt [1/2]

◆ My_mesh_pt [2/2]

◆ Radial_sample_point_pt

◆ Temp_pressure_pt

◆ Trace_file

◆ Volume_constraint_element_pt

◆ Volume_constraint_mesh_pt

template<class ELEMENT>
class UnstructuredFvKProblem< ELEMENT >