//
//  marching_cubes.h
//
//  Created by Olga Diamanti on 05/03/12.

#ifndef IGL_MARCHINGCUBES_H
#define IGL_MARCHINGCUBES_H
#include "igl_inline.h"

#include <Eigen/Core>
namespace igl 
{
  // marching_cubes( values, points, x_res, y_res, z_res, vertices, faces )
  // performs marching cubes reconstruction on the grid defined by values, and points, and generates vertices and faces
  //
  // Input:
  //  xres, yres, zres: resolutions of the grid in x,y,z dimensions
  //  
  //  values: #number_of_grid_points x 1 array -- the scalar values of an implicit function defined on the grid points (<0 in the inside of the surface, 0 on the border, >0 outside)
  //  
  //  points: #number_of_grid_points x 3 array -- 3-D positions of the grid points, ordered in x,y,z order:
  //  points[index] = the point at (x,y,z) where :
  //  x = (index % (xres -1), 
  //  y = (index / (xres-1)) %(yres-1),
  //  z = index / (xres -1) / (yres -1) ).
  //  where x,y,z index x, y, z dimensions
  //  i.e. index = x + y*xres + z*xres*yres
  //  
  // Output:
  // vertices,faces: mesh description (vertices and faces)
  
  template <typename DerivedV, typename DerivedF>
  IGL_INLINE void marching_cubes(const Eigen::Matrix<typename DerivedV::Scalar, Eigen::Dynamic, 1> &values,
                                 const Eigen::Matrix<typename DerivedV::Scalar, Eigen::Dynamic, 3> &points,
                                 const unsigned x_res,
                                 const unsigned y_res,
                                 const unsigned z_res,
                                 Eigen::PlainObjectBase<DerivedV> &vertices,
                                 Eigen::PlainObjectBase<DerivedF> &faces);
}

#ifdef IGL_HEADER_ONLY
#  include "marching_cubes.cpp"
#endif

#endif