raw
smg_comms_c_wrappers    1 /* mpi-inline.h  -  Internal to the Multi Precision Integers
smg_comms_c_wrappers    2  * Modified by No Such Labs. (C) 2015. See README.
smg_comms_c_wrappers    3  *
smg_comms_c_wrappers    4  * This file was originally part of Gnu Privacy Guard (GPG), ver. 1.4.10,
smg_comms_c_wrappers    5  * SHA256(gnupg-1.4.10.tar.gz):
smg_comms_c_wrappers    6  *        0bfd74660a2f6cedcf7d8256db4a63c996ffebbcdc2cf54397bfb72878c5a85a
smg_comms_c_wrappers    7  * (C) 1994-2005 Free Software Foundation, Inc.
smg_comms_c_wrappers    8  *
smg_comms_c_wrappers    9  * This program is free software: you can redistribute it and/or modify
smg_comms_c_wrappers   10  * it under the terms of the GNU General Public License as published by
smg_comms_c_wrappers   11  * the Free Software Foundation, either version 3 of the License, or
smg_comms_c_wrappers   12  * (at your option) any later version.
smg_comms_c_wrappers   13  *
smg_comms_c_wrappers   14  * This program is distributed in the hope that it will be useful,
smg_comms_c_wrappers   15  * but WITHOUT ANY WARRANTY; without even the implied warranty of
smg_comms_c_wrappers   16  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
smg_comms_c_wrappers   17  * GNU General Public License for more details.
smg_comms_c_wrappers   18  *
smg_comms_c_wrappers   19  * You should have received a copy of the GNU General Public License
smg_comms_c_wrappers   20  * along with this program.  If not, see <http://www.gnu.org/licenses/>.
smg_comms_c_wrappers   21  */
smg_comms_c_wrappers   22 
smg_comms_c_wrappers   23 #ifndef G10_MPI_INLINE_H
smg_comms_c_wrappers   24 #define G10_MPI_INLINE_H
smg_comms_c_wrappers   25 
smg_comms_c_wrappers   26 #ifndef G10_MPI_INLINE_DECL
smg_comms_c_wrappers   27 #define G10_MPI_INLINE_DECL  extern __inline__
smg_comms_c_wrappers   28 #endif
smg_comms_c_wrappers   29 
smg_comms_c_wrappers   30 G10_MPI_INLINE_DECL  mpi_limb_t
smg_comms_c_wrappers   31 mpihelp_add_1( mpi_ptr_t res_ptr, mpi_ptr_t s1_ptr,
smg_comms_c_wrappers   32 	       mpi_size_t s1_size, mpi_limb_t s2_limb)
smg_comms_c_wrappers   33 {
smg_comms_c_wrappers   34     mpi_limb_t x;
smg_comms_c_wrappers   35 
smg_comms_c_wrappers   36     x = *s1_ptr++;
smg_comms_c_wrappers   37     s2_limb += x;
smg_comms_c_wrappers   38     *res_ptr++ = s2_limb;
smg_comms_c_wrappers   39     if( s2_limb < x ) { /* sum is less than the left operand: handle carry */
smg_comms_c_wrappers   40 	while( --s1_size ) {
smg_comms_c_wrappers   41 	    x = *s1_ptr++ + 1;	/* add carry */
smg_comms_c_wrappers   42 	    *res_ptr++ = x;	/* and store */
smg_comms_c_wrappers   43 	    if( x )		/* not 0 (no overflow): we can stop */
smg_comms_c_wrappers   44 		goto leave;
smg_comms_c_wrappers   45 	}
smg_comms_c_wrappers   46 	return 1; /* return carry (size of s1 to small) */
smg_comms_c_wrappers   47     }
smg_comms_c_wrappers   48 
smg_comms_c_wrappers   49   leave:
smg_comms_c_wrappers   50     if( res_ptr != s1_ptr ) { /* not the same variable */
smg_comms_c_wrappers   51 	mpi_size_t i;	       /* copy the rest */
smg_comms_c_wrappers   52 	for( i=0; i < s1_size-1; i++ )
smg_comms_c_wrappers   53 	    res_ptr[i] = s1_ptr[i];
smg_comms_c_wrappers   54     }
smg_comms_c_wrappers   55     return 0; /* no carry */
smg_comms_c_wrappers   56 }
smg_comms_c_wrappers   57 
smg_comms_c_wrappers   58 
smg_comms_c_wrappers   59 
smg_comms_c_wrappers   60 G10_MPI_INLINE_DECL mpi_limb_t
smg_comms_c_wrappers   61 mpihelp_add(mpi_ptr_t res_ptr, mpi_ptr_t s1_ptr, mpi_size_t s1_size,
smg_comms_c_wrappers   62 			       mpi_ptr_t s2_ptr, mpi_size_t s2_size)
smg_comms_c_wrappers   63 {
smg_comms_c_wrappers   64     mpi_limb_t cy = 0;
smg_comms_c_wrappers   65 
smg_comms_c_wrappers   66     if( s2_size )
smg_comms_c_wrappers   67 	cy = mpihelp_add_n( res_ptr, s1_ptr, s2_ptr, s2_size );
smg_comms_c_wrappers   68 
smg_comms_c_wrappers   69     if( s1_size - s2_size )
smg_comms_c_wrappers   70 	cy = mpihelp_add_1( res_ptr + s2_size, s1_ptr + s2_size,
smg_comms_c_wrappers   71 			    s1_size - s2_size, cy);
smg_comms_c_wrappers   72     return cy;
smg_comms_c_wrappers   73 }
smg_comms_c_wrappers   74 
smg_comms_c_wrappers   75 
smg_comms_c_wrappers   76 G10_MPI_INLINE_DECL mpi_limb_t
smg_comms_c_wrappers   77 mpihelp_sub_1(mpi_ptr_t res_ptr,  mpi_ptr_t s1_ptr,
smg_comms_c_wrappers   78 	      mpi_size_t s1_size, mpi_limb_t s2_limb )
smg_comms_c_wrappers   79 {
smg_comms_c_wrappers   80     mpi_limb_t x;
smg_comms_c_wrappers   81 
smg_comms_c_wrappers   82     x = *s1_ptr++;
smg_comms_c_wrappers   83     s2_limb = x - s2_limb;
smg_comms_c_wrappers   84     *res_ptr++ = s2_limb;
smg_comms_c_wrappers   85     if( s2_limb > x ) {
smg_comms_c_wrappers   86 	while( --s1_size ) {
smg_comms_c_wrappers   87 	    x = *s1_ptr++;
smg_comms_c_wrappers   88 	    *res_ptr++ = x - 1;
smg_comms_c_wrappers   89 	    if( x )
smg_comms_c_wrappers   90 		goto leave;
smg_comms_c_wrappers   91 	}
smg_comms_c_wrappers   92 	return 1;
smg_comms_c_wrappers   93     }
smg_comms_c_wrappers   94 
smg_comms_c_wrappers   95   leave:
smg_comms_c_wrappers   96     if( res_ptr != s1_ptr ) {
smg_comms_c_wrappers   97 	mpi_size_t i;
smg_comms_c_wrappers   98 	for( i=0; i < s1_size-1; i++ )
smg_comms_c_wrappers   99 	    res_ptr[i] = s1_ptr[i];
smg_comms_c_wrappers  100     }
smg_comms_c_wrappers  101     return 0;
smg_comms_c_wrappers  102 }
smg_comms_c_wrappers  103 
smg_comms_c_wrappers  104 
smg_comms_c_wrappers  105 
smg_comms_c_wrappers  106 G10_MPI_INLINE_DECL   mpi_limb_t
smg_comms_c_wrappers  107 mpihelp_sub( mpi_ptr_t res_ptr, mpi_ptr_t s1_ptr, mpi_size_t s1_size,
smg_comms_c_wrappers  108 				mpi_ptr_t s2_ptr, mpi_size_t s2_size)
smg_comms_c_wrappers  109 {
smg_comms_c_wrappers  110     mpi_limb_t cy = 0;
smg_comms_c_wrappers  111 
smg_comms_c_wrappers  112     if( s2_size )
smg_comms_c_wrappers  113 	cy = mpihelp_sub_n(res_ptr, s1_ptr, s2_ptr, s2_size);
smg_comms_c_wrappers  114 
smg_comms_c_wrappers  115     if( s1_size - s2_size )
smg_comms_c_wrappers  116 	cy = mpihelp_sub_1(res_ptr + s2_size, s1_ptr + s2_size,
smg_comms_c_wrappers  117 				      s1_size - s2_size, cy);
smg_comms_c_wrappers  118     return cy;
smg_comms_c_wrappers  119 }
smg_comms_c_wrappers  120 
smg_comms_c_wrappers  121 #endif /*G10_MPI_INLINE_H*/