sparse_assemble_test.cc File Reference

#include <time.h>
#include "generic.h"
#include "navier_stokes.h"
#include "meshes/simple_rectangular_quadmesh.h"

Classes
class	RectangularDrivenCavityProblem< ELEMENT >
	Driven cavity problem in rectangular domain. More...

Namespaces
	Global_Variables
	Namespace for physical parameters.

Functions
int	main (int argc, char *argv[])
	Driver to compare different assembly strategies. More...

Variables
unsigned	Global_Variables::Nx =4
	Number of elements in x direction. More...
unsigned	Global_Variables::Ny =4
	Number of elements in x direction. More...
bool	Global_Variables::Dump_matrices =true
	By default we're dumping the matrices for comparison. More...
bool	Global_Variables::Halt_code =false

Function Documentation

main()

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

Driver to compare different assembly strategies.

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

{
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 cout << "Testing sparse assembly" << endl;
 cout << "-----------------------" << endl;
 
 // Loop over all assembly strategies or just do one?
 unsigned lo_method=1;
 unsigned hi_method=8;
 
 if (CommandLineArgs::Argc==1)
  {
   std::cout 
    << "No command line args specified -- running all assembly strategies"
    << std::endl;
  }
 else if (CommandLineArgs::Argc==6)
  {
   Global_Variables::Nx=atoi(argv[1]);
   Global_Variables::Ny=atoi(argv[1]);
   lo_method=atoi(argv[3]);
   hi_method=atoi(argv[3]);
   std::cout << "Assembly method: " << hi_method << std::endl;
   Global_Variables::Dump_matrices=atoi(argv[4]);
   Global_Variables::Halt_code=atoi(argv[5]);
  }
 else
  {
   std::cout 
    << "Wrong number of command line args specified!\n" 
    << "Specify none or five: \n"
    << "- Number of elements in x-direction\n"
    << "- Number of elements in y-direction\n"
    << "- assembly method [1-4]\n" 
    << "- output of matrices [0/1]\n" 
    << "- halt code to analyse memory usage [0/1]\n" 
    << std::endl;
   exit(1);
  }
 
 
 cout << "Mesh with " 
      << Global_Variables::Nx << " x " 
      << Global_Variables::Ny << " elements."
      << std::endl;
 
 for (unsigned method = lo_method; method<=hi_method; method++)
  {
   // Build the problem with QTaylorHoodElements
   RectangularDrivenCavityProblem<QTaylorHoodElement<2> > problem;
 
   // Perform actions before solve to get the problem ready
   problem.actions_before_newton_solve();
 
   // Run with assembly with specified method
   problem.compare_assembly_strategies(method);
 
  }
 
} // end_of_main

References oomph::CommandLineArgs::Argc, Global_Variables::Dump_matrices, Global_Variables::Halt_code, Global_Variables::Nx, Global_Variables::Ny, problem, and Flag_definition::setup().