Function relEps

test/pinverse.cpp:111–119

template<typename T>
double relEps(array in) {
    typedef typename af::dtype_traits<T>::base_type InBaseType;
    double fixed_eps = eps<T>();
    double calc_eps  = std::numeric_limits<InBaseType>::epsilon() *
                      std::max(in.dims(0), in.dims(1)) * af::max<double>(in);
    // Use the fixed values above if calculated error tolerance is unnecessarily
    // too small
    return std::max(fixed_eps, calc_eps);
}

typedef ::testing::Types<float, cfloat, double, cdouble> TestTypes;
TYPED_TEST_SUITE(Pinverse, TestTypes);