bridge/class-master/depend/pari/src/basemath/mellininv.c
/* Compute inverse Mellin at s. K from gammamellininv OR a Vga, in which
 * case the initialization data is computed. */
Compute inverse Mellin at s. K from gammamellininv OR a Vga, in which
case the initialization data is computed.
Compute inverse Mellin at s . K from gammamellininv OR a Vga in which case the initialization data is computed .
GEN
gammamellininv(GEN K, GEN s, long m, long bitprec)
{
  pari_sp av = avma;
  GEN z, s2d;
  long d;
  if (!is_vec_t(typ(K))) pari_err_TYPE("gammamellininv",K);
  if (lg(K) != 6 || !is_vec_t(typ(GMi_get_Vga(K))))
    K = gammamellininvinit(K, m, bitprec);
  d = lg(GMi_get_Vga(K))-1;
  s2d = gpow(s, gdivgs(gen_2, d), nbits2prec(bitprec));
  if (dblmodulus(s2d) < GMi_get_tmax(K, bitprec))
    z = Kderivsmall(K, s, s2d, bitprec);
  else
    z = Kderivlarge(K, s, s2d, bitprec);
  return gerepileupto(av, z);
}
GEN gammamellininv ( GEN K , GEN s , long m , long bitprec ) { pari_sp av = avma ; GEN z , s2d ; long d ; if ( ! is_vec_t ( typ ( K ) ) ) pari_err_TYPE ( " " , K ) ; if ( lg ( K ) != 6 || ! is_vec_t ( typ ( GMi_get_Vga ( K ) ) ) ) K = gammamellininvinit ( K , m , bitprec ) ; d = lg ( GMi_get_Vga ( K ) ) - 1 ; s2d = gpow ( s , gdivgs ( gen_2 , d ) , nbits2prec ( bitprec ) ) ; if ( dblmodulus ( s2d ) < GMi_get_tmax ( K , bitprec ) ) z = Kderivsmall ( K , s , s2d , bitprec ) ; else z = Kderivlarge ( K , s , s2d , bitprec ) ; return gerepileupto ( av , z ) ; }
Compute inverse Mellin at s. K from gammamellininv OR a Vga, in which
 case the initialization data is computed.
Compute inverse Mellin at s . K from gammamellininv OR a Vga in which case the initialization data is computed .