unstructured_two_d_act_inhibit.cc File Reference

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

Classes
class	RefineableActivatorInhibitorProblem< 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 2D ActivatorInhibitor 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...

Function Documentation

main()

int main

(

)

Driver code for 2D ActivatorInhibitor 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.25;
 
 //Construct the problem
 RefineableActivatorInhibitorProblem 
   <ProjectableAdvectionDiffusionReactionElement<TAdvectionDiffusionReactionElement<2,2,3> > > problem;
 
 //Timestep it
 problem.timestep(dt,2);
}

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