/* TomsFastMath, a fast ISO C bignum library. * * This project is meant to fill in where LibTomMath * falls short. That is speed ;-) * * This project is public domain and free for all purposes. * * Tom St Denis, tomstdenis@gmail.com */ #include "bignum_fast.h" /* b = a*a */ void fp_sqr(fp_int *A, fp_int *B) { int y; /* call generic if we're out of range */ if (A->used + A->used > FP_SIZE) { fp_sqr_comba(A, B); return ; } y = A->used; #if defined(TFM_SQR3) if (y <= 3) { fp_sqr_comba3(A,B); return; } #endif #if defined(TFM_SQR4) if (y == 4) { fp_sqr_comba4(A,B); return; } #endif #if defined(TFM_SQR6) if (y <= 6) { fp_sqr_comba6(A,B); return; } #endif #if defined(TFM_SQR7) if (y == 7) { fp_sqr_comba7(A,B); return; } #endif #if defined(TFM_SQR8) if (y == 8) { fp_sqr_comba8(A,B); return; } #endif #if defined(TFM_SQR9) if (y == 9) { fp_sqr_comba9(A,B); return; } #endif #if defined(TFM_SQR12) if (y <= 12) { fp_sqr_comba12(A,B); return; } #endif #if defined(TFM_SQR17) if (y <= 17) { fp_sqr_comba17(A,B); return; } #endif #if defined(TFM_SMALL_SET) if (y <= 16) { fp_sqr_comba_small(A,B); return; } #endif #if defined(TFM_SQR20) if (y <= 20) { fp_sqr_comba20(A,B); return; } #endif #if defined(TFM_SQR24) if (y <= 24) { fp_sqr_comba24(A,B); return; } #endif #if defined(TFM_SQR28) if (y <= 28) { fp_sqr_comba28(A,B); return; } #endif #if defined(TFM_SQR32) if (y <= 32) { fp_sqr_comba32(A,B); return; } #endif #if defined(TFM_SQR48) if (y <= 48) { fp_sqr_comba48(A,B); return; } #endif #if defined(TFM_SQR64) if (y <= 64) { fp_sqr_comba64(A,B); return; } #endif fp_sqr_comba(A, B); } /* $Source: /cvs/libtom/tomsfastmath/src/sqr/fp_sqr.c,v $ */ /* $Revision: 1.1 $ */ /* $Date: 2006/12/31 21:25:53 $ */