two_d_linear_elasticity_with_simple_block_diagonal_preconditioner.cc File Reference

#include "generic.h"
#include "linear_elasticity.h"
#include "meshes/rectangular_quadmesh.h"
#include "simple_block_preconditioners.h"

Classes
class	oomph::MyLinearElasticityElement< DIM >
class	oomph::FaceGeometry< MyLinearElasticityElement< DIM > >
class	PeriodicLoadProblem< ELEMENT >
	Periodic loading problem. More...

Namespaces
	oomph
	DRAIG: Change all instances of (SPATIAL_DIM) to (DIM-1).
	Global_Parameters
	Namespace for global parameters.

Functions
TimeHarmonicIsotropicElasticityTensor	Global_Parameters::E (Nu)
	The elasticity tensor for the solid. More...
void	Global_Parameters::exact_solution (const Vector< double > &x, Vector< double > &u)
	The exact solution for infinite depth case. More...
void	Global_Parameters::periodic_traction (const double &time, const Vector< double > &x, const Vector< double > &n, Vector< double > &result)
	The traction function. More...
int	main (int argc, char *argv[])
	Driver code for PeriodicLoad linearly elastic problem. More...

Function Documentation

main()

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

Driver code for PeriodicLoad linearly elastic problem.

{
#ifdef OOMPH_HAS_MPI
 // Initialise MPI
 MPI_Helpers::init(argc,argv);
#endif
 
 // Number of elements in x-direction
 unsigned nx=5;
 
 // Number of elements in y-direction (for (approximately) square elements)
 unsigned ny=unsigned(double(nx)*Global_Parameters::Ly/Global_Parameters::Lx);
 
 // Set up doc info
 DocInfo doc_info;
 
 // Set output directory
 doc_info.set_directory("RESLT");
 
 //Build the problem
 PeriodicLoadProblem<MyLinearElasticityElement<2> > 
  problem(nx,ny,Global_Parameters::Lx, Global_Parameters::Ly);
 
 // Solve
 problem.newton_solve();
 
 // Output the solution
 problem.doc_solution(doc_info);
 
#ifdef OOMPH_HAS_MPI
 // finalize MPI
 MPI_Helpers::finalize();
#endif
 
 return(EXIT_SUCCESS);  
} // end_of_main

References Global_Parameters::Lx, Global_Parameters::Ly, Mesh_Parameters::nx, Mesh_Parameters::ny, problem, and oomph::DocInfo::set_directory().