---

下位関数群
[ビット演算]

関数
size_t	ymp_rshiftabs (digit_t *result, const digit_t *source, size_t len, size_t width)
	多倍長の自然数に対し、右に向かってビットシフトを行なう。. より詳しく...
size_t	ymp_lshiftabs (digit_t *result, const digit_t *source, size_t len, size_t width)
	多倍長の自然数に対し、左に向かってビットシフトを行なう。. より詳しく...
size_t	ymp_modabs_2exp (digit_t *result, const digit_t *self, size_t len, size_t pow)
	多倍長自然数の2の冪乗による剰余. より詳しく...

---

関数の解説

size_t ymp_rshiftabs (  digit_t *    result, const digit_t *    source, size_t    len, size_t    width )

多倍長の自然数に対し、右に向かってビットシフトを行なう。. 溢れた分は切り捨てられる。 引数: result  結果を格納するdigit_tの配列の先頭へのポインタ source  元のdigit_tの列の先頭へのポインタ len  sourceがポイントする列の要素数 width  シフトのビット幅 事前条件: resultのポイントする領域は十分な大きさであること。 覚え書き: resultのポイントする領域は多くともMAX(0, len-width/DIGIT_BIT)個の 要素を格納できれば十分である。 resultとsourceは一致してもよい 戻り値: resultに格納された要素数 shift.c の 34 行で定義されています。 { register size_t i; size_t n_digits = width / DIGIT_BIT; size_t n_bits = width % DIGIT_BIT; digit_t prev; if (len <= n_digits) return 0; i=len; prev=0; while (i-- > n_digits) { double_digit_t tmp = LOW_DIGIT_TO_HIGH_DIGIT((double_digit_t)prev); prev = source[i]; tmp += source[i]; result[i-n_digits] = LOW_DIGIT(tmp >> n_bits); } /* 要素数の調整 */ len -= n_digits; /* 今、len > n_digits */ if (result[len-1]==0) --len; return len; }

size_t ymp_lshiftabs (  digit_t *    result, const digit_t *    source, size_t    len, size_t    width )

多倍長の自然数に対し、左に向かってビットシフトを行なう。. 右側は0で埋められる。 引数: result  結果を格納するdigit_tの配列の先頭へのポインタ source  元のdigit_tの列の先頭へのポインタ len  sourceがポイントする列の要素数 width  シフトのビット幅 事前条件: resultのポイントする領域は十分な大きさがあること 覚え書き: resultのポイントする領域は多くとも(len+width/DIGIT_BIT+1)個の要素を 格納できれば十分である。 resultとsourceは一致してもよい 戻り値: resultに格納された要素数 長さ0の列に対してはシフトを行わずに単に0を返す。 shift.c の 76 行で定義されています。 { register size_t i; size_t n_digits = width / DIGIT_BIT; size_t n_bits = width % DIGIT_BIT; double_digit_t tmp; if (len == 0) return 0; i = len; tmp = source[--i]; tmp <<= n_bits; if (HIGH_DIGIT(tmp)) result[len++ + n_digits] = HIGH_DIGIT(tmp); tmp = source[i]; while (i-- > 0) { tmp = LOW_DIGIT_TO_HIGH_DIGIT(tmp); tmp += source[i]; result[i+n_digits+1] = HIGH_DIGIT(tmp << n_bits); tmp = source[i]; } result[n_digits] = source[0] << n_bits; i = n_digits; while (i-- > 0) result[i] = 0; return len + n_digits; }

size_t ymp_modabs_2exp (  digit_t *    result, const digit_t *    self, size_t    len, size_t    pow )

多倍長自然数の2の冪乗による剰余. <self>を2のpow乗で割った剰余をresultに格納する。 引数: result  剰余を格納するdigit_tの列の先頭へのポインタ self  元の数を表すdigit_tの列の先頭へのポインタ len  selfのポイントする領域の長さ(要素数) pow  冪指数 事前条件: resultは pow / DIGIT_BIT + (powDIGIT_BIT)?1:0個の要素を格納できること 戻り値: resultに格納された要素数 覚え書き: selfとresultは一致してもよい 格納領域の長さは最大でも pow / DIGIT_BIT + 1 あれば十分 shift.c の 120 行で定義されています。 参照 DIGIT_BIT, と digit_t. 呼出 ymp_mod_2exp, と ymp_modpowabs_2exp_z. { size_t n_digits = pow / DIGIT_BIT; size_t n_bits = pow % DIGIT_BIT; if (n_digits < len) { len = n_digits; if (n_bits) { result[n_digits] = self[n_digits] & ( (1u << n_bits) - 1u ); ++n_digits; } } else if (n_digits > len) { n_digits = len; } if (result != self) memcpy(result, self, len*sizeof(digit_t)); while (n_digits-- > 0) { if (self[n_digits] != 0) break; } return n_digits+1; }

---