external/Cornucopia/Fresnel.cpp
void fresnel(const VectorXd &t, VectorXd *s, VectorXd *c)
{
    s->resize(t.size());
    c->resize(t.size());
    for(int i = 0; i < t.size(); ++i)
        fresnel(t[i], &((*s)[i]), &((*c)[i]));
}
//The double precision version is not vectorized because the scalar version is actually faster
The double precision version is not vectorized because the scalar version is actually faster
The double precision version is not vectorized because the scalar version is actually faster
void fresnel(const VectorXd &t, VectorXd *s, VectorXd *c)
{
    s->resize(t.size());
    c->resize(t.size());
    for(int i = 0; i < t.size(); ++i)
        fresnel(t[i], &((*s)[i]), &((*c)[i]));
}
void fresnel ( const VectorXd & t , VectorXd * s , VectorXd * c ) { s -> resize ( t . size ( ) ) ; c -> resize ( t . size ( ) ) ; for ( int i = 0 ; i < t . size ( ) ; ++ i ) fresnel ( t [ i ] , & ( ( * s ) [ i ] ) , & ( ( * c ) [ i ] ) ) ; }
The double precision version is not vectorized because the scalar version is actually faster
The double precision version is not vectorized because the scalar version is actually faster