#ifndef MESH_BOOLEAN_H
#define MESH_BOOLEAN_H

#include <igl/igl_inline.h>
#include "MeshBooleanType.h"
#include <Eigen/Core>
#include <functional>

namespace igl
{
  //  MESH_BOOLEAN Compute boolean csg operations on "solid", consistently oriented
  //  meshes.
  // 
  //  Inputs:
  //    V  #V by 3 list of vertex positions of first mesh
  //    F  #F by 3 list of triangle indices into V
  //    U  #U by 3 list of vertex positions of second mesh
  //    G  #G by 3 list of triangle indices into U
  //    type  type of boolean operation
  //  Outputs:
  //    W  #W by 3 list of vertex positions of boolean result mesh
  //    H  #H by 3 list of triangle indices into W
  //    J  #H list of indices into [FA;FB] revealing "birth" facet
  //  
  //  See also: self_intersect
  //     
  template <
    typename DerivedVA,
    typename DerivedFA,
    typename DerivedVB,
    typename DerivedFB,
    typename DerivedVC,
    typename DerivedFC,
    typename DerivedJ>
  IGL_INLINE void mesh_boolean(
    const Eigen::PlainObjectBase<DerivedVA > & VA,
    const Eigen::PlainObjectBase<DerivedFA > & FA,
    const Eigen::PlainObjectBase<DerivedVB > & VB,
    const Eigen::PlainObjectBase<DerivedFB > & FB,
    const MeshBooleanType & type,
    Eigen::PlainObjectBase<DerivedVC > & VC,
    Eigen::PlainObjectBase<DerivedFC > & FC,
    Eigen::PlainObjectBase<DerivedJ > & J);
  template <
    typename DerivedVA,
    typename DerivedFA,
    typename DerivedVB,
    typename DerivedFB,
    typename DerivedVC,
    typename DerivedFC>
  IGL_INLINE void mesh_boolean(
    const Eigen::PlainObjectBase<DerivedVA > & VA,
    const Eigen::PlainObjectBase<DerivedFA > & FA,
    const Eigen::PlainObjectBase<DerivedVB > & VB,
    const Eigen::PlainObjectBase<DerivedFB > & FB,
    const MeshBooleanType & type,
    Eigen::PlainObjectBase<DerivedVC > & VC,
    Eigen::PlainObjectBase<DerivedFC > & FC);
  // Inputs:
  //   resolve_fun  function handle for computing resolve of a
  //     self-intersections of a mesh and outputting the new mesh.
  template <
    typename DerivedVA,
    typename DerivedFA,
    typename DerivedVB,
    typename DerivedFB,
    typename DerivedVC,
    typename DerivedFC,
    typename DerivedJ>
  IGL_INLINE void mesh_boolean(
    const Eigen::PlainObjectBase<DerivedVA > & VA,
    const Eigen::PlainObjectBase<DerivedFA > & FA,
    const Eigen::PlainObjectBase<DerivedVB > & VB,
    const Eigen::PlainObjectBase<DerivedFB > & FB,
    const MeshBooleanType & type,
    const std::function<void(
      const Eigen::Matrix<typename DerivedVC::Scalar,Eigen::Dynamic,3> &,
      const Eigen::Matrix<typename DerivedFC::Scalar,Eigen::Dynamic,3> &,
            Eigen::Matrix<typename DerivedVC::Scalar,Eigen::Dynamic,3> &,
            Eigen::Matrix<typename DerivedFC::Scalar,Eigen::Dynamic,3> &,
            Eigen::Matrix<typename  DerivedJ::Scalar,Eigen::Dynamic,1>&)> 
      & resolve_fun,
    Eigen::PlainObjectBase<DerivedVC > & VC,
    Eigen::PlainObjectBase<DerivedFC > & FC,
    Eigen::PlainObjectBase<DerivedJ > & J);
}

#ifndef IGL_STATIC_LIBRARY
#  include "mesh_boolean.cpp"
#endif

#endif