WallyHackenslacker/PhotonMF
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

More photons. Still buggy as hell.

Browse Source
This commit is contained in:
Wally Hackenslacker
2017-02-22 21:14:08 -04:00
parent e361dc516e
commit c6060521d4
5 changed files with 96 additions and 161 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
kd_tree.cpp
View File

@@ -140,7 +140,7 @@ void kdTree::addPhoton(Photon p)
} }
void kdTree::createNodeKdTree(treeNode** node, std::vector<Photon> originalData , int* xyz, int* yzx, int* zxy, superKey key, int begin, int end, int* xyz_2, int* yzx_2, int* zxy_2) void kdTree::createNodeKdTree(treeNode** node, std::vector<Photon> & originalData , int* xyz, int* yzx, int* zxy, superKey key, int begin, int end, int* xyz_2, int* yzx_2, int* zxy_2)
{ {
if(end - begin < 2) if(end - begin < 2)
{ {
@@ -196,7 +196,7 @@ void kdTree::createNodeKdTree(treeNode** node, std::vector<Photon> originalData
createNodeKdTree((*node)->getRightChild(), originalData, xyz_2, yzx_2, zxy_2, key, mid + 1, end, xyz, yzx, zxy); createNodeKdTree((*node)->getRightChild(), originalData, xyz_2, yzx_2, zxy_2, key, mid + 1, end, xyz, yzx, zxy);
} }
void kdTree::reorderArrays(std::vector<Photon> originalData, int* A1, int* A2, int begin, int mid, int end, int orderIndex, superKey key, int* B1, int* B2) void kdTree::reorderArrays(std::vector<Photon> & originalData, int* A1, int* A2, int begin, int mid, int end, int orderIndex, superKey key, int* B1, int* B2)
{ {
int lowerindex1 = begin, higherindex1 = mid + 1, lowerindex2 = begin, higherindex2 = mid + 1; int lowerindex1 = begin, higherindex1 = mid + 1, lowerindex2 = begin, higherindex2 = mid + 1;

16
kd_tree.hpp
View File

@@ -43,14 +43,14 @@ struct Vec3
struct Photon struct Photon
{ {
Vec3 position; Vec3 position;
Vec3 direction;
float ref_index;
unsigned char radiance[4]; unsigned char radiance[4];
char phi, theta;
short unused_flag; // for 20 bytes struct.
Photon(Vec3 _p = Vec3(), float red = 0.0f, float green = 0.0f, float blue = 0.0f, char _phi = 0.0f, char _theta = 0.0f): Photon(Vec3 _p = Vec3(), Vec3 _d = Vec3(), float red = 0.0f, float green = 0.0f, float blue = 0.0f, float _r = 1.0f):
position(_p), position(_p),
phi(_phi), direction(_d),
theta(_theta) ref_index(_r)
{ {
float2rgbe(radiance, red, green, blue); float2rgbe(radiance, red, green, blue);
} }
@@ -159,7 +159,7 @@ private:
std::vector<Photon> Photons; std::vector<Photon> Photons;
void createNodeKdTree(treeNode** node, void createNodeKdTree(treeNode** node,
std::vector<Photon> originalData , std::vector<Photon> & originalData ,
int* xyz, int* xyz,
int* yzx, int* yzx,
int* zxy, int* zxy,
@@ -170,7 +170,7 @@ private:
int* yzx_2, int* yzx_2,
int* zxy_2); int* zxy_2);
void reorderArrays(std::vector<Photon> originalData, void reorderArrays(std::vector<Photon> & originalData,
int* A1, int* A1,
int* A2, int* A2,
int begin, int begin,
@@ -183,5 +183,5 @@ private:
void printNode(treeNode* node); void printNode(treeNode* node);
void findInRange (Vec3 min, Vec3 max, std::vector<Photon> &photons, treeNode *node) const; void findInRange (Vec3 min, Vec3 max, std::vector<Photon> & photons, treeNode *node) const;
}; };

233
photon_tracer.cpp
View File

@@ -22,17 +22,14 @@ using namespace glm;
PhotonTracer::~PhotonTracer() { } PhotonTracer::~PhotonTracer() { }
vec3 PhotonTracer::trace_ray(Ray & r, Scene * s, unsigned int rec_level) const { vec3 PhotonTracer::trace_ray(Ray & r, Scene * s, unsigned int rec_level) const {
float t, _t; float t, _t, radius, red, green, blue, kr;
Figure * _f; Figure * _f;
vec3 n, color, i_pos, ref, dir_diff_color, dir_spec_color, p_contrib; vec3 n, color, i_pos, ref, dir_diff_color, dir_spec_color, p_contrib;
Ray mv_r, sr, rr; Ray mv_r, sr, rr;
bool vis, is_area_light; bool vis, is_area_light;
float kr;
AreaLight * al; AreaLight * al;
Vec3 mn; Vec3 mn, mx;
Vec3 mx;
vector<Photon> photons; vector<Photon> photons;
float red, green, blue;
t = numeric_limits<float>::max(); t = numeric_limits<float>::max();
_f = NULL; _f = NULL;
@@ -102,7 +99,26 @@ vec3 PhotonTracer::trace_ray(Ray & r, Scene * s, unsigned int rec_level) const {
dir_spec_color += vis ? s->m_lights[l]->specular(n, r, i_pos, *_f->m_mat) : vec3(0.0f); dir_spec_color += vis ? s->m_lights[l]->specular(n, r, i_pos, *_f->m_mat) : vec3(0.0f);
} }
color += (dir_diff_color * (_f->m_mat->m_diffuse / pi<float>())) + (_f->m_mat->m_specular * dir_spec_color); // TODO: Change photon map search method for hemisphere search.
radius = m_h_radius;
mn = Vec3(i_pos.x - radius, i_pos.y - radius, i_pos.z - radius);
mx = Vec3(i_pos.x + radius, i_pos.y + radius, i_pos.z + radius);
while((photons = m_photon_map.findInRange(mn, mx)).size() == 0 && radius < 5.0) {
radius *= 2;
mn = Vec3(i_pos.x - radius, i_pos.y - radius, i_pos.z - radius);
mx = Vec3(i_pos.x + radius, i_pos.y + radius, i_pos.z + radius);
}
if (photons.size() > 0) {
for (vector<Photon>::iterator it = photons.begin(); it != photons.end(); it++) {
(*it).getColor(red, green, blue);
p_contrib += vec3(red, green, blue);
}
p_contrib /= pi<float>() * (radius * radius) * photons.size();
}
color += ((dir_diff_color + p_contrib) * (_f->m_mat->m_diffuse / pi<float>())) + (_f->m_mat->m_specular * dir_spec_color);
// Determine the specular reflection color. // Determine the specular reflection color.
if (_f->m_mat->m_rho > 0.0f && rec_level < m_max_depth) { if (_f->m_mat->m_rho > 0.0f && rec_level < m_max_depth) {
@@ -111,17 +127,6 @@ vec3 PhotonTracer::trace_ray(Ray & r, Scene * s, unsigned int rec_level) const {
} else if (_f->m_mat->m_rho > 0.0f && rec_level >= m_max_depth) } else if (_f->m_mat->m_rho > 0.0f && rec_level >= m_max_depth)
return vec3(0.0f); return vec3(0.0f);
// TODO: Change photon map search method for hemisphere search.
mn = Vec3(i_pos.x - m_h_radius, i_pos.y - m_h_radius, i_pos.z - m_h_radius);
mx = Vec3(i_pos.x + m_h_radius, i_pos.y + m_h_radius, i_pos.z + m_h_radius);
photons = m_photon_map.findInRange(mn, mx);
for (vector<Photon>::iterator it = photons.begin(); it != photons.end(); it++) {
(*it).getColor(red, green, blue);
p_contrib += vec3(red, green, blue);
}
p_contrib /= pi<float>() * (m_h_radius * m_h_radius);
color += p_contrib;
} else { } else {
// If the material has transmission enabled, calculate the Fresnel term. // If the material has transmission enabled, calculate the Fresnel term.
kr = fresnel(r.m_direction, n, r.m_ref_index, _f->m_mat->m_ref_index); kr = fresnel(r.m_direction, n, r.m_ref_index, _f->m_mat->m_ref_index);
@@ -153,8 +158,9 @@ void PhotonTracer::build_photon_map(Scene * s, const size_t n_photons_per_ligth,
Light * l; Light * l;
AreaLight * al; AreaLight * al;
vec3 l_sample, s_normal, h_sample; vec3 l_sample, s_normal, h_sample;
Vec3 ls, dir;
float r1, r2; float r1, r2;
Ray rr; Photon ph;
size_t total = 0, current = 0; size_t total = 0, current = 0;
for (vector<Light *>::iterator it = s->m_lights.begin(); it != s->m_lights.end(); it++) { for (vector<Light *>::iterator it = s->m_lights.begin(); it != s->m_lights.end(); it++) {
@@ -173,7 +179,7 @@ void PhotonTracer::build_photon_map(Scene * s, const size_t n_photons_per_ligth,
al = static_cast<AreaLight *>(l); al = static_cast<AreaLight *>(l);
#pragma omp parallel for schedule(dynamic, 1) private(l_sample, s_normal, h_sample, r1, r2, rr) shared(current) #pragma omp parallel for schedule(dynamic, 1) private(l_sample, s_normal, h_sample, r1, r2) shared(current)
for (size_t p = 0; p < n_photons_per_ligth; p++) { for (size_t p = 0; p < n_photons_per_ligth; p++) {
if (!specular) { if (!specular) {
l_sample = al->sample_at_surface(); l_sample = al->sample_at_surface();
@@ -181,15 +187,22 @@ void PhotonTracer::build_photon_map(Scene * s, const size_t n_photons_per_ligth,
r1 = random01(); r1 = random01();
r2 = random01(); r2 = random01();
h_sample = sample_hemisphere(r1, r2); h_sample = normalize(sample_hemisphere(r1, r2));
rotate_sample(h_sample, s_normal); rotate_sample(h_sample, s_normal);
rr = Ray(normalize(h_sample), l_sample + (h_sample * BIAS)); ls = Vec3(l_sample.x, l_sample.y, l_sample.z);
dir = Vec3(h_sample.x, h_sample.y, h_sample.z);
ph = Photon(ls, dir, al->m_figure->m_mat->m_emission.r, al->m_figure->m_mat->m_emission.g, al->m_figure->m_mat->m_emission.b);
} else { } else {
// TODO: Generate photon from light source in direction of specular reflective objects. // TODO: Generate photon from light source in direction of specular reflective objects.
} }
trace_photon(rr, s, 0, specular); #pragma omp critical
{
m_photon_map.addPhoton(ph);
}
trace_photon(ph, s, 0);
#pragma omp atomic #pragma omp atomic
current++; current++;
@@ -203,21 +216,20 @@ void PhotonTracer::build_photon_map(Scene * s, const size_t n_photons_per_ligth,
m_photon_map.buildKdTree(); m_photon_map.buildKdTree();
} }
vec3 PhotonTracer::trace_photon(Ray &r, Scene * s, const unsigned int rec_level, const bool specular) { void PhotonTracer::trace_photon(Photon & ph, Scene * s, const unsigned int rec_level) {
Photon photon; Photon photon;
float t, _t; float t, _t, red, green, blue;
Figure * _f; Figure * _f;
vec3 n, color, i_pos, ref, sample, dir_diff_color, dir_spec_color, ind_color, amb_color; vec3 n, color, i_pos, sample, ph_dir, ph_pos;
Vec3 p_pos; Vec3 p_pos, p_dir;
Ray mv_r, sr, rr; Ray r;
bool vis, is_area_light = false;
float kr, r1, r2; float kr, r1, r2;
AreaLight * al;
t = numeric_limits<float>::max(); t = numeric_limits<float>::max();
_f = NULL; _f = NULL;
// Find the closest intersecting surface. // Find the closest intersecting surface.
r = Ray(ph.direction.x, ph.direction.y, ph.direction.z, ph.position.x, ph.position.y, ph.position.z);
for (size_t f = 0; f < s->m_figures.size(); f++) { for (size_t f = 0; f < s->m_figures.size(); f++) {
if (s->m_figures[f]->intersect(r, _t) && _t < t) { if (s->m_figures[f]->intersect(r, _t) && _t < t) {
t = _t; t = _t;
@@ -231,142 +243,65 @@ vec3 PhotonTracer::trace_photon(Ray &r, Scene * s, const unsigned int rec_level,
i_pos = r.m_origin + (t * r.m_direction); i_pos = r.m_origin + (t * r.m_direction);
n = _f->normal_at_int(r, t); n = _f->normal_at_int(r, t);
is_area_light = false; if (!_f->m_mat->m_refract && rec_level < m_max_depth){
// Check if the object is an area light;
for (vector<Light *>::iterator it = s->m_lights.begin(); it != s->m_lights.end(); it++) {
if ((*it)->light_type() == Light::AREA && static_cast<AreaLight *>(*it)->m_figure == _f)
is_area_light = true;
}
// If the object is an area light, return it's emission value.
if (is_area_light) {
p_pos = Vec3(i_pos.x, i_pos.y, i_pos.z);
photon = Photon(p_pos, _f->m_mat->m_emission.r, _f->m_mat->m_emission.g, _f->m_mat->m_emission.b);
#pragma omp critical
{
m_photon_map.addPhoton(photon);
}
return _f->m_mat->m_emission;
// Check if the material is not reflective/refractive.
} else if (!_f->m_mat->m_refract) {
// Calculate the direct lighting.
for (size_t l = 0; l < s->m_lights.size(); l++) {
// For every light source
vis = true;
if (s->m_lights[l]->light_type() == Light::INFINITESIMAL) {
// Cast a shadow ray to determine visibility.
sr = Ray(s->m_lights[l]->direction(i_pos), i_pos + (n * BIAS));
for (size_t f = 0; f < s->m_figures.size(); f++) {
if (s->m_figures[f]->intersect(sr, _t) && _t < s->m_lights[l]->distance(i_pos)) {
vis = false;
break;
}
}
// Evaluate the shading model accounting for visibility.
dir_diff_color += vis ? s->m_lights[l]->diffuse(n, r, i_pos, *_f->m_mat) : vec3(0.0f);
dir_spec_color += vis ? s->m_lights[l]->specular(n, r, i_pos, *_f->m_mat) : vec3(0.0f);
} else if (s->m_lights[l]->light_type() == Light::AREA) {
// Cast a shadow ray towards a sample point on the surface of the light source.
al = static_cast<AreaLight *>(s->m_lights[l]);
al->sample_at_surface();
sr = Ray(al->direction(i_pos), i_pos + (n * BIAS));
for (size_t f = 0; f < s->m_figures.size(); f++) {
// Avoid self-intersection with the light source.
if (al->m_figure != s->m_figures[f]) {
if (s->m_figures[f]->intersect(sr, _t) && _t < al->distance(i_pos)) {
vis = false;
break;
}
}
}
// Evaluate the shading model accounting for visibility.
dir_diff_color += vis ? s->m_lights[l]->diffuse(n, r, i_pos, *_f->m_mat) : vec3(0.0f);
dir_spec_color += vis ? s->m_lights[l]->specular(n, r, i_pos, *_f->m_mat) : vec3(0.0f);
}
}
// Calculate indirect lighting contribution.
if (rec_level < m_max_depth) { if (rec_level < m_max_depth) {
r1 = random01(); r1 = random01();
r2 = random01(); r2 = random01();
sample = sample_hemisphere(r1, r2); sample = sample_hemisphere(r1, r2);
rotate_sample(sample, n); rotate_sample(sample, n);
rr = Ray(normalize(sample), i_pos + (sample * BIAS)); normalize(sample);
ind_color += r1 * trace_ray(rr, s, rec_level + 1) / PDF; } else
sample = vec3(0.0f);
ph.getColor(red, green, blue);
color = (1.0f - _f->m_mat->m_rho) * (vec3(red, green, blue) * (_f->m_mat->m_diffuse / pi<float>()));
p_pos = Vec3(i_pos.x, i_pos.y, i_pos.z);
p_dir = Vec3(sample.x, sample.y, sample.z);
photon = Photon(p_pos, p_dir, color.r, color.g, color.b);
#pragma omp critical
{
m_photon_map.addPhoton(photon);
} }
// Calculate environment light contribution trace_photon(photon, s, rec_level + 1);
vis = true; }
r1 = random01(); // Determine the specular reflection color.
r2 = random01(); if (!_f->m_mat->m_refract && _f->m_mat->m_rho > 0.0f && rec_level < m_max_depth) {
sample = sample_hemisphere(r1, r2); color = (_f->m_mat->m_rho) * vec3(red, green, blue);
rotate_sample(sample, n); i_pos += n * BIAS;
rr = Ray(normalize(sample), i_pos + (sample * BIAS)); p_pos = Vec3(i_pos.x, i_pos.y, i_pos.z);
ph_dir = normalize(reflect(vec3(ph.direction.x, ph.direction.y, ph.direction.z), n));
p_dir = Vec3(ph_dir.x, ph_dir.y, ph_dir.z);
photon = Photon(p_pos, p_dir, color.r, color.g, color.b);
trace_photon(photon, s, rec_level + 1);
// Cast a shadow ray to determine visibility. } else if (_f->m_mat->m_refract && rec_level >= m_max_depth) {
for (size_t f = 0; f < s->m_figures.size(); f++) {
if (s->m_figures[f]->intersect(rr, _t)) {
vis = false;
break;
}
}
amb_color = vis ? s->m_env->get_color(rr) * max(dot(n, rr.m_direction), 0.0f) / PDF : vec3(0.0f);
// Add lighting.
color += ((dir_diff_color + ind_color + amb_color) * (_f->m_mat->m_diffuse / pi<float>())) + (_f->m_mat->m_specular * dir_spec_color);
//if (specular) {
// Determine the specular reflection color.
if (_f->m_mat->m_rho > 0.0f && rec_level < m_max_depth) {
rr = Ray(normalize(reflect(r.m_direction, n)), i_pos + n * BIAS);
color += _f->m_mat->m_rho * trace_ray(rr, s, rec_level + 1);
} else if (_f->m_mat->m_rho > 0.0f && rec_level >= m_max_depth)
return vec3(0.0f);
//}
} else {
// If the material has transmission enabled, calculate the Fresnel term. // If the material has transmission enabled, calculate the Fresnel term.
kr = fresnel(r.m_direction, n, r.m_ref_index, _f->m_mat->m_ref_index); kr = fresnel(r.m_direction, n, r.m_ref_index, _f->m_mat->m_ref_index);
// Determine the specular reflection color. // Determine the specular reflection color.
if (kr > 0.0f && rec_level < m_max_depth) { if (kr > 0.0f && rec_level < m_max_depth) {
rr = Ray(normalize(reflect(r.m_direction, n)), i_pos + n * BIAS); color = kr * vec3(red, green, blue);
color += kr * trace_ray(rr, s, rec_level + 1); i_pos += n * BIAS;
} else if (rec_level >= m_max_depth) p_pos = Vec3(i_pos.x, i_pos.y, i_pos.z);
return vec3(0.0f); ph_dir = normalize(reflect(vec3(ph.direction.x, ph.direction.y, ph.direction.z), n));
p_dir = Vec3(ph_dir.x, ph_dir.y, ph_dir.z);
photon = Photon(p_pos, p_dir, color.r, color.g, color.b);
trace_photon(photon, s, rec_level + 1);
}
// Determine the transmission color. // Determine the transmission color.
if (_f->m_mat->m_refract && kr < 1.0f && rec_level < m_max_depth) { if (_f->m_mat->m_refract && kr < 1.0f && rec_level < m_max_depth) {
rr = Ray(normalize(refract(r.m_direction, n, r.m_ref_index / _f->m_mat->m_ref_index)), i_pos - n * BIAS, _f->m_mat->m_ref_index); color = (1.0f - kr) * vec3(red, green, blue);
color += (1.0f - kr) * trace_ray(rr, s, rec_level + 1); i_pos -= n * (2 * BIAS);
} else if (rec_level >= m_max_depth) p_pos = Vec3(i_pos.x, i_pos.y, i_pos.z);
return vec3(0.0f); ph_dir = normalize(refract(vec3(ph.direction.x, ph.direction.y, ph.direction.z), n, ph.ref_index / _f->m_mat->m_ref_index));
p_dir = Vec3(ph_dir.x, ph_dir.y, ph_dir.z);
photon = Photon(p_pos, p_dir, color.r, color.g, color.b, _f->m_mat->m_ref_index);
trace_photon(photon, s, rec_level + 1);
}
} }
}
color += _f->m_mat->m_emission;
p_pos = Vec3(i_pos.x, i_pos.y, i_pos.z);
photon = Photon(p_pos, color.r, color.g, color.b);
#pragma omp critical
{
m_photon_map.addPhoton(photon);
}
// Return final color.
return color;
} else
return s->m_env->get_color(r);
} }

2
photon_tracer.hpp
View File

@@ -18,7 +18,7 @@ public:
private: private:
float m_h_radius; float m_h_radius;
kdTree m_photon_map; kdTree m_photon_map;
vec3 trace_photon(Ray &r, Scene * s, const unsigned int rec_level, const bool specular = false); void trace_photon(Photon & ph, Scene * s, const unsigned int rec_level);
}; };
#endif #endif

2
whitted_tracer.cpp
View File

@@ -88,7 +88,7 @@ vec3 WhittedTracer::trace_ray(Ray & r, Scene * s, unsigned int rec_level) const
dir_spec_color += vis ? s->m_lights[l]->specular(n, r, i_pos, *_f->m_mat) : vec3(0.0f); dir_spec_color += vis ? s->m_lights[l]->specular(n, r, i_pos, *_f->m_mat) : vec3(0.0f);
} }
color += (dir_diff_color * (_f->m_mat->m_diffuse / pi<float>())) + (_f->m_mat->m_specular * dir_spec_color); color += (1.0f - _f->m_mat->m_rho) * ((dir_diff_color * (_f->m_mat->m_diffuse / pi<float>())) + (_f->m_mat->m_specular * dir_spec_color));
// Determine the specular reflection color. // Determine the specular reflection color.
if (_f->m_mat->m_rho > 0.0f && rec_level < m_max_depth) { if (_f->m_mat->m_rho > 0.0f && rec_level < m_max_depth) {