Lines Matching refs:dest
34 fp_fabs(struct fp_ext *dest, struct fp_ext *src) in fp_fabs() argument
38 fp_monadic_check(dest, src); in fp_fabs()
40 dest->sign = 0; in fp_fabs()
42 return dest; in fp_fabs()
46 fp_fneg(struct fp_ext *dest, struct fp_ext *src) in fp_fneg() argument
50 fp_monadic_check(dest, src); in fp_fneg()
52 dest->sign = !dest->sign; in fp_fneg()
54 return dest; in fp_fneg()
63 fp_fadd(struct fp_ext *dest, struct fp_ext *src) in fp_fadd() argument
69 fp_dyadic_check(dest, src); in fp_fadd()
71 if (IS_INF(dest)) { in fp_fadd()
73 if (IS_INF(src) && (src->sign != dest->sign)) in fp_fadd()
74 fp_set_nan(dest); in fp_fadd()
75 return dest; in fp_fadd()
78 fp_copy_ext(dest, src); in fp_fadd()
79 return dest; in fp_fadd()
82 if (IS_ZERO(dest)) { in fp_fadd()
84 if (src->sign != dest->sign) { in fp_fadd()
86 dest->sign = 1; in fp_fadd()
88 dest->sign = 0; in fp_fadd()
91 fp_copy_ext(dest, src); in fp_fadd()
92 return dest; in fp_fadd()
95 dest->lowmant = src->lowmant = 0; in fp_fadd()
97 if ((diff = dest->exp - src->exp) > 0) in fp_fadd()
100 fp_denormalize(dest, diff); in fp_fadd()
102 if (dest->sign == src->sign) { in fp_fadd()
103 if (fp_addmant(dest, src)) in fp_fadd()
104 if (!fp_addcarry(dest)) in fp_fadd()
105 return dest; in fp_fadd()
107 if (dest->mant.m64 < src->mant.m64) { in fp_fadd()
108 fp_submant(dest, src, dest); in fp_fadd()
109 dest->sign = !dest->sign; in fp_fadd()
111 fp_submant(dest, dest, src); in fp_fadd()
114 return dest; in fp_fadd()
123 fp_fsub(struct fp_ext *dest, struct fp_ext *src) in fp_fsub() argument
128 return fp_fadd(dest, src); in fp_fsub()
133 fp_fcmp(struct fp_ext *dest, struct fp_ext *src) in fp_fcmp() argument
137 FPDATA->temp[1] = *dest; in fp_fcmp()
143 fp_ftst(struct fp_ext *dest, struct fp_ext *src) in fp_ftst() argument
147 (void)dest; in fp_ftst()
153 fp_fmul(struct fp_ext *dest, struct fp_ext *src) in fp_fmul() argument
160 fp_dyadic_check(dest, src); in fp_fmul()
163 dest->sign = src->sign ^ dest->sign; in fp_fmul()
166 if (IS_INF(dest)) { in fp_fmul()
168 fp_set_nan(dest); in fp_fmul()
169 return dest; in fp_fmul()
172 if (IS_ZERO(dest)) in fp_fmul()
173 fp_set_nan(dest); in fp_fmul()
175 fp_copy_ext(dest, src); in fp_fmul()
176 return dest; in fp_fmul()
182 if (IS_ZERO(dest) || IS_ZERO(src)) { in fp_fmul()
183 dest->exp = 0; in fp_fmul()
184 dest->mant.m64 = 0; in fp_fmul()
185 dest->lowmant = 0; in fp_fmul()
187 return dest; in fp_fmul()
190 exp = dest->exp + src->exp - 0x3ffe; in fp_fmul()
195 if ((long)dest->mant.m32[0] >= 0) in fp_fmul()
196 exp -= fp_overnormalize(dest); in fp_fmul()
201 fp_multiplymant(&temp, dest, src); in fp_fmul()
207 fp_putmant128(dest, &temp, 1); in fp_fmul()
209 fp_putmant128(dest, &temp, 0); in fp_fmul()
212 fp_set_ovrflw(dest); in fp_fmul()
213 return dest; in fp_fmul()
215 dest->exp = exp; in fp_fmul()
218 fp_denormalize(dest, -exp); in fp_fmul()
221 return dest; in fp_fmul()
231 fp_fdiv(struct fp_ext *dest, struct fp_ext *src) in fp_fdiv() argument
238 fp_dyadic_check(dest, src); in fp_fdiv()
241 dest->sign = src->sign ^ dest->sign; in fp_fdiv()
244 if (IS_INF(dest)) { in fp_fdiv()
247 fp_set_nan(dest); in fp_fdiv()
249 return dest; in fp_fdiv()
253 dest->exp = 0; in fp_fdiv()
254 dest->mant.m64 = 0; in fp_fdiv()
255 dest->lowmant = 0; in fp_fdiv()
257 return dest; in fp_fdiv()
261 if (IS_ZERO(dest)) { in fp_fdiv()
264 fp_set_nan(dest); in fp_fdiv()
266 return dest; in fp_fdiv()
271 dest->exp = 0x7fff; in fp_fdiv()
272 dest->mant.m64 = 0; in fp_fdiv()
274 return dest; in fp_fdiv()
277 exp = dest->exp - src->exp + 0x3fff; in fp_fdiv()
282 if ((long)dest->mant.m32[0] >= 0) in fp_fdiv()
283 exp -= fp_overnormalize(dest); in fp_fdiv()
288 fp_dividemant(&temp, dest, src); in fp_fdiv()
294 fp_putmant128(dest, &temp, 32); in fp_fdiv()
296 fp_putmant128(dest, &temp, 31); in fp_fdiv()
299 fp_set_ovrflw(dest); in fp_fdiv()
300 return dest; in fp_fdiv()
302 dest->exp = exp; in fp_fdiv()
305 fp_denormalize(dest, -exp); in fp_fdiv()
308 return dest; in fp_fdiv()
312 fp_fsglmul(struct fp_ext *dest, struct fp_ext *src) in fp_fsglmul() argument
318 fp_dyadic_check(dest, src); in fp_fsglmul()
321 dest->sign = src->sign ^ dest->sign; in fp_fsglmul()
324 if (IS_INF(dest)) { in fp_fsglmul()
326 fp_set_nan(dest); in fp_fsglmul()
327 return dest; in fp_fsglmul()
330 if (IS_ZERO(dest)) in fp_fsglmul()
331 fp_set_nan(dest); in fp_fsglmul()
333 fp_copy_ext(dest, src); in fp_fsglmul()
334 return dest; in fp_fsglmul()
340 if (IS_ZERO(dest) || IS_ZERO(src)) { in fp_fsglmul()
341 dest->exp = 0; in fp_fsglmul()
342 dest->mant.m64 = 0; in fp_fsglmul()
343 dest->lowmant = 0; in fp_fsglmul()
345 return dest; in fp_fsglmul()
348 exp = dest->exp + src->exp - 0x3ffe; in fp_fsglmul()
351 fp_mul64(dest->mant.m32[0], dest->mant.m32[1], in fp_fsglmul()
352 dest->mant.m32[0] & 0xffffff00, in fp_fsglmul()
356 fp_set_ovrflw(dest); in fp_fsglmul()
357 return dest; in fp_fsglmul()
359 dest->exp = exp; in fp_fsglmul()
362 fp_denormalize(dest, -exp); in fp_fsglmul()
365 return dest; in fp_fsglmul()
369 fp_fsgldiv(struct fp_ext *dest, struct fp_ext *src) in fp_fsgldiv() argument
376 fp_dyadic_check(dest, src); in fp_fsgldiv()
379 dest->sign = src->sign ^ dest->sign; in fp_fsgldiv()
382 if (IS_INF(dest)) { in fp_fsgldiv()
385 fp_set_nan(dest); in fp_fsgldiv()
387 return dest; in fp_fsgldiv()
391 dest->exp = 0; in fp_fsgldiv()
392 dest->mant.m64 = 0; in fp_fsgldiv()
393 dest->lowmant = 0; in fp_fsgldiv()
395 return dest; in fp_fsgldiv()
399 if (IS_ZERO(dest)) { in fp_fsgldiv()
402 fp_set_nan(dest); in fp_fsgldiv()
404 return dest; in fp_fsgldiv()
409 dest->exp = 0x7fff; in fp_fsgldiv()
410 dest->mant.m64 = 0; in fp_fsgldiv()
412 return dest; in fp_fsgldiv()
415 exp = dest->exp - src->exp + 0x3fff; in fp_fsgldiv()
417 dest->mant.m32[0] &= 0xffffff00; in fp_fsgldiv()
421 if (dest->mant.m32[0] >= src->mant.m32[0]) { in fp_fsgldiv()
422 fp_sub64(dest->mant, src->mant); in fp_fsgldiv()
423 fp_div64(quot, rem, dest->mant.m32[0], 0, src->mant.m32[0]); in fp_fsgldiv()
424 dest->mant.m32[0] = 0x80000000 | (quot >> 1); in fp_fsgldiv()
425 dest->mant.m32[1] = (quot & 1) | rem; /* only for rounding */ in fp_fsgldiv()
427 fp_div64(quot, rem, dest->mant.m32[0], 0, src->mant.m32[0]); in fp_fsgldiv()
428 dest->mant.m32[0] = quot; in fp_fsgldiv()
429 dest->mant.m32[1] = rem; /* only for rounding */ in fp_fsgldiv()
434 fp_set_ovrflw(dest); in fp_fsgldiv()
435 return dest; in fp_fsgldiv()
437 dest->exp = exp; in fp_fsgldiv()
440 fp_denormalize(dest, -exp); in fp_fsgldiv()
443 return dest; in fp_fsgldiv()
452 static void fp_roundint(struct fp_ext *dest, int mode) in fp_roundint() argument
457 if (!fp_normalize_ext(dest)) in fp_roundint()
461 if (IS_INF(dest) || IS_ZERO(dest)) in fp_roundint()
465 oldmant = dest->mant; in fp_roundint()
466 switch (dest->exp) { in fp_roundint()
468 dest->mant.m64 = 0; in fp_roundint()
471 dest->mant.m32[0] &= 0xffffffffU << (0x401e - dest->exp); in fp_roundint()
472 dest->mant.m32[1] = 0; in fp_roundint()
473 if (oldmant.m64 == dest->mant.m64) in fp_roundint()
477 dest->mant.m32[1] &= 0xffffffffU << (0x403e - dest->exp); in fp_roundint()
478 if (oldmant.m32[1] == dest->mant.m32[1]) in fp_roundint()
500 switch (dest->exp) { in fp_roundint()
512 mask = 1 << (0x401d - dest->exp); in fp_roundint()
517 if (!(oldmant.m32[0] << (dest->exp - 0x3ffd)) && in fp_roundint()
530 mask = 1 << (0x403d - dest->exp); in fp_roundint()
535 if (!(oldmant.m32[1] << (dest->exp - 0x401d))) in fp_roundint()
545 if (dest->sign ^ (mode - FPCR_ROUND_RM)) in fp_roundint()
550 switch (dest->exp) { in fp_roundint()
552 dest->exp = 0x3fff; in fp_roundint()
553 dest->mant.m64 = 1ULL << 63; in fp_roundint()
556 mask = 1 << (0x401e - dest->exp); in fp_roundint()
557 if (dest->mant.m32[0] += mask) in fp_roundint()
559 dest->mant.m32[0] = 0x80000000; in fp_roundint()
560 dest->exp++; in fp_roundint()
563 mask = 1 << (0x403e - dest->exp); in fp_roundint()
564 if (dest->mant.m32[1] += mask) in fp_roundint()
566 if (dest->mant.m32[0] += 1) in fp_roundint()
568 dest->mant.m32[0] = 0x80000000; in fp_roundint()
569 dest->exp++; in fp_roundint()
579 modrem_kernel(struct fp_ext *dest, struct fp_ext *src, int mode) in modrem_kernel() argument
583 fp_dyadic_check(dest, src); in modrem_kernel()
586 if (IS_INF(dest) || IS_ZERO(src)) { in modrem_kernel()
587 fp_set_nan(dest); in modrem_kernel()
588 return dest; in modrem_kernel()
590 if (IS_ZERO(dest) || IS_INF(src)) in modrem_kernel()
591 return dest; in modrem_kernel()
594 fp_copy_ext(&tmp, dest); in modrem_kernel()
598 fp_fsub(dest, &tmp); in modrem_kernel()
601 fp_set_quotient((dest->mant.m64 & 0x7f) | (dest->sign << 7)); in modrem_kernel()
602 return dest; in modrem_kernel()
613 fp_fmod(struct fp_ext *dest, struct fp_ext *src) in fp_fmod() argument
616 return modrem_kernel(dest, src, FPCR_ROUND_RZ); in fp_fmod()
625 fp_frem(struct fp_ext *dest, struct fp_ext *src) in fp_frem() argument
628 return modrem_kernel(dest, src, FPCR_ROUND_RN); in fp_frem()
632 fp_fint(struct fp_ext *dest, struct fp_ext *src) in fp_fint() argument
636 fp_copy_ext(dest, src); in fp_fint()
638 fp_roundint(dest, FPDATA->rnd); in fp_fint()
640 return dest; in fp_fint()
644 fp_fintrz(struct fp_ext *dest, struct fp_ext *src) in fp_fintrz() argument
648 fp_copy_ext(dest, src); in fp_fintrz()
650 fp_roundint(dest, FPCR_ROUND_RZ); in fp_fintrz()
652 return dest; in fp_fintrz()
656 fp_fscale(struct fp_ext *dest, struct fp_ext *src) in fp_fscale() argument
662 fp_dyadic_check(dest, src); in fp_fscale()
666 fp_set_nan(dest); in fp_fscale()
667 return dest; in fp_fscale()
669 if (IS_INF(dest)) in fp_fscale()
670 return dest; in fp_fscale()
673 if (IS_ZERO(src) || IS_ZERO(dest)) in fp_fscale()
674 return dest; in fp_fscale()
678 fp_set_ovrflw(dest); in fp_fscale()
679 return dest; in fp_fscale()
689 scale += dest->exp; in fp_fscale()
692 fp_set_ovrflw(dest); in fp_fscale()
695 fp_denormalize(dest, -scale); in fp_fscale()
697 dest->exp = scale; in fp_fscale()
699 return dest; in fp_fscale()