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

Comparing ray/src/util/rmatrix.c (file contents):
Revision 2.24 by greg, Tue Aug 30 15:11:22 2016 UTC vs.
Revision 2.28 by greg, Wed Jun 13 18:36:38 2018 UTC

# Line 9 | Line 9 | static const char RCSid[] = "$Id$";
9   #include <stdlib.h>
10   #include <string.h>
11   #include <fcntl.h>
12 + #include <errno.h>
13   #include "rtio.h"
14   #include "platform.h"
15   #include "resolu.h"
# Line 145 | Line 146 | rmx_load_float(RMATRIX *rm, FILE *fp)
146   static int
147   rmx_load_double(RMATRIX *rm, FILE *fp)
148   {
149 <        int     i, j, k;
149 <        double  val[100];
149 >        int     i, j;
150  
151 –        if (rm->ncomp > 100) {
152 –                fputs("Unsupported # components in rmx_load_double()\n", stderr);
153 –                exit(1);
154 –        }
151          for (i = 0; i < rm->nrows; i++)
152 <            for (j = 0; j < rm->ncols; j++) {
153 <                if (getbinary(val, sizeof(val[0]), rm->ncomp, fp) != rm->ncomp)
152 >            for (j = 0; j < rm->ncols; j++)
153 >                if (getbinary(&rmx_lval(rm,i,j,0), sizeof(double), rm->ncomp, fp) != rm->ncomp)
154                      return(0);
159 –                for (k = rm->ncomp; k--; )
160 –                     rmx_lval(rm,i,j,k) = val[k];
161 –            }
155          return(1);
156   }
157  
# Line 333 | Line 326 | rmx_write_float(const RMATRIX *rm, FILE *fp)
326   static int
327   rmx_write_double(const RMATRIX *rm, FILE *fp)
328   {
329 <        int     i, j, k;
337 <        double  val[100];
329 >        int     i, j;
330  
339 –        if (rm->ncomp > 100) {
340 –                fputs("Unsupported # components in rmx_write_double()\n", stderr);
341 –                exit(1);
342 –        }
331          for (i = 0; i < rm->nrows; i++)
332 <            for (j = 0; j < rm->ncols; j++) {
333 <                for (k = rm->ncomp; k--; )
346 <                    val[k] = rmx_lval(rm,i,j,k);
347 <                if (putbinary(val, sizeof(val[0]), rm->ncomp, fp) != rm->ncomp)
332 >            for (j = 0; j < rm->ncols; j++)
333 >                if (putbinary(&rmx_lval(rm,i,j,0), sizeof(double), rm->ncomp, fp) != rm->ncomp)
334                          return(0);
349 –            }
335          return(1);
336   }
337  
338   static int
339   rmx_write_rgbe(const RMATRIX *rm, FILE *fp)
340   {
341 <        COLOR   *scan = (COLOR *)malloc(sizeof(COLOR)*rm->ncols);
341 >        COLR    *scan = (COLR *)malloc(sizeof(COLR)*rm->ncols);
342          int     i, j;
343  
344          if (scan == NULL)
345                  return(0);
346          for (i = 0; i < rm->nrows; i++) {
347              for (j = rm->ncols; j--; )
348 <                setcolor(scan[j],       rmx_lval(rm,i,j,0),
348 >                setcolr(scan[j],        rmx_lval(rm,i,j,0),
349                                          rmx_lval(rm,i,j,1),
350                                          rmx_lval(rm,i,j,2)      );
351 <            if (fwritescan(scan, rm->ncols, fp) < 0) {
351 >            if (fwritecolrs(scan, rm->ncols, fp) < 0) {
352                  free(scan);
353                  return(0);
354              }
# Line 522 | Line 507 | rmx_multiply(const RMATRIX *m1, const RMATRIX *m2)
507          return(mres);
508   }
509  
510 + /* Element-wise multiplication (or division) of m2 into m1 */
511 + int
512 + rmx_elemult(RMATRIX *m1, const RMATRIX *m2, int divide)
513 + {
514 +        int     zeroDivides = 0;
515 +        int     i, j, k;
516 +
517 +        if ((m1 == NULL) | (m2 == NULL) ||
518 +                        (m1->ncols != m2->ncols) | (m1->nrows != m2->nrows))
519 +                return(0);
520 +        if ((m2->ncomp > 1) & (m2->ncomp != m1->ncomp))
521 +                return(0);
522 +        i = rmx_newtype(m1->dtype, m2->dtype);
523 +        if (i)
524 +                m1->dtype = i;
525 +        else
526 +                rmx_addinfo(m1, rmx_mismatch_warn);
527 +        for (i = m1->nrows; i--; )
528 +            for (j = m1->ncols; j--; )
529 +                if (divide) {
530 +                    double      d;
531 +                    if (m2->ncomp == 1) {
532 +                        d = rmx_lval(m2,i,j,0);
533 +                        if (d == 0) {
534 +                            ++zeroDivides;
535 +                            for (k = m1->ncomp; k--; )
536 +                                rmx_lval(m1,i,j,k) = 0;
537 +                        } else {
538 +                            d = 1./d;
539 +                            for (k = m1->ncomp; k--; )
540 +                                rmx_lval(m1,i,j,k) *= d;
541 +                        }
542 +                    } else
543 +                        for (k = m1->ncomp; k--; ) {
544 +                            d = rmx_lval(m2,i,j,k);
545 +                            if (d == 0) {
546 +                                ++zeroDivides;
547 +                                rmx_lval(m1,i,j,k) = 0;
548 +                            } else
549 +                                rmx_lval(m1,i,j,k) /= d;
550 +                        }
551 +                } else {
552 +                    if (m2->ncomp == 1) {
553 +                        const double    d = rmx_lval(m2,i,j,0);
554 +                        for (k = m1->ncomp; k--; )
555 +                            rmx_lval(m1,i,j,k) *= d;
556 +                    } else
557 +                        for (k = m1->ncomp; k--; )
558 +                            rmx_lval(m1,i,j,k) *= rmx_lval(m2,i,j,k);
559 +                }
560 +        if (zeroDivides) {
561 +                fputs("Divide by zero in rmx_elemult()\n", stderr);
562 +                errno = ERANGE;
563 +        }
564 +        return(1);
565 + }
566 +
567   /* Sum second matrix into first, applying scale factor beforehand */
568   int
569   rmx_sum(RMATRIX *msum, const RMATRIX *madd, const double sf[])
# Line 569 | Line 611 | rmx_scale(RMATRIX *rm, const double sf[])
611                  for (k = rm->ncomp; k--; )
612                      rmx_lval(rm,i,j,k) *= sf[k];
613  
614 +        if (rm->info)
615 +                rmx_addinfo(rm, "Applied scalar\n");
616          return(1);
617   }
618  
# Line 584 | Line 628 | rmx_transform(const RMATRIX *msrc, int n, const double
628          dnew = rmx_alloc(msrc->nrows, msrc->ncols, n);
629          if (dnew == NULL)
630                  return(NULL);
631 +        if (msrc->info) {
632 +                char    buf[128];
633 +                sprintf(buf, "Applied %dx%d matrix transform\n",
634 +                                dnew->ncomp, msrc->ncomp);
635 +                rmx_addinfo(dnew, msrc->info);
636 +                rmx_addinfo(dnew, buf);
637 +        }
638          dnew->dtype = msrc->dtype;
639          for (i = dnew->nrows; i--; )
640              for (j = dnew->ncols; j--; )

Diff Legend

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