Lines Matching refs:dest
32 fp_fabs(struct fp_ext *dest, struct fp_ext *src)
36 fp_monadic_check(dest, src);
38 dest->sign = 0;
40 return dest;
44 fp_fneg(struct fp_ext *dest, struct fp_ext *src)
48 fp_monadic_check(dest, src);
50 dest->sign = !dest->sign;
52 return dest;
61 fp_fadd(struct fp_ext *dest, struct fp_ext *src)
67 fp_dyadic_check(dest, src);
69 if (IS_INF(dest)) {
71 if (IS_INF(src) && (src->sign != dest->sign))
72 fp_set_nan(dest);
73 return dest;
76 fp_copy_ext(dest, src);
77 return dest;
80 if (IS_ZERO(dest)) {
82 if (src->sign != dest->sign) {
84 dest->sign = 1;
86 dest->sign = 0;
89 fp_copy_ext(dest, src);
90 return dest;
93 dest->lowmant = src->lowmant = 0;
95 if ((diff = dest->exp - src->exp) > 0)
98 fp_denormalize(dest, diff);
100 if (dest->sign == src->sign) {
101 if (fp_addmant(dest, src))
102 if (!fp_addcarry(dest))
103 return dest;
105 if (dest->mant.m64 < src->mant.m64) {
106 fp_submant(dest, src, dest);
107 dest->sign = !dest->sign;
109 fp_submant(dest, dest, src);
112 return dest;
121 fp_fsub(struct fp_ext *dest, struct fp_ext *src)
126 return fp_fadd(dest, src);
131 fp_fcmp(struct fp_ext *dest, struct fp_ext *src)
135 FPDATA->temp[1] = *dest;
141 fp_ftst(struct fp_ext *dest, struct fp_ext *src)
145 (void)dest;
151 fp_fmul(struct fp_ext *dest, struct fp_ext *src)
158 fp_dyadic_check(dest, src);
161 dest->sign = src->sign ^ dest->sign;
164 if (IS_INF(dest)) {
166 fp_set_nan(dest);
167 return dest;
170 if (IS_ZERO(dest))
171 fp_set_nan(dest);
173 fp_copy_ext(dest, src);
174 return dest;
180 if (IS_ZERO(dest) || IS_ZERO(src)) {
181 dest->exp = 0;
182 dest->mant.m64 = 0;
183 dest->lowmant = 0;
185 return dest;
188 exp = dest->exp + src->exp - 0x3ffe;
193 if ((long)dest->mant.m32[0] >= 0)
194 exp -= fp_overnormalize(dest);
199 fp_multiplymant(&temp, dest, src);
205 fp_putmant128(dest, &temp, 1);
207 fp_putmant128(dest, &temp, 0);
210 fp_set_ovrflw(dest);
211 return dest;
213 dest->exp = exp;
216 fp_denormalize(dest, -exp);
219 return dest;
226 of src / dest, the result is actually dest / src. */
229 fp_fdiv(struct fp_ext *dest, struct fp_ext *src)
236 fp_dyadic_check(dest, src);
239 dest->sign = src->sign ^ dest->sign;
242 if (IS_INF(dest)) {
245 fp_set_nan(dest);
247 return dest;
251 dest->exp = 0;
252 dest->mant.m64 = 0;
253 dest->lowmant = 0;
255 return dest;
259 if (IS_ZERO(dest)) {
262 fp_set_nan(dest);
264 return dest;
269 dest->exp = 0x7fff;
270 dest->mant.m64 = 0;
272 return dest;
275 exp = dest->exp - src->exp + 0x3fff;
280 if ((long)dest->mant.m32[0] >= 0)
281 exp -= fp_overnormalize(dest);
286 fp_dividemant(&temp, dest, src);
292 fp_putmant128(dest, &temp, 32);
294 fp_putmant128(dest, &temp, 31);
297 fp_set_ovrflw(dest);
298 return dest;
300 dest->exp = exp;
303 fp_denormalize(dest, -exp);
306 return dest;
310 fp_fsglmul(struct fp_ext *dest, struct fp_ext *src)
316 fp_dyadic_check(dest, src);
319 dest->sign = src->sign ^ dest->sign;
322 if (IS_INF(dest)) {
324 fp_set_nan(dest);
325 return dest;
328 if (IS_ZERO(dest))
329 fp_set_nan(dest);
331 fp_copy_ext(dest, src);
332 return dest;
338 if (IS_ZERO(dest) || IS_ZERO(src)) {
339 dest->exp = 0;
340 dest->mant.m64 = 0;
341 dest->lowmant = 0;
343 return dest;
346 exp = dest->exp + src->exp - 0x3ffe;
349 fp_mul64(dest->mant.m32[0], dest->mant.m32[1],
350 dest->mant.m32[0] & 0xffffff00,
354 fp_set_ovrflw(dest);
355 return dest;
357 dest->exp = exp;
360 fp_denormalize(dest, -exp);
363 return dest;
367 fp_fsgldiv(struct fp_ext *dest, struct fp_ext *src)
374 fp_dyadic_check(dest, src);
377 dest->sign = src->sign ^ dest->sign;
380 if (IS_INF(dest)) {
383 fp_set_nan(dest);
385 return dest;
389 dest->exp = 0;
390 dest->mant.m64 = 0;
391 dest->lowmant = 0;
393 return dest;
397 if (IS_ZERO(dest)) {
400 fp_set_nan(dest);
402 return dest;
407 dest->exp = 0x7fff;
408 dest->mant.m64 = 0;
410 return dest;
413 exp = dest->exp - src->exp + 0x3fff;
415 dest->mant.m32[0] &= 0xffffff00;
419 if (dest->mant.m32[0] >= src->mant.m32[0]) {
420 fp_sub64(dest->mant, src->mant);
421 fp_div64(quot, rem, dest->mant.m32[0], 0, src->mant.m32[0]);
422 dest->mant.m32[0] = 0x80000000 | (quot >> 1);
423 dest->mant.m32[1] = (quot & 1) | rem; /* only for rounding */
425 fp_div64(quot, rem, dest->mant.m32[0], 0, src->mant.m32[0]);
426 dest->mant.m32[0] = quot;
427 dest->mant.m32[1] = rem; /* only for rounding */
432 fp_set_ovrflw(dest);
433 return dest;
435 dest->exp = exp;
438 fp_denormalize(dest, -exp);
441 return dest;
450 static void fp_roundint(struct fp_ext *dest, int mode)
455 if (!fp_normalize_ext(dest))
459 if (IS_INF(dest) || IS_ZERO(dest))
463 oldmant = dest->mant;
464 switch (dest->exp) {
466 dest->mant.m64 = 0;
469 dest->mant.m32[0] &= 0xffffffffU << (0x401e - dest->exp);
470 dest->mant.m32[1] = 0;
471 if (oldmant.m64 == dest->mant.m64)
475 dest->mant.m32[1] &= 0xffffffffU << (0x403e - dest->exp);
476 if (oldmant.m32[1] == dest->mant.m32[1])
498 switch (dest->exp) {
510 mask = 1 << (0x401d - dest->exp);
515 if (!(oldmant.m32[0] << (dest->exp - 0x3ffd)) &&
528 mask = 1 << (0x403d - dest->exp);
533 if (!(oldmant.m32[1] << (dest->exp - 0x401d)))
543 if (dest->sign ^ (mode - FPCR_ROUND_RM))
548 switch (dest->exp) {
550 dest->exp = 0x3fff;
551 dest->mant.m64 = 1ULL << 63;
554 mask = 1 << (0x401e - dest->exp);
555 if (dest->mant.m32[0] += mask)
557 dest->mant.m32[0] = 0x80000000;
558 dest->exp++;
561 mask = 1 << (0x403e - dest->exp);
562 if (dest->mant.m32[1] += mask)
564 if (dest->mant.m32[0] += 1)
566 dest->mant.m32[0] = 0x80000000;
567 dest->exp++;
577 modrem_kernel(struct fp_ext *dest, struct fp_ext *src, int mode)
581 fp_dyadic_check(dest, src);
584 if (IS_INF(dest) || IS_ZERO(src)) {
585 fp_set_nan(dest);
586 return dest;
588 if (IS_ZERO(dest) || IS_INF(src))
589 return dest;
592 fp_copy_ext(&tmp, dest);
596 fp_fsub(dest, &tmp);
599 fp_set_quotient((dest->mant.m64 & 0x7f) | (dest->sign << 7));
600 return dest;
608 fmod(src,dest) = (dest - (src * floor(dest / src))) */
611 fp_fmod(struct fp_ext *dest, struct fp_ext *src)
614 return modrem_kernel(dest, src, FPCR_ROUND_RZ);
619 frem(src,dest) = (dest - (src * round(dest / src)))
623 fp_frem(struct fp_ext *dest, struct fp_ext *src)
626 return modrem_kernel(dest, src, FPCR_ROUND_RN);
630 fp_fint(struct fp_ext *dest, struct fp_ext *src)
634 fp_copy_ext(dest, src);
636 fp_roundint(dest, FPDATA->rnd);
638 return dest;
642 fp_fintrz(struct fp_ext *dest, struct fp_ext *src)
646 fp_copy_ext(dest, src);
648 fp_roundint(dest, FPCR_ROUND_RZ);
650 return dest;
654 fp_fscale(struct fp_ext *dest, struct fp_ext *src)
660 fp_dyadic_check(dest, src);
664 fp_set_nan(dest);
665 return dest;
667 if (IS_INF(dest))
668 return dest;
671 if (IS_ZERO(src) || IS_ZERO(dest))
672 return dest;
676 fp_set_ovrflw(dest);
677 return dest;
687 scale += dest->exp;
690 fp_set_ovrflw(dest);
693 fp_denormalize(dest, -scale);
695 dest->exp = scale;
697 return dest;