init.c


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27

#include <stdio.h>
#include "gmp.h"
#include "gmp-impl.h"
#include "mpfr.h"

void
#if __STDC__
mpfr_init2 (mpfr_t x, unsigned long int p)
#else
mpfr_init2 (x, p)
     mpfr_t x;
     unsigned long int p;
#endif
{
  unsigned long xsize; 

  if (p==0) {
    printf("*** cannot initialize mpfr with precision 0\n"); exit(1);
  }

  xsize = (p - 1)/BITS_PER_MP_LIMB + 1; 

  x -> _mp_prec = p;
  x -> _mp_d = (mp_ptr) (*_mp_allocate_func) 
    (xsize * BYTES_PER_MP_LIMB);
  x -> _mp_size = xsize; 
}