/**
 * @file   interface_vtk.cpp
 * @author Julian Iseringhausen <isering@ins.uni-bonn.de>
 * @date   Mon Apr 18 12:58:39 2011
 *
 * @brief  VMG::InterfaceVTK
 *
 */

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

#ifdef HAVE_VTK

#include "interface_vtk.hpp"

void MGInterfaceVTK::ImportRightHandSide(MGMultigrid& multigrid)
{
  MGGrid& grid = multigrid(multigrid.GlobalMaxLevel());

  ptr<vtkStructuredPointsReader> source = ptr<vtkStructuredPointsReader>::New();

  source->SetFileName(filename.c_str());
  source->Update();

  for (int i=0; i<grid.Local().SizeTotal().X(); ++i)
    for (int j=0; j<grid.Local().SizeTotal().Y(); ++j)
      for (int k=0; k<grid.Local().SizeTotal().Z(); ++k)
        grid(i,j,k) = source->GetOutput()->GetScalarComponentAsDouble(i,j,k,0);
}

void MGInterfaceVTK::ExportSolution(MGGrid& grid)
{
}

#endif