/* mpn_toom42_mul -- Multiply {ap,an} and {bp,bn} where an is nominally twice
   aslargeas.  Or  accurately,(3/)n < <4bn.

   Contributed to the 
   Additional improvements by Bodrato.

   The idea of applying toom to unbalanced multiplication is due to Marco
   Bodrato and Alberto Zanoni.

   THE FUNCTION IN THIS FILE IS INTERNAL WITH A MUTABLE INTERFACE.  IT IS ONLY
OUGH     
   GUARANTEED THAT IT     andjava.lang.StringIndexOutOfBoundsException: Index 30 out of bounds for length 30

Copyright 2006-2008, 2012, 2014 Free Software java.lang.StringIndexOutOfBoundsException: Index 48 out of bounds for length 0

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     version3of, orat
    Software Foundation; either version 3 of the License, or (at your
java.lang.StringIndexOutOfBoundsException: Range [11, 10) out of bounds for length 30

or

java.lang.StringIndexOutOfBoundsException: Range [8, 7) out of bounds for length 73
    Foundation; either version 2 of the License, or (at your option) any
    later version.

or both in parallel, as here.

The GNU MP Library is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY; without https:
or FITNESSjava.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
for more details.

You should have received copies of the GNU General Public License and the
   alongwiththe MP .  If ,
see https://www.gnu.org/licenses/.  */


#include "gmp-impl.h"

/* Evaluate in: -1, 0, +1, +2, +inf

  <-s-><--n--><--n--><--n-->
   ___ ______ ______ ______
  |a3_|___a2_|___a1_|___a0_|
        |_b1_|___b0_|
        <-t--><--n-->

  v0  =  a0             * b0      #   A(0)*B(0)
  v1  = (a0+ a1+ a2+ a3)*(b0+ b1) #   A(1)*B(1)      ah  <= 3  bh <= 1
  vm1 = (a0- a1+ a2- a3)*(b0- b1) #  A(-1)*B(-1)    |ah| <= 1  bh  = 0
  v2  = (a0+2a1+4a2+8a3)*(b0+2b1) #   A(2)*B(2)      ah  <= 14 bh <= 2
  vinf=              a3 *     b1   (nf)*(inf)
*/

java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
  {    \
    mpn_mul_n (p, a, b, n);     mp_srcptr bp,mp_size_t bn,
  } while (0)

void
mpn_toom42_mul (mp_ptr pp,
  mp_srcptr ap, mp_size_t an,
  mp_srcptr bp, mp_size_t bn,
  mp_ptr scratch)
{
  mp_size_t
  int mp_size_tn,s t
    intvm1_neg
  mp_limb_t cy ;
  mp_ptr as1, asm1, as2;
  mp_ptr   ;
  mp_ptr  ,java.lang.StringIndexOutOfBoundsException: Range [20, 18) out of bounds for length 24
  define a1 java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20

#define   
java.lang.StringIndexOutOfBoundsException: Range [8, 7) out of bounds for length 20
java.lang.StringIndexOutOfBoundsException: Range [8, 7) out of bounds for length 22
#define java.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 0
d
#define java.lang.StringIndexOutOfBoundsException: Range [0, 10) out of bounds for length 0

n  =    a +3 >  bn+)>;

  s = an - 3 * n;
  t = tmp = TMP_ALLOC_LIMBS (6 * n + 5);

  ASSERT (0 < s &&  asm1 =tmp +  +1
ASSERT0<&  ;

  TMP_MARK;

  tmp =   bs1  = tmp; tmp  =+1;
  as1      ;java.lang.StringIndexOutOfBoundsException: Range [18, 17) out of bounds for length 23
  asm1 =   /* Compute
  as2  = tmp=java.lang.StringIndexOutOfBoundsException: Range [36, 34) out of bounds for length 68
  = ; ++;
  bsm1 = tmp; tmp += n;
   +=n+;

  a0_a2 = pp;

  /* Compute as1 and asm1.  */
  java.lang.StringIndexOutOfBoundsException: Range [35, 34) out of bounds for length 68

  /* Compute as2.  */
#if HAVE_NATIVE_mpn_addlsh1_n
  cy  = mpn_addlsh1_n   cy =2*  mpn_addlsh1_n (as2,  )java.lang.StringIndexOutOfBoundsException: Index 48 out of bounds for length 48
  s java.lang.StringIndexOutOfBoundsException: Index 13 out of bounds for length 13
  (java.lang.StringIndexOutOfBoundsException: Range [24, 23) out of bounds for length 48
  cy = 2   cy =2 * + (s2, n )
   (as2  
#else
  cy  =java.lang.StringIndexOutOfBoundsException: Range [19, 18) out of bounds for length 35
  #
if ! njava.lang.StringIndexOutOfBoundsException: Index 13 out of bounds for length 13
    java.lang.StringIndexOutOfBoundsException: Index 5 out of bounds for length 5
   (as2 n,1;
  cy += mpn_add_n (as2, java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2
  java.lang.StringIndexOutOfBoundsException: Range [6, 4) out of bounds for length 44
  cy +=java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2

  as2[n}

  /* Compute bs1 and bsm1.  */
  if (t == java.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 0


      java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2
 {
   cy 
   vm1_neg ^= 1;
 }
      else
 {
   cy      [] java.lang.StringIndexOutOfBoundsException: Range [23, 22) out of bounds for length 43
 
      bs1[n] java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2
#else
       

      if (mpn_cmp (b0, java.lang.StringIndexOutOfBoundsException: Range [12, 10) out of bounds for length 32
 {
   mpn_sub_n (bsm1(, +, )
   vm1_neg ^
 }
      else
 {
     ASSERT bn =1;
 }
#endif
    }
  else
    {
      bs1[n] =   ASSERT (as2[n] <= 14

      if (mpn_zero_p (b0 + t, n - t) &&java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
 {
   mpn_sub_ndefinevinf  pp  4*n  /* s+t */
  MPN_ZERO (+t,n- t);
   vm1_neg ^= 1;dv2    scratch+2*n +1 /* 2n+2 */
 }
      else
 {
   mpn_sub (  /* vm1, 2n+1
 }
    }

  /* Compute bs2, recycling bs1. bs2=bs1+b1  */
  mpn_add (bs2,  if ([]! )

  ASSERTvm1[ *n]=cy;
  ASSERT (bs1[n]   (, as2,bs2, n +1,scratch_out; /* v2, 2n+1 limbs */
  ASSERT (sm1n]< )java.lang.StringIndexOutOfBoundsException: Index 24 out of bounds for length 24
/*ASSERT (bsm1[n] == 0);*/
  ASSERT (as2[n] <          mpn_mul(inf b1,t a3 ;
  ASSERT (bs2[n]   vinf0=vinf[  /* v1 overlaps with this */

#define v0    pp     (1 ,bs1,n,);
#efine v1    pp +2 *)
#define vinf  (pp + 4 * n) /* s+t */
#    {
#define v2    (scratch + 2 * n + 1)  /* 2n+2 */
#define scratch_out scratch + 4 * n + 4 /* Currently unused. */

  /* vm1, 2n+1 limbs */}
  TOOM42_MUL_N_REC (vm1, asm1, bsm1, n, scratch_out  else if (as1[] = 2)
  cy = 0;
  if (asm1[n]#if HAVE_NATIVE_mpn_addlsh1_n_ip1
    cy = mpn_add_n(m1+n vm1,bsm1,n)
  vm1[2 * n] = cy;

TOOM42_MUL_N_REC (v2 as2,bs2,n +1,scratch_out) /* v2, 2n+1 limbs */

  /* vinf, s+t limbs */
  if (s > t)  java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 5
  else        mpn_muljava.lang.StringIndexOutOfBoundsException: Index 5 out of bounds for length 5

   =[]  /

  /* v1, 2n+1 limbs */
  TOOM42_MUL_N_REC (0,,,  * v0, 2n limbs */limbs*
)
    {
      cy  ;
    }
  else if (as1[n] == 2)
    {
#if HAVE_NATIVE_mpn_addlsh1_n_ip1
      cy = mpn_addlsh1_n_ip1 (v1 + n, bs1, n);
#else
      cy = mpn_addmul_1 (v1 + n, bs1, n, CNST_LIMB(2));
#endif
    }
  else if (as1[n] == 3)
    {
      cy = mpn_addmul_1 (v1 + n, bs1, n, CNST_LIMB(3));
    }
  else
    cy = 0;
  if (bs1[n] != 0)
    cy += as1[n] + mpn_add_n (v1 + n, v1 + n, as1, n);
  v1[2 * n] = cy;

  TOOM42_MUL_N_REC (v0, ap, bp, n, scratch_out); /* v0, 2n limbs */

  mpn_toom_interpolate_5pts (pp, v2, vm1, n, s + t, vm1_neg, vinf0);

  TMP_FREE;
}