21 |
|
#define ANG2R(r) (int)((r)*((1<<16)/M_PI)) |
22 |
|
#define R2ANG(c) (((c)+.5)*(M_PI/(1<<16))) |
23 |
|
|
24 |
+ |
/* moderated cosine factor */ |
25 |
+ |
#define COSF(z) (fabs(z)*0.98 + 0.02) |
26 |
+ |
|
27 |
|
typedef union { |
28 |
|
struct { |
29 |
|
float v; /* DSF sum */ |
34 |
|
|
35 |
|
typedef struct { |
36 |
|
float peak; /* lobe value at peak */ |
37 |
+ |
C_CHROMA chroma; /* encoded chromaticity */ |
38 |
|
unsigned short crad; /* radius (coded angle) */ |
39 |
< |
unsigned char gx, gy; /* grid position */ |
39 |
> |
unsigned short gx, gy; /* grid position */ |
40 |
|
} RBFVAL; /* radial basis function value */ |
41 |
|
|
42 |
|
struct s_rbfnode; /* forward declaration of RBF struct */ |
83 |
|
extern int input_orient; |
84 |
|
extern int output_orient; |
85 |
|
|
86 |
+ |
/* represented colorimetry */ |
87 |
+ |
typedef enum {RBCphotopic, RBCtristimulus, RBCspectral, RBCunknown} RBColor; |
88 |
+ |
|
89 |
+ |
extern RBColor rbf_colorimetry; |
90 |
+ |
|
91 |
+ |
extern const char *RBCident[]; |
92 |
+ |
|
93 |
|
/* log BSDF histogram */ |
94 |
|
#define HISTLEN 256 |
95 |
|
#define BSDF2BIG (1./M_PI) |
101 |
|
|
102 |
|
/* BSDF value for boundary regions */ |
103 |
|
extern double bsdf_min; |
104 |
+ |
extern double bsdf_spec_val; |
105 |
+ |
extern double bsdf_spec_rad; |
106 |
|
|
107 |
|
/* processed incident DSF measurements */ |
108 |
|
extern RBFNODE *dsf_list; |
122 |
|
|
123 |
|
#define BSDFREP_FMT "BSDF_RBFmesh" |
124 |
|
|
125 |
+ |
#define BSDFREP_MAGIC 0x5a3c |
126 |
+ |
|
127 |
|
/* global argv[0] */ |
128 |
|
extern char *progname; |
129 |
|
|
135 |
|
/* our loaded grid for this incident angle */ |
136 |
|
extern double theta_in_deg, phi_in_deg; |
137 |
|
extern GRIDVAL dsf_grid[GRIDRES][GRIDRES]; |
138 |
+ |
extern float (*spec_grid)[GRIDRES][GRIDRES]; |
139 |
+ |
extern int nspec_grid; |
140 |
|
|
141 |
|
/* Register new input direction */ |
142 |
|
extern int new_input_direction(double new_theta, double new_phi); |
165 |
|
/* Compute grid position from normalized input/output vector */ |
166 |
|
extern void pos_from_vec(int pos[2], const FVECT vec); |
167 |
|
|
168 |
< |
/* Evaluate RBF for DSF at the given normalized outgoing direction */ |
168 |
> |
/* Evaluate BSDF at the given normalized outgoing direction */ |
169 |
|
extern double eval_rbfrep(const RBFNODE *rp, const FVECT outvec); |
170 |
|
|
171 |
+ |
extern SDError eval_rbfcol(SDValue *sv, |
172 |
+ |
const RBFNODE *rp, const FVECT outvec); |
173 |
+ |
|
174 |
|
/* Insert a new directional scattering function in our global list */ |
175 |
|
extern int insert_dsf(RBFNODE *newrbf); |
176 |
|
|
197 |
|
/* Read a BSDF mesh interpolant from the given binary stream */ |
198 |
|
extern int load_bsdf_rep(FILE *ifp); |
199 |
|
|
200 |
+ |
/* Set up visible spectrum sampling */ |
201 |
+ |
extern void set_spectral_samples(int nspec); |
202 |
+ |
|
203 |
|
/* Start new DSF input grid */ |
204 |
|
extern void new_bsdf_data(double new_theta, double new_phi); |
205 |
|
|
206 |
|
/* Add BSDF data point */ |
207 |
|
extern void add_bsdf_data(double theta_out, double phi_out, |
208 |
< |
double val, int isDSF); |
208 |
> |
const double val[], int isDSF); |
209 |
|
|
210 |
|
/* Count up filled nodes and build RBF representation from current grid */ |
211 |
|
extern RBFNODE * make_rbfrep(void); |
215 |
|
|
216 |
|
/* Find edge(s) for interpolating the given vector, applying symmetry */ |
217 |
|
extern int get_interp(MIGRATION *miga[3], FVECT invec); |
218 |
+ |
|
219 |
+ |
/* Return single-lobe specular RBF for the given incident direction */ |
220 |
+ |
extern RBFNODE * def_rbf_spec(const FVECT invec); |
221 |
|
|
222 |
|
/* Advect and allocate new RBF along edge (internal call) */ |
223 |
|
extern RBFNODE * e_advect_rbf(const MIGRATION *mig, |