OLD | NEW |
| (Empty) |
1 /* mpz_fib2_ui -- calculate Fibonacci numbers. | |
2 | |
3 Copyright 2001 Free Software Foundation, Inc. | |
4 | |
5 This file is part of the GNU MP Library. | |
6 | |
7 The GNU MP Library is free software; you can redistribute it and/or modify | |
8 it under the terms of the GNU Lesser General Public License as published by | |
9 the Free Software Foundation; either version 3 of the License, or (at your | |
10 option) any later version. | |
11 | |
12 The GNU MP Library is distributed in the hope that it will be useful, but | |
13 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY | |
14 or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public | |
15 License for more details. | |
16 | |
17 You should have received a copy of the GNU Lesser General Public License | |
18 along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */ | |
19 | |
20 #include <stdio.h> | |
21 #include "gmp.h" | |
22 #include "gmp-impl.h" | |
23 | |
24 | |
25 void | |
26 mpz_fib2_ui (mpz_ptr fn, mpz_ptr fnsub1, unsigned long n) | |
27 { | |
28 mp_ptr fp, f1p; | |
29 mp_size_t size; | |
30 | |
31 size = MPN_FIB2_SIZE (n); | |
32 MPZ_REALLOC (fn, size); | |
33 MPZ_REALLOC (fnsub1, size); | |
34 fp = PTR (fn); | |
35 f1p = PTR (fnsub1); | |
36 | |
37 size = mpn_fib2_ui (fp, f1p, n); | |
38 | |
39 SIZ(fn) = size - (n == 0); | |
40 SIZ(fnsub1) = size - (f1p[size-1] == 0); | |
41 } | |
OLD | NEW |