sys/altivec_ast/vecast.c

142 	volatile vector float in1, in2, in3;
147 	/* in1 + in2 */
148 	in1 = VF(10, 0x10p-140, 0x20p-130, -0x2000p-134);
151 	ck_equal("vec_add", vec_add(in1, in2), ans);
153 	/* in1 - in2 */
154 	in1 = VF(0x4000p-134, 10, 0x10p-140,   0x3p-130);
157 	ck_equal("vec_sub", vec_sub(in1, in2), ans);
159 	/* in1 * in2 + in3 */
160 	in1 = VF( 0x6p-70,   0x6p-140, 6,   0x6p-100);
164 	ck_equal("vec_madd", vec_madd(in1, in2, in3), ans);
166 	/* in3 - in1 * in2 */
167 	in1 = VF( 0xbp-30,    0xbp-70,   0xbp44,  11);
171 	ck_equal("vec_nmsub", vec_nmsub(in1, in2, in3), ans);
173 	/* 1 / in1 */
174 	in1 = VF(      3,       0x3p126,       0x3p-126, 0x1p127);
176 	ck_estimate("vec_re", vec_re(in1), dan, REL_1_IN, 4096);
178 	/* 1 / sqrt(in1) */
179 	in1 = VF(1,       2,        0x1p-128,        0x5p-135);
181 	ck_estimate("vec_rsqrt", vec_rsqrte(in1), dan, REL_1_IN, 4096);
183 	/* log2(in1) */
184 	in1 = VF(0x1p-130, 0x1p-149, 32, 0x1p-10);
186 	ck_estimate("vec_loge", vec_loge(in1), dan, ABS_1_IN, 32);
187 	in1 = VF(     0x123p-139,       0xabcp-145,  1, 1);
189 	ck_estimate("vec_loge", vec_loge(in1), dan, ABS_1_IN, 32);
191 	/* 2**in1 */
192 	in1 = VF(    -149,     -138,     -127,   10);
194 	ck_equal("vec_expte", vec_expte(in1), ans);
195 	in1 = VF(    -10,      -145.3,       -136.9,       -127.1);
197 	ck_estimate("vec_expte", vec_expte(in1), dan, REL_1_IN, 16);
199 	/* (int)(in1 * 2**exponent) */
200 	in1 = VF(0x1p-127, 2.34, -0xfedp-140, -19.8);
202 	ck_equal_i("vec_cts", vec_cts(in1, 0), ian);
203 	in1 = VF(0x1p-113, -1, -0xabcp-143, 0x1fp-10);
205 	ck_equal_i("vec_cts", vec_cts(in1, 10), ian);
207 	/* (unsigned)(in1 * 2**exponent) */
208 	in1 = VF(0x1.ap-130, 0x1.ep-140, 24000012, 0);
210 	ck_equal_u("vec_ctu", vec_ctu(in1, 7), uan);