two_d_act_inhibit.cc File Reference

#include <iostream>
#include <fstream>
#include <cmath>
#include "generic.h"
#include "navier_stokes.h"
#include "advection_diffusion_reaction.h"
#include "meshes/rectangular_quadmesh.template.cc"

Classes
class	ActivatorInhibitorProblem< 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 ()

Variables
const double	GlobalVariables::pi = MathematicalConstants::Pi
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

(

)

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

  {
   double epsilon = 0.05;
   //Set the control parameters (should produce a spot)
   GlobalVariables::Tau[1] = 0.5;
   GlobalVariables::D[1] = epsilon*epsilon;
   GlobalVariables::A = -0.1;
 
   //Set the timestep
   double dt = 0.1;
 
   ActivatorInhibitorProblem 
    <QAdvectionDiffusionReactionElement<2,2,3> >problem(20,20);
   problem.timestep(dt,5);
  }

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