Changeset 364 for trunk/gui/vizservers/nanovis

Timestamp:

Mar 14, 2006, 10:39:24 PM (19 years ago)

Author:

mmc

Message:

• Added src2, which is the start of Rappture 2.0, with various objects for handling complex data and serialization.

• Fixed Wei's visualization server to load data from a DX file with a structured mesh, using FieldRect3D to interpolate values.

File:

• trunk/gui/vizservers/nanovis/nanovis.cpp (modified) (3 diffs)

trunk/gui/vizservers/nanovis/nanovis.cpp

@@ -13 +13 @@
  * ======================================================================
  */
+#include <stdio.h>
+#include <math.h>
+#include <fstream>
+#include <iostream>
+#include <string>
 
 #include "nanovis.h"
+#include "RpFieldRect3D.h"
 
 ParticleSystem* psys;
@@ -161 +167 @@
 }
 
+/* Load a 3D volume from a dx-format file
+ */
+void
+load_volume_file(int index, char *fname) {
+    int dummy, nx, ny, nz, npts;
+    double x0, y0, z0, dx, dy, dz, ddx, ddy, ddz;
+    char line[128], *start;
+    std::ifstream fin(fname);
+
+    do {
+        fin.getline(line,sizeof(line)-1);
+        for (start=&line[0]; *start == ' ' || *start == '\t'; start++)
+            ;  // skip leading blanks
+
+        if (*start != '#') {  // skip comment lines
+            if (sscanf(start, "object %d class gridpositions counts %d %d %d", &dummy, &nx, &ny, &nz) == 4) {
+                // found grid size
+            }
+            else if (sscanf(start, "origin %lg %lg %lg", &x0, &y0, &z0) == 3) {
+                // found origin
+            }
+            else if (sscanf(start, "delta %lg %lg %lg", &ddx, &ddy, &ddz) == 3) {
+                // found one of the delta lines
+                if (ddx != 0.0) { dx = ddx; }
+                else if (ddy != 0.0) { dy = ddy; }
+                else if (ddz != 0.0) { dz = ddz; }
+            }
+            else if (sscanf(start, "object %d class array type double rank 0 items %d data follows", &dummy, &npts) == 2) {
+                if (npts != nx*ny*nz) {
+                    std::cerr << "inconsistent data: expected " << nx*ny*nz << " points but found " << npts << " points" << std::endl;
+                    return;
+                }
+                break;
+            }
+            else if (sscanf(start, "object %d class gridconnections", &dummy) == 1) {
+                // ignore this
+            }
+        }
+    } while (!fin.eof());
+
+    // read data points
+    if (!fin.eof()) {
+        Rappture::Mesh1D xgrid(x0, x0+nx*dx, nx);
+        Rappture::Mesh1D ygrid(y0, y0+ny*dy, ny);
+        Rappture::Mesh1D zgrid(z0, z0+nz*dz, nz);
+        Rappture::FieldRect3D field(xgrid, ygrid, zgrid);
+
+        double dval;
+        int nread = 0;
+        while (!fin.eof()) {
+            if (!(fin >> dval).fail()) {
+                field.define(nread++, dval);
+            }
+        }
+
+        // make sure that we read all of the expected points
+        if (nread != nx*ny*nz) {
+            std::cerr << "inconsistent data: expected " << nx*ny*nz << " points but found " << nread << " points" << std::endl;
+            return;
+        }
+
+        // figure out a good mesh spacing
+        int nsample = 100;
+        dx = field.rangeMax(Rappture::xaxis) - field.rangeMin(Rappture::xaxis);
+        dy = field.rangeMax(Rappture::yaxis) - field.rangeMin(Rappture::yaxis);
+        dz = field.rangeMax(Rappture::zaxis) - field.rangeMin(Rappture::zaxis);
+        double dmin = pow((dx*dy*dz)/(nsample*nsample*nsample), 0.333);
+
+        nx = (int)ceil(dx/dmin);
+        ny = (int)ceil(dy/dmin);
+        nz = (int)ceil(dz/dmin);
+        float *data = new float[nx*ny*nz];
+        std::cout << "generating " << nx << "x" << ny << "x" << nz << " = " << nx*ny*nz << " points" << std::endl;
+
+        // generate the uniformly sampled data that we need for a volume
+        int ngen = 0;
+        for (int iz=0; iz < nz; iz++) {
+            double zval = z0 + iz*dmin;
+            for (int iy=0; iy < ny; iy++) {
+                double yval = y0 + iy*dmin;
+                for (int ix=0; ix < nx; ix++) {
+                    double xval = x0 + ix*dmin;
+                    data[ngen++] = field.value(xval,yval,zval);
+                }
+            }
+        }
+        load_volume(index, nx, ny, nz, 1, data);
+        delete [] data;
+    } else {
+        std::cerr << "WARNING: data not found in file " << fname << std::endl;
+    }
+}
 
 /* Load a 3D volume
@@ -416 +514 @@
    }
 
-   init_vector_field(); //3d vector field
+   load_volume_file(0, "./data/A-apbs-2-out-potential-PE0.dx");
    init_fbo();  //frame buffer objects
    init_cg();   //init cg shaders
 
-   psys = new ParticleSystem(NMESH, NMESH, g_context, volume[0]->id);
-   init_particles();    //fill initial particles
-
-   get_slice_vectors();
+   //init_vector_field();       //3d vector field
+   //psys = new ParticleSystem(NMESH, NMESH, g_context, volume[0]->id);
+   //init_particles();  //fill initial particles
+
+   //get_slice_vectors();
 }