# mpfr_root is deprecated in favour of mpfr_rootn_ui in mpfr-4.0.0
# The 2 variants should produce identical results, except for the
# nth root of -0 when n is zero.
use strict;
use warnings;
use Math::MPFR qw(:mpfr);
my $old = 0;
if(MPFR_VERSION_MAJOR > 3) {
  *_ROOT = \&Rmpfr_rootn_ui;
}
else {
  *_ROOT = \&Rmpfr_root;
  $old = 1;
}
print "1..57\n";
my($inex1, $inex2, $check);
my($rop1, $rop2) = (Rmpfr_init(), Rmpfr_init());
my $root = 2;
my $op = Math::MPFR->new(17);
$inex1 = Rmpfr_rec_sqrt($rop1, $op,        MPFR_RNDN);
$inex2 = Rmpfr_rec_root($rop2, $op, $root, MPFR_RNDN);
if($inex1 == $inex2) {print "ok 1\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 1\n";
}
if($rop1 == $rop2) {print "ok 2\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\n";
  print "not ok 2\n";
}
$op *= -1; # -17
Rmpfr_clear_nanflag();
$inex1 = Rmpfr_rec_sqrt($rop1, $op,        MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 3\n"}
else {print "not ok 3\n"}
Rmpfr_clear_nanflag();
$inex2 = Rmpfr_rec_root($rop2, $op, $root, MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 4\n"}
else {print "not ok 4\n"}
if($inex1 == $inex2) {print "ok 5\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 5\n";
}
if(Rmpfr_nan_p($rop1) && Rmpfr_nan_p($rop2)) {print "ok 6\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\n";
  print "not ok 6\n";
}
$root = 5;
$inex1 = _ROOT         ($rop1, $op, $root, MPFR_RNDN);
$inex2 = Rmpfr_rec_root($rop2, $op, $root, MPFR_RNDN);
$check = $rop1 * $rop2;
if($check < 1.0000001 && $check > 0.9999999) {print "ok 7\n"}
else {
  warn "\n \$check: $check\n";
  print "not ok 7\n";
}
# see tests 50-57 for root = 0.
## $op is +/- 0 ##
## root is even, root is odd, root is 0
my $pzero = Math::MPFR->new(0);
my $nzero = $pzero * -1;
# $root is 5
Rmpfr_clear_divby0();
$inex1 = _ROOT         ($rop1, $pzero, $root, MPFR_RNDN);
$inex2 = Rmpfr_rec_root($rop2, $pzero, $root, MPFR_RNDN);
if($inex1 == $inex2 && Rmpfr_divby0_p()) {print "ok 8\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\ndivby0: ", Rmpfr_divby0_p(), "\n";
  print "not ok 8\n";
}
Rmpfr_clear_divby0();
if($rop1 == 1 / $rop2 && !Rmpfr_divby0_p() && Rmpfr_signbit($rop1) == Rmpfr_signbit($rop2)) {print "ok 9\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\ndivby0: ", Rmpfr_divby0_p(), "\n",
       "signbits \$rop1: ", Rmpfr_signbit($rop1), " \$rop2: ", Rmpfr_signbit($rop2), "\n";
  print "not ok 9\n";
}
$inex1 = _ROOT         ($rop1, $nzero, $root, MPFR_RNDN);
$inex2 = Rmpfr_rec_root($rop2, $nzero, $root, MPFR_RNDN);
if($inex1 == $inex2 && Rmpfr_divby0_p()) {print "ok 10\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\ndivby0: ", Rmpfr_divby0_p(), "\n";
  print "not ok 10\n";
}
Rmpfr_clear_divby0();
if($rop1 == 1 / $rop2 && !Rmpfr_divby0_p() && Rmpfr_signbit($rop1) == Rmpfr_signbit($rop2)) {print "ok 11\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\ndivby0: ", Rmpfr_divby0_p(), "\n",
       "signbits \$rop1: ", Rmpfr_signbit($rop1), " \$rop2: ", Rmpfr_signbit($rop2), "\n";
  print "not ok 11\n";
}
$inex1 = _ROOT         ($rop1, $pzero, $root - 1, MPFR_RNDN);
$inex2 = Rmpfr_rec_root($rop2, $pzero, $root - 1, MPFR_RNDN);
if($inex1 == $inex2 && Rmpfr_divby0_p()) {print "ok 12\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\ndivby0: ", Rmpfr_divby0_p(), "\n";
  print "not ok 12\n";
}
Rmpfr_clear_divby0();
if($rop1 == 1 / $rop2 && !Rmpfr_divby0_p() && Rmpfr_signbit($rop1) == Rmpfr_signbit($rop2)) {print "ok 13\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\ndivby0: ", Rmpfr_divby0_p(), "\n",
       "signbits \$rop1: ", Rmpfr_signbit($rop1), " \$rop2: ", Rmpfr_signbit($rop2), "\n";
  print "not ok 13\n";
}
$inex1 = _ROOT         ($rop1, $nzero, $root - 1, MPFR_RNDN);
$inex2 = Rmpfr_rec_root($rop2, $nzero, $root - 1, MPFR_RNDN);
if($inex1 == $inex2 && Rmpfr_divby0_p()) {print "ok 14\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\ndivby0: ", Rmpfr_divby0_p(), "\n";
  print "not ok 14\n";
}
Rmpfr_clear_divby0();
if($old) {
  if($rop1 == 1 / $rop2 && !Rmpfr_divby0_p() && Rmpfr_signbit($rop1) != Rmpfr_signbit($rop2)) {print "ok 15\n"}
  else {
    warn "\n \$rop1: $rop1\n \$rop2: $rop2\ndivby0: ", Rmpfr_divby0_p(), "\n",
         "signbits \$rop1: ", Rmpfr_signbit($rop1), " \$rop2: ", Rmpfr_signbit($rop2), "\n";
    print "not ok 15\n";
  }
}
else {
  if($rop1 == 1 / $rop2 && !Rmpfr_divby0_p() && Rmpfr_signbit($rop1) == Rmpfr_signbit($rop2)) {print "ok 15\n"}
  else {
    warn "\n \$rop1: $rop1\n \$rop2: $rop2\ndivby0: ", Rmpfr_divby0_p(), "\n",
         "signbits \$rop1: ", Rmpfr_signbit($rop1), " \$rop2: ", Rmpfr_signbit($rop2), "\n";
    print "not ok 15\n";
  }
}
Rmpfr_clear_divby0();
Rmpfr_clear_nanflag();
$inex1 = _ROOT($rop1, $pzero, 0, MPFR_RNDN);
if(Rmpfr_nanflag_p() && !Rmpfr_divby0_p()) {print "ok 16\n"}
else {print "not ok 16\n"}
Rmpfr_clear_divby0();
Rmpfr_clear_nanflag();
$inex2 = Rmpfr_rec_root($rop2, $nzero, 0, MPFR_RNDN);
if(Rmpfr_nanflag_p() && !Rmpfr_divby0_p()) {print "ok 17\n"}
else {print "not ok 17\n"}
Rmpfr_clear_divby0();
Rmpfr_clear_nanflag();
if($inex1 == $inex2) {print "ok 18\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 18\n";
}
if(Rmpfr_nan_p($rop1) && Rmpfr_nan_p($rop2)) {print "ok 19\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\n";
  print "not ok 19\n";
}
## $op is +/- Inf ##
## root is even, root is odd, root is 0.
$root = 2;
Rmpfr_set_d($op, 999**(999**999), MPFR_RNDN);
$inex1 = Rmpfr_rec_sqrt($rop1, $op,        MPFR_RNDN);
$inex2 = Rmpfr_rec_root($rop2, $op, $root, MPFR_RNDN);
if($inex1 == $inex2) {print "ok 20\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 20\n";
}
if($rop1 == $rop2) {print "ok 21\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\n";
  print "not ok 21\n";
}
$op *= -1; # -Inf
$inex1 = Rmpfr_rec_sqrt($rop1, $op,        MPFR_RNDN);
$inex2 = Rmpfr_rec_root($rop2, $op, $root, MPFR_RNDN);
if($inex1 == $inex2) {print "ok 22\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 22\n";
}
if(Rmpfr_nan_p($rop1) && Rmpfr_nan_p($rop2)) {print "ok 23\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\n";
  print "not ok 23\n";
}
$root = 5;
$inex1 = _ROOT         ($rop1, $op, $root, MPFR_RNDN);
$inex2 = Rmpfr_rec_root($rop2, $op, $root, MPFR_RNDN);
if($inex1 == $inex2) {print "ok 24\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 24\n";
}
if($rop1 == 1 / $rop2) {print "ok 25\n"}
else {
  warn "\n \$check: $check\n";
  print "not ok 25\n";
}
$op *= -1; # +Inf
$inex1 = _ROOT         ($rop1, $op, $root, MPFR_RNDN);
$inex2 = Rmpfr_rec_root($rop2, $op, $root, MPFR_RNDN);
if($inex1 == $inex2) {print "ok 26\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 26\n";
}
if($rop1 == 1 / $rop2) {print "ok 27\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\n";
  print "not ok 27\n";
}
$inex1 = _ROOT         ($rop1, $op, $root - 1, MPFR_RNDN);
$inex2 = Rmpfr_rec_root($rop2, $op, $root - 1, MPFR_RNDN);
if($inex1 == $inex2) {print "ok 28\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 28\n";
}
if($rop1 == 1 / $rop2) {print "ok 29\n"}
else {
  warn "\n \$check: $check\n";
  print "not ok 29\n";
}
Rmpfr_clear_nanflag();
$inex1 = _ROOT($rop1, $op, 0, MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 30\n"}
else {print "not ok 30\n"}
Rmpfr_clear_nanflag();
$inex2 = Rmpfr_rec_root($rop2, $op, 0, MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 31\n"}
else {print "not ok 31\n"}
if($inex1 == $inex2) {print "ok 32\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 32\n";
}
if(Rmpfr_nan_p($rop1) && Rmpfr_nan_p($rop2) && Rmpfr_nanflag_p()) {print "ok 33\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\n ", Rmpfr_nanflag_p(), "\n";
  print "not ok 33\n";
}
$op *= -1; # -Inf
Rmpfr_clear_nanflag();
$inex1 = _ROOT($rop1, $op, 0, MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 34\n"}
else {print "not ok 34\n"}
Rmpfr_clear_nanflag();
$inex2 = Rmpfr_rec_root($rop2, $op, 0, MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 35\n"}
else {print "not ok 35\n"}
if($inex1 == $inex2) {print "ok 36\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 36\n";
}
if(Rmpfr_nan_p($rop1) && Rmpfr_nan_p($rop2) && Rmpfr_nanflag_p()) {print "ok 37\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\n ", Rmpfr_nanflag_p(), "\n";
  print "not ok 37\n";
}
## $op is NaN
## root is even, root is odd, root is 0
Rmpfr_set_nan($op);
$root = 2;
Rmpfr_clear_nanflag();
$inex1 = Rmpfr_rec_sqrt($rop1, $op,        MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 38\n"}
else {print "not ok 38\n"}
Rmpfr_clear_nanflag();
$inex2 = Rmpfr_rec_root($rop2, $op, $root, MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 39\n"}
else {print "not ok 39\n"}
Rmpfr_clear_nanflag();
if($inex1 == $inex2) {print "ok 40\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 40\n";
}
if(Rmpfr_nan_p($rop1) && Rmpfr_nan_p($rop2)) {print "ok 41\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\n ", Rmpfr_nanflag_p(), "\n";
  print "not ok 41\n";
}
$root = 5;
Rmpfr_clear_nanflag();
$inex1 = _ROOT($rop1, $op, $root, MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 42\n"}
else {print "not ok 42\n"}
Rmpfr_clear_nanflag();
$inex2 = Rmpfr_rec_root($rop2, $op, $root, MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 43\n"}
else {print "not ok 43\n"}
Rmpfr_clear_nanflag();
if($inex1 == $inex2) {print "ok 44\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 44\n";
}
if(Rmpfr_nan_p($rop1) && Rmpfr_nan_p($rop2)) {print "ok 45\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\n ", Rmpfr_nanflag_p(), "\n";
  print "not ok 45\n";
}
$root = 0;
Rmpfr_clear_nanflag();
$inex1 = _ROOT($rop1, $op, $root, MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 46\n"}
else {print "not ok 46\n"}
Rmpfr_clear_nanflag();
$inex2 = Rmpfr_rec_root($rop2, $op, $root, MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 47\n"}
else {print "not ok 47\n"}
Rmpfr_clear_nanflag();
if($inex1 == $inex2) {print "ok 48\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 48\n";
}
if(Rmpfr_nan_p($rop1) && Rmpfr_nan_p($rop2)) {print "ok 49\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\n ", Rmpfr_nanflag_p(), "\n";
  print "not ok 49\n";
}
Rmpfr_set_ui($op, 42, MPFR_RNDN);
# Now check regular values for root = 0 (which was missed earlier).
Rmpfr_clear_nanflag();
$inex1 = _ROOT($rop1, $op, 0, MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 50\n"}
else {print "not ok 50\n"}
Rmpfr_clear_nanflag();
$inex2 = Rmpfr_rec_root($rop2, $op, 0, MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 51\n"}
else {print "not ok 51\n"}
if($inex1 == $inex2) {print "ok 52\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 52\n";
}
if(Rmpfr_nan_p($rop1) && Rmpfr_nan_p($rop2)) {print "ok 53\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\n ", Rmpfr_nanflag_p(), "\n";
  print "not ok 53\n";
}
$op *= -1; # -42
Rmpfr_clear_nanflag();
$inex1 = _ROOT($rop1, $op, 0, MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 54\n"}
else {print "not ok 54\n"}
Rmpfr_clear_nanflag();
$inex2 = Rmpfr_rec_root($rop2, $op, 0, MPFR_RNDN);
if(Rmpfr_nanflag_p()) {print "ok 55\n"}
else {print "not ok 55\n"}
if($inex1 == $inex2) {print "ok 56\n"}
else {
  warn "\n \$inex1: $inex1\n \$inex2: $inex2\n";
  print "not ok 56\n";
}
if(Rmpfr_nan_p($rop1) && Rmpfr_nan_p($rop2)) {print "ok 57\n"}
else {
  warn "\n \$rop1: $rop1\n \$rop2: $rop2\n ", Rmpfr_nanflag_p(), "\n";
  print "not ok 57\n";
}