rev |
line source |
nuclear@0
|
1 /*
|
nuclear@0
|
2 libvmath - a vector math library
|
nuclear@0
|
3 Copyright (C) 2004-2011 John Tsiombikas <nuclear@member.fsf.org>
|
nuclear@0
|
4
|
nuclear@0
|
5 This program is free software: you can redistribute it and/or modify
|
nuclear@0
|
6 it under the terms of the GNU Lesser General Public License as published
|
nuclear@0
|
7 by the Free Software Foundation, either version 3 of the License, or
|
nuclear@0
|
8 (at your option) any later version.
|
nuclear@0
|
9
|
nuclear@0
|
10 This program is distributed in the hope that it will be useful,
|
nuclear@0
|
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
|
nuclear@0
|
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
nuclear@0
|
13 GNU Lesser General Public License for more details.
|
nuclear@0
|
14
|
nuclear@0
|
15 You should have received a copy of the GNU Lesser General Public License
|
nuclear@0
|
16 along with this program. If not, see <http://www.gnu.org/licenses/>.
|
nuclear@0
|
17 */
|
nuclear@0
|
18
|
nuclear@0
|
19 #ifdef __cplusplus
|
nuclear@0
|
20 extern "C" {
|
nuclear@0
|
21 #endif /* __cplusplus */
|
nuclear@0
|
22
|
nuclear@0
|
23 static inline ray_t ray_cons(vec3_t origin, vec3_t dir)
|
nuclear@0
|
24 {
|
nuclear@0
|
25 ray_t r;
|
nuclear@0
|
26 r.origin = origin;
|
nuclear@0
|
27 r.dir = dir;
|
nuclear@0
|
28 return r;
|
nuclear@0
|
29 }
|
nuclear@0
|
30
|
nuclear@0
|
31 #ifdef __cplusplus
|
nuclear@0
|
32 }
|
nuclear@0
|
33
|
nuclear@0
|
34 inline Ray reflect_ray(const Ray &inray, const Vector3 &norm)
|
nuclear@0
|
35 {
|
nuclear@0
|
36 Ray ray = inray;
|
nuclear@0
|
37 ray.iter--;
|
nuclear@0
|
38 ray.dir = ray.dir.reflection(norm);
|
nuclear@0
|
39 return ray;
|
nuclear@0
|
40 }
|
nuclear@0
|
41
|
nuclear@0
|
42 inline Ray refract_ray(const Ray &inray, const Vector3 &norm, scalar_t mat_ior, bool entering, scalar_t ray_mag)
|
nuclear@0
|
43 {
|
nuclear@0
|
44 Ray ray = inray;
|
nuclear@0
|
45 ray.iter--;
|
nuclear@0
|
46
|
nuclear@0
|
47 scalar_t ior = ray.calc_ior(entering, mat_ior);
|
nuclear@0
|
48
|
nuclear@0
|
49 if(entering) {
|
nuclear@0
|
50 ray.enter(mat_ior);
|
nuclear@0
|
51 } else {
|
nuclear@0
|
52 ray.leave();
|
nuclear@0
|
53 }
|
nuclear@0
|
54
|
nuclear@0
|
55 if(ray_mag < 0.0) {
|
nuclear@0
|
56 ray_mag = ray.dir.length();
|
nuclear@0
|
57 }
|
nuclear@0
|
58 ray.dir = (ray.dir / ray_mag).refraction(norm, ior) * ray_mag;
|
nuclear@0
|
59
|
nuclear@0
|
60 /* check TIR */
|
nuclear@0
|
61 if(dot_product(ray.dir, norm) > 0.0) {
|
nuclear@0
|
62 if(entering) {
|
nuclear@0
|
63 ray.leave();
|
nuclear@0
|
64 } else {
|
nuclear@0
|
65 ray.enter(mat_ior);
|
nuclear@0
|
66 }
|
nuclear@0
|
67 }
|
nuclear@0
|
68 return ray;
|
nuclear@0
|
69 }
|
nuclear@0
|
70 #endif /* __cplusplus */
|