//=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*
//  "isosurf.h"
//  3次元密度分布の等高面を描くライブラリ isosurf.h の利用例
//=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include "gluthack.h"
#include "isosurf.h"
#include "stdconf.h"

// 1辺あたりのデータ格子点の数
const int     N = 32;

const double dL = L/N;
inline double itol( int i ){
  return i*dL - L/2;
}
inline double sq( double d ){
  return d*d;
}

int main( int argc, char* argv[] )
{
  static double data[N][N][N];
  static double iso_value = 0.125;

  // 等高面の作成を管理するオブジェクトの宣言
  static Isosurf<N> Iso;

  // OpenGLの典型的なセットアップ
  stdSetup( argc, argv );

  for( int ix=0; ix<N; ix++ ){
    for( int iy=0; iy<N; iy++ ){
      for( int iz=0; iz<N; iz++ ){
	double x = itol(ix);
	double y = itol(iy);
	double z = itol(iz);
	double r = sqrt(x*x+y*y+z*z);
	data[ix][iy][iz] = sq((x*y)*exp(-r/3));
      }
    }
  }

  // 等高面作成オブジェクトにデータと等高値を渡す
  Iso.Assign( data, iso_value );

  // 等高面と空間の枠を display listに作成する
  glNewList( 1, GL_COMPILE );{
    const GLfloat  OBJ_COL[]   = { 0.5, 0.5, 0.8, 0.5 };  // 物体の拡散色
    const GLfloat  OBJ_SPE[]   = { 1.0, 1.0, 1.0, 1.0 };  // 物体の反射色
    const GLfloat  OBJ_SHIN    = 256.0;                   // 物体の反射強度
    glMaterialfv( GL_FRONT_AND_BACK, GL_SPECULAR,            OBJ_SPE  );
    glMaterialfv( GL_FRONT_AND_BACK, GL_AMBIENT_AND_DIFFUSE, OBJ_COL  );
    glMaterialf(  GL_FRONT_AND_BACK, GL_SHININESS,           OBJ_SHIN );

    glPushMatrix();{
      glTranslated( -L/2, -L/2, -L/2 );
      glScaled( L, L, L );

      // 等高面作成オブジェクトによる等高面の作成
      Iso.Make();

      // 空間の枠の作成
      glDisable( GL_LIGHTING );
      glutBoxWire( 1.0, 1.0, 1.0 );
      glEnable( GL_LIGHTING );
    }glPopMatrix();

  }glEndList();

  glutMainLoop();

  return 0;
}