52 #define	test(func, x, y, z, result, exceptmask, excepts) do {		\  argument
53 	volatile long double _vx = (x), _vy = (y), _vz = (z);		\
57 	    #func, #x);							\
60 #define	testall(x, y, z, result, exceptmask, excepts)	do {		\  argument
61 	test(fma, (double)(x), (double)(y), (double)(z),		\
63 	test(fmaf, (float)(x), (float)(y), (float)(z),			\
65 	test(fmal, (x), (y), (z), (result), (exceptmask), (excepts));	\
69 #define	testrnd(func, x, y, z, rn, ru, rd, rz, exceptmask, excepts)	do { \  argument
71 	test((func), (x), (y), (z), (rn), (exceptmask), (excepts));	\
73 	test((func), (x), (y), (z), (ru), (exceptmask), (excepts));	\
75 	test((func), (x), (y), (z), (rd), (exceptmask), (excepts));	\
77 	test((func), (x), (y), (z), (rz), (exceptmask), (excepts));	\
166 	/* x*y should not raise an inexact/overflow/underflow if z is NaN. */  in test_nans()
177  * Tests for cases where z is very small compared to x*y.
182 	/* x*y positive, z positive */  in test_small_z()
195 	/* x*y negative, z negative */  in test_small_z()
208 	/* x*y positive, z negative */  in test_small_z()
221 	/* x*y negative, z positive */  in test_small_z()
236  * Tests for cases where z is very large compared to x*y.
241 	/* z positive, x*y positive */  in test_big_z()
254 	/* z negative, x*y negative */  in test_big_z()
267 	/* z negative, x*y positive */  in test_big_z()
280 	/* z positive, x*y negative */  in test_big_z()
298 	/* ilogb(x*y) - ilogb(z) = 20 */  in test_accuracy()
299 	testrnd(fmaf, -0x1.c139d8p-51, -0x1.600e7ap32, 0x1.26558cp-38,  in test_accuracy()
303 		0x1.26558cac31580p-38, 0x1.34e48a78aae97p-18,  in test_accuracy()
309 		0x1.26558cac315807eb07e448042101p-38L,  in test_accuracy()
317 		0x1.26558cac315807ecp-38L, 0x1.34e48a78aae96c78p-18L,  in test_accuracy()
322 		0x1.26558cac31580p-38L, 0x1.34e48a78aae97p-18L,  in test_accuracy()
327 	/* ilogb(x*y) - ilogb(z) = -40 */  in test_accuracy()
356 	/* ilogb(x*y) - ilogb(z) = 0 */  in test_accuracy()
385 	/* x*y (rounded) ~= -z */  in test_accuracy()
544 	 * - Cancellation tests (e.g., z = (double)x*y, but x*y is inexact)  in ATF_TP_ADD_TCS()