33 lines
		
	
	
		
			656 B
		
	
	
	
		
			C
		
	
	
	
	
	
			
		
		
	
	
			33 lines
		
	
	
		
			656 B
		
	
	
	
		
			C
		
	
	
	
	
	
| #include "tommath_private.h"
 | |
| #ifdef BN_MP_REDUCE_2K_SETUP_C
 | |
| /* LibTomMath, multiple-precision integer library -- Tom St Denis */
 | |
| /* SPDX-License-Identifier: Unlicense */
 | |
| 
 | |
| /* determines the setup value */
 | |
| mp_err mp_reduce_2k_setup(const mp_int *a, mp_digit *d)
 | |
| {
 | |
|    mp_err err;
 | |
|    mp_int tmp;
 | |
|    int    p;
 | |
| 
 | |
|    if ((err = mp_init(&tmp)) != MP_OKAY) {
 | |
|       return err;
 | |
|    }
 | |
| 
 | |
|    p = mp_count_bits(a);
 | |
|    if ((err = mp_2expt(&tmp, p)) != MP_OKAY) {
 | |
|       mp_clear(&tmp);
 | |
|       return err;
 | |
|    }
 | |
| 
 | |
|    if ((err = s_mp_sub(&tmp, a, &tmp)) != MP_OKAY) {
 | |
|       mp_clear(&tmp);
 | |
|       return err;
 | |
|    }
 | |
| 
 | |
|    *d = tmp.dp[0];
 | |
|    mp_clear(&tmp);
 | |
|    return MP_OKAY;
 | |
| }
 | |
| #endif
 | 
