#include "AnT.hpp"

#define kappa parameters[0]

bool heat_conduction ( const CellularState& Temperature,
		       const Array<real_t>& parameters,
		       int cellIndex,
		       real_t deltaX,
		       StateCell& rhs)
{
  rhs[0] = kappa * D<0,0> (Temperature,cellIndex,0,deltaX);
  return true;
}

#undef kappa

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