ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cv/bsdf2rad.c
(Generate patch)

Comparing ray/src/cv/bsdf2rad.c (file contents):
Revision 2.24 by greg, Tue Apr 11 18:26:55 2017 UTC vs.
Revision 2.25 by greg, Wed Apr 12 00:31:40 2017 UTC

# Line 270 | Line 270 | build_wRBF(void)
270  
271   /* Put out mirror arrow for the given incident vector */
272   static void
273 < put_mirror_arrow(const FVECT ivec, int inc_side)
273 > put_mirror_arrow(const FVECT origin, const FVECT nrm)
274   {
275          const double    arrow_len = 1.2*bsdf_rad;
276          const double    tip_len = 0.2*bsdf_rad;
277 <        FVECT           origin, refl;
277 >        FVECT           refl;
278          int             i;
279  
280 <        cvt_sposition(origin, ivec, inc_side);
280 >        refl[0] = 2.*nrm[2]*nrm[0];
281 >        refl[1] = 2.*nrm[2]*nrm[1];
282 >        refl[2] = 2.*nrm[2]*nrm[2] - 1.;
283  
282 –        refl[0] = -2.*ivec[2]*ivec[0];
283 –        refl[1] = -2.*ivec[2]*ivec[1];
284 –        refl[2] = 2.*ivec[2]*ivec[2] - 1.;
285 –
284          printf("\n# Mirror arrow\n");
285          printf("\narrow_mat cylinder inc_dir\n0\n0\n7");
286          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n",
# Line 309 | Line 307 | put_mirror_arrow(const FVECT ivec, int inc_side)
307  
308   /* Put out transmitted direction arrow for the given incident vector */
309   static void
310 < put_trans_arrow(const FVECT ivec, int inc_side)
310 > put_trans_arrow(const FVECT origin)
311   {
312          const double    arrow_len = 1.2*bsdf_rad;
313          const double    tip_len = 0.2*bsdf_rad;
316 –        FVECT           origin;
314          int             i;
315  
319 –        cvt_sposition(origin, ivec, inc_side);
320 –
316          printf("\n# Transmission arrow\n");
317          printf("\narrow_mat cylinder trans_dir\n0\n0\n7");
318          printf("\n\t%f %f %f\n\t%f %f %f\n\t%f\n",
# Line 372 | Line 367 | static int
367   put_BSDFs(void)
368   {
369          const double    scalef = bsdf_rad/(log10(overall_max) - min_log10);
370 <        FVECT           ivec, sorg, upv;
370 >        FVECT           ivec, sorg, nrm, upv;
371          RREAL           vMtx[3][3];
372          char            *fname;
373          char            cmdbuf[256];
374 <        char            xfargs[128];
375 <        int             nxfa;
374 >        char            rotargs[64];
375 >        int             nrota;
376          int             i;
377  
378          printf("\n# Gensurf output corresponding to %d incident directions\n",
# Line 387 | Line 382 | put_BSDFs(void)
382          printf("\nvoid mixfunc arrow_mat\n4 arrow_glow void 0.25 .\n0\n0\n");
383  
384          if (front_comp & SDsampR)                       /* front reflection */
385 <                for (i = 0; i < NINCIDENT; i++) {
386 <                        get_ivector(ivec, i);
387 <                        put_mirror_arrow(ivec, 1);
385 >        for (i = 0; i < NINCIDENT; i++) {
386 >                get_ivector(ivec, i);
387 >                nrm[0] = -ivec[0]; nrm[1] = -ivec[1]; nrm[2] = ivec[2];
388 >                upv[0] = nrm[0]*nrm[1]*(nrm[2] - 1.);
389 >                upv[1] = nrm[0]*nrm[0] + nrm[1]*nrm[1]*nrm[2];
390 >                upv[2] = -nrm[1]*(nrm[0]*nrm[0] + nrm[1]*nrm[1]);
391 >                if (SDcompXform(vMtx, nrm, upv) != SDEnone)
392 >                        continue;
393 >                nrota = addrot(rotargs, vMtx[0], vMtx[1], vMtx[2]);
394 >                if (front_comp) {
395                          cvt_sposition(sorg, ivec, 1);
394 –                        ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */
395 –                        upv[0] = ivec[0]*ivec[1]*(ivec[2] - 1.);
396 –                        upv[1] = ivec[0]*ivec[0] + ivec[1]*ivec[1]*ivec[2];
397 –                        upv[2] = -ivec[1]*(ivec[0]*ivec[0] + ivec[1]*ivec[1]);
398 –                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
399 –                                        sorg[0], sorg[1], sorg[2]);
400 –                        nxfa = 6;
396                          printf("\nvoid colorfunc scale_pat\n");
397 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
398 <                                        4+nxfa, xfargs);
397 >                        printf("10 bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n");
398 >                        printf("\t-s %f -t %f %f %f\n0\n0\n",
399 >                                        bsdf_rad, sorg[0], sorg[1], sorg[2]);
400                          printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
401 <                        if (SDcompXform(vMtx, ivec, upv) != SDEnone)
402 <                                continue;
403 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
408 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
409 <                                        scalef, sorg[0], sorg[1], sorg[2]);
410 <                        nxfa += 6;
401 >                }
402 >                if (front_comp & SDsampR) {
403 >                        put_mirror_arrow(sorg, nrm);
404                          fname = tfile_name(frpref, dsuffix, i);
405 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform %s",
405 >                        sprintf(cmdbuf,
406 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform %s -s %f -t %f %f %f",
407                                          frpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
408 <                                        xfargs);
408 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
409                          if (!run_cmd(cmdbuf))
410                                  return(0);
411                  }
412 <        if (front_comp & SDsampT)                       /* front transmission */
413 <                for (i = 0; i < NINCIDENT; i++) {
420 <                        get_ivector(ivec, i);
421 <                        put_trans_arrow(ivec, 1);
422 <                        cvt_sposition(sorg, ivec, 1);
423 <                        ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */
424 <                        upv[0] = ivec[0]*ivec[1]*(ivec[2] - 1.);
425 <                        upv[1] = ivec[0]*ivec[0] + ivec[1]*ivec[1]*ivec[2];
426 <                        upv[2] = -ivec[1]*(ivec[0]*ivec[0] + ivec[1]*ivec[1]);
427 <                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
428 <                                        sorg[0], sorg[1], sorg[2]);
429 <                        nxfa = 6;
430 <                        printf("\nvoid colorfunc scale_pat\n");
431 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
432 <                                        4+nxfa, xfargs);
433 <                        printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
434 <                        if (SDcompXform(vMtx, ivec, upv) != SDEnone)
435 <                                continue;
436 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
437 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
438 <                                        scalef, sorg[0], sorg[1], sorg[2]);
439 <                        nxfa += 6;
412 >                if (front_comp & SDsampT) {
413 >                        put_trans_arrow(sorg);
414                          fname = tfile_name(ftpref, dsuffix, i);
415 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I %s",
415 >                        sprintf(cmdbuf,
416 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform -I %s -s %f -t %f %f %f",
417                                          ftpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
418 <                                        xfargs);
418 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
419                          if (!run_cmd(cmdbuf))
420                                  return(0);
421                  }
422 <        if (back_comp & SDsampR)                        /* rear reflection */
448 <                for (i = 0; i < NINCIDENT; i++) {
449 <                        get_ivector(ivec, i);
450 <                        put_mirror_arrow(ivec, -1);
422 >                if (back_comp) {
423                          cvt_sposition(sorg, ivec, -1);
452 –                        ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */
453 –                        upv[0] = ivec[0]*ivec[1]*(ivec[2] - 1.);
454 –                        upv[1] = ivec[0]*ivec[0] + ivec[1]*ivec[1]*ivec[2];
455 –                        upv[2] = -ivec[1]*(ivec[0]*ivec[0] + ivec[1]*ivec[1]);
456 –                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
457 –                                        sorg[0], sorg[1], sorg[2]);
458 –                        nxfa = 6;
424                          printf("\nvoid colorfunc scale_pat\n");
425 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
426 <                                        4+nxfa, xfargs);
425 >                        printf("10 bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n");
426 >                        printf("\t-s %f -t %f %f %f\n0\n0\n",
427 >                                        bsdf_rad, sorg[0], sorg[1], sorg[2]);
428                          printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
429 <                        if (SDcompXform(vMtx, ivec, upv) != SDEnone)
430 <                                continue;
431 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
466 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
467 <                                        scalef, sorg[0], sorg[1], sorg[2]);
468 <                        nxfa += 6;
429 >                }
430 >                if (back_comp & SDsampR) {
431 >                        put_mirror_arrow(sorg, nrm);
432                          fname = tfile_name(brpref, dsuffix, i);
433 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -ry 180 %s",
433 >                        sprintf(cmdbuf,
434 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform -I -ry 180 %s -s %f -t %f %f %f",
435                                          brpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
436 <                                        xfargs);
436 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
437                          if (!run_cmd(cmdbuf))
438                                  return(0);
439                  }
440 <        if (back_comp & SDsampT)                        /* rear transmission */
441 <                for (i = 0; i < NINCIDENT; i++) {
478 <                        get_ivector(ivec, i);
479 <                        put_trans_arrow(ivec, -1);
480 <                        cvt_sposition(sorg, ivec, -1);
481 <                        ivec[0] = -ivec[0]; ivec[1] = -ivec[1]; /* normal */
482 <                        upv[0] = ivec[0]*ivec[1]*(ivec[2] - 1.);
483 <                        upv[1] = ivec[0]*ivec[0] + ivec[1]*ivec[1]*ivec[2];
484 <                        upv[2] = -ivec[1]*(ivec[0]*ivec[0] + ivec[1]*ivec[1]);
485 <                        sprintf(xfargs, "-s %f -t %f %f %f", bsdf_rad,
486 <                                        sorg[0], sorg[1], sorg[2]);
487 <                        nxfa = 6;
488 <                        printf("\nvoid colorfunc scale_pat\n");
489 <                        printf("%d bsdf_red bsdf_grn bsdf_blu bsdf2rad.cal\n\t%s\n0\n0\n",
490 <                                        4+nxfa, xfargs);
491 <                        printf("\nscale_pat glow scale_mat\n0\n0\n4 1 1 1 0\n");
492 <                        if (SDcompXform(vMtx, ivec, upv) != SDEnone)
493 <                                continue;
494 <                        nxfa = addrot(xfargs, vMtx[0], vMtx[1], vMtx[2]);
495 <                        sprintf(xfargs+strlen(xfargs), " -s %f -t %f %f %f",
496 <                                        scalef, sorg[0], sorg[1], sorg[2]);
497 <                        nxfa += 6;
440 >                if (back_comp & SDsampT) {
441 >                        put_trans_arrow(sorg);
442                          fname = tfile_name(btpref, dsuffix, i);
443 <                        sprintf(cmdbuf, "gensurf scale_mat %s%d %s %s %s %d %d | xform -ry 180 %s",
443 >                        sprintf(cmdbuf,
444 >                "gensurf scale_mat %s%d %s %s %s %d %d | xform -ry 180 %s -s %f -t %f %f %f",
445                                          btpref, i, fname, fname, fname, SAMPRES-1, SAMPRES-1,
446 <                                        xfargs);
446 >                                        rotargs, scalef, sorg[0], sorg[1], sorg[2]);
447                          if (!run_cmd(cmdbuf))
448                                  return(0);
449                  }
450 +        }
451          return(1);
452   }
453  

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines