one_d_act_inhibit.cc File Reference

#include "generic.h"
#include "advection_diffusion_reaction.h"
#include "meshes/one_d_mesh.h"

Classes
class	RefineableOneDAdvectionDiffusionReactionProblem< ELEMENT >

Namespaces
	GlobalVariables

Functions
void	GlobalVariables::activator_inhibitor_reaction (const Vector< double > &C, Vector< double > &R)
void	GlobalVariables::activator_inhibitor_reaction_derivative (const Vector< double > &C, DenseMatrix< double > &dRdC)
	Derivative of simple reaction kinetics above. More...
int	main ()
	Driver code for 1D AdvectionDiffusionReaction problem. More...

Variables
Vector< double >	GlobalVariables::Tau (2, 1.0)
	The vector of timescales. More...
Vector< double >	GlobalVariables::D (2, 1.0)
	The vector of diffusion coefficients. More...
double	GlobalVariables::A = -0.1
	Bifurcation (loss) parameter. More...

Function Documentation

main()

int main

(

)

Driver code for 1D AdvectionDiffusionReaction problem.

{
 //Diffusive length-scale
 double epsilon = 0.05;
 //Set the control parameters
 GlobalVariables::Tau[1] = 0.2;
 GlobalVariables::D[1] = epsilon*epsilon;
 GlobalVariables::A = -0.4;
 
 //Set the (largeish) timestep
 double dt = 0.1;
 
 //Set up the problem
 //------------------
 // DREIGIAU: There's an inherent problem with the 1D elements. Doesn't
 // seem to be a problem with 2D elements but we can't even initialise
 // a new object of the type used to template the problem below.
 
 // Create the problem with 1D three-node refineable elements from the
 // RefineableLineAdvectionDiffusionReactionElement family. 
 RefineableOneDAdvectionDiffusionReactionProblem<
  RefineableQAdvectionDiffusionReactionElement<2,1,3> > problem;
 
 //Take four levels of uniform refinement to start things off
 for(unsigned i=0;i<4;i++) { problem.refine_uniformly(); }
 
 //Now timestep the problem
 problem.timestep(dt,10);
 
} // End of main

References GlobalVariables::A, GlobalVariables::D, oomph::SarahBL::epsilon, i, problem, and GlobalVariables::Tau.