Changeset 6518 in orxonox.OLD for trunk/src/world_entities/environments/water.cc

Timestamp:

Jan 18, 2006, 12:09:11 AM (18 years ago)

Author:

bensch

Message:

orxonox/trunk: some waves

File:

• trunk/src/world_entities/environments/water.cc (modified) (5 diffs)

trunk/src/world_entities/environments/water.cc

@@ -44 +44 @@
   this->grid = NULL;
 
+  this->velocities = NULL;
+  this->viscosity = 5;
+  this->cohesion = .00002;
+
   if (root != NULL)
     this->loadParams(root);
@@ -51 +55 @@
   this->waterShader = (Shader*)ResourceManager::getInstance()->load("shaders/water.vert", SHADER, RP_GAME, (void*)"shaders/water.frag");
 
+  this->grid->height(this->grid->columns()/2,this->grid->rows()/2) = 1000;
 }
 
@@ -78 +83 @@
 void Water::rebuildGrid()
 {
+  if (this->velocities != NULL)
+  {
+    assert (this->grid != NULL);
+    for (unsigned int i = 0; i < this->grid->rows(); i++)
+      delete[] this->velocities[i];
+    delete[] this->velocities;
+  }
   if (this->grid != NULL)
     delete this->grid;
   this->grid = new Grid(this->sizeX, this->sizeY, this->resX, this->resY);
+  this->velocities = new float*[this->resX];
+  for (unsigned int i = 0; i < this->grid->rows(); i++)
+    this->velocities[i] = new float[this->resY];
 
   this->setModel(this->grid, 0);
@@ -117 +132 @@
 void Water::tick(float dt)
 {
+/*
+    THE OLD USELESS ALGORITHM
   phase += dt *.1;
   for (unsigned int i = 0; i < this->grid->rows(); i++)
@@ -126 +143 @@
     }
   }
-  this->grid->rebuildNormals(this->height);
-}
+  this->grid->rebuildNormals(this->height);*/
+
+
+  unsigned int i, j;
+  float u;
+
+  // wave/advection
+  // calc movement
+  for(j = 1; j < this->grid->rows() - 1; j++) {
+    for(i = 1; i < this->grid->columns() - 1; i++) {
+      u =  this->grid->height(i+1,j)+ this->grid->height(i-1, j) +
+          this->grid->height(i, j+1) + this->grid->height(i, j-1) -
+          4 * this->grid->height(i, j);
+      this->velocities[i][j] += dt * this->viscosity * this->viscosity * u / this->height;
+    }
+  }
+  // advect
+  for(j = 1; j < this->grid->rows() - 1; j++) {
+    for(i = 1; i < this->grid->columns() - 1; i++) {
+      this->grid->height(i, j) += dt * this->velocities[i][j];
+    }
+  }
+  // bound
+//   unsigned int w = this->grid->columns - 1;
+//   for(i = 0; i < this->grid->columns; i++) {
+//     _map[i][0].u[1] = _map[i][1  ].u[1];
+//     _map[i][w].u[1] = _map[i][w-1].u[1];
+//     _map[0][i].u[1] = _map[1  ][i].u[1];
+//     _map[w][i].u[1] = _map[w-1][i].u[1];
+//   }
+
+  // diffusion
+  for(j = 1; j < this->grid->rows() - 1; j++) {
+    for(i = 1; i < this->grid->columns() - 1 ; i++) {
+      u = this->grid->height(i+1, j) + this->grid->height(i-1, j) +
+          this->grid->height(i, j+1) + this->grid->height(i, j-1) -
+          4* this->grid->height(i, j);
+      this->grid->height(i,j) += dt * this->cohesion * u / this->height;
+    }
+  }
+
+  // calc normals
+/*  float l[3];
+  float m[3];
+  for(j = 1; j < _width-1; j++) {
+    for(i = 1; i < _width - 1; i++) {
+      l[0] = _map[i][j-1].u[0] - _map[i][j+1].u[0];
+      l[1] = _map[i][j-1].u[1] - _map[i][j+1].u[1];
+      l[2] = _map[i][j-1].u[2] - _map[i][j+1].u[2];
+      m[0] = _map[i-1][j].u[0] - _map[i+1][j].u[0];
+      m[1] = _map[i-1][j].u[1] - _map[i+1][j].u[1];
+      m[2] = _map[i-1][j].u[2] - _map[i+1][j].u[2];
+      _map[i][j].n[0] = l[1]*m[2]-l[2]*m[1];
+      _map[i][j].n[1] = l[2]*m[0]-l[0]*m[2];
+      _map[i][j].n[2] = l[0]*m[1]-l[1]*m[0];
+    }
+  }*/
+}