ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/util/rmatrix.h
Revision: 2.14
Committed: Fri Mar 4 01:27:12 2022 UTC (2 years ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.13: +13 -5 lines
Log Message:
perf(rmtxop): Added memory-mapping for double matrices

File Contents

# User Rev Content
1 greg 2.14 /* RCSid $Id: rmatrix.h,v 2.13 2021/01/19 23:32:00 greg Exp $ */
2 greg 2.1 /*
3     * Header file for general matrix routines.
4     */
5    
6     #ifndef _RAD_RMATRIX_H_
7     #define _RAD_RMATRIX_H_
8    
9     #include "cmatrix.h"
10    
11     #ifdef __cplusplus
12     extern "C" {
13     #endif
14    
15 greg 2.13 /* Preferred BSDF component:
16     transmission, reflection front (normal side), reflection back */
17     typedef enum {RMPtrans=0, RMPreflF, RMPreflB} RMPref;
18    
19 greg 2.14 /* General [row][col][cmp] component matrix */
20 greg 2.1 typedef struct {
21 greg 2.14 char *info;
22     void *mapped;
23     double *mtx;
24     COLOR cexp;
25 greg 2.11 int nrows, ncols;
26     short ncomp;
27     uby8 dtype;
28     uby8 swapin;
29 greg 2.1 } RMATRIX;
30    
31 greg 2.9 #define rmx_lval(rm,r,c,i) (rm)->mtx[(i)+(rm)->ncomp*((c)+(size_t)(rm)->ncols*(r))]
32 greg 2.1
33 greg 2.14 /* Initialize a RMATRIX struct but don't allocate array space */
34     extern RMATRIX *rmx_new(int nr, int nc, int n);
35    
36     /* Prepare a RMATRIX for writing (allocate array if needed) */
37     extern int rmx_prepare(RMATRIX *rm);
38    
39     /* Call rmx_new() and rmx_prepare() */
40 greg 2.1 extern RMATRIX *rmx_alloc(int nr, int nc, int n);
41    
42 greg 2.4 /* Free a RMATRIX array */
43     extern void rmx_free(RMATRIX *rm);
44    
45     /* Resolve data type based on two input types (returns 0 for mismatch) */
46     extern int rmx_newtype(int dtyp1, int dtyp2);
47    
48 greg 2.7 /* Load matrix from supported file type (NULL for stdin, '!' with command) */
49 greg 2.13 extern RMATRIX *rmx_load(const char *inspec, RMPref rmp);
50 greg 2.1
51 greg 2.3 /* Append header information associated with matrix data */
52     extern int rmx_addinfo(RMATRIX *rm, const char *info);
53    
54 greg 2.4 /* Write matrix to file type indicated by dtype */
55 greg 2.6 extern int rmx_write(const RMATRIX *rm, int dtype, FILE *fp);
56 greg 2.1
57     /* Allocate and assign square identity matrix with n components */
58     extern RMATRIX *rmx_identity(int dim, int n);
59    
60     /* Duplicate the given matrix */
61     extern RMATRIX *rmx_copy(const RMATRIX *rm);
62    
63 greg 2.2 /* Allocate and assign transposed matrix */
64     extern RMATRIX *rmx_transpose(const RMATRIX *rm);
65 greg 2.1
66     /* Multiply (concatenate) two matrices and allocate the result */
67     extern RMATRIX *rmx_multiply(const RMATRIX *m1, const RMATRIX *m2);
68    
69 greg 2.8 /* Element-wise multiplication (or division) of m2 into m1 */
70     extern int rmx_elemult(RMATRIX *m1, const RMATRIX *m2, int divide);
71    
72 greg 2.1 /* Sum second matrix into first, applying scale factor beforehand */
73     extern int rmx_sum(RMATRIX *msum, const RMATRIX *madd, const double sf[]);
74    
75     /* Scale the given matrix by the indicated scalar component vector */
76     extern int rmx_scale(RMATRIX *rm, const double sf[]);
77    
78     /* Allocate new matrix and apply component transformation */
79     extern RMATRIX *rmx_transform(const RMATRIX *msrc, int n, const double cmat[]);
80    
81     /* Convert a color matrix to newly allocated RMATRIX buffer */
82     extern RMATRIX *rmx_from_cmatrix(const CMATRIX *cm);
83    
84     /* Convert general matrix to newly allocated CMATRIX buffer */
85     extern CMATRIX *cm_from_rmatrix(const RMATRIX *rm);
86    
87     #ifdef __cplusplus
88     }
89     #endif
90     #endif /* _RAD_RMATRIX_H_ */