//! @tparam Point2D 二次元座標を一つ表現する。operator[](int)が必要
//! @brief 多角形を作成する
//! @param [out] points 作成した多角形の頂点
//! @param [in] N 頂点数
//! @param [in] radius 中心から頂点までの距離
//! @return なし
template<typename Point2D>
void createNgon(
  std::vector< Point2D >& points,
  const int N,
  const double radius = 1.0
) {


  for (size_t i = 0; i < N; i++) {
    double t = 2 * 3.14 / N * i;
    double x = sin(t);
    double y = cos(t);

    Point2D p;
    p[0] = x * radius;
    p[1] = y * radius;
    points.push_back(p);
  }

}

  //多角形作成
  std::vector<std::array<double, 2>> ngon;
  createNgon(ngon,3);
  
  //表示
  glPointSize(1);
  glBegin(GL_LINE_LOOP);
  for (size_t i = 0; i < ngon.size(); i++) {
    glVertex2dv(ngon[i].data());
  }
  glEnd();

#pragma once
#include <vector>
#include <array>
#include <limits>
#include <valarray>

template<typename real_t>
real_t Length3(const real_t * const vec) {
  const int X = 0;
  const int Y = 1;
  const int Z = 2;
  return sqrt(vec[X] * vec[X] + vec[Y] * vec[Y] + vec[Z] * vec[Z]);
}
template<typename real_t>
bool Normalize3(real_t* const pV)
{

  const int X = 0;
  const int Y = 1;
  const int Z = 2;

  real_t len;
  real_t& x = pV[X];
  real_t& y = pV[Y];
  real_t& z = pV[Z];

  len = static_cast<real_t>(sqrt(x * x + y * y + z * z));

  if (len < static_cast<real_t>(1e-6)) return false;

  len = static_cast<real_t>(1.0) / len;
  x *= len;
  y *= len;
  z *= len;

  return true;
}

template<typename real_t>
real_t Inner3(real_t const* const vec1, real_t const* const vec2) {
  const int X = 0;
  const int Y = 1;
  const int Z = 2;
  return ((vec1[X]) * (vec2[X]) + (vec1[Y]) * (vec2[Y]) + (vec1[Z]) * (vec2[Z]));
}