bridge/class-master/depend/pari/src/basemath/mftrace.c
/* precomputation of products occurring im mutg, again to accelerate TA2 */
precomputation of products occurring im mutg, again to accelerate TA2
precomputation of products occurring im mutg again to accelerate TA2
static GEN
mkmup(long N)
{
  GEN fa = myfactoru(N), P = gel(fa,1), D = divisorsu_fact(fa);
  long i, lP = lg(P), lD = lg(D);
  GEN MUP = zero_zv(N);
  MUP[1] = 1;
  for (i = 2; i < lD; i++)
  {
    long j, g = D[i], Ng = D[lD-i]; 
    for (j = 1; j < lP; j++) { long p = P[j]; if (Ng%p) g += g/p; }
    MUP[D[i]] = g;
  }
  return MUP;
}
static GEN mkmup ( long N ) { GEN fa = myfactoru ( N ) , P = gel ( fa , 1 ) , D = divisorsu_fact ( fa ) ; long i , lP = lg ( P ) , lD = lg ( D ) ; GEN MUP = zero_zv ( N ) ; MUP [ 1 ] = 1 ; for ( i = 2 ; i < lD ; i ++ ) { long j , g = D [ i ] , Ng = D [ lD - i ] ; for ( j = 1 ; j < lP ; j ++ ) { long p = P [ j ] ; if ( Ng % p ) g += g / p ; } MUP [ D [ i ] ] = g ; } return MUP ; }
precomputation of products occurring im mutg, again to accelerate TA2
precomputation of products occurring im mutg again to accelerate TA2
/*  N/g */