| 1 |
< |
/* RCSid: $Id$ */ |
| 1 |
> |
/* RCSid $Id$ */ |
| 2 |
|
/* |
| 3 |
|
* ray.h - header file for routines using rays. |
| 4 |
|
*/ |
| 5 |
+ |
#ifndef _RAD_RAY_H_ |
| 6 |
+ |
#define _RAD_RAY_H_ |
| 7 |
|
|
| 6 |
– |
/* ==================================================================== |
| 7 |
– |
* The Radiance Software License, Version 1.0 |
| 8 |
– |
* |
| 9 |
– |
* Copyright (c) 1990 - 2002 The Regents of the University of California, |
| 10 |
– |
* through Lawrence Berkeley National Laboratory. All rights reserved. |
| 11 |
– |
* |
| 12 |
– |
* Redistribution and use in source and binary forms, with or without |
| 13 |
– |
* modification, are permitted provided that the following conditions |
| 14 |
– |
* are met: |
| 15 |
– |
* |
| 16 |
– |
* 1. Redistributions of source code must retain the above copyright |
| 17 |
– |
* notice, this list of conditions and the following disclaimer. |
| 18 |
– |
* |
| 19 |
– |
* 2. Redistributions in binary form must reproduce the above copyright |
| 20 |
– |
* notice, this list of conditions and the following disclaimer in |
| 21 |
– |
* the documentation and/or other materials provided with the |
| 22 |
– |
* distribution. |
| 23 |
– |
* |
| 24 |
– |
* 3. The end-user documentation included with the redistribution, |
| 25 |
– |
* if any, must include the following acknowledgment: |
| 26 |
– |
* "This product includes Radiance software |
| 27 |
– |
* (http://radsite.lbl.gov/) |
| 28 |
– |
* developed by the Lawrence Berkeley National Laboratory |
| 29 |
– |
* (http://www.lbl.gov/)." |
| 30 |
– |
* Alternately, this acknowledgment may appear in the software itself, |
| 31 |
– |
* if and wherever such third-party acknowledgments normally appear. |
| 32 |
– |
* |
| 33 |
– |
* 4. The names "Radiance," "Lawrence Berkeley National Laboratory" |
| 34 |
– |
* and "The Regents of the University of California" must |
| 35 |
– |
* not be used to endorse or promote products derived from this |
| 36 |
– |
* software without prior written permission. For written |
| 37 |
– |
* permission, please contact [email protected]. |
| 38 |
– |
* |
| 39 |
– |
* 5. Products derived from this software may not be called "Radiance", |
| 40 |
– |
* nor may "Radiance" appear in their name, without prior written |
| 41 |
– |
* permission of Lawrence Berkeley National Laboratory. |
| 42 |
– |
* |
| 43 |
– |
* THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED |
| 44 |
– |
* WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES |
| 45 |
– |
* OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE |
| 46 |
– |
* DISCLAIMED. IN NO EVENT SHALL Lawrence Berkeley National Laboratory OR |
| 47 |
– |
* ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, |
| 48 |
– |
* SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT |
| 49 |
– |
* LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF |
| 50 |
– |
* USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND |
| 51 |
– |
* ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, |
| 52 |
– |
* OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT |
| 53 |
– |
* OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF |
| 54 |
– |
* SUCH DAMAGE. |
| 55 |
– |
* ==================================================================== |
| 56 |
– |
* |
| 57 |
– |
* This software consists of voluntary contributions made by many |
| 58 |
– |
* individuals on behalf of Lawrence Berkeley National Laboratory. For more |
| 59 |
– |
* information on Lawrence Berkeley National Laboratory, please see |
| 60 |
– |
* <http://www.lbl.gov/>. |
| 61 |
– |
*/ |
| 62 |
– |
|
| 8 |
|
#include "standard.h" |
| 64 |
– |
|
| 9 |
|
#include "octree.h" |
| 66 |
– |
|
| 10 |
|
#include "object.h" |
| 68 |
– |
|
| 11 |
|
#include "color.h" |
| 12 |
+ |
#include "pmapparm.h" |
| 13 |
|
|
| 14 |
+ |
#ifdef __cplusplus |
| 15 |
+ |
extern "C" { |
| 16 |
+ |
#endif |
| 17 |
+ |
|
| 18 |
+ |
#ifndef RNUMBER |
| 19 |
+ |
#define RNUMBER unsigned long /* ray counter (>= sizeof pointer) */ |
| 20 |
+ |
#endif |
| 21 |
+ |
|
| 22 |
|
#define MAXDIM 32 /* maximum number of dimensions */ |
| 23 |
|
|
| 24 |
|
/* ray type flags */ |
| 33 |
|
/* reflected ray types */ |
| 34 |
|
#define RAYREFL (SHADOW|REFLECTED|AMBIENT|SPECULAR) |
| 35 |
|
|
| 36 |
+ |
/* Arrange so double's come first for optimal alignment */ |
| 37 |
+ |
/* Pointers and long's come second for 64-bit mode */ |
| 38 |
+ |
/* Int's next (unknown length), then floats, followed by short's & char's */ |
| 39 |
|
typedef struct ray { |
| 86 |
– |
unsigned long rno; /* unique ray number */ |
| 87 |
– |
int rlvl; /* number of reflections for this ray */ |
| 88 |
– |
float rweight; /* cumulative weight of this ray */ |
| 89 |
– |
short rtype; /* ray type */ |
| 90 |
– |
short crtype; /* cumulative ray type */ |
| 91 |
– |
struct ray *parent; /* ray this originated from */ |
| 40 |
|
FVECT rorg; /* origin of ray */ |
| 41 |
|
FVECT rdir; /* normalized direction of ray */ |
| 42 |
< |
double rmax; /* maximum distance (aft clipping plane) */ |
| 43 |
< |
int rsrc; /* source we're aiming for */ |
| 42 |
> |
RREAL rmax; /* maximum distance (aft clipping plane) */ |
| 43 |
> |
RREAL rot; /* distance to object */ |
| 44 |
> |
FVECT rop; /* intersection point */ |
| 45 |
> |
FVECT ron; /* intersection surface normal */ |
| 46 |
> |
RREAL rod; /* -DOT(rdir, ron) */ |
| 47 |
> |
RREAL uv[2]; /* local coordinates */ |
| 48 |
> |
FVECT pert; /* surface normal perturbation */ |
| 49 |
> |
RREAL rmt; /* returned mirrored ray length */ |
| 50 |
> |
RREAL rxt; /* returned unmirrored ray length */ |
| 51 |
> |
const struct ray *parent; /* ray this originated from */ |
| 52 |
|
OBJECT *clipset; /* set of objects currently clipped */ |
| 53 |
|
OBJECT *newcset; /* next clipset, used for transmission */ |
| 54 |
< |
void (*revf)(); /* evaluation function for this ray */ |
| 55 |
< |
OBJECT robj; /* intersected object number */ |
| 54 |
> |
void (*revf)(struct ray *); /* ray evaluation function */ |
| 55 |
> |
void (*hitf)(OBJECT *, struct ray *); /* custom hit test */ |
| 56 |
|
OBJREC *ro; /* intersected object (one with material) */ |
| 101 |
– |
double rot; /* distance to object */ |
| 102 |
– |
FVECT rop; /* intersection point */ |
| 103 |
– |
FVECT ron; /* intersection surface normal */ |
| 104 |
– |
double rod; /* -DOT(rdir, ron) */ |
| 57 |
|
FULLXF *rox; /* object transformation */ |
| 58 |
< |
FVECT pert; /* surface normal perturbation */ |
| 58 |
> |
int *slights; /* list of lights to test for scattering */ |
| 59 |
> |
RNUMBER rno; /* unique ray number */ |
| 60 |
> |
OBJECT robj; /* intersected object number */ |
| 61 |
> |
int rsrc; /* source we're aiming for */ |
| 62 |
> |
float rweight; /* cumulative weight (for termination) */ |
| 63 |
> |
float gecc; /* scattering eccentricity coefficient */ |
| 64 |
> |
COLOR rcoef; /* contribution coefficient w.r.t. parent */ |
| 65 |
|
COLOR pcol; /* pattern color */ |
| 66 |
< |
COLOR rcol; /* returned ray value */ |
| 67 |
< |
double rt; /* returned effective ray length */ |
| 66 |
> |
COLOR mcol; /* mirrored color contribution */ |
| 67 |
> |
COLOR rcol; /* returned radiance value */ |
| 68 |
|
COLOR cext; /* medium extinction coefficient */ |
| 69 |
|
COLOR albedo; /* medium scattering albedo */ |
| 70 |
< |
float gecc; /* scattering eccentricity coefficient */ |
| 71 |
< |
int *slights; /* list of lights to test for scattering */ |
| 70 |
> |
short rflips; /* surface orientation has been reversed */ |
| 71 |
> |
short rlvl; /* number of reflections for this ray */ |
| 72 |
> |
short rtype; /* ray type */ |
| 73 |
> |
short crtype; /* cumulative ray type */ |
| 74 |
|
} RAY; |
| 75 |
|
|
| 76 |
|
#define rayvalue(r) (*(r)->revf)(r) |
| 77 |
|
|
| 78 |
+ |
#define raydistance(r) (bright((r)->mcol) > 0.5*bright((r)->rcol) ? \ |
| 79 |
+ |
(r)->rmt : (r)->rxt) |
| 80 |
+ |
|
| 81 |
+ |
#define rayreorient(r) if ((r)->rflips & 1) flipsurface(r); else |
| 82 |
+ |
|
| 83 |
|
extern char VersionID[]; /* Radiance version ID string */ |
| 84 |
|
|
| 85 |
|
extern CUBE thescene; /* our scene */ |
| 86 |
|
extern OBJECT nsceneobjs; /* number of objects in our scene */ |
| 87 |
|
|
| 88 |
< |
extern unsigned long raynum; /* next ray ID */ |
| 89 |
< |
extern unsigned long nrays; /* total rays traced so far */ |
| 88 |
> |
extern RNUMBER raynum; /* next ray ID */ |
| 89 |
> |
extern RNUMBER nrays; /* total rays traced so far */ |
| 90 |
|
|
| 91 |
|
extern OBJREC Lamb; /* a Lambertian surface */ |
| 92 |
|
extern OBJREC Aftplane; /* aft clipping object */ |
| 101 |
|
|
| 102 |
|
extern int do_irrad; /* compute irradiance? */ |
| 103 |
|
|
| 104 |
+ |
extern int rand_samp; /* pure Monte Carlo sampling? */ |
| 105 |
+ |
|
| 106 |
|
extern double dstrsrc; /* square source distribution */ |
| 107 |
|
extern double shadthresh; /* shadow threshold */ |
| 108 |
|
extern double shadcert; /* shadow testing certainty */ |
| 138 |
|
extern int ray_pnprocs; /* number of child processes */ |
| 139 |
|
extern int ray_pnidle; /* number of idle processes */ |
| 140 |
|
|
| 141 |
< |
#define AMBLLEN 128 /* max. ambient list length */ |
| 142 |
< |
#define AMBWORD 8 /* average word length */ |
| 141 |
> |
#ifndef AMBLLEN |
| 142 |
> |
#define AMBLLEN 512 /* max. ambient list length */ |
| 143 |
> |
#endif |
| 144 |
> |
#define AMBWORD 12 /* average word length */ |
| 145 |
|
|
| 146 |
|
typedef struct { /* rendering parameter holder */ |
| 147 |
|
int do_irrad; |
| 148 |
+ |
int rand_samp; |
| 149 |
|
double dstrsrc; |
| 150 |
|
double shadthresh; |
| 151 |
|
double shadcert; |
| 173 |
|
int ambincl; |
| 174 |
|
short amblndx[AMBLLEN+1]; |
| 175 |
|
char amblval[AMBLLEN*AMBWORD]; |
| 176 |
+ |
|
| 177 |
+ |
/* PMAP: photon mapping parameters */ |
| 178 |
+ |
PhotonMapParams pmapParams [NUM_PMAP_TYPES]; |
| 179 |
|
} RAYPARAMS; |
| 180 |
|
|
| 181 |
|
#define rpambmod(p,i) ( (i)>=AMBLLEN||(p)->amblndx[i]<0 ? \ |
| 182 |
|
(char *)NULL : (p)->amblval+(p)->amblndx[i] ) |
| 183 |
|
|
| 211 |
– |
#ifdef NOPROTO |
| 212 |
– |
|
| 213 |
– |
extern void headclean(); |
| 214 |
– |
extern void openheader(); |
| 215 |
– |
extern void dupheader(); |
| 216 |
– |
extern void pfdetach(); |
| 217 |
– |
extern void pfclean(); |
| 218 |
– |
extern void pflock(); |
| 219 |
– |
extern void pfhold(); |
| 220 |
– |
extern void io_process(); |
| 221 |
– |
extern int free_objs(); |
| 222 |
– |
extern void free_objmem(); |
| 223 |
– |
extern int load_os(); |
| 224 |
– |
extern void preload_objs(); |
| 225 |
– |
extern void ray_init(); |
| 226 |
– |
extern void ray_trace(); |
| 227 |
– |
extern void ray_done(); |
| 228 |
– |
extern void ray_save(); |
| 229 |
– |
extern void ray_restore(); |
| 230 |
– |
extern void ray_defaults(); |
| 231 |
– |
extern void ray_pinit(); |
| 232 |
– |
extern void ray_psend(); |
| 233 |
– |
extern int ray_pqueue(); |
| 234 |
– |
extern int ray_presult(); |
| 235 |
– |
extern void ray_pdone(); |
| 236 |
– |
extern void ray_popen(); |
| 237 |
– |
extern void ray_pclose(); |
| 238 |
– |
extern int rayorigin(); |
| 239 |
– |
extern void rayclear(); |
| 240 |
– |
extern void raytrace(); |
| 241 |
– |
extern void raycont(); |
| 242 |
– |
extern void raytrans(); |
| 243 |
– |
extern int rayshade(); |
| 244 |
– |
extern void rayparticipate(); |
| 245 |
– |
extern int raymixture(); |
| 246 |
– |
extern double raydist(); |
| 247 |
– |
extern double raynormal(); |
| 248 |
– |
extern void newrayxf(); |
| 249 |
– |
extern void flipsurface(); |
| 250 |
– |
extern int localhit(); |
| 251 |
– |
extern int getrenderopt(); |
| 252 |
– |
extern void print_rdefaults(); |
| 253 |
– |
extern void drawsources(); |
| 254 |
– |
extern void rtrace(); |
| 255 |
– |
extern void rview(); |
| 256 |
– |
extern void rpict(); |
| 257 |
– |
|
| 258 |
– |
#else |
| 184 |
|
/* defined in duphead.c */ |
| 185 |
|
extern void headclean(void); |
| 186 |
|
extern void openheader(void); |
| 187 |
|
extern void dupheader(void); |
| 188 |
+ |
/* defined in persist.c */ |
| 189 |
+ |
extern void persistfile(char *pfn); |
| 190 |
|
extern void pfdetach(void); |
| 191 |
|
extern void pfclean(void); |
| 192 |
|
extern void pflock(int lf); |
| 207 |
|
extern void ray_defaults(RAYPARAMS *rp); |
| 208 |
|
/* defined in raypcalls.c */ |
| 209 |
|
extern void ray_pinit(char *otnm, int nproc); |
| 210 |
< |
extern void ray_psend(RAY *r); |
| 210 |
> |
extern int ray_psend(RAY *r); |
| 211 |
|
extern int ray_pqueue(RAY *r); |
| 212 |
|
extern int ray_presult(RAY *r, int poll); |
| 213 |
|
extern void ray_pdone(int freall); |
| 214 |
|
extern void ray_popen(int nadd); |
| 215 |
|
extern void ray_pclose(int nsub); |
| 216 |
+ |
/* defined in ray_fifo.c */ |
| 217 |
+ |
extern int (*ray_fifo_out)(RAY *r); |
| 218 |
+ |
extern int ray_fifo_in(RAY *r); |
| 219 |
+ |
extern int ray_fifo_flush(void); |
| 220 |
|
/* defined in raytrace.c */ |
| 221 |
< |
extern int rayorigin(RAY *r, RAY *ro, int rt, double rw); |
| 221 |
> |
extern int rayorigin(RAY *r, int rt, const RAY *ro, const COLOR rc); |
| 222 |
|
extern void rayclear(RAY *r); |
| 223 |
|
extern void raytrace(RAY *r); |
| 224 |
+ |
extern int rayreject(OBJREC *o, RAY *r, double t); |
| 225 |
+ |
extern void rayhit(OBJECT *oset, RAY *r); |
| 226 |
|
extern void raycont(RAY *r); |
| 227 |
|
extern void raytrans(RAY *r); |
| 228 |
+ |
extern int raytirrad(OBJREC *m, RAY *r); |
| 229 |
|
extern int rayshade(RAY *r, int mod); |
| 230 |
|
extern void rayparticipate(RAY *r); |
| 231 |
+ |
extern void raytexture(RAY *r, OBJECT mod); |
| 232 |
|
extern int raymixture(RAY *r, OBJECT fore, OBJECT back, double coef); |
| 233 |
< |
extern double raydist(RAY *r, int flags); |
| 233 |
> |
extern void raycontrib(RREAL rc[3], const RAY *r, int flags); |
| 234 |
> |
extern double raydist(const RAY *r, int flags); |
| 235 |
|
extern double raynormal(FVECT norm, RAY *r); |
| 236 |
|
extern void newrayxf(RAY *r); |
| 237 |
|
extern void flipsurface(RAY *r); |
| 242 |
|
/* defined in srcdraw.c */ |
| 243 |
|
extern void drawsources(COLOR *pic[], float *zbf[], |
| 244 |
|
int x0, int xsiz, int y0, int ysiz); |
| 245 |
+ |
extern void init_drawsources(int rad); |
| 246 |
+ |
/* defined in rt/initotypes.c */ |
| 247 |
+ |
extern void initotypes(void); |
| 248 |
|
/* module main procedures */ |
| 249 |
< |
extern void rtrace(char *fname); |
| 249 |
> |
extern void rtrace(char *fname, int nproc); |
| 250 |
> |
extern char *formstr(int f); |
| 251 |
|
extern void rview(void); |
| 252 |
|
extern void rpict(int seq, char *pout, char *zout, char *prvr); |
| 253 |
|
|
| 254 |
+ |
#ifdef __FAST_MATH__ |
| 255 |
+ |
#define checknorm(vn) (void)normalize(vn) |
| 256 |
+ |
#else |
| 257 |
+ |
#define checknorm(vn) |
| 258 |
|
#endif |
| 259 |
+ |
|
| 260 |
+ |
#ifdef __cplusplus |
| 261 |
+ |
} |
| 262 |
+ |
#endif |
| 263 |
+ |
#endif /* _RAD_RAY_H_ */ |
| 264 |
+ |
|
