...

[phys/latlib.git] / culooks_cube.cpp

diff --git a/culooks_cube.cpp b/culooks_cube.cpp
new file mode 100644 (file)
index 0000000..2498fef
--- /dev/null
+++ b/culooks_cube.cpp
@@ -0,0 +1,97 @@
+#include "culooks_cube.h"
+
+#include <GL/glut.h>
+
+int cube::allid = 0;
+
+void cube::draw()
+{
+  glMatrixMode(GL_MODELVIEW);
+  glPushMatrix();
+
+  glScalef(zoom, zoom, zoom);
+  glRotatef(az, 0, 1, 0);
+  glRotatef(alt, 1, 0, 0);
+
+  drawAll();
+
+  glPopMatrix();
+}
+
+void cube::drawAll()
+{
+  glMatrixMode(GL_MODELVIEW);
+  glPushMatrix();
+
+  glTranslatef(-1,-1,-1);
+
+  for (int iz=0; iz<l; iz++) {
+    for (int iy=0; iy<l; iy++) {
+      for (int ix=0; ix<l; ix++) {
+       glPushMatrix();
+       glTranslatef(2.0*ix/l,2.0*iy/l,2.0*iz/l);
+       glScalef(2.0/l,2.0/l,2.0/l);
+
+       /* draw links */
+       glBegin(GL_LINES);
+       glColor4f(link[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 0*4 + 0],
+                 link[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 0*4 + 1],
+                 link[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 0*4 + 2],
+                 link[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 0*4 + 3]);
+       glVertex3f(0,0,0); glVertex3f(1,0,0);
+       glColor4f(link[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 1*4 + 0],
+                 link[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 1*4 + 1],
+                 link[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 1*4 + 2],
+                 link[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 1*4 + 3]);
+       glVertex3f(0,0,0); glVertex3f(0,1,0);
+       glColor4f(link[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 2*4 + 0],
+                 link[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 2*4 + 1],
+                 link[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 2*4 + 2],
+                 link[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 2*4 + 3]);
+       glVertex3f(0,0,0); glVertex3f(0,0,1);
+       glEnd();
+
+       /* draw plaquettes */
+       glColor4f(1,1,1,0.5);
+       glBegin(GL_QUADS);
+       glColor4f(plaq[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 2*4 + 0],
+                 plaq[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 2*4 + 1],
+                 plaq[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 2*4 + 2],
+                 plaq[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 2*4 + 3]);
+       glVertex3f(0, 0, 0); glVertex3f(1, 0, 0); glVertex3f(1, 1, 0); glVertex3f(0, 1, 0);
+       glColor4f(plaq[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 1*4 + 0],
+                 plaq[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 1*4 + 1],
+                 plaq[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 1*4 + 2],
+                 plaq[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 1*4 + 3]);
+       glVertex3f(0, 0, 0); glVertex3f(1, 0, 0); glVertex3f(1, 0, 1); glVertex3f(0, 0, 1);
+       glColor4f(plaq[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 0*4 + 0],
+                 plaq[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 0*4 + 1],
+                 plaq[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 0*4 + 2],
+                 plaq[iz*l*l*3*4 + iy*l*3*4 + ix*3*4 + 0*4 + 3]);
+       glVertex3f(0, 0, 0); glVertex3f(0, 1, 0); glVertex3f(0, 1, 1); glVertex3f(0, 0, 1);
+       glEnd();
+
+       glPopMatrix();
+      }
+    }
+  }
+    
+  glPopMatrix();
+}
+
+cube::cube(int _l)
+{
+  l = _l;
+  link = new float[l*l*l*3 * 4];
+  plaq = new float[l*l*l*3 * 4];
+  az = 30;
+  alt = -20;
+  zoom = 0.6;
+  id = allid;
+  allid++;
+
+  for (int i=0; i<l*l*l*3; i++) {
+    link[i*4+0]=1; link[i*4+3]=0.6;
+    plaq[i*4+1]=1; plaq[i*4+3]=0.2;
+  }
+}