/**
 * @file   stencil.cpp
 * @author Julian Iseringhausen <isering@ins.uni-bonn.de>
 * @date   Mon Apr 18 12:24:43 2011
 *
 * @brief  VMG::Stencil.
 *
 */


#ifdef HAVE_CONFIG_H
#include <config.h>
#endif

#include "base/stencil.hpp"
#include "grid/grid.hpp"

using namespace VMG;

vmg_float Stencil::Apply(const Grid& grid, const Index& index) const
{
  vmg_float result = diag * grid.GetVal(index);

  for (Stencil::iterator iter=disp.begin(); iter!=disp.end(); ++iter)
    result += iter->val * grid.GetVal(index + iter);

  return result;
}