mpi/distribution/clamped_shell/clamped_shell_with_arclength_cont.cc File Reference

#include <iostream>
#include <fstream>
#include <cmath>
#include <typeinfo>
#include <algorithm>
#include <cstdio>
#include "generic.h"
#include "shell.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	ShellMesh< ELEMENT >
class	ShellProblem< ELEMENT >

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
double &	Global_Physical_Variables::external_pressure ()
	Access function to value of external pressure. More...
void	Global_Physical_Variables::press_load (const Vector< double > &xi, const Vector< double > &x, const Vector< double > &N, Vector< double > &load)
	Load function: Perturbation pressure to force non-axisymmetric deformation. More...
int	main (int argc, char *argv[])
	Driver. More...

Variables
double	Global_Physical_Variables::Prescribed_y = 1.0
	Prescribed position of control point. More...

Function Documentation

main()

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

Driver.

{
#ifdef OOMPH_HAS_MPI
 MPI_Helpers::init(argc,argv);
#endif
 { 
  //SuperLUSolver::Suppress_incorrect_rhs_distribution_warning_in_resolve=true;;
 
 //Length of domain
 double L = 10.0;
 double L_phi=0.5*MathematicalConstants::Pi;
 
 //Set up the problem
 ShellProblem<StorableShapeSolidElement<DiagHermiteShellElement> > 
  problem(5,5,L,L_phi);
 
 //Let's just be crazy and distribut it
#ifdef OOMPH_HAS_MPI
  //Set up a dummy partition
  unsigned n_element = problem.mesh_pt()->nelement();
  Vector<unsigned> element_partition(n_element);
  for(unsigned e=0;e<n_element/2;e++) {element_partition[e]=0;}
  for(unsigned e=n_element/2;e<n_element;e++) {element_partition[e]=1;}
  
  DocInfo mesh_doc_info;
  bool report_stats=true;
  mesh_doc_info.set_directory("RESLT_MESH");
  problem.distribute(element_partition,mesh_doc_info,report_stats);
  problem.check_halo_schemes(mesh_doc_info);
#endif
 
 //Solve the problem
 problem.solve();
 }
#ifdef OOMPH_HAS_MPI
 MPI_Helpers::finalize();
#endif
}

References e(), L, BiharmonicTestFunctions2::Pi, problem, and oomph::DocInfo::set_directory().