lin_unsteady_ring.cc File Reference

#include "generic.h"
#include "beam.h"
#include "meshes/one_d_lagrangian_mesh.h"

Classes
class	ElasticRingProblem< ELEMENT >
	Ring problem. More...

Namespaces
	Global_Physical_Variables
	Global variables.

Functions
int	main (int argc, char *argv[])
	Driver for ring that performs small-amplitude oscillations. More...

Variables
unsigned	Global_Physical_Variables::Long_run_flag =1
	Flag for long/short run: Default = perform long run. More...
unsigned	Global_Physical_Variables::Fixed_timestep_flag =1
	Flag for fixed timestep: Default = fixed timestep. More...
bool	Global_Physical_Variables::Consistent_newmark_ic

Function Documentation

main()

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

Driver for ring that performs small-amplitude oscillations.

Convert command line arguments (if any) into flags:

{
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 if (argc==2)
  {
   // Nontrivial command line input: Setup Newmark IC directly
   // (rather than consistently with PVD)
   if (atoi(argv[1])==1) 
    {
     Global_Physical_Variables::Consistent_newmark_ic=true;
     cout << "Setting Newmark IC consistently" << std::endl;
    }
   else
    {
     Global_Physical_Variables::Consistent_newmark_ic=false;
     cout << "Setting Newmark IC directly" << std::endl;
    }
   
   cout << "Not enough command line arguments specified -- using defaults." 
        << std::endl;
  } // end of 1 argument
 else if (argc==4)
  {
   cout << "Three command line arguments specified:" << std::endl;
   // Nontrivial command line input: Setup Newmark IC directly
   // (rather than consistently with PVD)
   if (atoi(argv[1])==1) 
    {
     Global_Physical_Variables::Consistent_newmark_ic=true;
     cout << "Setting Newmark IC consistently" << std::endl;
    }
   else
    {
     Global_Physical_Variables::Consistent_newmark_ic=false;
     cout << "Setting Newmark IC directly" << std::endl;
    }
   // Flag for long run
   Global_Physical_Variables::Long_run_flag=atoi(argv[2]);
   // Flag for fixed timestep
   Global_Physical_Variables::Fixed_timestep_flag=atoi(argv[3]);
  } // end of 3 arguments
 else
  {
   std::string error_message =
    "Wrong number of command line arguments. Specify one or three.\n";
   error_message += "Arg1: Long_run_flag [0/1]\n";
   error_message += "Arg2: Impulsive_start_flag [0/1]\n";
   error_message += "Arg3: Restart_flag [restart_file] (optional)\n";
 
   throw OomphLibError(error_message,
                       OOMPH_CURRENT_FUNCTION,
                       OOMPH_EXCEPTION_LOCATION);
  } // too many arguments
 cout << "Setting Newmark IC consistently: "
      <<  Global_Physical_Variables::Consistent_newmark_ic << std::endl;
 cout << "Long run flag: " 
      <<  Global_Physical_Variables::Long_run_flag << std::endl;
 cout << "Fixed timestep flag: " 
      <<  Global_Physical_Variables::Fixed_timestep_flag << std::endl;
 
 //Length of domain
 double L = MathematicalConstants::Pi/2.0;
  
 // Number of elements
 unsigned nelem = 13;
 
 //Set up the problem
 ElasticRingProblem<HermiteBeamElement,Newmark<3> > 
  problem(nelem,L);
 
 // Do unsteady run
 problem.unsteady_run();
 
} // end of main

References Global_Physical_Variables::Consistent_newmark_ic, Global_Physical_Variables::Fixed_timestep_flag, L, Global_Physical_Variables::Long_run_flag, OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, BiharmonicTestFunctions2::Pi, problem, Flag_definition::setup(), and oomph::Global_string_for_annotation::string().