ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/util/rmatrix.h
Revision: 2.8
Committed: Mon Aug 28 15:59:46 2017 UTC (6 years, 7 months ago) by greg
Content type: text/plain
Branch: MAIN
CVS Tags: rad5R2
Changes since 2.7: +4 -1 lines
Log Message:
Added element-wise multiplication and division to rmtxop command

File Contents

# Content
1 /* RCSid $Id: rmatrix.h,v 2.7 2015/05/04 20:53:21 greg Exp $ */
2 /*
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 /* General plane-ordered component matrix */
16 typedef struct {
17 int nrows, ncols, ncomp;
18 int dtype;
19 char *info;
20 double mtx[1]; /* extends struct */
21 } RMATRIX;
22
23 #define rmx_lval(rm,r,c,i) (rm)->mtx[(i)+(rm)->ncomp*((c)+(rm)->ncols*(r))]
24
25 /* Allocate a nr x nc matrix with n components */
26 extern RMATRIX *rmx_alloc(int nr, int nc, int n);
27
28 /* Free a RMATRIX array */
29 extern void rmx_free(RMATRIX *rm);
30
31 /* Resolve data type based on two input types (returns 0 for mismatch) */
32 extern int rmx_newtype(int dtyp1, int dtyp2);
33
34 /* Load matrix from supported file type (NULL for stdin, '!' with command) */
35 extern RMATRIX *rmx_load(const char *inspec);
36
37 /* Append header information associated with matrix data */
38 extern int rmx_addinfo(RMATRIX *rm, const char *info);
39
40 /* Write matrix to file type indicated by dtype */
41 extern int rmx_write(const RMATRIX *rm, int dtype, FILE *fp);
42
43 /* Allocate and assign square identity matrix with n components */
44 extern RMATRIX *rmx_identity(int dim, int n);
45
46 /* Duplicate the given matrix */
47 extern RMATRIX *rmx_copy(const RMATRIX *rm);
48
49 /* Allocate and assign transposed matrix */
50 extern RMATRIX *rmx_transpose(const RMATRIX *rm);
51
52 /* Multiply (concatenate) two matrices and allocate the result */
53 extern RMATRIX *rmx_multiply(const RMATRIX *m1, const RMATRIX *m2);
54
55 /* Element-wise multiplication (or division) of m2 into m1 */
56 extern int rmx_elemult(RMATRIX *m1, const RMATRIX *m2, int divide);
57
58 /* Sum second matrix into first, applying scale factor beforehand */
59 extern int rmx_sum(RMATRIX *msum, const RMATRIX *madd, const double sf[]);
60
61 /* Scale the given matrix by the indicated scalar component vector */
62 extern int rmx_scale(RMATRIX *rm, const double sf[]);
63
64 /* Allocate new matrix and apply component transformation */
65 extern RMATRIX *rmx_transform(const RMATRIX *msrc, int n, const double cmat[]);
66
67 /* Convert a color matrix to newly allocated RMATRIX buffer */
68 extern RMATRIX *rmx_from_cmatrix(const CMATRIX *cm);
69
70 /* Convert general matrix to newly allocated CMATRIX buffer */
71 extern CMATRIX *cm_from_rmatrix(const RMATRIX *rm);
72
73 #ifdef __cplusplus
74 }
75 #endif
76 #endif /* _RAD_RMATRIX_H_ */