Method unitvector

library/modules/Random.cpp:127–146

template<class T>
void MersenneRNG::unitvector(T *p, int size)
{
    for (;;)
    {
        T rsqr = 0;
        for (int i = 0; i < size; i++)
        {
            p[i] = (T)unitrandom();
            rsqr += p[i]*p[i];
        }

        if (rsqr > 0 && rsqr <= 1)
        {
            rsqr = std::sqrt(rsqr);
            for (int i = 0; i < size; i++)
                p[i] /= rsqr;
            break;
        }
    }
}

template void DFHACK_EXPORT MersenneRNG::unitvector<float>(float *p, int size);
template void DFHACK_EXPORT MersenneRNG::unitvector<double>(double *p, int size);