unstructured_cylinder.cc File Reference

#include "generic.h"
#include "time_harmonic_fourier_decomposed_linear_elasticity.h"
#include "meshes/triangle_mesh.h"

Classes
class	FourierDecomposedTimeHarmonicLinearElasticityProblem< ELEMENT >

Namespaces
	Global_Parameters
	Namespace for global parameters.

Functions
std::complex< double >	Global_Parameters::Nu (0.3, 0.05)
	Define Poisson's ratio Nu. More...
std::complex< double >	Global_Parameters::E (1.0, 0.01)
	Define the non-dimensional Young's modulus. More...
std::complex< double >	Global_Parameters::Omega_sq (10.0, 5.0)
const std::complex< double >	Global_Parameters::I (0.0, 1.0)
	Define the imaginary unit. More...
void	Global_Parameters::boundary_traction (const Vector< double > &x, const Vector< double > &n, Vector< std::complex< double > > &result)
	The traction function at r=rmin: (t_r, t_z, t_theta) More...
void	Global_Parameters::body_force (const Vector< double > &x, Vector< std::complex< double > > &b)
	Body force. More...
void	Global_Parameters::exact_solution (const Vector< double > &x, Vector< double > &u)
	The exact solution for infinite depth case. More...
int	main (int argc, char *argv[])
	Driver code. More...

Function Documentation

main()

int main	(	int argc	,
		char *	argv[]
	)

Driver code.

{
 
 // Set up doc info
 DocInfo doc_info;
 
 // Set output directory
 doc_info.set_directory("RESLT");
 
#ifdef ADAPTIVE
 
 // Set up problem
 FourierDecomposedTimeHarmonicLinearElasticityProblem
  <ProjectableTimeHarmonicFourierDecomposedLinearElasticityElement
   <TTimeHarmonicFourierDecomposedLinearElasticityElement<3> > >
  problem(Global_Parameters::rmin,Global_Parameters::rmax,
          Global_Parameters::zmin,Global_Parameters::zmax);
 
 // Solve
 unsigned max_adapt=3;
 problem.newton_solve(max_adapt);
 
#else
 
 // Set up problem
 FourierDecomposedTimeHarmonicLinearElasticityProblem
  <TTimeHarmonicFourierDecomposedLinearElasticityElement<3> > 
  problem(Global_Parameters::rmin,Global_Parameters::rmax,
          Global_Parameters::zmin,Global_Parameters::zmax);
 
 // Solve
 problem.newton_solve();
 
#endif
 
 
 // Output the solution
 problem.doc_solution(doc_info);
  
} // end_of_main

References problem, Global_Parameters::rmax, Global_Parameters::rmin, oomph::DocInfo::set_directory(), Global_Parameters::zmax, and Global_Parameters::zmin.