#include "AnT.hpp"

#define a parameters[0]
#define u uState[cellIndex][0]

bool KdV ( const CellularState& uState,
	   const Array<real_t>& parameters,
	   int cellIndex,
	   real_t deltaX,
	   StateCell& rhs)
{

  rhs[0] = a * u * D<0> (uState,cellIndex,0,deltaX) - 
    D<0,0,0> (uState,cellIndex,0,deltaX);

  return true;
}


#undef a
#undef u

extern "C" 
{ 
  void connectSystem () 
  { 
    PDE_1d_Proxy::systemFunction = KdV; 
  }
}