SCM Repository

[diderot] Diff of /benchmarks/programs/vr-lite-cam/bmark-teem.c

[diderot] / benchmarks / programs / vr-lite-cam / bmark-teem.c

Diff of /benchmarks/programs/vr-lite-cam/bmark-teem.c

-revision 1574, Sat Oct 29 14:49:18 2011 UTC
+revision 1575, Sat Oct 29 16:50:18 2011 UTC
 Line 143
    Nrrd *nin;
    int status;
    nin = nrrdNew();
-   if(nin == NULL)
+   if(nin == NULL) {
-   {
      err = biffGetDone(NRRD);
      fprintf(stderr, "Trouble allocating nrrd struct:\n%s", err);
      free(err);
-Line 153
+Line 152
    status = nrrdLoad(nin, infile, NULL);
-   if (status)
+   if (status) {
-   {
      err = biffGetDone(NRRD);
      fprintf(stderr, "Trouble reading \"%s\":\n%s", infile, err);
      free(err);
-Line 170
+Line 168
    double kparm[NRRD_KERNEL_PARMS_NUM];
    ctx = gageContextNew();
-   if(ctx == NULL)
+   if(ctx == NULL) {
-   {
      err = biffGetDone(GAGE);
      fprintf(stderr, "Trouble allocating gage context:\n%s", err);
      free(err);
-Line 180
+Line 177
    }
    pvl = gagePerVolumeNew(ctx, nin, gageKindScl);
-   if(pvl == NULL)
+   if(pvl == NULL) {
-   {
      err = biffGetDone(GAGE);
      fprintf(stderr, "Trouble allocating gage PerVolume:\n%s", err);
      free(err);
-Line 191
+Line 187
    }
    status = gagePerVolumeAttach(ctx, pvl);
-   if(status)
+   if(status) {
-   {
      err = biffGetDone(GAGE);
      fprintf(stderr, "Trouble attaching gage PerVolume:\n%s", err);
      free(err);
-Line 203
+Line 198
    // It's okay to use kparm unitialized here: gageKernelSet ignores it.
    status = gageKernelSet(ctx, gageKernel00, nrrdKernelBSpline3, kparm);
-   if(status)
+   if(status) {
-   {
      err = biffGetDone(GAGE);
      fprintf(stderr, "Trouble setting kernel:\n%s", err);
      free(err);
-Line 215
+Line 209
    // It's okay to use kparm unitialized here: gageKernelSet ignores it.
    status = gageKernelSet(ctx, gageKernel11, nrrdKernelBSpline3D, kparm);
-   if(status)
+   if(status) {
-   {
      err = biffGetDone(GAGE);
      fprintf(stderr, "Trouble setting kernel:\n%s", err);
      free(err);
-Line 226
+Line 219
    }
    status = gageQueryItemOn(ctx, pvl, gageSclValue);
-   if(status)
+   if(status) {
-   {
      err = biffGetDone(GAGE);
      fprintf(stderr, "Trouble with gageQueryItemOn:\n%s", err);
      free(err);
-Line 237
+Line 229
    }
    status = gageQueryItemOn(ctx, pvl, gageSclGradVec);
-   if(status)
+   if(status) {
-   {
      err = biffGetDone(GAGE);
      fprintf(stderr, "Trouble with gageQueryItemOn:\n%s", err);
      free(err);
-Line 248
+Line 239
    }
    status = gageUpdate(ctx);
-   if(status)
+   if(status) {
-   {
      err = biffGetDone(GAGE);
      fprintf(stderr, "Trouble with gageUpdate:\n%s", err);
      free(err);
-Line 259
+Line 249
    }
    val = gageAnswerPointer(ctx, pvl, gageSclValue);
-   if(val == NULL)
+   if(val == NULL) {
-   {
      err = biffGetDone(GAGE);
      fprintf(stderr, "Trouble getting answer pointer:\n%s", err);
      free(err);
-Line 270
+Line 259
    }
    grad = gageAnswerPointer(ctx, pvl, gageSclGradVec);
-   if(grad == NULL)
+   if(grad == NULL) {
-   {
      err = biffGetDone(GAGE);
      fprintf(stderr, "Trouble getting answer pointer:\n%s", err);
      free(err);
-Line 280
+Line 268
      return -1;
    }
    int ui, vi;
    double rayU, rayV, rayN;
    double rayVec[3];
-Line 297
+Line 284
    double max_gray = 0;
    double min_gray = 0;
    double *out_data = malloc(sizeof(double) * 4 * imgResU * imgResV);
-   if(out_data == NULL)
+   if(out_data == NULL) {
-   {
      fprintf(stderr, "Trouble with malloc.\n");
      gageContextNix(ctx);
      nrrdNuke(nin);
      return -1;
    }
-   double t0 = GetTime();
+   double t0 = GetTime(); // start timing
-   for(vi = 0; vi < imgResV; vi++)
+   for(vi = 0; vi < imgResV; vi++) {
-   {
      rayV = lerp(-camVmax, camVmax, -0.5, (double)(vi), (double)(imgResV) - 0.5);
      scale_vec3(rayV, camV, temp);
-     for(ui = 0; ui < imgResU; ui++)
+     for(ui = 0; ui < imgResU; ui++) {
-     {
        rayU = lerp(-camUmax, camUmax, -0.5, (double)(ui), (double)(imgResU) - 0.5);
        scale_vec3(rayU, camU, temp2);
        add3(temp, temp2, temp2);
-Line 325
+Line 309
        output[1] = 0;
        output[2] = 0;
        output[3] = 0;
-       for(rayN = camVspNear; rayN <= camVspFar; rayN += rayStep)
+       for(rayN = camVspNear; rayN <= camVspFar; rayN += rayStep) {
-       {
          scale_vec3(rayN, rayDir, temp2);
          add3(temp2, camEye, rayPos);
-         if(inside(rayPos))
+         if(inside(rayPos)) {
-         {
            status = gageProbe(ctx, rayPos[0], rayPos[1], rayPos[2]);
-           if(status)
+           if(status) {
-           {
              err = ctx->errStr;
              fprintf(stderr, "Trouble with gageProbe:\n%s", err);
              free(err);
-Line 342
+Line 323
              nrrdNuke(nin);
              return -1;
            }
-           if(*val > valOpacMin)
+           if(*val > valOpacMin) {
-           {
              temp2[0] = *(grad + 0);
              temp2[1] = *(grad + 1);
              temp2[2] = *(grad + 2);
-Line 366
+Line 346
              transp *= 1 - alpha;
            }
          }
-         if(transp < 0.01)
+         if(transp < 0.01) {
-         {
            transp = 0;
            break;
          }
-Line 376
+Line 355
        *(out_data + ui * imgResV * 4 + vi * 4 + 1) = gray;
        *(out_data + ui * imgResV * 4 + vi * 4 + 2) = gray;
        *(out_data + ui * imgResV * 4 + vi * 4 + 3) = 1.0f - transp;
-       if(!set)
+       if(!set) {
-       {
          max_gray = gray;
          min_gray = gray;
          set = true;
-Line 388
+Line 366
    }
    unsigned char *uc_out_data = malloc(sizeof(unsigned char) * imgResU * imgResV);
-   for(vi = 0; vi < imgResV; vi++)
+   for(vi = 0; vi < imgResV; vi++) {
-   {
+     for(ui = 0; ui < imgResU; ui++) {
-     for(ui = 0; ui < imgResU; ui++)
-     {
        *(uc_out_data + vi * imgResU + ui) = ulerp(min_gray, *(out_data + ui * imgResV * 4 + vi * 4 + 0), max_gray);
      }
    }
-   double totalTime = GetTime() - t0;
+   double totalTime = GetTime() - t0;
    printf("usr=%f\n", totalTime);
    // Write out data
    Nrrd *nout;
    nout = nrrdNew();
-   if(nout == NULL)
+   if(nout == NULL) {
-   {
      err = biffGetDone(NRRD);
      fprintf(stderr, "Trouble allocating nrrd struct:\n%s", err);
      free(err);
-Line 415
+Line 390
    }
    status = nrrdWrap_va(nout, uc_out_data, nrrdTypeUChar, 2, imgResU, imgResV);
-   if(status)
+   if(status) {
-   {
      err = biffGetDone(NRRD);
      fprintf(stderr, "Trouble wrapping nrrd struct:\n%s", err);
      free(err);
-Line 428
+Line 402
    }
    status = nrrdSave(outfile, nout, NULL);
-   if(status)
+   if(status) {
-   {
      err = biffGetDone(NRRD);
      fprintf(stderr, "Trouble saving nrrd struct:\n%s", err);
      free(err);

 Legend:



Removed from v.1574
 


changed lines


 
Added in v.1575
 Legend:



Removed from v.1574
 


changed lines


 
Added in v.1575
-Removed from v.1574
+Added in v.1575

root@smlnj-gforge.cs.uchicago.edu	ViewVC Help
Powered by ViewVC 1.0.0