#ifndef IGL_PSO_H
#define IGL_PSO_H
#include "igl_inline.h"
#include <Eigen/Core>
#include <functional>

namespace igl
{
  /// Global optimization with the particle swarm algorithm.
  ///
  /// Solve the problem:
  ///
  ///   minimize f(x)
  ///   subject to lb ≤ x ≤ ub 
  /// 
  /// by particle swarm optimization (PSO).
  ///
  /// @param[in] f  function that evaluates the objective for a given "particle" location
  /// @param[in] LB  #X vector of lower bounds 
  /// @param[in] UB  #X vector of upper bounds 
  /// @param[in] max_iters  maximum number of iterations
  /// @param[in] population  number of particles in swarm
  /// @param[out] X  best particle seen so far
  /// @return objective corresponding to best particle seen so far
  template <
    typename Scalar, 
    typename DerivedX,
    typename DerivedLB, 
    typename DerivedUB>
  IGL_INLINE Scalar pso(
    const std::function< Scalar (DerivedX &) > f,
    const Eigen::MatrixBase<DerivedLB> & LB,
    const Eigen::MatrixBase<DerivedUB> & UB,
    const int max_iters,
    const int population,
    DerivedX & X);
  /// \overload
  /// @param[out] P  whether each DOF is periodic
  ///
  /// \bug `P` appears to be unused
  template <
    typename Scalar, 
    typename DerivedX,
    typename DerivedLB, 
    typename DerivedUB,
    typename DerivedP>
  IGL_INLINE Scalar pso(
    const std::function< Scalar (DerivedX &) > f,
    const Eigen::MatrixBase<DerivedLB> & LB,
    const Eigen::MatrixBase<DerivedUB> & UB,
    const Eigen::DenseBase<DerivedP> & P,
    const int max_iters,
    const int population,
    DerivedX & X);
}

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

#endif