t_convergence_3d.cc File Reference

#include "generic.h"
#include "poisson.h"
#include "meshes/simple_cubic_tet_mesh.h"

Classes
class	PoissonProblem< ELEMENT >
	Micky mouse Poisson problem. More...

Namespaces
	TanhSolnForPoisson
	Namespace for exact solution for Poisson equation with "sharp step".

Functions
void	TanhSolnForPoisson::get_exact_u (const Vector< double > &x, Vector< double > &u)
	Exact solution as a Vector. More...
void	TanhSolnForPoisson::get_exact_u (const Vector< double > &x, double &u)
	Exact solution as a scalar. More...
void	TanhSolnForPoisson::get_source (const Vector< double > &x, double &source)
	Source function to make it an exact solution. More...
int	main (int argc, char *argv[])
	Demonstrate how to solve Poisson problem. More...

Function Documentation

main()

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

Demonstrate how to solve Poisson problem.

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

Label for output

{
 
  // Store command line arguments
  CommandLineArgs::setup(argc,argv);
 
  DocInfo doc_info;
 
  // Output directory
  doc_info.set_directory("RESLT");
 
  // Step number
  doc_info.number()=0;
 
  unsigned nstep;  
 
  // Output for convergence data
  ofstream convergence_file("RESLT/convergence.dat");
 
  convergence_file << "VARIABLES=\"log h\",\"log e\",\"theory\",\"log u\" " 
                   << std::endl;
 
 
  //Tecplot output separator for convergence results
  convergence_file << "ZONE T=\"Linear element\"" << std::endl;
 
  cout << "Doing linear elements" << std::endl;
  cout << "=====================" << std::endl;
 
  // Number of steps (reduced for validation)
  nstep = 5;   
  if (CommandLineArgs::Argc>1)
   {
    nstep=2;
   }
  for (unsigned h_power=0;h_power<nstep;h_power++)
   {
 
    //Set up the problem
    PoissonProblem<TPoissonElement<3,2> > 
     problem(&TanhSolnForPoisson::get_source,h_power);
 
    cout << "Self test: " << problem.self_test() << std::endl;
 
    // Solve the problem
    problem.newton_solve();
 
    //Output solution & convergence data
    problem.doc_solution(doc_info,convergence_file);
 
    //Increment counter for solutions 
    doc_info.number()++;
   }
 
 
  cout << "Doing quadratic elements" << std::endl;
  cout << "========================" << std::endl;
  
  //Tecplot output separator for convergence results
  convergence_file << "ZONE T=\"Quadratic element\"" << std::endl;
 
  // Number of steps (reduced for validation)
  nstep = 4;  
  if (CommandLineArgs::Argc>1)
   {
    nstep=2;
   }  
  for (unsigned h_power=0;h_power<nstep;h_power++)
   {
    
    //Set up the problem
    PoissonProblem<TPoissonElement<3,3> > 
     problem(&TanhSolnForPoisson::get_source,h_power);
    
    // Solve the problem
    problem.newton_solve();
    
    //Output solution & convergence data
    problem.doc_solution(doc_info,convergence_file);
    
    //Increment counter for solutions 
    doc_info.number()++;
   }
  
 
  convergence_file.close();
 
}

References oomph::CommandLineArgs::Argc, TanhSolnForPoisson::get_source(), oomph::DocInfo::number(), problem, oomph::DocInfo::set_directory(), and Flag_definition::setup().