config/libbid/bid_round.c

123   UINT64 Cstar;
152   // Cstar = P128 >> Ex
155   Cstar = P128.w[1] >> shift;
189     if (Cstar & 0x01) {	// Cstar is odd; MP in [EVEN, ODD]
191       Cstar--;	// Cstar is now even
201   // check for rounding overflow, which occurs if Cstar = 10^(q-x)
203   if (Cstar == ten2k64[ind]) {	// if  Cstar = 10^(q-x)
204     Cstar = ten2k64[ind - 1];	// Cstar = 10^(q-x-1)
206   } else {	// 10^33 <= Cstar <= 10^34 - 1
209   *ptr_Cstar = Cstar;
226   UINT128 Cstar;
275   // Cstar = P256 >> Ex
279     Cstar.w[0] = (P256.w[2] >> shift) | (P256.w[3] << (64 - shift));
280     Cstar.w[1] = (P256.w[3] >> shift);
286     Cstar.w[0] = P256.w[3] >> shift;
287     Cstar.w[1] = 0x0ULL;
342     if (Cstar.w[0] & 0x01) {	// Cstar is odd; MP in [EVEN, ODD]
344       Cstar.w[0]--;	// Cstar is now even
345       if (Cstar.w[0] == 0xffffffffffffffffULL) {
346 	Cstar.w[1]--;
357   // check for rounding overflow, which occurs if Cstar = 10^(q-x)
360     if (Cstar.w[1] == 0x0ULL && Cstar.w[0] == ten2k64[ind]) {
361       // if  Cstar = 10^(q-x)
362       Cstar.w[0] = ten2k64[ind - 1];	// Cstar = 10^(q-x-1)
369     if (Cstar.w[1] == ten2k128[0].w[1]
370 	&& Cstar.w[0] == ten2k128[0].w[0]) {
371       // if  Cstar = 10^(q-x)
372       Cstar.w[0] = ten2k64[19];	// Cstar = 10^(q-x-1)
373       Cstar.w[1] = 0x0ULL;
379     if (Cstar.w[1] == ten2k128[ind - 20].w[1] &&
380 	Cstar.w[0] == ten2k128[ind - 20].w[0]) {
381       // if  Cstar = 10^(q-x)
382       Cstar.w[0] = ten2k128[ind - 21].w[0];	// Cstar = 10^(q-x-1)
383       Cstar.w[1] = ten2k128[ind - 21].w[1];
389   ptr_Cstar->w[1] = Cstar.w[1];
390   ptr_Cstar->w[0] = Cstar.w[0];
407   UINT192 Cstar;
478   // Cstar = P384 >> Ex
482     Cstar.w[2] = (P384.w[5] >> shift);
483     Cstar.w[1] = (P384.w[5] << (64 - shift)) | (P384.w[4] >> shift);
484     Cstar.w[0] = (P384.w[4] << (64 - shift)) | (P384.w[3] >> shift);
492     Cstar.w[2] = 0x0ULL;
493     Cstar.w[1] = P384.w[5] >> shift;
494     Cstar.w[0] = (P384.w[5] << (64 - shift)) | (P384.w[4] >> shift);
502     Cstar.w[2] = 0x0ULL;
503     Cstar.w[1] = 0x0ULL;
504     Cstar.w[0] = P384.w[5] >> shift;
580     if (Cstar.w[0] & 0x01) {	// Cstar is odd; MP in [EVEN, ODD]
582       Cstar.w[0]--;	// Cstar is now even
583       if (Cstar.w[0] == 0xffffffffffffffffULL) {
584 	Cstar.w[1]--;
585 	if (Cstar.w[1] == 0xffffffffffffffffULL) {
586 	  Cstar.w[2]--;
598   // check for rounding overflow, which occurs if Cstar = 10^(q-x)
601     if (Cstar.w[2] == 0x0ULL && Cstar.w[1] == 0x0ULL &&
602 	Cstar.w[0] == ten2k64[ind]) {
603       // if  Cstar = 10^(q-x)
604       Cstar.w[0] = ten2k64[ind - 1];	// Cstar = 10^(q-x-1)
611     if (Cstar.w[2] == 0x0ULL && Cstar.w[1] == ten2k128[0].w[1] &&
612 	Cstar.w[0] == ten2k128[0].w[0]) {
613       // if  Cstar = 10^(q-x)
614       Cstar.w[0] = ten2k64[19];	// Cstar = 10^(q-x-1)
615       Cstar.w[1] = 0x0ULL;
621     if (Cstar.w[2] == 0x0ULL && Cstar.w[1] == ten2k128[ind - 20].w[1] &&
622 	Cstar.w[0] == ten2k128[ind - 20].w[0]) {
623       // if  Cstar = 10^(q-x)
624       Cstar.w[0] = ten2k128[ind - 21].w[0];	// Cstar = 10^(q-x-1)
625       Cstar.w[1] = ten2k128[ind - 21].w[1];
631     if (Cstar.w[2] == ten2k256[0].w[2] && Cstar.w[1] == ten2k256[0].w[1]
632 	&& Cstar.w[0] == ten2k256[0].w[0]) {
633       // if  Cstar = 10^(q-x)
634       Cstar.w[0] = ten2k128[18].w[0];	// Cstar = 10^(q-x-1)
635       Cstar.w[1] = ten2k128[18].w[1];
636       Cstar.w[2] = 0x0ULL;
642     if (Cstar.w[2] == ten2k256[ind - 39].w[2] &&
643 	Cstar.w[1] == ten2k256[ind - 39].w[1] &&
644 	Cstar.w[0] == ten2k256[ind - 39].w[0]) {
645       // if  Cstar = 10^(q-x)
646       Cstar.w[0] = ten2k256[ind - 40].w[0];	// Cstar = 10^(q-x-1)
647       Cstar.w[1] = ten2k256[ind - 40].w[1];
648       Cstar.w[2] = ten2k256[ind - 40].w[2];
654   ptr_Cstar->w[2] = Cstar.w[2];
655   ptr_Cstar->w[1] = Cstar.w[1];
656   ptr_Cstar->w[0] = Cstar.w[0];
673   UINT256 Cstar;
790   // Cstar = P512 >> Ex
794     Cstar.w[3] = (P512.w[7] >> shift);
795     Cstar.w[2] = (P512.w[7] << (64 - shift)) | (P512.w[6] >> shift);
796     Cstar.w[1] = (P512.w[6] << (64 - shift)) | (P512.w[5] >> shift);
797     Cstar.w[0] = (P512.w[5] << (64 - shift)) | (P512.w[4] >> shift);
807     Cstar.w[3] = 0x0ULL;
808     Cstar.w[2] = P512.w[7] >> shift;
809     Cstar.w[1] = (P512.w[7] << (64 - shift)) | (P512.w[6] >> shift);
810     Cstar.w[0] = (P512.w[6] << (64 - shift)) | (P512.w[5] >> shift);
820     Cstar.w[3] = 0x0ULL;
821     Cstar.w[2] = 0x0ULL;
822     Cstar.w[1] = P512.w[7] >> shift;
823     Cstar.w[0] = (P512.w[7] << (64 - shift)) | (P512.w[6] >> shift);
833     Cstar.w[3] = 0x0ULL;
834     Cstar.w[2] = 0x0ULL;
835     Cstar.w[1] = 0x0ULL;
836     Cstar.w[0] = P512.w[7];
846     Cstar.w[3] = 0x0ULL;
847     Cstar.w[2] = 0x0ULL;
848     Cstar.w[1] = 0x0ULL;
849     Cstar.w[0] = P512.w[7] >> shift;
949     if (Cstar.w[0] & 0x01) {	// Cstar is odd; MP in [EVEN, ODD]
951       Cstar.w[0]--;	// Cstar is now even
952       if (Cstar.w[0] == 0xffffffffffffffffULL) {
953 	Cstar.w[1]--;
954 	if (Cstar.w[1] == 0xffffffffffffffffULL) {
955 	  Cstar.w[2]--;
956 	  if (Cstar.w[2] == 0xffffffffffffffffULL) {
957 	    Cstar.w[3]--;
970   // check for rounding overflow, which occurs if Cstar = 10^(q-x)
973     if (Cstar.w[3] == 0x0ULL && Cstar.w[2] == 0x0ULL &&
974 	Cstar.w[1] == 0x0ULL && Cstar.w[0] == ten2k64[ind]) {
975       // if  Cstar = 10^(q-x)
976       Cstar.w[0] = ten2k64[ind - 1];	// Cstar = 10^(q-x-1)
983     if (Cstar.w[3] == 0x0ULL && Cstar.w[2] == 0x0ULL &&
984 	Cstar.w[1] == ten2k128[0].w[1]
985 	&& Cstar.w[0] == ten2k128[0].w[0]) {
986       // if  Cstar = 10^(q-x)
987       Cstar.w[0] = ten2k64[19];	// Cstar = 10^(q-x-1)
988       Cstar.w[1] = 0x0ULL;
994     if (Cstar.w[3] == 0x0ULL && Cstar.w[2] == 0x0ULL &&
995 	Cstar.w[1] == ten2k128[ind - 20].w[1] &&
996 	Cstar.w[0] == ten2k128[ind - 20].w[0]) {
997       // if  Cstar = 10^(q-x)
998       Cstar.w[0] = ten2k128[ind - 21].w[0];	// Cstar = 10^(q-x-1)
999       Cstar.w[1] = ten2k128[ind - 21].w[1];
1005     if (Cstar.w[3] == 0x0ULL && Cstar.w[2] == ten2k256[0].w[2] &&
1006 	Cstar.w[1] == ten2k256[0].w[1]
1007 	&& Cstar.w[0] == ten2k256[0].w[0]) {
1008       // if  Cstar = 10^(q-x)
1009       Cstar.w[0] = ten2k128[18].w[0];	// Cstar = 10^(q-x-1)
1010       Cstar.w[1] = ten2k128[18].w[1];
1011       Cstar.w[2] = 0x0ULL;
1017     if (Cstar.w[3] == 0x0ULL && Cstar.w[2] == ten2k256[ind - 39].w[2] &&
1018 	Cstar.w[1] == ten2k256[ind - 39].w[1] &&
1019 	Cstar.w[0] == ten2k256[ind - 39].w[0]) {
1020       // if  Cstar = 10^(q-x)
1021       Cstar.w[0] = ten2k256[ind - 40].w[0];	// Cstar = 10^(q-x-1)
1022       Cstar.w[1] = ten2k256[ind - 40].w[1];
1023       Cstar.w[2] = ten2k256[ind - 40].w[2];
1030     if (Cstar.w[3] == ten2k256[ind - 39].w[3] &&
1031 	Cstar.w[2] == ten2k256[ind - 39].w[2] &&
1032 	Cstar.w[1] == ten2k256[ind - 39].w[1] &&
1033 	Cstar.w[0] == ten2k256[ind - 39].w[0]) {
1034       // if  Cstar = 10^(q-x)
1035       Cstar.w[0] = ten2k256[ind - 40].w[0];	// Cstar = 10^(q-x-1)
1036       Cstar.w[1] = ten2k256[ind - 40].w[1];
1037       Cstar.w[2] = ten2k256[ind - 40].w[2];
1038       Cstar.w[3] = ten2k256[ind - 40].w[3];
1044   ptr_Cstar->w[3] = Cstar.w[3];
1045   ptr_Cstar->w[2] = Cstar.w[2];
1046   ptr_Cstar->w[1] = Cstar.w[1];
1047   ptr_Cstar->w[0] = Cstar.w[0];