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

Comparing ray/src/cv/bsdf2rado.c (file contents):
Revision 2.1 by greg, Thu Jun 7 15:39:36 2018 UTC vs.
Revision 2.4 by greg, Thu Sep 9 15:17:45 2021 UTC

# Line 75 | Line 75 | main(int argc, char *argv[])
75                          bsdf_min = myBSDF.rLambFront.cieY/M_PI;
76                  if (myBSDF.rb != NULL && myBSDF.rLambBack.cieY < bsdf_min*M_PI)
77                          bsdf_min = myBSDF.rLambBack.cieY/M_PI;
78 <                if ((myBSDF.tf != NULL) | (myBSDF.tb != NULL) &&
79 <                                myBSDF.tLamb.cieY < bsdf_min*M_PI)
80 <                        bsdf_min = myBSDF.tLamb.cieY/M_PI;
78 >                if (myBSDF.tf != NULL && myBSDF.tLambFront.cieY < bsdf_min*M_PI)
79 >                        bsdf_min = myBSDF.tLambFront.cieY/M_PI;
80 >                if (myBSDF.tb != NULL && myBSDF.tLambBack.cieY < bsdf_min*M_PI)
81 >                        bsdf_min = myBSDF.tLambBack.cieY/M_PI;
82                  if (doTrans && (myBSDF.tf == NULL) & (myBSDF.tb == NULL)) {
83                          fprintf(stderr, "%s: no transmitted component in '%s'\n",
84                                          progname, argv[1]);
# Line 95 | Line 96 | main(int argc, char *argv[])
96                  fclose(fp);
97          }
98   #ifdef DEBUG
99 <        fprintf(stderr, "Minimum BSDF set to %.4f\n", bsdf_min);
99 >        fprintf(stderr, "Minimum BSDF set to %.4g\n", bsdf_min);
100   #endif
101 <        min_log = log(bsdf_min*.5 + 1e-5);
101 >        min_log = log(bsdf_min*.5 + 1e-6);
102                                                  /* output BSDF rep. */
103          for (n = 0; (n < 6) & (2*n+3 < argc); n++) {
104                  double  theta = (M_PI/180.)*atof(argv[2*n+2]);
# Line 116 | Line 117 | main(int argc, char *argv[])
117                  idir[1] = idir[2] * sin(phi);
118                  idir[2] = input_orient * sqrt(1. - idir[2]*idir[2]);
119   #ifdef DEBUG
120 <                fprintf(stderr, "Computing BSDF for incident direction (%.1f,%.1f)\n",
120 >                fprintf(stderr, "Computing BSDF for incident direction (%.2f,%.2f)\n",
121                                  get_theta180(idir), get_phi360(idir));
122   #endif
123                  if (!inpXML)
124                          rbf = advect_rbf(idir, 15000);
125   #ifdef DEBUG
126                  if (inpXML)
127 <                        fprintf(stderr, "Hemispherical %s: %.3f\n",
127 >                        fprintf(stderr, "Hemispherical %s: %.4f\n",
128                                  (output_orient > 0 ^ input_orient > 0 ?
129                                          "transmission" : "reflection"),
130                                  SDdirectHemi(idir, SDsampSp|SDsampDf |
# Line 132 | Line 133 | main(int argc, char *argv[])
133                  else if (rbf == NULL)
134                          fputs("Empty RBF\n", stderr);
135                  else
136 <                        fprintf(stderr, "Hemispherical %s: %.3f\n",
136 >                        fprintf(stderr, "Hemispherical %s: %.4f\n",
137                                  (output_orient > 0 ^ input_orient > 0 ?
138                                          "transmission" : "reflection"),
139                                  rbf->vtotal);
# Line 155 | Line 156 | main(int argc, char *argv[])
156                  }
157                  fflush(stdout);
158                  sprintf(buf, "gensurf tmat bsdf%d - - - %d %d %s", n+1,
159 <                                                GRIDRES-1, GRIDRES-1, validf);
159 >                                                grid_res-1, grid_res-1, validf);
160                  fp = popen(buf, "w");
161                  if (fp == NULL) {
162                          fprintf(stderr, "%s: cannot open '| %s'\n", progname, buf);
163                          return(1);
164                  }
165 <                for (i = 0; i < GRIDRES; i++)
166 <                    for (j = 0; j < GRIDRES; j++) {
165 >                for (i = 0; i < grid_res; i++)
166 >                    for (j = 0; j < grid_res; j++) {
167                          ovec_from_pos(odir, i, j);
168                          if (inpXML) {
169                                  SDValue sval;
# Line 172 | Line 173 | main(int argc, char *argv[])
173                                  bsdf = sval.cieY;
174                          } else
175                                  bsdf = eval_rbfrep(rbf, odir);
176 <                        bsdf = log(bsdf + 1e-5) - min_log;
176 >                        bsdf = log(bsdf + 1e-6) - min_log;
177                          fprintf(fp, "%.8e %.8e %.8e\n",
178                                          odir[0]*bsdf, odir[1]*bsdf, odir[2]*bsdf);
179                      }

Diff Legend

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