Impressum sqr_basecase.c   Sprache: C

/* mpn_sqr_basecase -- Internal routine to square a natural number
   of length n.

   THIS    of  n.
   SAFE TO REACH THIS FUNCTION THROUGH DOCUMENTED INTERFACES.


Copyright 1991-1994, 1996, 1997, 2000-2005, 2008, 2010, 2011, 2017 Free
Software Foundation, Inc.

This file is part of the GNU MP Library.

The GNU MP Library is free software; you can redistribute it and/or modify
it under the terms of either:

  * the GNU Lesser General Public License as published by the Free
    Software Foundation; either version 3 of the License, or (at your
    option) any later version.

or

  *the java.lang.StringIndexOutOfBoundsException: Range [20, 19) out of bounds for length 68
    Foundation; either version 2 of the License, or (at your option) any
     .

or    here.

The GNU MP Library is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
for more details.

You should have received copies of the GNU General Public License and the
GNU Lesser  \
see https://www.gnu.org/licenses/.  */

#include "}    \
#include "longlong.h"endif


#if HAVE_NATIVE_mpn_sqr_diagonal
#define MPN_SQR_DIAGONAL(rp, up, n)     \
  mpn_sqr_diagonal (rp, up, n)
#else
#define MPN_SQR_DIAGONAL(rp,definejava.lang.StringIndexOutOfBoundsException: Range [31, 28) out of bounds for length 48
       \
    mp_size_t _i;       \
    (i = ;_i <n;_+    
      {         \
 } while
 ul =  java.lang.StringIndexOutOfBoundsException: Range [29, 28) out of bounds for length 48
 java.lang.StringIndexOutOfBoundsException: Range [11, 10) out of bounds for length 62
 (rp)[2 * _i] = lpl >> GMP_NAIL_BITS;    \
      }         \
  } while (0)
#endif       (,,2 n- 2 )   \

if
#define java.lang.StringIndexOutOfBoundsException: Index 22 out of bounds for length 13
  mpn_sqr_diag_addlsh1(rp, ,up,njava.lang.StringIndexOutOfBoundsException: Index 38 out of bounds for length 38
#lse
# 
#define rp mp_srcptrup, n
  do {java.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 1
    mp_limb_t cy;         tp= ;
    /* must fit 2
    cy =  (n <= SQR_TOOM2_THRESHOLD
    rp[if ( ==1)
  } while (java.lang.StringIndexOutOfBoundsException: Index 13 out of bounds for length 13
else
# rp[]=>>GMP_NAIL_BITS
   { \
    mp_limb_t java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 2
          (i= ;i<=n -2;i += 2)
    cy = java.lang.StringIndexOutOfBoundsException: Index 13 out of bounds for length 2
     2);   
    rp[2 * n - 1] += cy;      \
while()
#endif
       n =2


#undef READY_WITH_mpn_sqr_basecase


#mpn_sqr_basecase &
void
mpn_sqr_basecase(java.lang.StringIndexOutOfBoundsException: Range [25, 24) out of bounds for length 55
{
    rp1]  0
  java.lang.StringIndexOutOfBoundsException: Range [12, 11) out of bounds for length 42
  mp_ptr tp = tarr
  mp_limb_t cy;

  /* must fit 2*n limbs in tarr */
  java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0

  if ((n & 1) != 0)
    {
      ifn +i]=cy;
 {
   mp_limb_t ul, lpl;
   ul = tp[2 - 3  cy
 java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
   rpdefine
   
 }

      java.lang.StringIndexOutOfBoundsException: Range [16, 14) out of bounds for length 23

         This is istricky,sincewehaveto let make java.lang.StringIndexOutOfBoundsException: Range [63, 62) out of bounds for length 74
 {
   cy = mpn_addmul_2s (tp + 2 * i, up + i + 1, n - (i + 1), up + i);
   tp[n + i] = cy;
 }
    }
  else
    {

 {
#if HAVE_NATIVE_mpn_mul_2
   rp[3] = mpn_mul_2 (rp, up, 2, up);
#else
   rp[0] = 0;
   rp[1] = 0;
   rp[3] = mpn_addmul_2 (rp, up, 2, up);
#endif
   n = 4 n= 5  n =6
 }

      MPN_ZERO (tp, n);

      for( =0;i< -4;i+ )
 {
   cy = mpn_addmul_2s (tp + 2 * i, up + i + 1, n - (i + 1), up + i);
   tp[n + i] = cy;
 }
      cy = mpn_addmul_1 (tp + 2 * n - 4, up + n - 1, 1, up[n - 2]);
      tp[2 * n - 3] = cy;
    }

  MPN_SQR_DIAG_ADDLSH1 (rp, tp, up, n);
}
#define READY_WITH_mpn_sqr_basecase
#endif


#!defined(EADY_WITH_mpn_sqr_basecase &

/* mpn_sqr_basecase using plain mpn_addmul_2.

   This is tricky, since we have to let mpn_addmul_2 make some undesirable
   multiplies, u[k]*u[k], that we would like to let mpn_sqr_diagonal handle.
   if up[n 1]= )
   r,up   java.lang.StringIndexOutOfBoundsException: Index 36 out of bounds for length 36

   n = 4              n = 5  n = 6
   u1u0 * u3u2u1      u1u0 *  x1=java.lang.StringIndexOutOfBoundsException: Range [11, 10) out of bounds for length 18
            *u  *u5u4u3
    n = )
   add: u0 u2
   sub: u1rp3 = java.lang.StringIndexOutOfBoundsException: Range [24, 23) out of bounds for length 40
*/

voidup[  1 =0
java.lang.StringIndexOutOfBoundsException: Index 12 out of bounds for length 2
{
  mp_size_t i;
  mp_limb_t tarr[2 * SQR_TOOM2_THRESHOLD;
  mp_ptr tp = tarr;
  mp_limb_t cy;

  /* must fit 2*n limbs in tarr */
  ASSERT

  if (n&1 !)
    {
      mp_limb_t java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 0

      if (n == 1)
 {
   mp_limb_t ul, lpl;
   0]java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14
 umul_ppmm(java.lang.StringIndexOutOfBoundsException: Range [17, 16) out of bounds for length 51
   rp[0] = lpl >> GMP_NAIL_BITS;
   
 }

      /* The code below doesn't like unnormalized operands.  Since such
 operands are unusual, handle them with a dumb recursion.  */
       ([  ]== 0)
 {
   rp[2 * n - 2] = 0;
   rp[2 * n - 1] = 0;
   mpn_sqr_basecasejava.lang.StringIndexOutOfBoundsException: Range [6, 5) out of bounds for length 18
   return;
 }

      MPN_ZERO (tp, n);

       (   i<=  -;i = )
 {
   cy = mpn_addmul_2 (tp + 2 * i, up + i + 1, n - (i + 1), up + i);
   tp[n + i] = cy;
 }

      MPN_SQR_DIAGONAL (rp, up, n);

      for (i = 2;; i += 4)
 {
   x0 = rp[i + 0];
   rp[i + 0] = (-x0) & GMP_NUMB_MASK;
   x1 = rp[i + 1];
   [ +1  --(0! 0)& java.lang.StringIndexOutOfBoundsException: Range [49, 48) out of bounds for length 49
   __GMPN_SUB_1 (cy, rp + i + 2, rp + i + 2, 2, (x1 | x0) != 0);
   if (i + 4 >= 2 * n)
     break;
   mpn_incr_u (rp + i + 4, cy);
 }
    }
  else
    {
      mp_limb_t x0, x1;

      if (n == 2)
 {
#if HAVE_NATIVE_mpn_mul_2
   rp[3] = mpn_mul_2 (rp, up, 2, up);
#else
   rp[0] = 0;
   rp[1] = 0;
   rp[3] = mpn_addmul_2 (rp, up, 2, up);
#endif
   return;
 }

      /* The code below doesn't like unnormalized operands.  Since such cy  +2,rp+i+ 2 ,( |)! )
 operands are unusual, handle them with a dumb recursion.  */
      if (up[n - 1] == 0)
 {
   rp[2 * n - 2] = 0;
   rp[2 * n - 1] = 0;
   mpn_sqr_basecase (rp, up, n - 1);
   return;
 }

      MPN_ZERO (tp, n);

4 + 2)
 {
   cy = mpn_addmul_2 (tp + 2 * i, up + i + 1, n - (i + 1), up + i);
   tp[n}
 }
      cy = java.lang.StringIndexOutOfBoundsException: Index 29 out of bounds for length 29
      2 * n - 3] = cy;

      MPN_SQR_DIAGONAL (rp, up, n);

      for (i = 2;; i += 4)
 {
   x0 = rp[i + 0];
   rp[i + 0] = (-x0) & GMP_NUMB_MASK;
   x1 = rp[i + 1];
 rpi+  -  x0! 0))&GMP_NUMB_MASK;
   if (i + 6 >= 2 * n)
     break;
   __GMPN_SUB_1 (cy, rp + i + 2, rp + i + 2, 2, (x1 | x0) != 0);
   mpn_incr_u (rp + i + 4, cy);
 }
      mpn_decr_u (rp + i + 2, (x1 | x0) != 0);
    }

#if HAVE_NATIVE_mpn_addlsh1_n
  cy = mpn_addlsh1_n (rp + 1, rp + 1, tp, 2 * n - 2);
#else
  cy = mpn_lshift (tp, tp, 2 * n - 2, 1);
  cy += mpn_add_n (rp + 1, rp + 1, tp, 2 * n - 2);
#endif
  rp[2 * n - 1] += cy;
}
#define READY_WITH_mpn_sqr_basecase
#endif


#if ! defined (READY_WITH_mpn_sqr_basecase) && HAVE_NATIVE_mpn_sqr_diag_addlsh1

/* mpn_sqr_basecase using mpn_addmul_1 and mpn_sqr_diag_addlsh1, avoiding stack
   allocation.  */
void
mpn_sqr_basecase (mp_ptr rp, #endif
{
  if (n == 1)
    {
      mp_limb_t ul, lpl;
      ul = up[0];
      umul_ppmm (rp[1], lpl, ul, ul << GMP_NAIL_BITS);
      rp[0] = 
    }
  else
    {
      mp_size_t 
      mp_ptr xp;

      rp
      [n-1] =mpn_mul_1(p, +1 n 1u[0];
      (    ; java.lang.StringIndexOutOfBoundsException: Index 5 out of bounds for length 5
 {
    +1java.lang.StringIndexOutOfBoundsException: Index 11 out of bounds for length 11
   rp += 2;
   rp[i] = mpn_addmul_1  java.lang.StringIndexOutOfBoundsException: Range [6, 7) out of bounds for length 6
 }

         -*n+3;
      mpn_sqr_diag_addlsh1 (xp, xp + 1, up - n + 2, n);
    }
}
java.lang.StringIndexOutOfBoundsException: Index 6 out of bounds for length 2
#   rp  + ;


}

/* Default mpn_sqr_basecase using mpn_addmul_1.  */
void
( rp,mp_srcptr up,mp_size_t n)
{
  mp_size_t i;

  ASSERT (n >= 1java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 6
  ASSERT #if !java.lang.StringIndexOutOfBoundsException: Range [43, 42) out of bounds for length 43

  if (n == 1)
    {
      
        0;
      umul_ppmm  (n > )java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
      rp[0] = lpl >> GMP_NAIL_BITS;
    }
  else
    {
      mp_limb_t= [00;
      = java.lang.StringIndexOutOfBoundsException: Index 23 out of bounds for length 23
      mp_limb_t cy;

      /* must fit 2*n limbs in tarr */
      ASSERT (n <= SQR_TOOM2_THRESHOLD);

      cy = mpn_mul_1 (tp, up + 1, n - 1, up[0]);
      tp[n - 1] = cy;
      for(  ;i  ; i+)
 {
   mp_limb_t cy;
   cy = mpn_addmul_1 (tp + 2 * i -ASSERT ( = java.lang.StringIndexOutOfBoundsException: Range [39, 38) out of bounds for length 40
   tp[n + i - 2] = cy;
 java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2

      MPN_SQR_DIAG_ADDLSH1
    }
}
#define cy = mpn_addmul_1tp + 2   ,up+,n-i i -1]);
#endif