arm/vfp/vfpdouble.c

43 	.exponent	= 2047,
50 	pr_debug("VFP: %s: sign=%d exponent=%d significand=%016llx\n",
51 		 str, d->sign != 0, d->exponent, d->significand);
63 		vd->exponent -= bits - 1;
73 	int exponent, shift, underflow;
81 	if (vd->exponent == 2047 && (vd->significand == 0 || exceptions))
88 		vd->exponent = 0;
92 	exponent = vd->exponent;
99 		exponent -= shift;
104 	vd->exponent = exponent;
112 	underflow = exponent < 0;
114 		significand = vfp_shiftright64jamming(significand, -exponent);
115 		exponent = 0;
117 		vd->exponent = exponent;
146 		exponent += 1;
150 		vd->exponent = exponent;
171 	if (exponent >= 2046) {
174 			vd->exponent = 2045;
177 			vd->exponent = 2047;		/* infinity */
182 			exponent = 0;
183 		if (exponent || significand > 0x8000000000000000ULL)
187 		vd->exponent = exponent;
317 	vdd.exponent = ((vdm.exponent - 1023) >> 1) + 1023;
318 	vdd.significand = (u64)vfp_estimate_sqrt_significand(vdm.exponent, vdm.significand >> 32) << 31;
322 	vdm.significand >>= 1 + (vdm.exponent & 1);
463 	 * If we have an infinity or a NaN, the exponent must be 255
466 		vsd.exponent = 255;
471 		vsd.exponent = 0;
473 		vsd.exponent = vdm.exponent - (1023 - 127);
488 	vdm.exponent = 1023 + 63 - 1;
500 	vdm.exponent = 1023 + 63 - 1;
525 	if (vdm.exponent >= 1023 + 32) {
528 	} else if (vdm.exponent >= 1023 - 1) {
529 		int shift = 1023 + 63 - vdm.exponent;
562 		if (vdm.exponent | vdm.significand) {
605 	} else if (vdm.exponent >= 1023 + 32) {
610 	} else if (vdm.exponent >= 1023 - 1) {
611 		int shift = 1023 + 63 - vdm.exponent;	/* 58 */
639 		if (vdm.exponent | vdm.significand) {
743 	if (vdn->exponent < vdm->exponent) {
753 	if (vdn->exponent == 2047)
766 	exp_diff = vdn->exponent - vdm->exponent;
801 	if (vdn->exponent < vdm->exponent) {
813 	if (vdn->exponent == 2047) {
814 		if (vdn->significand || (vdm->exponent == 2047 && vdm->significand))
816 		if ((vdm->exponent | vdm->significand) == 0) {
820 		vdd->exponent = vdn->exponent;
829 	if ((vdm->exponent | vdm->significand) == 0) {
830 		vdd->exponent = 0;
836 	 * We add 2 to the destination exponent for the same reason
840 	vdd->exponent = vdn->exponent + vdm->exponent - 1023 + 2;
857 	if (vdn.exponent == 0 && vdn.significand)
861 	if (vdm.exponent == 0 && vdm.significand)
869 	if (vdn.exponent == 0 && vdn.significand)
924 	if (vdn.exponent == 0 && vdn.significand)
928 	if (vdm.exponent == 0 && vdm.significand)
944 	if (vdn.exponent == 0 && vdn.significand)
948 	if (vdm.exponent == 0 && vdm.significand)
966 	if (vdn.exponent == 0 && vdn.significand)
970 	if (vdm.exponent == 0 && vdm.significand)
987 	if (vdn.exponent == 0 && vdn.significand)
991 	if (vdm.exponent == 0 && vdm.significand)
1066 	vdd.exponent = vdn.exponent - vdm.exponent + 1023 - 1;
1070 		vdd.exponent++;
1096 	vdd.exponent = 0;
1103 	vdd.exponent = 2047;