12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394 |
- /*
- * linux/arch/arm/vfp/vfpdouble.c
- *
- * This code is derived in part from John R. Housers softfloat library, which
- * carries the following notice:
- *
- * ===========================================================================
- * This C source file is part of the SoftFloat IEC/IEEE Floating-point
- * Arithmetic Package, Release 2.
- *
- * Written by John R. Hauser. This work was made possible in part by the
- * International Computer Science Institute, located at Suite 600, 1947 Center
- * Street, Berkeley, California 94704. Funding was partially provided by the
- * National Science Foundation under grant MIP-9311980. The original version
- * of this code was written as part of a project to build a fixed-point vector
- * processor in collaboration with the University of California at Berkeley,
- * overseen by Profs. Nelson Morgan and John Wawrzynek. More information
- * is available through the web page `http://HTTP.CS.Berkeley.EDU/~jhauser/
- * arithmetic/softfloat.html'.
- *
- * THIS SOFTWARE IS DISTRIBUTED AS IS, FOR FREE. Although reasonable effort
- * has been made to avoid it, THIS SOFTWARE MAY CONTAIN FAULTS THAT WILL AT
- * TIMES RESULT IN INCORRECT BEHAVIOR. USE OF THIS SOFTWARE IS RESTRICTED TO
- * PERSONS AND ORGANIZATIONS WHO CAN AND WILL TAKE FULL RESPONSIBILITY FOR ANY
- * AND ALL LOSSES, COSTS, OR OTHER PROBLEMS ARISING FROM ITS USE.
- *
- * Derivative works are acceptable, even for commercial purposes, so long as
- * (1) they include prominent notice that the work is derivative, and (2) they
- * include prominent notice akin to these three paragraphs for those parts of
- * this code that are retained.
- * ===========================================================================
- */
- #include <linux/kernel.h>
- #include <linux/bitops.h>
- #include <asm/div64.h>
- #include <asm/vfp.h>
- #include "vfpinstr.h"
- #include "vfp.h"
- static struct vfp_double vfp_double_default_qnan = {
- .exponent = 2047,
- .sign = 0,
- .significand = VFP_DOUBLE_SIGNIFICAND_QNAN,
- };
- static void vfp_double_dump(const char *str, struct vfp_double *d)
- {
- pr_debug("VFP: %s: sign=%d exponent=%d significand=%016llx\n",
- str, d->sign != 0, d->exponent, d->significand);
- }
- static void vfp_double_normalise_denormal(struct vfp_double *vd)
- {
- int bits = 31 - fls(vd->significand >> 32);
- if (bits == 31)
- bits = 63 - fls(vd->significand);
- vfp_double_dump("normalise_denormal: in", vd);
- if (bits) {
- vd->exponent -= bits - 1;
- vd->significand <<= bits;
- }
- vfp_double_dump("normalise_denormal: out", vd);
- }
- u32 vfp_double_normaliseround(int dd, struct vfp_double *vd, u32 fpscr, u32 exceptions, const char *func)
- {
- u64 significand, incr;
- int exponent, shift, underflow;
- u32 rmode;
- vfp_double_dump("pack: in", vd);
- /*
- * Infinities and NaNs are a special case.
- */
- if (vd->exponent == 2047 && (vd->significand == 0 || exceptions))
- goto pack;
- /*
- * Special-case zero.
- */
- if (vd->significand == 0) {
- vd->exponent = 0;
- goto pack;
- }
- exponent = vd->exponent;
- significand = vd->significand;
|