Subversion Repositories gelsvn

Rev

Rev 349 | Only display areas with differences | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 349 Rev 353
1 
#include "glass.hpp"
 1 
#include "glass.h"
2 
 
 2 
 
3 
using namespace CGLA;
 3 
using namespace CGLA;
4 
 
 4 
 
5 
glass::glass(const CGLA::Vec3f& c, float ior)
 5 
glass::glass(const CGLA::Vec3f& c, float ior)
6 
: material(ior, 0.f), color_(c)
 6 
: material(ior, 0.f), color_(c)
7 
{
 7 
{
8 
}
 8 
}
9 
 
 9 
 
10 
void glass::sample(const ray& r, hit_info& hi) const
 10 
void glass::sample(const ray& r, hit_info& hi) const
11 
{
 11 
{
12 
    hi.ior = hi.inside ? 1.f / ior_ : ior_;
 12 
    hi.ior = hi.inside ? 1.f / ior_ : ior_;
13 
 
 13 
 
14 
    //compute fresnel
 14 
    //compute fresnel
15 
    float f;
 15 
    float f;
16 
    float cosi = dot(-r.direction, hi.shading_normal);
 16 
    float cosi = dot(-r.direction, hi.shading_normal);
17 
    cosi = clamp(cosi, -1.f, 1.f);
 17 
    cosi = clamp(cosi, -1.f, 1.f);
18 
    float sint = 1.f / hi.ior * std::sqrt(std::max(0.f, 1.f - cosi*cosi));
 18 
    float sint = 1.f / hi.ior * std::sqrt(std::max(0.f, 1.f - cosi*cosi));
19 
 
 19 
 
20 
    if (sint > 1.f)
 20 
    if (sint > 1.f)
21 
    {
 21 
    {
22 
        //total internal reflection
 22 
        //total internal reflection
23 
        f = 1.f;
 23 
        f = 1.f;
24 
    }
 24 
    }
25 
    else if (cosi < 0.)
 25 
    else if (cosi < 0.)
26 
    {
 26 
    {
27 
        f = 0.f;
 27 
        f = 0.f;
28 
    }
 28 
    }
29 
    else
 29 
    else
30 
    {
 30 
    {
31 
        float cost = std::sqrt(std::max(0.f, 1.f - sint * sint));
 31 
        float cost = std::sqrt(std::max(0.f, 1.f - sint * sint));
32 
        cost = clamp(cost, -1.f, 1.f);
 32 
        cost = clamp(cost, -1.f, 1.f);
33 
        f = fresnel_dielectric(cosi, cost, hi.ior);
 33 
        f = fresnel_dielectric(cosi, cost, hi.ior);
34 
 
 34 
 
35 
        if (f < 0.f || f > 1.f)
 35 
        if (f < 0.f || f > 1.f)
36 
            assert(false);
 36 
            assert(false);
37 
    }
 37 
    }
38 
 
 38 
 
39 
    hi.reflection = f * color_;
 39 
    hi.reflection = f * color_;
40 
    hi.refraction = (1.f - f) * color_;
 40 
    hi.refraction = (1.f - f) * color_;
41 
}
 41 
}
42 
 
 42 
 
43 
//02566 framework, Anders Wang Kristensen, awk@imm.dtu.dk, 2007
 43 
//02566 framework, Anders Wang Kristensen, awk@imm.dtu.dk, 2007
44 
 
 44 
 
45 
 
 45 
 
46 
Generated by GNU Enscript 1.6.6.
 46 
Generated by GNU Enscript 1.6.6.
47 
 
 47 
 
48 
 
 48 
 
49 
 
 49