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

Comparing src/util/rmtxop.c (file contents):
Revision 2.30 by greg, Fri Dec 8 00:12:31 2023 UTC vs.
Revision 2.40 by greg, Fri Apr 4 22:47:56 2025 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  
14 – #define MAXCOMP         MAXCSAMP        /* #components we support */
15 –
13   /* Unary matrix operation(s) */
14   typedef struct {
15          double          cmat[MAXCOMP*MAXCOMP];  /* component transformation */
# Line 35 | Line 32 | typedef struct {
32   int     verbose = 0;                    /* verbose reporting? */
33  
34   /* Load matrix */
35 < static int
35 > int
36   loadmatrix(ROPMAT *rop)
37   {
38          if (rop->mtx != NULL)           /* already loaded? */
# Line 46 | Line 43 | loadmatrix(ROPMAT *rop)
43          return(!rop->mtx ? -1 : 1);
44   }
45  
46 < static int      checksymbolic(ROPMAT *rop);
46 > extern int      checksymbolic(ROPMAT *rop);
47  
48   /* Check/set transform based on a reference input file */
49 < static int
49 > int
50   checkreffile(ROPMAT *rop)
51   {
52          static const char       *curRF = NULL;
# Line 108 | Line 105 | checkreffile(ROPMAT *rop)
105   }
106  
107   /* Compute conversion row from spectrum to one channel of RGB */
108 < static void
108 > void
109   rgbrow(ROPMAT *rop, int r, int p)
110   {
111          const int       nc = rop->mtx->ncomp;
# Line 125 | Line 122 | rgbrow(ROPMAT *rop, int r, int p)
122   }
123  
124   /* Compute conversion row from spectrum to one channel of XYZ */
125 < static void
125 > void
126   xyzrow(ROPMAT *rop, int r, int p)
127   {
128          const int       nc = rop->mtx->ncomp;
# Line 142 | Line 139 | xyzrow(ROPMAT *rop, int r, int p)
139   }
140  
141   /* Use the spectral sensitivity function to compute matrix coefficients */
142 < static void
143 < sensrow(ROPMAT *rop, int r, double (*sf)(SCOLOR sc, int ncs, const float wlpt[4]))
142 > void
143 > sensrow(ROPMAT *rop, int r, double (*sf)(const SCOLOR sc, int ncs, const float wlpt[4]))
144   {
145          const int       nc = rop->mtx->ncomp;
146          int             i;
# Line 157 | Line 154 | sensrow(ROPMAT *rop, int r, double (*sf)(SCOLOR sc, in
154   }
155  
156   /* Check/set symbolic transform */
157 < static int
157 > int
158   checksymbolic(ROPMAT *rop)
159   {
160          const int       nc = rop->mtx->ncomp;
161          const int       dt = rop->mtx->dtype;
162 +        double          cf = 1;
163          int             i, j;
164                                          /* check suffix => reference file */
165          if (strchr(rop->preop.csym, '.') > rop->preop.csym)
# Line 182 | Line 180 | checksymbolic(ROPMAT *rop)
180                  int     comp = 0;
181                  switch (rop->preop.csym[j]) {
182                  case 'B':
183 +                case 'b':
184                          ++comp;
185                          /* fall through */
186                  case 'G':
187 +                case 'g':
188                          ++comp;
189                          /* fall through */
190                  case 'R':
191 +                case 'r':
192 +                        if (rop->preop.csym[j] <= 'Z')
193 +                                cf = 1./WHTEFFICACY;
194                          if (dt == DTxyze) {
195                                  for (i = 3; i--; )
196 <                                        rop->preop.cmat[j*nc+i] = 1./WHTEFFICACY *
194 <                                                        xyz2rgbmat[comp][i];
196 >                                        rop->preop.cmat[j*nc+i] = cf*xyz2rgbmat[comp][i];
197                          } else if (nc == 3)
198                                  rop->preop.cmat[j*nc+comp] = 1.;
199                          else
200                                  rgbrow(rop, j, comp);
201                          break;
202                  case 'Z':
203 +                case 'z':
204                          ++comp;
205                          /* fall through */
206                  case 'Y':
207 +                case 'y':
208                          ++comp;
209                          /* fall through */
210                  case 'X':
211 +                case 'x':
212 +                        if ((rop->preop.csym[j] <= 'Z') & (dt != DTxyze))
213 +                                cf = WHTEFFICACY;
214                          if (dt == DTxyze) {
215                                  rop->preop.cmat[j*nc+comp] = 1.;
216                          } else if (nc == 3) {
# Line 215 | Line 222 | checksymbolic(ROPMAT *rop)
222                          else
223                                  xyzrow(rop, j, comp);
224  
225 <                        for (i = nc*(dt != DTxyze); i--; )
226 <                                rop->preop.cmat[j*nc+i] *= WHTEFFICACY;
225 >                        for (i = nc*(cf != 1); i--; )
226 >                                rop->preop.cmat[j*nc+i] *= cf;
227                          break;
228                  case 'S':               /* scotopic (il)luminance */
229 +                        cf = WHTSCOTOPIC;
230 +                        /* fall through */
231 +                case 's':
232                          sensrow(rop, j, scolor2scotopic);
233 <                        for (i = nc; i--; )
234 <                                rop->preop.cmat[j*nc+i] *= WHTSCOTOPIC;
233 >                        for (i = nc*(cf != 1); i--; )
234 >                                rop->preop.cmat[j*nc+i] *= cf;
235                          break;
236                  case 'M':               /* melanopic (il)luminance */
237 +                        cf = WHTMELANOPIC;
238 +                        /* fall through */
239 +                case 'm':
240                          sensrow(rop, j, scolor2melanopic);
241 <                        for (i = nc; i--; )
242 <                                rop->preop.cmat[j*nc+i] *= WHTMELANOPIC;
241 >                        for (i = nc*(cf != 1); i--; )
242 >                                rop->preop.cmat[j*nc+i] *= cf;
243                          break;
244                  case 'A':               /* average component */
245 +                case 'a':
246                          for (i = nc; i--; )
247                                  rop->preop.cmat[j*nc+i] = 1./(double)nc;
248                          break;
# Line 239 | Line 253 | checksymbolic(ROPMAT *rop)
253                  }
254          }
255                                          /* return recommended output type */
256 <        if (!strcmp(rop->preop.csym, "XYZ")) {
256 >        if (!strcasecmp(rop->preop.csym, "XYZ")) {
257                  if (dt <= DTspec)
258                          return(DTxyze);
259 <        } else if (!strcmp(rop->preop.csym, "RGB")) {
259 >        } else if (!strcasecmp(rop->preop.csym, "RGB")) {
260                  if (dt <= DTspec)
261                          return(DTrgbe);
262          } else if (dt == DTspec)
# Line 251 | Line 265 | checksymbolic(ROPMAT *rop)
265   }
266  
267   /* Get matrix and perform unary operations */
268 < static RMATRIX *
268 > RMATRIX *
269   loadop(ROPMAT *rop)
270   {
271          int     outtype = 0;
# Line 327 | Line 341 | loadop(ROPMAT *rop)
341                  }
342          }
343          if (rop->preop.transpose) {             /* transpose matrix? */
344 <                mres = rmx_transpose(rop->mtx);
331 <                if (mres == NULL) {
344 >                if (!rmx_transpose(rop->mtx)) {
345                          fputs(rop->inspec, stderr);
346                          fputs(": transpose failed\n", stderr);
347                          goto failure;
# Line 337 | Line 350 | loadop(ROPMAT *rop)
350                          fputs(rop->inspec, stderr);
351                          fputs(": transposed rows and columns\n", stderr);
352                  }
340 –                rmx_free(rop->mtx);
341 –                rop->mtx = mres;
353          }
354          mres = rop->mtx;
355          rop->mtx = NULL;
# Line 351 | Line 362 | failure:
362   }
363  
364   /* Execute binary operation, free matrix arguments and return new result */
365 < static RMATRIX *
365 > RMATRIX *
366   binaryop(const char *inspec, RMATRIX *mleft, int op, RMATRIX *mright)
367   {
368          RMATRIX *mres = NULL;
# Line 428 | Line 439 | binaryop(const char *inspec, RMATRIX *mleft, int op, R
439   }
440  
441   /* Perform matrix operations from left to right */
442 < static RMATRIX *
442 > RMATRIX *
443   op_left2right(ROPMAT *mop)
444   {
445          RMATRIX *mleft = loadop(mop);
# Line 444 | Line 455 | op_left2right(ROPMAT *mop)
455   }
456  
457   /* Perform matrix operations from right to left */
458 < static RMATRIX *
458 > RMATRIX *
459   op_right2left(ROPMAT *mop)
460   {
461          RMATRIX *mright;
# Line 473 | Line 484 | op_right2left(ROPMAT *mop)
484                                                  : (mop)->mtx->ncols)
485  
486   /* Should we prefer concatenating from rightmost matrix towards left? */
487 < static int
487 > int
488   prefer_right2left(ROPMAT *mop)
489   {
490          int     mri = 0;
# Line 500 | Line 511 | prefer_right2left(ROPMAT *mop)
511          return(t_ncols(mop+mri) < t_nrows(mop));
512   }
513  
514 < static int
514 > int
515   get_factors(double da[], int n, char *av[])
516   {
517          int     ac;
# Line 510 | Line 521 | get_factors(double da[], int n, char *av[])
521          return(ac);
522   }
523  
524 < static ROPMAT *
524 > ROPMAT *
525   resize_moparr(ROPMAT *mop, int n2alloc)
526   {
527          int     nmats = 0;
# Line 518 | Line 529 | resize_moparr(ROPMAT *mop, int n2alloc)
529  
530          while (mop[nmats++].binop)
531                  ;
532 <        for (i = nmats; i > n2alloc; i--)
532 >        for (i = nmats; i >= n2alloc; i--)
533                  rmx_free(mop[i].mtx);
534          mop = (ROPMAT *)realloc(mop, n2alloc*sizeof(ROPMAT));
535          if (mop == NULL) {
# Line 665 | Line 676 | main(int argc, char *argv[])
676          mres = loadop(mop+nmats);
677          if (mres == NULL)
678                  return(1);
679 <        if (outfmt == DTfromHeader)     /* check data type */
679 >        if ((outfmt == DTfromHeader) & (mres->dtype < DTspec))
680                  outfmt = mres->dtype;
681 <        if (outfmt == DTrgbe) {
681 >        if (outfmt == DTrgbe) {         /* check data type */
682                  if (mres->ncomp > 3)
683                          outfmt = DTspec;
684                  else if (mres->dtype == DTxyze)
685                          outfmt = DTxyze;
686          }
687 + #if DTrmx_native==DTfloat
688 +        if (outfmt == DTdouble)
689 +                fprintf(stderr,
690 +                        "%s: warning - writing float result as double\n",
691 +                                argv[0]);
692 + #endif
693          newheader("RADIANCE", stdout);  /* write result to stdout */
694          printargs(argc, argv, stdout);
695          return(rmx_write(mres, outfmt, stdout) ? 0 : 1);

Diff Legend

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