bridge/class-master/depend/pari/src/basemath/mftrace.c
/* initialization for mfgaexpansion: what does not depend on cusp */
initialization for mfgaexpansion: what does not depend on cusp
initialization for mfgaexpansion : what does not depend on cusp
static GEN
mf_eisendec(GEN mf, GEN F, long prec)
{
  GEN B = liftpol_shallow(mfeisensteindec(mf, F)), v = variables_vecsmall(B);
  GEN Mvecj = obj_check(mf, MF_EISENSPACE);
  long l = lg(v), i, ord;
  if (lg(Mvecj) < 5) Mvecj = gel(Mvecj,1);
  ord = itou(gel(Mvecj,4));
  for (i = 1; i < l; i++)
    if (v[i] != 1) { B = gsubst(B, v[i], rootsof1u_cx(ord, prec)); break; }
  return B;
}
static GEN mf_eisendec ( GEN mf , GEN F , long prec ) { GEN B = liftpol_shallow ( mfeisensteindec ( mf , F ) ) , v = variables_vecsmall ( B ) ; GEN Mvecj = obj_check ( mf , MF_EISENSPACE ) ; long l = lg ( v ) , i , ord ; if ( lg ( Mvecj ) < 5 ) Mvecj = gel ( Mvecj , 1 ) ; ord = itou ( gel ( Mvecj , 4 ) ) ; for ( i = 1 ; i < l ; i ++ ) if ( v [ i ] != 1 ) { B = gsubst ( B , v [ i ] , rootsof1u_cx ( ord , prec ) ) ; break ; } return B ; }
initialization for mfgaexpansion: what does not depend on cusp
initialization for mfgaexpansion : what does not depend on cusp