uint64_t ulps(vec<T, N> x, vec<T, N> y)
{
uint64_t u = 0;
for (size_t i = 0; i < N; i++)
{
u = std::max(u, ulps(x[i], y[i]));
}
return u;
}
NT2_TEST_UNIT_DECL double max_ulps(float a, float b ) { return ulps(a,b); }
NT2_TEST_UNIT_DECL double max_ulps(double a, double b) { return ulps(a,b); }