#ifndef ELEMENTUPDATE_2D_H
#define ELEMENTUPDATE_2D_H

#include "FixVec.h"
#include "ElementUpdate.h"
#include "VelocityExt.h"

using namespace AMDiS;

class ElementUpdate_2d : public ElementUpdate
{
public:
  ElementUpdate_2d(VelocityExt *velExt_ = NULL)
    : ElementUpdate(velExt_)
  {}

  /** \brief
   * Realization of ElementUpdate::calcElementUpdate. Calculates the Bornemann 
   * update for an element.
   */
  double calcElementUpdate(const FixVec<WorldVector<double> *, VERTEX> &vert,
			   FixVec<double, VERTEX> &uhVal);

private:
  WorldVector<double> xhminusYh, zhminusYh, xhminusZh;
};

#endif // ELEMENTUPDATE_2D_H