Function greatestPrimeFactor

src/backend/common/dispatch.hpp:37–48

template<typename T>
inline T greatestPrimeFactor(T n) {
    T v{2};

    while (v <= n) {
        if (n % v == 0 && isPrime(v))
            n /= v;
        else
            v += 1;
    }

    return v;
}
// Empty columns (dim==1) in refDims are removed from dims & strides.
// INPUT: refDims, refNdims
// UPDATE: dims, strides