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

Comparing ray/src/util/rmtxop.c (file contents):
Revision 2.27 by greg, Sat Dec 2 00:42:21 2023 UTC vs.
Revision 2.33 by greg, Thu May 16 18:59:19 2024 UTC

# Line 7 | Line 7 | static const char RCSid[] = "$Id$";
7  
8   #include <errno.h>
9   #include "rtio.h"
10 – #include "resolu.h"
10   #include "rmatrix.h"
11   #include "platform.h"
12  
# Line 69 | Line 68 | checkreffile(ROPMAT *rop)
68                  fclose(fp);
69                  curRF = rop->preop.csym;
70          }
71 <        if ((refm.ncomp == 3) & (refm.dtype != DTspec)) {
71 >        if (refm.ncomp == 3) {
72                  rop->preop.csym = (refm.dtype == DTxyze) ? "XYZ" : "RGB";
73                  return(checksymbolic(rop));
74          }
# Line 103 | Line 102 | checkreffile(ROPMAT *rop)
102                  for (j = refm.ncomp; j-- > 0; )
103                          rop->preop.cmat[j*nc + i] = scresp[j];
104          }
105 +        memcpy(rop->mtx->wlpart, refm.wlpart, sizeof(rop->mtx->wlpart));
106          return(0);
107   }
108  
# Line 161 | Line 161 | checksymbolic(ROPMAT *rop)
161   {
162          const int       nc = rop->mtx->ncomp;
163          const int       dt = rop->mtx->dtype;
164 +        double          cf = 1;
165          int             i, j;
166                                          /* check suffix => reference file */
167          if (strchr(rop->preop.csym, '.') > rop->preop.csym)
# Line 181 | Line 182 | checksymbolic(ROPMAT *rop)
182                  int     comp = 0;
183                  switch (rop->preop.csym[j]) {
184                  case 'B':
185 +                case 'b':
186                          ++comp;
187                          /* fall through */
188                  case 'G':
189 +                case 'g':
190                          ++comp;
191                          /* fall through */
192                  case 'R':
193 +                case 'r':
194 +                        if (rop->preop.csym[j] <= 'Z')
195 +                                cf = 1./WHTEFFICACY;
196                          if (dt == DTxyze) {
197                                  for (i = 3; i--; )
198 <                                        rop->preop.cmat[j*nc+i] = 1./WHTEFFICACY *
193 <                                                        xyz2rgbmat[comp][i];
198 >                                        rop->preop.cmat[j*nc+i] = cf*xyz2rgbmat[comp][i];
199                          } else if (nc == 3)
200                                  rop->preop.cmat[j*nc+comp] = 1.;
201                          else
202                                  rgbrow(rop, j, comp);
203                          break;
204                  case 'Z':
205 +                case 'z':
206                          ++comp;
207                          /* fall through */
208                  case 'Y':
209 +                case 'y':
210                          ++comp;
211                          /* fall through */
212                  case 'X':
213 +                case 'x':
214 +                        if ((rop->preop.csym[j] <= 'Z') & (dt != DTxyze))
215 +                                cf = WHTEFFICACY;
216                          if (dt == DTxyze) {
217                                  rop->preop.cmat[j*nc+comp] = 1.;
218                          } else if (nc == 3) {
# Line 214 | Line 224 | checksymbolic(ROPMAT *rop)
224                          else
225                                  xyzrow(rop, j, comp);
226  
227 <                        for (i = nc*(dt != DTxyze); i--; )
228 <                                rop->preop.cmat[j*nc+i] *= WHTEFFICACY;
227 >                        for (i = nc*(cf != 1); i--; )
228 >                                rop->preop.cmat[j*nc+i] *= cf;
229                          break;
230                  case 'S':               /* scotopic (il)luminance */
231 +                        cf = WHTSCOTOPIC;
232 +                        /* fall through */
233 +                case 's':
234                          sensrow(rop, j, scolor2scotopic);
235 <                        for (i = nc; i--; )
236 <                                rop->preop.cmat[j*nc+i] *= WHTSCOTOPIC;
235 >                        for (i = nc*(cf != 1); i--; )
236 >                                rop->preop.cmat[j*nc+i] *= cf;
237                          break;
238                  case 'M':               /* melanopic (il)luminance */
239 +                        cf = WHTMELANOPIC;
240 +                        /* fall through */
241 +                case 'm':
242                          sensrow(rop, j, scolor2melanopic);
243 <                        for (i = nc; i--; )
244 <                                rop->preop.cmat[j*nc+i] *= WHTMELANOPIC;
243 >                        for (i = nc*(cf != 1); i--; )
244 >                                rop->preop.cmat[j*nc+i] *= cf;
245                          break;
246                  case 'A':               /* average component */
247 +                case 'a':
248                          for (i = nc; i--; )
249                                  rop->preop.cmat[j*nc+i] = 1./(double)nc;
250                          break;
# Line 238 | Line 255 | checksymbolic(ROPMAT *rop)
255                  }
256          }
257                                          /* return recommended output type */
258 <        if (!strcmp(rop->preop.csym, "XYZ")) {
258 >        if (!strcasecmp(rop->preop.csym, "XYZ")) {
259                  if (dt <= DTspec)
260                          return(DTxyze);
261 <        } else if (!strcmp(rop->preop.csym, "RGB")) {
261 >        } else if (!strcasecmp(rop->preop.csym, "RGB")) {
262                  if (dt <= DTspec)
263                          return(DTrgbe);
264 <        }
248 <        if ((nc > 3) & (dt <= DTspec))
264 >        } else if (dt == DTspec)
265                  return(DTfloat);        /* probably not actual spectrum */
266          return(0);
267   }
# Line 274 | Line 290 | loadop(ROPMAT *rop)
290                          if (rop->preop.nsf == 1) {
291                                  for (i = rop->preop.clen; i--; )
292                                          rop->preop.cmat[i] *= rop->preop.sca[0];
293 <                        } else if (rop->preop.nsf != rop->mtx->ncomp) {
293 >                        } else if (rop->preop.nsf*rop->mtx->ncomp != rop->preop.clen) {
294                                  fprintf(stderr, "%s: -s must have one or %d factors\n",
295 <                                                rop->inspec, rop->mtx->ncomp);
295 >                                                rop->inspec,
296 >                                                rop->preop.clen/rop->mtx->ncomp);
297                                  goto failure;
298                          } else {
299 <                                for (j = rop->preop.clen/rop->preop.nsf; j--; )
300 <                                        for (i = rop->preop.nsf; i--; )
301 <                                                rop->preop.cmat[j*rop->preop.nsf+i] *=
302 <                                                                rop->preop.sca[i];
299 >                                for (i = rop->preop.nsf; i--; )
300 >                                        for (j = rop->mtx->ncomp; j--; )
301 >                                                rop->preop.cmat[i*rop->mtx->ncomp+j]
302 >                                                                *= rop->preop.sca[i];
303                          }
304                  }
305                  mres = rmx_transform(rop->mtx, rop->preop.clen/rop->mtx->ncomp,

Diff Legend

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