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
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
|
use strict;
use warnings;
use Config;
use Math::MPFR qw(:mpfr);
use Test::More;
my($have_gmp, $have_mpz, $have_mpq, $have_mpf) = (0, 0, 0, 0);
my $long_double_formats_ok = 1;
if($^O =~ /^MSWin/ && $Config{libc} !~ /ucrt/) {
$long_double_formats_ok = 0 if WIN32_FMT_BUG;
}
eval {require Math::GMP;};
$have_gmp = 1 unless $@;
eval {require Math::GMPz;};
$have_mpz = 1 unless $@;
eval {require Math::GMPq;};
$have_mpq = 1 unless $@;
eval {require Math::GMPf;};
$have_mpf = 1 unless $@;
my $buflen = 16;
my $buf;
my $nv = sqrt(2);
if($Config{nvsize} == 8) {
Rmpfr_snprintf($buf, 7, "%.14g", $nv, $buflen);
cmp_ok($buf, 'eq', '1.4142', "sqrt 2 ok for 'double'");
Rmpfr_snprintf($buf, 8, "%a", $nv, $buflen);
like($buf, qr/^0x1\.6a0$|^0xb\.504$|^0x2\.d4$|^0x5\.a8$/, 'sqrt 2 ok for "%a" formatting');
Rmpfr_snprintf($buf, 8, "%A", $nv, $buflen);
like($buf, qr/^0X1\.6A0$|^0XB\.504$|^0X2\.D4$|^0x5\.A8$/, 'sqrt 2 ok for "%A" formatting');
}
if($Config{nvtype} eq 'long double' && $long_double_formats_ok) {
Rmpfr_snprintf($buf, 7, "%.14Lg", $nv, $buflen * 2);
cmp_ok($buf, 'eq', '1.4142', "sqrt 2 ok for 'long double'");
if(length(sqrt(2.0)) > 25) {
# IEEE 754 long double
Rmpfr_snprintf($buf, 8, "%La", $nv, $buflen * 2);
cmp_ok($buf, 'eq', '0x1.6a0', 'sqrt 2 ok for 128-bit "%La" formatting');
Rmpfr_snprintf($buf, 8, "%LA", $nv, $buflen * 2);
cmp_ok($buf, 'eq', '0X1.6A0', 'sqrt 2 ok for 128-bit "%La" formatting');
}
else {
# 80-bit extended precision long double
my($buf1, $buf2);
Rmpfr_snprintf($buf1, 8, "%La", $nv, $buflen * 2);
like($buf1, qr/^0x1.6a0$|^0x2.d41$|^0x5.a82$|^0xb.504$/, 'sqrt 2 ok for 80-bit "%La" formatting');
Rmpfr_snprintf($buf2, 8, "%LA", $nv, $buflen * 2);
cmp_ok($buf2, 'eq', uc($buf1), 'sqrt 2 ok for 80-bit "%LA" formatting');
}
}
Rmpfr_snprintf($buf, 8, "%s", 'hello world', $buflen);
cmp_ok($buf, 'eq', 'hello w', "'hello world' ok for PV");
if($have_gmp) {
Rmpfr_snprintf($buf, 7, "%Zd", Math::GMP->new(12345678), $buflen);
cmp_ok($buf, 'eq', '123456', "Math::GMP: 12345678 ok");
}
if($have_mpz) {
Rmpfr_snprintf($buf, 7, "%Zd", Math::GMPz->new(12345678), $buflen);
cmp_ok($buf, 'eq', 123456, "Math::GMPz: 12345678 ok");
}
if($have_mpq) {
Rmpfr_snprintf($buf, 4, "%Qd", Math::GMPq->new('19/21'), $buflen);
cmp_ok($buf, 'eq', '19/', "Math::GMPq: 19/21 ok");
}
if($have_mpf) {
Rmpfr_snprintf($buf, 7, "%.14Fg", sqrt(Math::GMPf->new(2)), $buflen);
cmp_ok($buf, 'eq', '1.4142', "Math::GMPf: sqrt 2 ok");
}
my $fr = Math::MPFR->new($nv);
Rmpfr_snprintf($buf, 7, "%.14RDg", $fr, $buflen);
cmp_ok($buf, 'eq', '1.4142', "Math::MPFR: sqrt 2 ok");
Rmpfr_snprintf($buf, 3, "%Pd", prec_cast(Rmpfr_get_prec($fr)), $buflen);
cmp_ok($buf, 'eq', '53', "Math::MPFR precision is '53'");
done_testing();
|
