Fixed ticket #147 refract for scalar

This commit is contained in:
Christophe Riccio 2011-11-14 18:13:18 +00:00
parent 390498ef87
commit 369e929910
2 changed files with 28 additions and 2 deletions

View File

@ -283,6 +283,25 @@ namespace glm
} }
// refract // refract
template <typename genType>
GLM_FUNC_QUALIFIER genType refract
(
genType const & I,
genType const & N,
genType const & eta
)
{
//It could be a vector
//GLM_STATIC_ASSERT(detail::type<genType>::is_float);
genType dotValue = dot(N, I);
genType k = genType(1) - eta * eta * (genType(1) - dotValue * dotValue);
if(k < genType(0))
return genType(0);
else
return eta * I - (eta * dotValue + sqrt(k)) * N;
}
template <typename genType> template <typename genType>
GLM_FUNC_QUALIFIER genType refract GLM_FUNC_QUALIFIER genType refract
( (

View File

@ -34,17 +34,24 @@ int test_refract()
{ {
int Error = 0; int Error = 0;
{
float A(1.0f);
float B(1.0f);
float C = glm::refract(A, B, 0.5f);
Error += C == 1.0f ? 0 : 1;
}
{ {
glm::vec2 A(1.0f, 0.0f); glm::vec2 A(1.0f, 0.0f);
glm::vec2 B(0.0f, 1.0f); glm::vec2 B(0.0f, 1.0f);
glm::vec2 C = glm::reflect(A, B); glm::vec2 C = glm::refract(A, B, 0.5f);
Error += C == glm::vec2(-1.0, 0.0) ? 0 : 1; Error += C == glm::vec2(-1.0, 0.0) ? 0 : 1;
} }
{ {
glm::dvec2 A(1.0f, 0.0f); glm::dvec2 A(1.0f, 0.0f);
glm::dvec2 B(0.0f, 1.0f); glm::dvec2 B(0.0f, 1.0f);
glm::dvec2 C = glm::reflect(A, B); glm::dvec2 C = glm::refract(A, B, 0.5);
Error += C == glm::dvec2(-1.0, 0.0) ? 0 : 1; Error += C == glm::dvec2(-1.0, 0.0) ? 0 : 1;
} }