#include "bone_parents.h"

template <typename DerivedBE, typename DerivedP>
IGL_INLINE void igl::bone_parents(
  const Eigen::PlainObjectBase<DerivedBE>& BE,
  Eigen::PlainObjectBase<DerivedP>& P)
{
  P.resize(BE.rows(),1);
  // Stupid O(n²) version
  for(int e = 0;e<BE.rows();e++)
  {
    P(e) = -1;
    for(int f = 0;f<BE.rows();f++)
    {
      if(BE(e,0) == BE(f,1))
      {
        P(e) = f;
      }
    }
  }
}

#ifdef IGL_STATIC_LIBRARY
template void igl::bone_parents<Eigen::Matrix<int, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, 1, 0, -1, 1> >(Eigen::PlainObjectBase<Eigen::Matrix<int, -1, -1, 0, -1, -1> > const&, Eigen::PlainObjectBase<Eigen::Matrix<int, -1, 1, 0, -1, 1> >&);
#endif