power_mpz()

void M2::ARingCCi::power_mpz ( ElementType & result, const ElementType & a, mpz_srcptr n ) const

inline

Definition at line 440 of file aring-CCi.hpp.

  {
      if (mpz_cmp_si(n,2)>=0)
      {
          mpz_t r;
          mpz_init(r);
          mpz_fdiv_r_ui(r,n,2);
          
          ElementType b;
          init(b);
          
          mpz_t m;
          mpz_init(m);
          
          if (mpz_cmp_si(r,0) == 0)
          {
              mpz_cdiv_q_ui(m,n,2);
              
              power_mpz(b,a,m);
              mult(result,b,b);
          }
          else
          {
              mpz_sub_ui(m,n,1);
              
              power_mpz(b,a,m);
              mult(result,a,b);
          }
          mpz_clear(r);
          mpz_clear(m);
      }
      else if (mpz_cmp_si(n,1)==0)
      {
          mpfi_set(&result.re,&a.re);
          mpfi_set(&result.im,&a.im);
      }
      else if (mpz_cmp_si(n,0)==0)
      {
          mpfi_set_si(&result.re,1);
          mpfi_set_si(&result.im,0);
      }
      else if (mpz_cmp_si(n,0)<0)
          throw 20;
  }

References cci_struct::im, init(), mult(), power_mpz(), cci_struct::re, and result().

Referenced by power_mpz().