axisym_foeppl_von_karman/axisym_fvk.cc File Reference

#include <fenv.h>
#include "generic.h"
#include "axisym_foeppl_von_karman.h"
#include "meshes/one_d_mesh.h"

Classes
class	AxisymFvKProblem< ELEMENT >
	Problem class. More...

Namespaces
	AxisymFvKParameters
	Namespace for pressure solution.

Functions
void	AxisymFvKParameters::pressure_function (const double &r, double &pressure)
	Function that computes the pressure at radius r. More...
void	AxisymFvKParameters::get_exact_u (const Vector< double > &r, Vector< double > &u)
	Function to get the exact solution for the pure bending model. More...
int	main (int argc, char **argv)
	Driver. More...

Function Documentation

main()

int main	(	int argc	,
		char **	argv
	)

Driver.

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

{
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Define possible command line arguments and parse the ones that
 // were actually specified
 
 // Value of eta
 CommandLineArgs::specify_command_line_flag("--eta",
                                            &AxisymFvKParameters::Eta);
 
 // Choose between pure bending or fvk model
 CommandLineArgs::specify_command_line_flag("--use_linear_eqns");
 
 // Output directory
 CommandLineArgs::specify_command_line_flag("--dir",
                                            &AxisymFvKParameters::Directory);
 
 // Number of elements in the mesh
 unsigned n_element=100;
 CommandLineArgs::specify_command_line_flag("--n_element", &n_element);
 
 // Increment displacements
 unsigned n_step=10; // hierher 10;
 CommandLineArgs::specify_command_line_flag("--n_step", &n_step);
 double dp=10;
 CommandLineArgs::specify_command_line_flag("--dp", &dp);
 
 // Parse command line
 CommandLineArgs::parse_and_assign(); 
 
 // Doc what has actually been specified on the command line
 CommandLineArgs::doc_specified_flags();
 
 // Set up the problem: 
 AxisymFvKProblem<AxisymFoepplvonKarmanElement<3> >  
  problem(n_element);
 
 // Set initial value for pressure 
 AxisymFvKParameters::Pressure=0.0;
 
 for (unsigned i=0;i<n_step;i++)
  {
   // Solve the problem
   problem.newton_solve();
   
   //Output solution
   problem.doc_solution();
   
   // Increment pressure
   AxisymFvKParameters::Pressure+=dp; 
  }
 
 cout << "\n Done \n" << std::endl;
 
} // end of main

References AxisymFvKParameters::Directory, oomph::CommandLineArgs::doc_specified_flags(), AxisymFvKParameters::Eta, i, oomph::CommandLineArgs::parse_and_assign(), AxisymFvKParameters::Pressure, problem, Flag_definition::setup(), and oomph::CommandLineArgs::specify_command_line_flag().