Page MenuHome
Differential D13466 Diff 46665 intern/cycles/test/util_avxf_test.h

Changeset View

Standalone View

View Options

intern/cycles/test/util_avxf_test.h

Show All 26 Lines#ifdef __KERNEL_AVX2__
return system_cpu_support_avx2(); return system_cpu_support_avx2();
#else #else
# ifdef __KERNEL_AVX__ # ifdef __KERNEL_AVX__
return system_cpu_support_avx(); return system_cpu_support_avx();
# endif # endif
#endif #endif
} }
#define VALIDATECPU \ #define INIT_AVX_TEST \
if (!validate_cpu_capabilities()) \ if (!validate_cpu_capabilities()) \
return; return; \
\
const avxf avxf_a(0.1f, 0.2f, 0.3f, 0.4f, 0.5f, 0.6f, 0.7f, 0.8f); \
const avxf avxf_b(1.0f, 2.0f, 3.0f, 4.0f, 5.0f, 6.0f, 7.0f, 8.0f); \
const avxf avxf_c(1.1f, 2.2f, 3.3f, 4.4f, 5.5f, 6.6f, 7.7f, 8.8f);
#define compare_vector_scalar(a, b) \ #define compare_vector_scalar(a, b) \
for (size_t index = 0; index < a.size; index++) \ for (size_t index = 0; index < a.size; index++) \
EXPECT_FLOAT_EQ(a[index], b); EXPECT_FLOAT_EQ(a[index], b);
#define compare_vector_vector(a, b) \ #define compare_vector_vector(a, b) \
for (size_t index = 0; index < a.size; index++) \ for (size_t index = 0; index < a.size; index++) \
EXPECT_FLOAT_EQ(a[index], b[index]); EXPECT_FLOAT_EQ(a[index], b[index]);
#define compare_vector_vector_near(a, b, abserror) \ #define compare_vector_vector_near(a, b, abserror) \
for (size_t index = 0; index < a.size; index++) \ for (size_t index = 0; index < a.size; index++) \
EXPECT_NEAR(a[index], b[index], abserror); EXPECT_NEAR(a[index], b[index], abserror);
#define basic_test_vv(a, b, op) \ #define basic_test_vv(a, b, op) \
VALIDATECPU \ INIT_AVX_TEST \
avxf c = a op b; \ avxf c = a op b; \
for (size_t i = 0; i < a.size; i++) \ for (size_t i = 0; i < a.size; i++) \
EXPECT_FLOAT_EQ(c[i], a[i] op b[i]); EXPECT_FLOAT_EQ(c[i], a[i] op b[i]);
/* vector op float tests */ /* vector op float tests */
#define basic_test_vf(a, b, op) \ #define basic_test_vf(a, b, op) \
VALIDATECPU \ INIT_AVX_TEST \
avxf c = a op b; \ avxf c = a op b; \
for (size_t i = 0; i < a.size; i++) \ for (size_t i = 0; i < a.size; i++) \
EXPECT_FLOAT_EQ(c[i], a[i] op b); EXPECT_FLOAT_EQ(c[i], a[i] op b);
static const avxf avxf_a(0.1f, 0.2f, 0.3f, 0.4f, 0.5f, 0.6f, 0.7f, 0.8f);
static const avxf avxf_b(1.0f, 2.0f, 3.0f, 4.0f, 5.0f, 6.0f, 7.0f, 8.0f);
static const avxf avxf_c(1.1f, 2.2f, 3.3f, 4.4f, 5.5f, 6.6f, 7.7f, 8.8f);
static const float float_b = 1.5f; static const float float_b = 1.5f;
TEST(TEST_CATEGORY_NAME, avxf_add_vv){basic_test_vv(avxf_a, avxf_b, +)} TEST(TEST_CATEGORY_NAME, TEST(TEST_CATEGORY_NAME, avxf_add_vv){basic_test_vv(avxf_a, avxf_b, +)} TEST(TEST_CATEGORY_NAME,
avxf_sub_vv){ avxf_sub_vv){
basic_test_vv(avxf_a, avxf_b, -)} TEST(TEST_CATEGORY_NAME, avxf_mul_vv){ basic_test_vv(avxf_a, avxf_b, -)} TEST(TEST_CATEGORY_NAME, avxf_mul_vv){
basic_test_vv(avxf_a, avxf_b, *)} TEST(TEST_CATEGORY_NAME, avxf_div_vv){ basic_test_vv(avxf_a, avxf_b, *)} TEST(TEST_CATEGORY_NAME, avxf_div_vv){
basic_test_vv(avxf_a, avxf_b, /)} TEST(TEST_CATEGORY_NAME, avxf_add_vf){ basic_test_vv(avxf_a, avxf_b, /)} TEST(TEST_CATEGORY_NAME, avxf_add_vf){
basic_test_vf(avxf_a, float_b, +)} TEST(TEST_CATEGORY_NAME, avxf_sub_vf){ basic_test_vf(avxf_a, float_b, +)} TEST(TEST_CATEGORY_NAME, avxf_sub_vf){
basic_test_vf(avxf_a, float_b, -)} TEST(TEST_CATEGORY_NAME, avxf_mul_vf){ basic_test_vf(avxf_a, float_b, -)} TEST(TEST_CATEGORY_NAME, avxf_mul_vf){
basic_test_vf(avxf_a, float_b, *)} TEST(TEST_CATEGORY_NAME, basic_test_vf(avxf_a, float_b, *)} TEST(TEST_CATEGORY_NAME,
avxf_div_vf){basic_test_vf(avxf_a, float_b, /)} avxf_div_vf){basic_test_vf(avxf_a, float_b, /)}
TEST(TEST_CATEGORY_NAME, avxf_ctor) TEST(TEST_CATEGORY_NAME, avxf_ctor)
{ {
VALIDATECPU INIT_AVX_TEST
compare_vector_scalar(avxf(7.0f, 6.0f, 5.0f, 4.0f, 3.0f, 2.0f, 1.0f, 0.0f), compare_vector_scalar(avxf(7.0f, 6.0f, 5.0f, 4.0f, 3.0f, 2.0f, 1.0f, 0.0f),
static_cast<float>(index)); static_cast<float>(index));
compare_vector_scalar(avxf(1.0f), 1.0f); compare_vector_scalar(avxf(1.0f), 1.0f);
compare_vector_vector(avxf(1.0f, 2.0f), avxf(1.0f, 1.0f, 1.0f, 1.0f, 2.0f, 2.0f, 2.0f, 2.0f)); compare_vector_vector(avxf(1.0f, 2.0f), avxf(1.0f, 1.0f, 1.0f, 1.0f, 2.0f, 2.0f, 2.0f, 2.0f));
compare_vector_vector(avxf(1.0f, 2.0f, 3.0f, 4.0f), compare_vector_vector(avxf(1.0f, 2.0f, 3.0f, 4.0f),
avxf(1.0f, 2.0f, 3.0f, 4.0f, 1.0f, 2.0f, 3.0f, 4.0f)); avxf(1.0f, 2.0f, 3.0f, 4.0f, 1.0f, 2.0f, 3.0f, 4.0f));
compare_vector_vector(avxf(make_float3(1.0f, 2.0f, 3.0f)), compare_vector_vector(avxf(make_float3(1.0f, 2.0f, 3.0f)),
avxf(0.0f, 3.0f, 2.0f, 1.0f, 0.0f, 3.0f, 2.0f, 1.0f)); avxf(0.0f, 3.0f, 2.0f, 1.0f, 0.0f, 3.0f, 2.0f, 1.0f));
} }
TEST(TEST_CATEGORY_NAME, avxf_sqrt) TEST(TEST_CATEGORY_NAME, avxf_sqrt)
{ {
VALIDATECPU INIT_AVX_TEST
compare_vector_vector(mm256_sqrt(avxf(1.0f, 4.0f, 9.0f, 16.0f, 25.0f, 36.0f, 49.0f, 64.0f)), compare_vector_vector(mm256_sqrt(avxf(1.0f, 4.0f, 9.0f, 16.0f, 25.0f, 36.0f, 49.0f, 64.0f)),
avxf(1.0f, 2.0f, 3.0f, 4.0f, 5.0f, 6.0f, 7.0f, 8.0f)); avxf(1.0f, 2.0f, 3.0f, 4.0f, 5.0f, 6.0f, 7.0f, 8.0f));
} }
TEST(TEST_CATEGORY_NAME, avxf_min_max) TEST(TEST_CATEGORY_NAME, avxf_min_max)
{ {
VALIDATECPU INIT_AVX_TEST
compare_vector_vector(min(avxf_a, avxf_b), avxf_a); compare_vector_vector(min(avxf_a, avxf_b), avxf_a);
compare_vector_vector(max(avxf_a, avxf_b), avxf_b); compare_vector_vector(max(avxf_a, avxf_b), avxf_b);
} }
TEST(TEST_CATEGORY_NAME, avxf_set_sign) TEST(TEST_CATEGORY_NAME, avxf_set_sign)
{ {
VALIDATECPU INIT_AVX_TEST
avxf res = set_sign_bit<1, 0, 0, 0, 0, 0, 0, 0>(avxf_a); avxf res = set_sign_bit<1, 0, 0, 0, 0, 0, 0, 0>(avxf_a);
compare_vector_vector(res, avxf(0.1f, 0.2f, 0.3f, 0.4f, 0.5f, 0.6f, 0.7f, -0.8f)); compare_vector_vector(res, avxf(0.1f, 0.2f, 0.3f, 0.4f, 0.5f, 0.6f, 0.7f, -0.8f));
} }
TEST(TEST_CATEGORY_NAME, avxf_msub) TEST(TEST_CATEGORY_NAME, avxf_msub)
{ {
VALIDATECPU INIT_AVX_TEST
avxf res = msub(avxf_a, avxf_b, avxf_c); avxf res = msub(avxf_a, avxf_b, avxf_c);
avxf exp = avxf((avxf_a[7] * avxf_b[7]) - avxf_c[7], avxf exp = avxf((avxf_a[7] * avxf_b[7]) - avxf_c[7],
(avxf_a[6] * avxf_b[6]) - avxf_c[6], (avxf_a[6] * avxf_b[6]) - avxf_c[6],
(avxf_a[5] * avxf_b[5]) - avxf_c[5], (avxf_a[5] * avxf_b[5]) - avxf_c[5],
(avxf_a[4] * avxf_b[4]) - avxf_c[4], (avxf_a[4] * avxf_b[4]) - avxf_c[4],
(avxf_a[3] * avxf_b[3]) - avxf_c[3], (avxf_a[3] * avxf_b[3]) - avxf_c[3],
(avxf_a[2] * avxf_b[2]) - avxf_c[2], (avxf_a[2] * avxf_b[2]) - avxf_c[2],
(avxf_a[1] * avxf_b[1]) - avxf_c[1], (avxf_a[1] * avxf_b[1]) - avxf_c[1],
(avxf_a[0] * avxf_b[0]) - avxf_c[0]); (avxf_a[0] * avxf_b[0]) - avxf_c[0]);
compare_vector_vector(res, exp); compare_vector_vector(res, exp);
} }
TEST(TEST_CATEGORY_NAME, avxf_madd) TEST(TEST_CATEGORY_NAME, avxf_madd)
{ {
VALIDATECPU INIT_AVX_TEST
avxf res = madd(avxf_a, avxf_b, avxf_c); avxf res = madd(avxf_a, avxf_b, avxf_c);
avxf exp = avxf((avxf_a[7] * avxf_b[7]) + avxf_c[7], avxf exp = avxf((avxf_a[7] * avxf_b[7]) + avxf_c[7],
(avxf_a[6] * avxf_b[6]) + avxf_c[6], (avxf_a[6] * avxf_b[6]) + avxf_c[6],
(avxf_a[5] * avxf_b[5]) + avxf_c[5], (avxf_a[5] * avxf_b[5]) + avxf_c[5],
(avxf_a[4] * avxf_b[4]) + avxf_c[4], (avxf_a[4] * avxf_b[4]) + avxf_c[4],
(avxf_a[3] * avxf_b[3]) + avxf_c[3], (avxf_a[3] * avxf_b[3]) + avxf_c[3],
(avxf_a[2] * avxf_b[2]) + avxf_c[2], (avxf_a[2] * avxf_b[2]) + avxf_c[2],
(avxf_a[1] * avxf_b[1]) + avxf_c[1], (avxf_a[1] * avxf_b[1]) + avxf_c[1],
(avxf_a[0] * avxf_b[0]) + avxf_c[0]); (avxf_a[0] * avxf_b[0]) + avxf_c[0]);
compare_vector_vector(res, exp); compare_vector_vector(res, exp);
} }
TEST(TEST_CATEGORY_NAME, avxf_nmadd) TEST(TEST_CATEGORY_NAME, avxf_nmadd)
{ {
VALIDATECPU INIT_AVX_TEST
avxf res = nmadd(avxf_a, avxf_b, avxf_c); avxf res = nmadd(avxf_a, avxf_b, avxf_c);
avxf exp = avxf(avxf_c[7] - (avxf_a[7] * avxf_b[7]), avxf exp = avxf(avxf_c[7] - (avxf_a[7] * avxf_b[7]),
avxf_c[6] - (avxf_a[6] * avxf_b[6]), avxf_c[6] - (avxf_a[6] * avxf_b[6]),
avxf_c[5] - (avxf_a[5] * avxf_b[5]), avxf_c[5] - (avxf_a[5] * avxf_b[5]),
avxf_c[4] - (avxf_a[4] * avxf_b[4]), avxf_c[4] - (avxf_a[4] * avxf_b[4]),
avxf_c[3] - (avxf_a[3] * avxf_b[3]), avxf_c[3] - (avxf_a[3] * avxf_b[3]),
avxf_c[2] - (avxf_a[2] * avxf_b[2]), avxf_c[2] - (avxf_a[2] * avxf_b[2]),
avxf_c[1] - (avxf_a[1] * avxf_b[1]), avxf_c[1] - (avxf_a[1] * avxf_b[1]),
avxf_c[0] - (avxf_a[0] * avxf_b[0])); avxf_c[0] - (avxf_a[0] * avxf_b[0]));
compare_vector_vector(res, exp); compare_vector_vector(res, exp);
} }
TEST(TEST_CATEGORY_NAME, avxf_compare) TEST(TEST_CATEGORY_NAME, avxf_compare)
{ {
VALIDATECPU INIT_AVX_TEST
avxf a(0.0f, 1.0f, 2.0f, 3.0f, 4.0f, 5.0f, 6.0f, 7.0f); avxf a(0.0f, 1.0f, 2.0f, 3.0f, 4.0f, 5.0f, 6.0f, 7.0f);
avxf b(7.0f, 6.0f, 5.0f, 4.0f, 3.0f, 2.0f, 1.0f, 0.0f); avxf b(7.0f, 6.0f, 5.0f, 4.0f, 3.0f, 2.0f, 1.0f, 0.0f);
avxb res = a <= b; avxb res = a <= b;
int exp[8] = { int exp[8] = {
a[0] <= b[0] ? -1 : 0, a[0] <= b[0] ? -1 : 0,
a[1] <= b[1] ? -1 : 0, a[1] <= b[1] ? -1 : 0,
a[2] <= b[2] ? -1 : 0, a[2] <= b[2] ? -1 : 0,
a[3] <= b[3] ? -1 : 0, a[3] <= b[3] ? -1 : 0,
a[4] <= b[4] ? -1 : 0, a[4] <= b[4] ? -1 : 0,
a[5] <= b[5] ? -1 : 0, a[5] <= b[5] ? -1 : 0,
a[6] <= b[6] ? -1 : 0, a[6] <= b[6] ? -1 : 0,
a[7] <= b[7] ? -1 : 0, a[7] <= b[7] ? -1 : 0,
}; };
compare_vector_vector(res, exp); compare_vector_vector(res, exp);
} }
TEST(TEST_CATEGORY_NAME, avxf_permute) TEST(TEST_CATEGORY_NAME, avxf_permute)
{ {
VALIDATECPU INIT_AVX_TEST
avxf res = permute<3, 0, 1, 7, 6, 5, 2, 4>(avxf_b); avxf res = permute<3, 0, 1, 7, 6, 5, 2, 4>(avxf_b);
compare_vector_vector(res, avxf(4.0f, 6.0f, 3.0f, 2.0f, 1.0f, 7.0f, 8.0f, 5.0f)); compare_vector_vector(res, avxf(4.0f, 6.0f, 3.0f, 2.0f, 1.0f, 7.0f, 8.0f, 5.0f));
} }
TEST(TEST_CATEGORY_NAME, avxf_blend) TEST(TEST_CATEGORY_NAME, avxf_blend)
{ {
VALIDATECPU INIT_AVX_TEST
avxf res = blend<0, 0, 1, 0, 1, 0, 1, 0>(avxf_a, avxf_b); avxf res = blend<0, 0, 1, 0, 1, 0, 1, 0>(avxf_a, avxf_b);
compare_vector_vector(res, avxf(0.1f, 0.2f, 3.0f, 0.4f, 5.0f, 0.6f, 7.0f, 0.8f)); compare_vector_vector(res, avxf(0.1f, 0.2f, 3.0f, 0.4f, 5.0f, 0.6f, 7.0f, 0.8f));
} }
TEST(TEST_CATEGORY_NAME, avxf_shuffle) TEST(TEST_CATEGORY_NAME, avxf_shuffle)
{ {
VALIDATECPU INIT_AVX_TEST
avxf res = shuffle<0, 1, 2, 3, 1, 3, 2, 0>(avxf_a); avxf res = shuffle<0, 1, 2, 3, 1, 3, 2, 0>(avxf_a);
compare_vector_vector(res, avxf(0.4f, 0.2f, 0.1f, 0.3f, 0.5f, 0.6f, 0.7f, 0.8f)); compare_vector_vector(res, avxf(0.4f, 0.2f, 0.1f, 0.3f, 0.5f, 0.6f, 0.7f, 0.8f));
} }
TEST(TEST_CATEGORY_NAME, avxf_cross) TEST(TEST_CATEGORY_NAME, avxf_cross)
{ {
VALIDATECPU INIT_AVX_TEST
avxf res = cross(avxf_b, avxf_c); avxf res = cross(avxf_b, avxf_c);
compare_vector_vector_near(res, compare_vector_vector_near(res,
avxf(0.0f, avxf(0.0f,
-9.5367432e-07f, -9.5367432e-07f,
0.0f, 0.0f,
4.7683716e-07f, 4.7683716e-07f,
0.0f, 0.0f,
-3.8146973e-06f, -3.8146973e-06f,
3.8146973e-06f, 3.8146973e-06f,
3.8146973e-06f), 3.8146973e-06f),
0.000002000f); 0.000002000f);
} }
TEST(TEST_CATEGORY_NAME, avxf_dot3) TEST(TEST_CATEGORY_NAME, avxf_dot3)
{ {
VALIDATECPU INIT_AVX_TEST
float den, den2; float den, den2;
dot3(avxf_a, avxf_b, den, den2); dot3(avxf_a, avxf_b, den, den2);
EXPECT_FLOAT_EQ(den, 14.9f); EXPECT_FLOAT_EQ(den, 14.9f);
EXPECT_FLOAT_EQ(den2, 2.9f); EXPECT_FLOAT_EQ(den2, 2.9f);
} }
CCL_NAMESPACE_END CCL_NAMESPACE_END