eluting_sphere_adapt.cc File Reference

#include "generic.h"
#include "spherical_advection_diffusion.h"
#include "meshes/rectangular_quadmesh.h"

Classes
class	RefineableElutingSphereProblem< ELEMENT >

Namespaces
	Global_Physical_Variables
	Global variables.
	Wind_Function

Functions
void	Wind_Function::get_wind (const Vector< double > &x, Vector< double > &wind)
	Wind. More...
int	main ()
	Driver for RefineableElutingSphere test problem using a. More...

Function Documentation

main()

int main

(

)

Driver for RefineableElutingSphere test problem using a.

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

{
 // Label for output
 DocInfo doc_info;
 
 // Set output directory
 doc_info.set_directory("RESLT_adapt");
 
 // ---------------
 // end of Set up doc info
 
 // Build the problem with QCrouzeixRaviartElements
 RefineableElutingSphereProblem<
 RefineableQSphericalAdvectionDiffusionElement<3> > 
  problem;
 
 // Over-ride the maximum and minimum permitted errors
 problem.mesh_pt()->max_permitted_error() = 1.0e-3; //Default = 1.0e-3
 problem.mesh_pt()->min_permitted_error() = 1.0e-5; //Default = 1.0e-5
 
 // Over-ride the maximum and minimum permitted refinement levels
 problem.mesh_pt()->max_refinement_level() = 4;//maximum_ref_level;
 problem.mesh_pt()->min_refinement_level() = 1;//minimum_ref_level;
 
  //Set the boundary conditions
 problem.set_boundary_conditions();
 
 //Set the maximum adaptation
 unsigned max_adapt = 10;
  
 // Solve the problem
 problem.steady_newton_solve(max_adapt);
 problem.doc_solution(doc_info);
 
 
} // end_of_main

References problem, and oomph::DocInfo::set_directory().