Commit ad0523a8 authored by Laurent Belcour's avatar Laurent Belcour
Browse files

Updating the retro-fresnel.

parent 5b06ca04
......@@ -26,7 +26,7 @@ vec retro_schlick::value(const vec& x) const
params::convert(&x[0], input_parametrization(), params::SCHLICK_VK, xp);
params::convert(&x[0], input_parametrization(), params::CARTESIAN, cart);
const double dotRK = -xp[0]*cart[0] -xp[1]*cart[1] + xp[2]*cart[2];
const double dotRK = xp[0]*cart[0] + xp[1]*cart[1] + xp[2]*cart[2];
vec res(_nY);
for(int i=0; i<_nY; ++i)
......@@ -66,7 +66,7 @@ vec retro_schlick::parametersJacobian(const vec& x) const
params::convert(&x[0], input_parametrization(), params::SCHLICK_VK, xp);
params::convert(&x[0], input_parametrization(), params::CARTESIAN, cart);
const double dotRK = -xp[0]*cart[0] -xp[1]*cart[1] + xp[2]*cart[2];
const double dotRK = xp[0]*cart[0] + xp[1]*cart[1] + xp[2]*cart[2];
vec jac(nbParameters()*nY);
for(int i=0; i<nY; ++i)
......
......@@ -20,7 +20,7 @@ int main(int argc, char** argv)
arguments args(argc, argv) ;
#ifdef __linux__
feenableexcept(FE_DIVBYZERO | FE_OVERFLOW | FE_INVALID);
//feenableexcept(FE_DIVBYZERO | FE_OVERFLOW | FE_INVALID);
#endif
fitter* fit = plugins_manager::get_fitter(args["fitter"]) ;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment