eighth_sphere_poisson_hp_adapt.cc File Reference

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

Classes
class	EighthSpherePoissonProblem< ELEMENT >
	Poisson problem in refineable eighth of a sphere mesh. 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[])

Function Documentation

main()

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

///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// Driver for 3D Poisson problem in eighth of a sphere. Solution has a sharp step. If there are any command line arguments, we regard this as a validation run and perform only a single adaptation.

{ 
 
 // Store command line arguments
 CommandLineArgs::setup(argc,argv);
 
 // Set up the problem with 27-node brick elements, pass pointer to 
 // source function
 EighthSpherePoissonProblem<PRefineableQPoissonElement<3> >
  problem(&TanhSolnForPoisson::get_source);
 
 // Setup labels for output
  DocInfo doc_info;
 
 // Output directory
 doc_info.set_directory("RESLT");
 
 // Step number
 doc_info.number()=0;
 
 // Check if we're ready to go
 std::cout << "Self test: " << problem.self_test() << std::endl;
 
 // Solve the problem
 problem.newton_solve();
 
 //Output solution
 problem.doc_solution(doc_info);
 
 //Increment counter for solutions 
 doc_info.number()++;
 
 // Perform initial mesh refinement and solve
 problem.refine_uniformly();
 problem.refine_uniformly();
 problem.p_refine_uniformly();
 
 //std::cout << "Self test: " << problem.self_test() << std::endl;
 problem.newton_solve();
 problem.doc_solution(doc_info);
 doc_info.number()++;
 
 // Adapt and solve
 problem.adapt();
 //std::cout << "Self test: " << problem.self_test() << std::endl;
 problem.newton_solve();
 problem.doc_solution(doc_info);
 doc_info.number()++;
 
 problem.p_adapt();
 //std::cout << "Self test: " << problem.self_test() << std::endl;
 problem.newton_solve();
 problem.doc_solution(doc_info);
 doc_info.number()++;
 
} // end of main

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