ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/source.h
Revision: 2.26
Committed: Fri Nov 15 20:47:42 2024 UTC (6 months, 1 week ago) by greg
Content type: text/plain
Branch: MAIN
Changes since 2.25: +25 -1 lines
Log Message:
feat(rpict): Experimental source skipping option with -DSSKIPOPT compile flag

File Contents

# Content
1 /* RCSid $Id: source.h,v 2.25 2024/11/09 15:21:32 greg Exp $ */
2 /*
3 * source.h - header file for ray tracing sources.
4 *
5 * Include after ray.h
6 */
7 #ifndef _RAD_SOURCE_H_
8 #define _RAD_SOURCE_H_
9
10 #include <string.h>
11
12 #ifdef __cplusplus
13 extern "C" {
14 #endif
15
16 #ifndef AIMREQT
17 #define AIMREQT 100 /* required aim success/failure */
18 #endif
19 #ifndef SHADCACHE
20 #define SHADCACHE 20 /* shadow cache resolution */
21 #endif
22 #ifndef MINSHADCNT
23 #define MINSHADCNT 2 /* test at least this many shadows */
24 #endif
25
26 #define SDISTANT 01 /* source distant flag */
27 #define SSKIP 02 /* source skip flag */
28 #define SPROX 04 /* source proximity flag */
29 #define SSPOT 010 /* source spotlight flag */
30 #define SVIRTUAL 020 /* source virtual flag */
31 #define SFLAT 040 /* source flat flag */
32 #define SCIR 0100 /* source circular flag */
33 #define SCYL 0200 /* source cylindrical flag */
34 #define SFOLLOW 0400 /* source follow path flag */
35
36 typedef struct {
37 FVECT aim; /* aim direction or center */
38 float siz; /* output solid angle or area */
39 float flen; /* focal length (negative if distant source) */
40 } SPOT; /* spotlight */
41
42 typedef struct {
43 union {
44 struct {
45 FVECT u, v; /* unit vectors */
46 } f; /* flat source indexing */
47 struct {
48 FVECT o; /* origin position */
49 RREAL e1, e2; /* 1/extent */
50 int ax; /* major direction */
51 } d; /* distant source indexing */
52 } p; /* indexing parameters */
53 OBJECT obs[1]; /* cache obstructors (extends struct) */
54 } OBSCACHE; /* obstructor cache */
55
56 typedef struct {
57 FVECT sloc; /* direction or position of source */
58 FVECT ss[3]; /* source dimension vectors, U, V, and W */
59 float srad; /* maximum source radius */
60 float ss2; /* solid angle or projected area */
61 OBJREC *so; /* source destination object */
62 struct {
63 float prox; /* proximity */
64 SPOT *s; /* spot */
65 } sl; /* localized source information */
66 union {
67 long success; /* successes - AIMREQT*failures */
68 struct {
69 short pn; /* projection number */
70 int sn; /* next source to aim for */
71 } sv; /* virtual source */
72 } sa; /* source aiming information */
73 unsigned long
74 ntests, nhits; /* shadow tests and hits */
75 #if SHADCACHE
76 OBSCACHE *obscache; /* obstructor cache */
77 #endif
78 int sflags; /* source flags */
79 } SRCREC; /* light source */
80
81 #define MAXSPART 64 /* maximum partitions per source */
82
83 #define SU 0 /* U vector or partition */
84 #define SV 1 /* V vector or partition */
85 #define SW 2 /* W vector or partition */
86 #define S0 3 /* leaf partition */
87
88 #define snorm ss[SW] /* normal vector for flat source */
89
90 typedef struct {
91 double dom; /* solid angle of partition */
92 int sn; /* source number */
93 short np; /* number of partitions */
94 short sp; /* this partition number */
95 unsigned char spt[MAXSPART/2]; /* source partitioning */
96 } SRCINDEX; /* source index structure */
97
98 #define initsrcindex(s) ((s)->sn = (s)->sp = -1, (s)->np = 0)
99
100 #define clrpart(pt) memset((char *)(pt), '\0', MAXSPART/2)
101 #define setpart(pt,i,v) ((pt)[(i)>>2] |= (v)<<(((i)&3)<<1))
102 #define spart(pt,pi) ((pt)[(pi)>>2] >> (((pi)&3)<<1) & 3)
103
104 /*
105 * Special support functions for sources
106 */
107
108 /*
109 * Virtual source materials must define the following.
110 *
111 * vproj(pm, op, sp, i) Compute i'th virtual projection
112 * of source sp in object op and assign
113 * the 4x4 transformation matrix pm.
114 * Return 1 on success, 0 if no i'th projection.
115 *
116 * nproj The number of projections. The value of
117 * i passed to vproj runs from 0 to nproj-1.
118 */
119
120 typedef struct {
121 int (*vproj)(); /* project virtual sources */
122 int nproj; /* number of possible projections */
123 } VSMATERIAL; /* virtual source material functions */
124
125 typedef struct {
126 void (*setsrc)(); /* set light source for object */
127 void (*partit)(); /* partition light source object */
128 double (*getpleq)(); /* plane equation for surface */
129 double (*getdisk)(); /* maximum disk for surface */
130 } SOBJECT; /* source object functions */
131
132 typedef union {
133 VSMATERIAL *mf; /* material functions */
134 SOBJECT *of; /* object functions */
135 } SRCFUNC; /* source functions */
136
137 extern SRCFUNC sfun[]; /* source dispatch table */
138
139 extern SRCREC *source; /* our source list */
140 extern int nsources; /* the number of sources */
141
142 #define sflatform(sn,dir) -DOT(source[sn].snorm, dir)
143
144 #define getplaneq(c,o) (*sfun[(o)->otype].of->getpleq)(c,o)
145 #define getmaxdisk(c,o) (*sfun[(o)->otype].of->getdisk)(c,o)
146 #define setsource(s,o) (*sfun[(o)->otype].of->setsrc)(s,o)
147
148 #define SSKIPFLSIZ ((nsources+7) >> 3)
149
150 extern uby8 *ssf_select; /* sources we may skip */
151
152 #define sskip_new() ((uby8 *)ecalloc(1,SSKIPFLSIZ))
153 #define sskip_free(fl) (efree(fl), fl=NULL)
154 #define sskip_eq(fl1,fl2) !memcmp(fl1, fl2, SSKIPFLSIZ)
155 #define sskip_cpy(dfl,sfl) memcpy(dfl,sfl,SSKIPFLSIZ)
156 #define sskip_op(fl,op,sn) ((fl)[(sn)>>3] op (1<<((sn)&7)))
157 #define sskip_chk(fl,sn) sskip_op(fl,&,sn)
158 #define sskip_set(fl,sn) sskip_op(fl,|=,sn)
159 #define sskip_clr(fl,sn) sskip_op(fl,&=~,sn)
160
161 /* defined in source.c */
162 extern void marksources(void);
163 extern void distantsources(void);
164 extern void freesources(void);
165 extern int srcray(RAY *sr, RAY *r, SRCINDEX *si);
166 extern void srcvalue(RAY *r);
167 extern int sourcehit(RAY *r);
168 typedef void srcdirf_t(SCOLOR cv, void *np, FVECT ldir, double omega);
169 extern void direct(RAY *r, srcdirf_t *f, void *p);
170 extern void srcscatter(RAY *r);
171 extern int m_light(OBJREC *m, RAY *r);
172 /* defined in srcobstr.c */
173 extern void initobscache(int sn);
174 extern int srcblocker(RAY *r);
175 extern int srcblocked(RAY *r);
176 extern void freeobscache(SRCREC *s);
177 extern void markclip(OBJREC *m);
178 /* defined in srcsamp.c */
179 extern int sskip_rsi(uby8 *flags);
180 extern uby8 *sskip_flags(int rsi);
181 extern void sskip_addflags(uby8 *dfl, const uby8 *sfl);
182 extern int srcskip(int sn, RAY *r);
183 extern double nextssamp(RAY *r, SRCINDEX *si);
184 extern int skipparts(int ct[3], int sz[3], int pp[2], unsigned char *pt);
185 extern void nopart(SRCINDEX *si, RAY *r);
186 extern void cylpart(SRCINDEX *si, RAY *r);
187 extern void flatpart(SRCINDEX *si, RAY *r);
188 extern double scylform(int sn, FVECT dir);
189 /* defined in srcskipload.c */
190 extern int sskip_dim[2]; /* source skip image size */
191 extern int srcskip_open(char *bmpspec, char *scorrimg);
192 extern int srcskip_getrow(int row, int *sndx, float *scorr);
193 extern int *srcskip_ndxmap(void);
194 extern float *srcskip_corrmap(void);
195 extern void srcskip_close(void);
196 extern void srcskip_free_maps(void);
197 /* defined in srcsupp.c */
198 extern void initstypes(void);
199 extern int newsource(void);
200 extern void setflatss(SRCREC *src);
201 extern void fsetsrc(SRCREC *src, OBJREC *so);
202 extern void ssetsrc(SRCREC *src, OBJREC *so);
203 extern void sphsetsrc(SRCREC *src, OBJREC *so);
204 extern void rsetsrc(SRCREC *src, OBJREC *so);
205 extern void cylsetsrc(SRCREC *src, OBJREC *so);
206 extern SPOT *makespot(OBJREC *m);
207 extern int spotout(RAY *r, SPOT *s);
208 extern double fgetmaxdisk(FVECT ocent, OBJREC *op);
209 extern double rgetmaxdisk(FVECT ocent, OBJREC *op);
210 extern double fgetplaneq(FVECT nvec, OBJREC *op);
211 extern double rgetplaneq(FVECT nvec, OBJREC *op);
212 extern int commonspot(SPOT *sp1, SPOT *sp2, FVECT org);
213 extern int commonbeam(SPOT *sp1, SPOT *sp2, FVECT org);
214 extern int checkspot(SPOT *sp, FVECT nrm);
215 extern double spotdisk(FVECT oc, OBJREC *op, SPOT *sp, FVECT pos);
216 extern double beamdisk(FVECT oc, OBJREC *op, SPOT *sp, FVECT dir);
217 extern double intercircle(FVECT cc, FVECT c1, FVECT c2,
218 double r1s, double r2s);
219 /* defined in virtuals.c */
220 extern void markvirtuals(void);
221 extern void addvirtuals(int sn, int nr);
222 extern void vproject(OBJREC *o, int sn, int n);
223 extern OBJREC *vsmaterial(OBJREC *o);
224 extern int makevsrc(OBJREC *op, int sn, MAT4 pm);
225 extern double getdisk(FVECT oc, OBJREC *op, int sn);
226 extern int vstestvis(int f, OBJREC *o, FVECT oc, double or2, int sn);
227 extern void virtverb(int sn, FILE *fp);
228
229
230 #ifdef __cplusplus
231 }
232 #endif
233 #endif /* _RAD_SOURCE_H_ */
234