Lines Matching refs:dest
31 struct fp_ext *fp_fabs(struct fp_ext *dest, struct fp_ext *src) in fp_fabs() argument
35 fp_monadic_check(dest, src); in fp_fabs()
37 dest->sign = 0; in fp_fabs()
39 return dest; in fp_fabs()
42 struct fp_ext *fp_fneg(struct fp_ext *dest, struct fp_ext *src) in fp_fneg() argument
46 fp_monadic_check(dest, src); in fp_fneg()
48 dest->sign = !dest->sign; in fp_fneg()
50 return dest; in fp_fneg()
58 struct fp_ext *fp_fadd(struct fp_ext *dest, struct fp_ext *src) in fp_fadd() argument
64 fp_dyadic_check(dest, src); in fp_fadd()
66 if (IS_INF(dest)) { in fp_fadd()
68 if (IS_INF(src) && (src->sign != dest->sign)) in fp_fadd()
69 fp_set_nan(dest); in fp_fadd()
70 return dest; in fp_fadd()
73 fp_copy_ext(dest, src); in fp_fadd()
74 return dest; in fp_fadd()
77 if (IS_ZERO(dest)) { in fp_fadd()
79 if (src->sign != dest->sign) { in fp_fadd()
81 dest->sign = 1; in fp_fadd()
83 dest->sign = 0; in fp_fadd()
86 fp_copy_ext(dest, src); in fp_fadd()
87 return dest; in fp_fadd()
90 dest->lowmant = src->lowmant = 0; in fp_fadd()
92 if ((diff = dest->exp - src->exp) > 0) in fp_fadd()
95 fp_denormalize(dest, diff); in fp_fadd()
97 if (dest->sign == src->sign) { in fp_fadd()
98 if (fp_addmant(dest, src)) in fp_fadd()
99 if (!fp_addcarry(dest)) in fp_fadd()
100 return dest; in fp_fadd()
102 if (dest->mant.m64 < src->mant.m64) { in fp_fadd()
103 fp_submant(dest, src, dest); in fp_fadd()
104 dest->sign = !dest->sign; in fp_fadd()
106 fp_submant(dest, dest, src); in fp_fadd()
109 return dest; in fp_fadd()
117 struct fp_ext *fp_fsub(struct fp_ext *dest, struct fp_ext *src) in fp_fsub() argument
122 return fp_fadd(dest, src); in fp_fsub()
126 struct fp_ext *fp_fcmp(struct fp_ext *dest, struct fp_ext *src) in fp_fcmp() argument
130 FPDATA->temp[1] = *dest; in fp_fcmp()
135 struct fp_ext *fp_ftst(struct fp_ext *dest, struct fp_ext *src) in fp_ftst() argument
139 (void)dest; in fp_ftst()
144 struct fp_ext *fp_fmul(struct fp_ext *dest, struct fp_ext *src) in fp_fmul() argument
151 fp_dyadic_check(dest, src); in fp_fmul()
154 dest->sign = src->sign ^ dest->sign; in fp_fmul()
157 if (IS_INF(dest)) { in fp_fmul()
159 fp_set_nan(dest); in fp_fmul()
160 return dest; in fp_fmul()
163 if (IS_ZERO(dest)) in fp_fmul()
164 fp_set_nan(dest); in fp_fmul()
166 fp_copy_ext(dest, src); in fp_fmul()
167 return dest; in fp_fmul()
173 if (IS_ZERO(dest) || IS_ZERO(src)) { in fp_fmul()
174 dest->exp = 0; in fp_fmul()
175 dest->mant.m64 = 0; in fp_fmul()
176 dest->lowmant = 0; in fp_fmul()
178 return dest; in fp_fmul()
181 exp = dest->exp + src->exp - 0x3ffe; in fp_fmul()
186 if ((long)dest->mant.m32[0] >= 0) in fp_fmul()
187 exp -= fp_overnormalize(dest); in fp_fmul()
192 fp_multiplymant(&temp, dest, src); in fp_fmul()
198 fp_putmant128(dest, &temp, 1); in fp_fmul()
200 fp_putmant128(dest, &temp, 0); in fp_fmul()
203 fp_set_ovrflw(dest); in fp_fmul()
204 return dest; in fp_fmul()
206 dest->exp = exp; in fp_fmul()
209 fp_denormalize(dest, -exp); in fp_fmul()
212 return dest; in fp_fmul()
221 struct fp_ext *fp_fdiv(struct fp_ext *dest, struct fp_ext *src) in fp_fdiv() argument
228 fp_dyadic_check(dest, src); in fp_fdiv()
231 dest->sign = src->sign ^ dest->sign; in fp_fdiv()
234 if (IS_INF(dest)) { in fp_fdiv()
237 fp_set_nan(dest); in fp_fdiv()
239 return dest; in fp_fdiv()
243 dest->exp = 0; in fp_fdiv()
244 dest->mant.m64 = 0; in fp_fdiv()
245 dest->lowmant = 0; in fp_fdiv()
247 return dest; in fp_fdiv()
251 if (IS_ZERO(dest)) { in fp_fdiv()
254 fp_set_nan(dest); in fp_fdiv()
256 return dest; in fp_fdiv()
261 dest->exp = 0x7fff; in fp_fdiv()
262 dest->mant.m64 = 0; in fp_fdiv()
264 return dest; in fp_fdiv()
267 exp = dest->exp - src->exp + 0x3fff; in fp_fdiv()
272 if ((long)dest->mant.m32[0] >= 0) in fp_fdiv()
273 exp -= fp_overnormalize(dest); in fp_fdiv()
278 fp_dividemant(&temp, dest, src); in fp_fdiv()
284 fp_putmant128(dest, &temp, 32); in fp_fdiv()
286 fp_putmant128(dest, &temp, 31); in fp_fdiv()
289 fp_set_ovrflw(dest); in fp_fdiv()
290 return dest; in fp_fdiv()
292 dest->exp = exp; in fp_fdiv()
295 fp_denormalize(dest, -exp); in fp_fdiv()
298 return dest; in fp_fdiv()
301 struct fp_ext *fp_fsglmul(struct fp_ext *dest, struct fp_ext *src) in fp_fsglmul() argument
307 fp_dyadic_check(dest, src); in fp_fsglmul()
310 dest->sign = src->sign ^ dest->sign; in fp_fsglmul()
313 if (IS_INF(dest)) { in fp_fsglmul()
315 fp_set_nan(dest); in fp_fsglmul()
316 return dest; in fp_fsglmul()
319 if (IS_ZERO(dest)) in fp_fsglmul()
320 fp_set_nan(dest); in fp_fsglmul()
322 fp_copy_ext(dest, src); in fp_fsglmul()
323 return dest; in fp_fsglmul()
329 if (IS_ZERO(dest) || IS_ZERO(src)) { in fp_fsglmul()
330 dest->exp = 0; in fp_fsglmul()
331 dest->mant.m64 = 0; in fp_fsglmul()
332 dest->lowmant = 0; in fp_fsglmul()
334 return dest; in fp_fsglmul()
337 exp = dest->exp + src->exp - 0x3ffe; in fp_fsglmul()
340 fp_mul64(dest->mant.m32[0], dest->mant.m32[1], in fp_fsglmul()
341 dest->mant.m32[0] & 0xffffff00, in fp_fsglmul()
345 fp_set_ovrflw(dest); in fp_fsglmul()
346 return dest; in fp_fsglmul()
348 dest->exp = exp; in fp_fsglmul()
351 fp_denormalize(dest, -exp); in fp_fsglmul()
354 return dest; in fp_fsglmul()
357 struct fp_ext *fp_fsgldiv(struct fp_ext *dest, struct fp_ext *src) in fp_fsgldiv() argument
364 fp_dyadic_check(dest, src); in fp_fsgldiv()
367 dest->sign = src->sign ^ dest->sign; in fp_fsgldiv()
370 if (IS_INF(dest)) { in fp_fsgldiv()
373 fp_set_nan(dest); in fp_fsgldiv()
375 return dest; in fp_fsgldiv()
379 dest->exp = 0; in fp_fsgldiv()
380 dest->mant.m64 = 0; in fp_fsgldiv()
381 dest->lowmant = 0; in fp_fsgldiv()
383 return dest; in fp_fsgldiv()
387 if (IS_ZERO(dest)) { in fp_fsgldiv()
390 fp_set_nan(dest); in fp_fsgldiv()
392 return dest; in fp_fsgldiv()
397 dest->exp = 0x7fff; in fp_fsgldiv()
398 dest->mant.m64 = 0; in fp_fsgldiv()
400 return dest; in fp_fsgldiv()
403 exp = dest->exp - src->exp + 0x3fff; in fp_fsgldiv()
405 dest->mant.m32[0] &= 0xffffff00; in fp_fsgldiv()
409 if (dest->mant.m32[0] >= src->mant.m32[0]) { in fp_fsgldiv()
410 fp_sub64(dest->mant, src->mant); in fp_fsgldiv()
411 fp_div64(quot, rem, dest->mant.m32[0], 0, src->mant.m32[0]); in fp_fsgldiv()
412 dest->mant.m32[0] = 0x80000000 | (quot >> 1); in fp_fsgldiv()
413 dest->mant.m32[1] = (quot & 1) | rem; /* only for rounding */ in fp_fsgldiv()
415 fp_div64(quot, rem, dest->mant.m32[0], 0, src->mant.m32[0]); in fp_fsgldiv()
416 dest->mant.m32[0] = quot; in fp_fsgldiv()
417 dest->mant.m32[1] = rem; /* only for rounding */ in fp_fsgldiv()
422 fp_set_ovrflw(dest); in fp_fsgldiv()
423 return dest; in fp_fsgldiv()
425 dest->exp = exp; in fp_fsgldiv()
428 fp_denormalize(dest, -exp); in fp_fsgldiv()
431 return dest; in fp_fsgldiv()
440 static void fp_roundint(struct fp_ext *dest, int mode) in fp_roundint() argument
445 if (!fp_normalize_ext(dest)) in fp_roundint()
449 if (IS_INF(dest) || IS_ZERO(dest)) in fp_roundint()
453 oldmant = dest->mant; in fp_roundint()
454 switch (dest->exp) { in fp_roundint()
456 dest->mant.m64 = 0; in fp_roundint()
459 dest->mant.m32[0] &= 0xffffffffU << (0x401e - dest->exp); in fp_roundint()
460 dest->mant.m32[1] = 0; in fp_roundint()
461 if (oldmant.m64 == dest->mant.m64) in fp_roundint()
465 dest->mant.m32[1] &= 0xffffffffU << (0x403e - dest->exp); in fp_roundint()
466 if (oldmant.m32[1] == dest->mant.m32[1]) in fp_roundint()
488 switch (dest->exp) { in fp_roundint()
500 mask = 1 << (0x401d - dest->exp); in fp_roundint()
505 if (!(oldmant.m32[0] << (dest->exp - 0x3ffd)) && in fp_roundint()
518 mask = 1 << (0x403d - dest->exp); in fp_roundint()
523 if (!(oldmant.m32[1] << (dest->exp - 0x401d))) in fp_roundint()
533 if (dest->sign ^ (mode - FPCR_ROUND_RM)) in fp_roundint()
538 switch (dest->exp) { in fp_roundint()
540 dest->exp = 0x3fff; in fp_roundint()
541 dest->mant.m64 = 1ULL << 63; in fp_roundint()
544 mask = 1 << (0x401e - dest->exp); in fp_roundint()
545 if (dest->mant.m32[0] += mask) in fp_roundint()
547 dest->mant.m32[0] = 0x80000000; in fp_roundint()
548 dest->exp++; in fp_roundint()
551 mask = 1 << (0x403e - dest->exp); in fp_roundint()
552 if (dest->mant.m32[1] += mask) in fp_roundint()
554 if (dest->mant.m32[0] += 1) in fp_roundint()
556 dest->mant.m32[0] = 0x80000000; in fp_roundint()
557 dest->exp++; in fp_roundint()
566 static struct fp_ext *modrem_kernel(struct fp_ext *dest, struct fp_ext *src, in modrem_kernel() argument
571 fp_dyadic_check(dest, src); in modrem_kernel()
574 if (IS_INF(dest) || IS_ZERO(src)) { in modrem_kernel()
575 fp_set_nan(dest); in modrem_kernel()
576 return dest; in modrem_kernel()
578 if (IS_ZERO(dest) || IS_INF(src)) in modrem_kernel()
579 return dest; in modrem_kernel()
582 fp_copy_ext(&tmp, dest); in modrem_kernel()
586 fp_fsub(dest, &tmp); in modrem_kernel()
589 fp_set_quotient((dest->mant.m64 & 0x7f) | (dest->sign << 7)); in modrem_kernel()
590 return dest; in modrem_kernel()
600 struct fp_ext *fp_fmod(struct fp_ext *dest, struct fp_ext *src) in fp_fmod() argument
603 return modrem_kernel(dest, src, FPCR_ROUND_RZ); in fp_fmod()
611 struct fp_ext *fp_frem(struct fp_ext *dest, struct fp_ext *src) in fp_frem() argument
614 return modrem_kernel(dest, src, FPCR_ROUND_RN); in fp_frem()
617 struct fp_ext *fp_fint(struct fp_ext *dest, struct fp_ext *src) in fp_fint() argument
621 fp_copy_ext(dest, src); in fp_fint()
623 fp_roundint(dest, FPDATA->rnd); in fp_fint()
625 return dest; in fp_fint()
628 struct fp_ext *fp_fintrz(struct fp_ext *dest, struct fp_ext *src) in fp_fintrz() argument
632 fp_copy_ext(dest, src); in fp_fintrz()
634 fp_roundint(dest, FPCR_ROUND_RZ); in fp_fintrz()
636 return dest; in fp_fintrz()
639 struct fp_ext *fp_fscale(struct fp_ext *dest, struct fp_ext *src) in fp_fscale() argument
645 fp_dyadic_check(dest, src); in fp_fscale()
649 fp_set_nan(dest); in fp_fscale()
650 return dest; in fp_fscale()
652 if (IS_INF(dest)) in fp_fscale()
653 return dest; in fp_fscale()
656 if (IS_ZERO(src) || IS_ZERO(dest)) in fp_fscale()
657 return dest; in fp_fscale()
661 fp_set_ovrflw(dest); in fp_fscale()
662 return dest; in fp_fscale()
672 scale += dest->exp; in fp_fscale()
675 fp_set_ovrflw(dest); in fp_fscale()
678 fp_denormalize(dest, -scale); in fp_fscale()
680 dest->exp = scale; in fp_fscale()
682 return dest; in fp_fscale()