コード例 #1
0
  /** this test compares result of MT and ST distance calculations */
  void devTestMTvsSTprecision() throws Exception {

    if (DEBUG) printf("%s.testPrecisionMTvsST()\n", this.getClass().getName());
    double vs = 0.5 * MM;
    double w = 10 * MM; // half width
    // double firstLayerThickness = 0.7;
    // double firstLayerThickness = 1.7;
    double firstLayerThickness = 2.7;
    // double xmin = -w, xmax = w, ymin = -w, ymax = w, zmin = -vs, zmax = 2*vs;
    double xmin = -w, xmax = w, ymin = -w, ymax = w, zmin = -w, zmax = w;
    int densityBitCount = 8;
    int distanceBitCount = 16;

    int subvoxelResolution = (1 << densityBitCount) - 1;

    int voxelSquareSize = 20; // for visualization
    boolean snapToGrid = false;
    int iterationsCount = 0;
    int threadCount = 4;
    boolean writeViz = false;
    boolean compareDistances = false;

    Bounds bounds = new Bounds(xmin, xmax, ymin, ymax, zmin, zmax);
    double pnts[];
    pnts = makeCircleZ(0.5 * vs, 0.5 * vs, 0.5 * vs, 0.9 * w, 128);
    // double pnts[] = makeCircleX(0.5*vs, 0.5*vs, 0.5*vs, 0.9*w, 128);
    // double pnts[]  = makeCircleY(0.5*vs, 0.5*vs, 0.5*vs, 0.9*w, 128);
    // pnts = makeUnion(makeUnion(makeCircleX(0.5*vs, 0.5*vs, 0.5*vs, 0.9*w,
    // 128),makeCircleY(0.5*vs, 0.5*vs, 0.5*vs, 0.9*w, 128)),makeCircleZ(0.5*vs, 0.5*vs, 0.5*vs,
    // 0.9*w, 128));
    // double pnts[] = makeUnion(makeUnion(makeCircleZ(0.5*vs, 0.25*vs, 0.5*w, 0.9*w,
    // 128),makeCircleZ(0.5*vs, 0.5*vs, 0.5*w, 0.9*w, 128)),makeCircleZ(0.5*vs, 0.5*vs, 0.5*vs,
    // 0.9*w, 128));

    ArrayAttributeGridInt indexGrid1 = new ArrayAttributeGridInt(bounds, vs, vs);
    ArrayAttributeGridInt indexGrid2 = new ArrayAttributeGridInt(bounds, vs, vs);

    int pcount = pnts.length / 3;
    double pntx[] = new double[pcount];
    double pnty[] = new double[pcount];
    double pntz[] = new double[pcount];

    ClosestPointIndexer.getPointsInGridUnits(indexGrid1, pnts, pntx, pnty, pntz);
    if (snapToGrid) {
      ClosestPointIndexer.snapToVoxels(pntx);
      ClosestPointIndexer.snapToVoxels(pnty);
      ClosestPointIndexer.snapToVoxels(pntz);
    }

    ClosestPointIndexer.initFirstLayer(indexGrid1, pntx, pnty, pntz, firstLayerThickness);
    ClosestPointIndexer.initFirstLayer(indexGrid2, pntx, pnty, pntz, firstLayerThickness);

    for (int z = 0; z < indexGrid1.getDepth(); z++) {
      renderDiff(
          indexGrid1,
          z,
          pntx,
          pnty,
          pntz,
          voxelSquareSize,
          fmt("/tmp/dist/dist00_%02d.png", z),
          true);
    }

    int usedCount = ClosestPointIndexer.removeUnusedPoints(indexGrid1, pntx, pnty, pntz);
    printf(
        "grid: [%d x %d x %d] threads: %d points: %d usedPoints: %d \n",
        indexGrid1.getWidth(),
        indexGrid1.getHeight(),
        indexGrid1.getDepth(),
        threadCount,
        pcount,
        usedCount);

    // distribute distances to the whole grid
    long t0 = time();
    ClosestPointIndexer.PI3_sorted(pntx, pnty, pntz, indexGrid1);
    printf("ClosestPointIndexer done: %d\n", (time() - t0));
    t0 = time();
    ClosestPointIndexerMT.PI3_MT(pntx, pnty, pntz, indexGrid2, threadCount);
    printf("ClosestPointIndexerMT done: %d\n ", (time() - t0));
    long diff = GridUtil.compareGrids(indexGrid1, indexGrid2);
    printf("index difference count: %d\n", diff);

    printf("maxDiff1: %6.4f\n", getMaxDiff(indexGrid1, pntx, pnty, pntz));
    printf("maxDiff2: %6.4f\n", getMaxDiff(indexGrid1, pntx, pnty, pntz));

    if (writeViz) {
      for (int z = 0; z < indexGrid1.getDepth(); z++) {
        // renderDiff(indexGrid1, z, pntx, pnty, pntz, voxelSquareSize,
        // fmt("/tmp/dist/dist01_%02d.png",z), true);
        renderDiff(
            indexGrid2,
            z,
            pntx,
            pnty,
            pntz,
            voxelSquareSize,
            fmt("/tmp/dist/distDiff2_%02d.png", z),
            true);
      }
    }

    if (compareDistances) {
      AttributeGrid distGrid1 = makeDistanceGrid(bounds, -w, w, distanceBitCount, vs);
      AttributeGrid distGrid2 = makeDistanceGrid(bounds, -w, w, distanceBitCount, vs);
      ClosestPointIndexer.getPointsInWorldUnits(indexGrid1, pntx, pnty, pntz);

      ClosestPointIndexer.makeDistanceGrid(indexGrid1, pntx, pnty, pntz, null, -w, w, distGrid1);
      ClosestPointIndexer.makeDistanceGrid(indexGrid2, pntx, pnty, pntz, null, -w, w, distGrid2);
      long diffDist = GridUtil.compareGrids(distGrid1, distGrid2);
      printf("distance difference count: %d\n", diffDist);
    }
  }
コード例 #2
0
  void makeTestPoint() throws Exception {

    if (DEBUG) printf("%s.makeTestPoint()\n", this.getClass().getName());
    double vs = 0.05 * MM;
    double w = 10 * MM; // half width
    double firstLayerThickness = 2.5; // 1.5;//3.5;//2.4;//1.8;//0.8
    double xmin = -w, xmax = w, ymin = -w, ymax = w, zmin = -w, zmax = w; // 10*vs;
    int subvoxelResolution = 100;
    double INF = 1.e5; // point to ignore
    int voxelSquareSize = 25; // for visualization
    boolean snapToGrid = false;
    int iterationsCount = 0;
    int threadCount = 1;

    ArrayAttributeGridInt indexGrid =
        new ArrayAttributeGridInt(new Bounds(xmin, xmax, ymin, ymax, zmin, zmax), vs, vs);
    double pnts[] =
        new double[] {
          0,
          0,
          0,
          // 1*vs, 0*vs, 0.5*vs,  0.2*vs, -2.3*vs,  0.5*vs,
          // 1*vs, 0*vs, 0.5*vs,  0.2*vs, -2.4*vs,  0.5*vs,
          // 1*vs, 0*vs, 0.5*vs,  0.5*vs, -2.4*vs,  0.5*vs, / bad
          // 7*vs, 0*vs, 0.5*vs,  6.5*vs, -2.4*vs,  0.5*vs, //6.5*vs, 2.4*vs,  0.5*vs,

          7 * vs,
          0 * vs,
          0.5 * vs,
          6.5 * vs,
          2.4 * vs,
          0.5 * vs,
          6.5 * vs,
          -2.4 * vs,
          0.5 * vs,
          -7 * vs,
          0 * vs,
          0.5 * vs,
          -6.5 * vs,
          2.4 * vs,
          0.5 * vs,
          -6.5 * vs,
          -2.4 * vs,
          0.5 * vs,
          0,
          7 * vs,
          0.5 * vs,
          2.4 * vs,
          6.5 * vs,
          0.5 * vs,
          -2.4 * vs,
          6.5 * vs,
          0.5 * vs,
          0,
          -7 * vs,
          0.5 * vs,
          2.4 * vs,
          -6.5 * vs,
          0.5 * vs,
          -2.4 * vs,
          -6.5 * vs,
          0.5 * vs,
        };
    // pnts = makeCircle(0.2*vs, 0.6*vs, 0, 7*vs, 32);
    // pnts = makeCircleZ(0.2*vs, 0.2*vs, 0.2*vs, 7*vs, 128);
    // pnts = makeCircleZ(0.5*vs, 0.5*vs, 0.5*vs, 7*vs, 64);
    // pnts = makeCircleZ(0.5*vs, 0.5*vs, 0.5*vs, 0.8*w, 64);
    // pnts = makeCircleZ(0.5*vs, 0.5*vs, 0.5*vs, 0.7*w, 80);
    // pnts = makeCircleZ(0.5*vs, 0.5*vs, 0.5*vs, 0.7*w, 256);
    // pnts = makeCircleX(0.5*vs, 0.5*vs, 0.5*vs, 7*vs, 32);
    // pnts = makeCircleY(0.5*vs, 0.5*vs, 0.5*vs, 7*vs, 32);
    // pnts = makeUnion(makeUnion(makeCircleX(0.5*vs, 0.5*vs, 0.5*vs, 7*vs, 32),makeCircleY(0.5*vs,
    // 0.5*vs, 0.5*vs, 7*vs, 32)),makeCircleZ(0.5*vs, 0.5*vs, 0.5*vs, 7*vs, 32));
    pnts =
        makeUnion(
            makeUnion(
                makeCircleX(0.5 * vs, 0.5 * vs, 0.5 * vs, 0.9 * w, 128),
                makeCircleY(0.5 * vs, 0.5 * vs, 0.5 * vs, 0.9 * w, 128)),
            makeCircleZ(0.5 * vs, 0.5 * vs, 0.5 * vs, 0.9 * w, 128));
    // pnts = makeUnion(makeUnion(makeCircleX(0.3*vs, 0.3*vs, 0.3*vs, 7*vs, 32),makeCircleY(0.3*vs,
    // 0.3*vs, 0.3*vs, 7*vs, 32)),makeCircleZ(0.3*vs, 0.3*vs, 0.3*vs, 7*vs, 32));
    // pnts = makeUnion(makeUnion(makeCircleX(0.3*vs, 0.3*vs, 0.3*vs, 7*vs, 64),makeCircleY(0.3*vs,
    // 0.3*vs, 0.3*vs, 7*vs, 64)),makeCircleZ(0.3*vs, 0.3*vs, 0.3*vs, 7*vs, 64));
    // pnts = makeUnion(makeUnion(makeCircleX(0.5*vs, 0.5*vs, 0.5*vs, 7*vs, 64),makeCircleY(0.5*vs,
    // 0.5*vs, 0.5*vs, 7*vs, 64)),makeCircleZ(0.5*vs, 0.5*vs, 0.5*vs, 7*vs, 64));
    // pnts = makeUnion(makeUnion(makeCircleX(0.5*vs, 0.5*vs, 0.5*vs, 7*vs, 128),makeCircleY(0.5*vs,
    // 0.5*vs, 0.5*vs, 7*vs, 128)),makeCircleZ(0.5*vs, 0.5*vs, 0.5*vs, 7*vs, 128));
    // pnts = makeUnion(makeCircle(0.5*vs, 0.5*vs, 0.5*vs, 7*vs, 32),makeCircle(2.5*vs, 0.5*vs,
    // 3.5*vs, 7*vs, 32));
    // pnts = makeUnion(makeCircle(-9.5*vs, 0.5*vs, 0, 7*vs, 32),makeCircle(9.5*vs, 0.5*vs, 0, 7*vs,
    // 32));
    // pnts = makeUnion(makeUnion(makeCircle(-10.5*vs, 0.5*vs, 0.5*vs, 7*vs, 32),makeCircle(11.5*vs,
    // 0.5*vs, 0.5*vs, 7*vs, 32)),
    //                 makeUnion(makeCircle(0.5*vs, -10.5*vs, 0.5*vs, 7*vs, 32),makeCircle(0.5*vs,
    // 11.5*vs, 0.5*vs, 7*vxs, 32)));
    // pnts = makeRandomPoints();
    int pcount = pnts.length / 3;
    // printf("pcount: %d\n", (pcount-1));
    double pntx[] = new double[pcount];
    double pnty[] = new double[pcount];
    double pntz[] = new double[pcount];

    ClosestPointIndexer.getPointsInGridUnits(indexGrid, pnts, pntx, pnty, pntz);
    if (snapToGrid) {
      ClosestPointIndexer.snapToVoxels(pntx);
      ClosestPointIndexer.snapToVoxels(pnty);
      ClosestPointIndexer.snapToVoxels(pntz);
    }

    ClosestPointIndexer.initFirstLayer(indexGrid, pntx, pnty, pntz, firstLayerThickness);
    int usedCount = ClosestPointIndexer.removeUnusedPoints(indexGrid, pntx, pnty, pntz);
    printf(
        "grid: [%d x %d x %d] threads: %d points: %d\n",
        indexGrid.getWidth(), indexGrid.getHeight(), indexGrid.getDepth(), threadCount, pcount);

    int nz = indexGrid.getDepth();

    // if(true) printIndices(indexGrid);
    // for(int z = 0; z < nz; z++)renderDiff(indexGrid, z, pntx, pnty, pntz, voxelSquareSize,
    // fmt("/tmp/dist/distDiff1_%02d.png",z), true);

    // distribute distances to the whole grid
    long t0 = time();
    ClosestPointIndexerMT.PI3_MT(pntx, pnty, pntz, indexGrid, threadCount);
    // ClosestPointIndexer.PI3(pntx, pnty, pntz, indexGrid);
    printf("ClosestPointIndexer time: %d\n ", (time() - t0));
    // printIndices(indexGrid);
    // printDiff(indexGrid, pntx, pnty, pntz, true);
    // printDiff(indexGrid, pntx, pnty, pntz, false);
    // printf("maxDiff: %10.3e\n", getMaxDiff(indexGrid, pntx, pnty, pntz));
    for (int z = 0; z < nz; z++) {
      // renderDiff(indexGrid, z, pntx, pnty, pntz, voxelSquareSize,
      // fmt("/tmp/dist/distDiff2_%02d.png",z), true);
    }
  }