+ double t = fabs(src0[i]) + fabs(src1[i]) + fabs(src0[i] * src1[i]) + 1.0;
+ if (!float_near_abs_eps(cdst[i], odst[i], t * 2 * FLT_EPSILON)) {
+ fprintf(stderr, "%d: %- .12f - %- .12f = % .12g\n",
+ i, cdst[i], odst[i], cdst[i] - odst[i]);
+ fail();
+ break;
+ }
+ }
+ bench_new(odst, src0, src1, LEN);
+}
+
+static void test_vector_dmul(const double *src0, const double *src1)
+{
+ LOCAL_ALIGNED_32(double, cdst, [LEN]);
+ LOCAL_ALIGNED_32(double, odst, [LEN]);
+ int i;
+
+ declare_func(void, double *dst, const double *src0, const double *src1,
+ int len);
+
+ call_ref(cdst, src0, src1, LEN);
+ call_new(odst, src0, src1, LEN);
+ for (i = 0; i < LEN; i++) {
+ double t = fabs(src0[i]) + fabs(src1[i]) + fabs(src0[i] * src1[i]) + 1.0;
+ if (!double_near_abs_eps(cdst[i], odst[i], t * 2 * DBL_EPSILON)) {