namespace nu {

  typedef float real_t;

  inline real_t to_radian(real_t degree) {
    return degree * (real_t)3.14159265358979 / (real_t)180.0;
  }

  void mytranslate(
    real_t* m,
    real_t x,
    real_t y,
    real_t z
  ) {
    m[0] = 1.0;
    m[1] = 0.0;
    m[2] = 0.0;
    m[3] = 0.0;

    m[4] = 0.0;
    m[5] = 1.0;
    m[6] = 0.0;
    m[7] = 0.0;

    m[8] = 0.0;
    m[9] = 0.0;
    m[10] = 1.0;
    m[11] = 0.0;

    m[12] = x;
    m[13] = y;
    m[14] = z;
    m[15] = 1.0;

  }



  void myscale(
    real_t* m,
    real_t sx,
    real_t sy,
    real_t sz
  ) {
    m[0] = sx;
    m[1] = 0.0;
    m[2] = 0.0;
    m[3] = 0.0;

    m[4] = 0.0;
    m[5] = sy;
    m[6] = 0.0;
    m[7] = 0.0;

    m[8] = 0.0;
    m[9] = 0.0;
    m[10] = sz;
    m[11] = 0.0;

    m[12] = 0.0;
    m[13] = 0.0;
    m[14] = 0.0;
    m[15] = 1.0;

  }

}

#include <windows.h>
#include <gl/GLU.h>
#include <gl/freeglut.h>
#include <cstdio>
#include <cmath>

//検証用
#include <cassert>

//自前実装版glTranslate* / glScale*
#include "mytranslatescale.h"

#pragma comment(lib,"opengl32.lib")
#pragma comment(lib,"glu32.lib")
#pragma comment(lib,"freeglut.lib")

void display(void)
{
  
  GLfloat glmat[16]; //gl関数で生成
  GLfloat mymat[16]; //自分関数で生成

  //gltranslate,mytranslateへの入力の値を色々変える
  static int i = 0;
  i++;

  GLfloat x = 1.0+i;
  GLfloat y = 1.0 * (i * 0.1);
  GLfloat z = 1.0 * (i / 0.2);

  glMatrixMode(GL_MODELVIEW);
  glLoadIdentity();

  //openglで計算
  glTranslatef( x, y, z);
  glGetFloatv(GL_MODELVIEW_MATRIX, glmat);

  //自前計算
  nu::mytranslate(mymat,x,y,z);

  puts("----------------------------------------------");
  printf("** %d :: %lf %lf %lf\n", i,x,y,z);
  for (int k = 0; k < 16; k++) {
    //比較
    printf("m[%2d] %+5.10lf %+5.10lf \n", k, glmat[k], mymat[k]);

    //自前計算とgltranslate計算の値があまりにかけ離れていたらassertで落ちる
    assert(std::abs(glmat[k] - mymat[k]) < 0.0001);

  }
  puts("----------------------------------------------");

}
void mouse(int button, int state, int x, int y) {
  display();
}

int main(int argc, char* argv[])
{

  glutInit(&argc, argv);
  glutInitDisplayMode(GLUT_RGBA);
  glutCreateWindow(argv[0]);
  glutInitWindowPosition(0, 0);
  glutInitWindowSize(100,100);
  glutDisplayFunc(display);
  glutMouseFunc(mouse);
  glutMainLoop();

  return 0;

}

#include <windows.h>
#include <gl/GLU.h>
#include <gl/freeglut.h>
#include <cstdio>
#include <cmath>

//検証用
#include <cassert>

//自前実装版glTranslate* / glScale*
#include "mytranslatescale.h"

#pragma comment(lib,"opengl32.lib")
#pragma comment(lib,"glu32.lib")
#pragma comment(lib,"freeglut.lib")

void display(void)
{
  
  GLfloat glmat[16]; //gl関数で生成
  GLfloat mymat[16]; //自分関数で生成

  //glscalef,myscaleへの入力の値を色々変える
  static int i = 0;
  i++;

  GLfloat x = 1.0+i;
  GLfloat y = 1.0 * (i * 0.1);
  GLfloat z = 1.0 * (i / 0.2);

  glMatrixMode(GL_MODELVIEW);
  glLoadIdentity();

  //openglで計算
  glScalef( x, y, z);
  glGetFloatv(GL_MODELVIEW_MATRIX, glmat);

  //自前計算
  nu::myscale(mymat,x,y,z);

  puts("----------------------------------------------");
  printf("** %d :: %lf %lf %lf\n", i,x,y,z);
  for (int k = 0; k < 16; k++) {
    //比較
    printf("m[%2d] %+5.10lf %+5.10lf \n", k, glmat[k], mymat[k]);

    //自前計算とglScalef計算の値があまりにかけ離れていたらassertで落ちる
    assert(std::abs(glmat[k] - mymat[k]) < 0.0001);

  }
  puts("----------------------------------------------");

}
void mouse(int button, int state, int x, int y) {
  display();
}

int main(int argc, char* argv[])
{

  glutInit(&argc, argv);
  glutInitDisplayMode(GLUT_RGBA);
  glutCreateWindow(argv[0]);
  glutInitWindowPosition(0, 0);
  glutInitWindowSize(100,100);
  glutDisplayFunc(display);
  glutMouseFunc(mouse);
  glutMainLoop();

  return 0;

}

#include <cmath>

namespace nu {

  typedef float real_t;

  inline real_t to_radian(real_t degree) {
    return degree * (real_t)3.14159265358979 / (real_t)180.0;
  }

  //! @brief 回転行列作成
  //! @param [out] m 結果を格納する要素数16の配列
  //! @param [in] angle_degree 回転角を度で指定
  //! @param [in] x 回転軸のX成分
  //! @param [in] y 回転軸のY成分
  //! @param [in] z 回転軸のZ成分
  //! @return なし
  void myrotate(
    real_t* m,
    real_t angle_degree,
    real_t x,
    real_t y,
    real_t z) {

    //len(x y z) != 0ならnormalizeする
    real_t len = sqrt(x * x + y * y + z * z);
    if (abs(1.0 - len) > 0.000001) {
      x = x / len;
      y = y / len;
      z = z / len;
    }


    real_t angle_rad = to_radian(angle_degree);

    real_t c = std::cos(angle_rad);
    real_t s = std::sin(angle_rad);

    m[0] = x * x * (1 - c) + c;
    m[1] = y * x * (1 - c) + z * s;
    m[2] = x * z * (1 - c) - y * s;
    m[3] = 0.0;

    m[4] = x * y * (1 - c) - z * s;
    m[5] = y * y * (1 - c) + c;
    m[6] = y * z * (1 - c) + x * s;
    m[7] = 0.0;


    m[8] = x * z * (1 - c) + y * s;
    m[9] = y * z * (1 - c) - x * s;
    m[10] = z * z * (1 - c) + c;
    m[11] = 0.0;


    m[12] = 0.0;
    m[13] = 0.0;
    m[14] = 0.0;
    m[15] = 1.0;

  }

}