gcclibs/libdecnumber/decNumber.c

994   decContext workset;		/* working context */
1083       workset = *set;		/* copy all fields */
1085       workset.digits = reqdigits + (inrhs->digits + inrhs->exponent) + 1;
1087       if (workset.digits > DECNUMMAXP)
1093       /* workset.digits is the count of digits for the accumulator we need */
1096 	sizeof (decNumber) + (D2U (workset.digits) - 1) * sizeof (Unit);
1136 	      decDivideOp (dac, &dnOne, lhs, &workset, DIVIDE, &status);
1180 	      decMultiplyOp (dac, dac, lhs, &workset, &status);	/* dac=dac*x */
1186 	  decMultiplyOp (dac, dac, dac, &workset, &status);	/* dac=dac*dac [square] */
1207 	  /* round subnormals [to set.digits rather than workset.digits] */
1219 	  decDivideOp (dac, &dnOne, dac, &workset, DIVIDE, &status);
1457   decContext workset, approxset;	/* work contexts */
1583       decContextDefault (&workset, DEC_INIT_DECIMAL64);
1591       workset.digits = set->digits;	/* p for initial calculation */
1642       decMultiplyOp (a, a, f, &workset, &ignore);	/* a=a*f */
1643       decAddOp (a, a, t, &workset, 0, &ignore);	/* ..+t */
1652       workset.digits = 3;	/* initial p */
1656 	  workset.digits = workset.digits * 2 - 2;
1657 	  if (workset.digits > maxp)
1658 	    workset.digits = maxp;
1661 	  decDivideOp (b, f, a, &workset, DIVIDE, &ignore);	/* b=f/a */
1662 	  decAddOp (b, b, a, &workset, 0, &ignore);	/* b=b+a */
1663 	  decMultiplyOp (a, b, t, &workset, &ignore);	/* a=b*0.5 */
1666 	  if (workset.digits == maxp)
1673       /* Here workset.digits=maxp and t=0.5 */
1674       workset.digits--;		/* maxp-1 is OK now */
1677       decAddOp (b, b, t, &workset, DECNEG, &ignore);	/* b = a - 0.5 ulp */
1678       workset.round = DEC_ROUND_UP;
1679       decMultiplyOp (b, b, b, &workset, &ignore);	/* b = mulru(b, b) */
1680       decCompareOp (b, f, b, &workset, COMPARE, &ignore);	/* b ? f, reversed */
1686 	  decAddOp (a, a, t, &workset, DECNEG, &ignore);	/* a = a - 1 ulp */
1693 	  decAddOp (b, b, t, &workset, 0, &ignore);	/* b = a + 0.5 ulp */
1694 	  workset.round = DEC_ROUND_DOWN;
1695 	  decMultiplyOp (b, b, b, &workset, &ignore);	/* b = mulrd(b, b) */
1696 	  decCompareOp (b, b, f, &workset, COMPARE, &ignore);	/* b ? f */
1701 	      decAddOp (a, a, t, &workset, 0, &ignore);	/* a = a + 1 ulp */
1730 	  decMultiplyOp (b, b, b, &workset, &mstatus);	/* try the multiply */
1737 	      decCompareOp (t, b, rhs, &workset, COMPARE, &mstatus);	/* b ? rhs */
1828   decContext workset;		/* working context */
1852   workset = *set;		/* clone rounding, etc. */
1853   workset.digits = rhs->digits;	/* no length rounding */
1854   workset.traps = 0;		/* no traps */
1856   return decNumberQuantize (res, rhs, &dn, &workset);
3719 	      decContext workset;	/* work */
3720 	      workset = *set;	/* clone rounding, etc. */
3721 	      workset.digits = lhs->digits - adjust;	/* set requested length */
3723 	      decCopyFit (res, lhs, &workset, &residue, status);	/* fit to result */
3724 	      decApplyRound (res, &workset, residue, status);	/* .. and round */
5269   decContext workset;		/* work */
5327   workset = *set;		/* clone rounding, etc. */
5328   workset.digits = dn->digits - adjust;	/* set requested length */
5329   workset.emin -= adjust;	/* and adjust emin to match */
5331   decSetCoeff (dn, &workset, dn->lsu, dn->digits, residue, status);
5332   decApplyRound (dn, &workset, *residue, status);