gnu_mpc 0.8.0

data/spec/proj_spec.rb

@@ -0,0 +1,20 @@
+require File.join(File.dirname(__FILE__), 'spec_helper')
+
+describe MPC, "#proj" do
+  it "should calculate the projection correctly" do
+    MPC.new([GMP::F(-1), -GMP::F(0)]).proj.real.should eq -1
+    MPC.new([GMP::F(-1), -GMP::F(0)]).proj.imag.should eq 0
+
+    MPC.new([-GMP::F(0),  GMP::F(-1)]).proj.real.should eq 0
+    MPC.new([-GMP::F(0),  GMP::F(-1)]).proj.imag.should eq -1
+
+    MPC.new([-GMP::F(0), -GMP::F(0)]).proj.real.should eq 0
+    MPC.new([-GMP::F(0), -GMP::F(0)]).proj.imag.should eq 0
+
+    MPC.new([-GMP::F(0),  GMP::F(1)]).proj.real.should eq 0
+    MPC.new([-GMP::F(0),  GMP::F(1)]).proj.imag.should eq 1
+
+    MPC.new([ GMP::F(1),  GMP::F(0)]).proj.real.should eq 1
+    MPC.new([ GMP::F(1),  GMP::F(0)]).proj.imag.should eq 0
+  end
+end

data/spec/real_spec.rb

@@ -0,0 +1,23 @@
+require File.join(File.dirname(__FILE__), 'spec_helper')
+
+describe MPC, "#real" do
+  it "should return correctly for an MPC generated with a Fixnum" do
+    MPC.new(0).real.should eq GMP::F(0)
+    MPC.new(1).real.should eq GMP::F(1)
+    MPC.new([1, 1]).real.should eq GMP::F(1)
+  end
+
+  it "should return correctly for an MPC generated with a GMP::Z" do
+    MPC.new(GMP::Z(0)).real.should eq GMP::F(0)
+    MPC.new(GMP::Z(1)).real.should eq GMP::F(1)
+    MPC.new(GMP::Z(2**64)).real.should eq GMP::F(2**64)
+    MPC.new([GMP::Z(1), GMP::Z(1)]).real.should eq GMP::F(1)
+  end
+
+  it "should return correctly for an MPC generated with a Float" do
+    MPC.new(0.0).real.should eq GMP::F(0)
+    MPC.new(1.0).real.should eq GMP::F(1)
+    MPC.new([1.0, 1.0]).real.should eq GMP::F(1)
+    MPC.new([2*Math::PI, Math::PI]).real.should eq GMP::F(2*Math::PI)
+  end
+end

data/spec/rounding_spec.rb

@@ -0,0 +1,199 @@
+require File.join(File.dirname(__FILE__), 'spec_helper')
+
+describe MPC::Rnd, '#mode' do
+  it 'should return MPC::MPC_RNDNN.mode correctly' do
+    MPC::MPC_RNDNN.mode.should eq 0
+  end
+
+  it 'should return MPC::MPC_RNDZN.mode correctly' do
+    MPC::MPC_RNDZN.mode.should eq 1
+  end
+
+  it 'should return MPC::MPC_RNDUN.mode correctly' do
+    MPC::MPC_RNDUN.mode.should eq 2
+  end
+
+  it 'should return MPC::MPC_RNDDN.mode correctly' do
+    MPC::MPC_RNDDN.mode.should eq 3
+  end
+
+  it 'should return MPC::MPC_RNDNZ.mode correctly' do
+    MPC::MPC_RNDNZ.mode.should eq 16
+  end
+
+  it 'should return MPC::MPC_RNDZZ.mode correctly' do
+    MPC::MPC_RNDZZ.mode.should eq 17
+  end
+
+  it 'should return MPC::MPC_RNDUZ.mode correctly' do
+    MPC::MPC_RNDUZ.mode.should eq 18
+  end
+
+  it 'should return MPC::MPC_RNDDZ.mode correctly' do
+    MPC::MPC_RNDDZ.mode.should eq 19
+  end
+
+  it 'should return MPC::MPC_RNDNU.mode correctly' do
+    MPC::MPC_RNDNU.mode.should eq 32
+  end
+
+  it 'should return MPC::MPC_RNDNU.mode correctly' do
+    MPC::MPC_RNDZU.mode.should eq 33
+  end
+
+  it 'should return MPC::MPC_RNDNU.mode correctly' do
+    MPC::MPC_RNDUU.mode.should eq 34
+  end
+
+  it 'should return MPC::MPC_RNDNU.mode correctly' do
+    MPC::MPC_RNDDU.mode.should eq 35
+  end
+
+  it 'should return MPC::MPC_RNDND.mode correctly' do
+    MPC::MPC_RNDND.mode.should eq 48
+  end
+
+  it 'should return MPC::MPC_RNDND.mode correctly' do
+    MPC::MPC_RNDZD.mode.should eq 49
+  end
+
+  it 'should return MPC::MPC_RNDND.mode correctly' do
+    MPC::MPC_RNDUD.mode.should eq 50
+  end
+
+  it 'should return MPC::MPC_RNDND.mode correctly' do
+    MPC::MPC_RNDDD.mode.should eq 51
+  end
+end
+
+describe MPC::Rnd, '#name' do
+  it 'should return MPC::MPC_RNDNN.name correctly' do
+    MPC::MPC_RNDNN.name.should eq 'MPC_RNDNN'
+  end
+
+  it 'should return MPC::MPC_RNDZN.name correctly' do
+    MPC::MPC_RNDZN.name.should eq 'MPC_RNDZN'
+  end
+
+  it 'should return MPC::MPC_RNDUN.name correctly' do
+    MPC::MPC_RNDUN.name.should eq 'MPC_RNDUN'
+  end
+
+  it 'should return MPC::MPC_RNDDN.name correctly' do
+    MPC::MPC_RNDDN.name.should eq 'MPC_RNDDN'
+  end
+
+  it 'should return MPC::MPC_RNDNZ.name correctly' do
+    MPC::MPC_RNDNZ.name.should eq 'MPC_RNDNZ'
+  end
+
+  it 'should return MPC::MPC_RNDZZ.name correctly' do
+    MPC::MPC_RNDZZ.name.should eq 'MPC_RNDZZ'
+  end
+
+  it 'should return MPC::MPC_RNDUZ.name correctly' do
+    MPC::MPC_RNDUZ.name.should eq 'MPC_RNDUZ'
+  end
+
+  it 'should return MPC::MPC_RNDDZ.name correctly' do
+    MPC::MPC_RNDDZ.name.should eq 'MPC_RNDDZ'
+  end
+
+  it 'should return MPC::MPC_RNDNU.name correctly' do
+    MPC::MPC_RNDNU.name.should eq 'MPC_RNDNU'
+  end
+
+  it 'should return MPC::MPC_RNDZU.name correctly' do
+    MPC::MPC_RNDZU.name.should eq 'MPC_RNDZU'
+  end
+
+  it 'should return MPC::MPC_RNDUU.name correctly' do
+    MPC::MPC_RNDUU.name.should eq 'MPC_RNDUU'
+  end
+
+  it 'should return MPC::MPC_RNDDU.name correctly' do
+    MPC::MPC_RNDDU.name.should eq 'MPC_RNDDU'
+  end
+
+  it 'should return MPC::MPC_RNDND.name correctly' do
+    MPC::MPC_RNDND.name.should eq 'MPC_RNDND'
+  end
+
+  it 'should return MPC::MPC_RNDZD.name correctly' do
+    MPC::MPC_RNDZD.name.should eq 'MPC_RNDZD'
+  end
+
+  it 'should return MPC::MPC_RNDUD.name correctly' do
+    MPC::MPC_RNDUD.name.should eq 'MPC_RNDUD'
+  end
+
+  it 'should return MPC::MPC_RNDDD.name correctly' do
+    MPC::MPC_RNDDD.name.should eq 'MPC_RNDDD'
+  end
+end
+
+describe MPC::Rnd, '#ieee754' do
+  it 'should return MPC::MPC_RNDNN.ieee754 correctly' do
+    MPC::MPC_RNDNN.ieee754.should eq '(roundTiesToEven,roundTiesToEven)'
+  end
+
+  it 'should return MPC::MPC_RNDZN.ieee754 correctly' do
+    MPC::MPC_RNDZN.ieee754.should eq '(roundTowardZero,roundTiesToEven)'
+  end
+
+  it 'should return MPC::MPC_RNDUN.ieee754 correctly' do
+    MPC::MPC_RNDUN.ieee754.should eq '(roundTowardPositive,roundTiesToEven)'
+  end
+
+  it 'should return MPC::MPC_RNDDN.ieee754 correctly' do
+    MPC::MPC_RNDDN.ieee754.should eq '(roundTowardNegative,roundTiesToEven)'
+  end
+
+  it 'should return MPC::MPC_RNDNZ.ieee754 correctly' do
+    MPC::MPC_RNDNZ.ieee754.should eq '(roundTiesToEven,roundTowardZero)'
+  end
+
+  it 'should return MPC::MPC_RNDZZ.ieee754 correctly' do
+    MPC::MPC_RNDZZ.ieee754.should eq '(roundTowardZero,roundTowardZero)'
+  end
+
+  it 'should return MPC::MPC_RNDUZ.ieee754 correctly' do
+    MPC::MPC_RNDUZ.ieee754.should eq '(roundTowardPositive,roundTowardZero)'
+  end
+
+  it 'should return MPC::MPC_RNDDZ.ieee754 correctly' do
+    MPC::MPC_RNDDZ.ieee754.should eq '(roundTowardNegative,roundTowardZero)'
+  end
+
+  it 'should return MPC::MPC_RNDNU.ieee754 correctly' do
+    MPC::MPC_RNDNU.ieee754.should eq '(roundTiesToEven,roundTowardPositive)'
+  end
+
+  it 'should return MPC::MPC_RNDZU.ieee754 correctly' do
+    MPC::MPC_RNDZU.ieee754.should eq '(roundTowardZero,roundTowardPositive)'
+  end
+
+  it 'should return MPC::MPC_RNDUU.ieee754 correctly' do
+    MPC::MPC_RNDUU.ieee754.should eq '(roundTowardPositive,roundTowardPositive)'
+  end
+
+  it 'should return MPC::MPC_RNDDU.ieee754 correctly' do
+    MPC::MPC_RNDDU.ieee754.should eq '(roundTowardNegative,roundTowardPositive)'
+  end
+
+  it 'should return MPC::MPC_RNDND.ieee754 correctly' do
+    MPC::MPC_RNDND.ieee754.should eq '(roundTiesToEven,roundTowardNegative)'
+  end
+
+  it 'should return MPC::MPC_RNDZD.ieee754 correctly' do
+    MPC::MPC_RNDZD.ieee754.should eq '(roundTowardZero,roundTowardNegative)'
+  end
+
+  it 'should return MPC::MPC_RNDUD.ieee754 correctly' do
+    MPC::MPC_RNDUD.ieee754.should eq '(roundTowardPositive,roundTowardNegative)'
+  end
+
+  it 'should return MPC::MPC_RNDDD.ieee754 correctly' do
+    MPC::MPC_RNDDD.ieee754.should eq '(roundTowardNegative,roundTowardNegative)'
+  end
+end

data/spec/sin_spec.rb

@@ -0,0 +1,66 @@
+require File.join(File.dirname(__FILE__), 'spec_helper')
+
+# All tests adapted from MPC 1.0.1's tests/sin.dat
+describe MPC, '#sin' do
+  it 'should calculate the sine of a pure real argument' do
+    data = [
+      ["0x4787C62AC28Bp-48",   MPC.new([GMP::F(-6), GMP::F(0)])],
+      ["0xC1BDCEEEE0F57p-52",  MPC.new([GMP::F(-4), GMP::F(0)])],
+      ["-0xE8C7B7568DA23p-52", MPC.new([GMP::F(-2), GMP::F(0)])],
+      ["-0xD76AA47848677p-52", MPC.new([GMP::F(-1), GMP::F(0)])],
+      ["0xD76AA47848677p-52",  MPC.new([GMP::F( 1), GMP::F(0)])],
+      ["0xE8C7B7568DA23p-52",  MPC.new([GMP::F( 2), GMP::F(0)])],
+      ["-0xC1BDCEEEE0F57p-52", MPC.new([GMP::F( 4), GMP::F(0)])],
+      ["-0x4787C62AC28Bp-48",  MPC.new([GMP::F( 6), GMP::F(0)])]
+    ]
+    data.each do |expected, input|
+      actual = input.sin
+      actual.real.should eq GMP::F.new(expected, 53, 16)
+      actual.imag.should eq GMP::F(0)
+    end
+  end
+
+  it 'should calculate the sine of a pure imaginary argument' do
+    data = [
+      ["-0x1936D22F67C805p-45", MPC.new([GMP::F(0), GMP::F(-6)])],
+      ["-0x1B4A3803703631p-48", MPC.new([GMP::F(0), GMP::F(-4)])],
+      ["-0x1D03CF63B6E19Fp-51", MPC.new([GMP::F(0), GMP::F(-2)])],
+      ["-0x966CFE2275CC1p-51",  MPC.new([GMP::F(0), GMP::F(-1)])],
+      ["0x966CFE2275CC1p-51",   MPC.new([GMP::F(0), GMP::F( 1)])],
+      ["0x1D03CF63B6E19Fp-51",  MPC.new([GMP::F(0), GMP::F( 2)])],
+      ["0x1B4A3803703631p-48",  MPC.new([GMP::F(0), GMP::F( 4)])],
+      ["0x1936D22F67C805p-45",  MPC.new([GMP::F(0), GMP::F( 6)])]
+    ]
+    data.each do |expected, input|
+      actual = input.sin
+      actual.imag.should eq GMP::F.new(expected, 53, 16)
+      actual.real.should eq GMP::F(0)
+    end
+  end
+
+  it 'should calculate the sine of ieee-754 double precision numbers' do
+    data = [
+      [[514],                            ["-0x11B7CA26B51951p-97", 53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDNN],
+      [[514],                            ["-0x11B7CA26B51951p-97", 53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDNZ],
+      [[514],                            ["-0x11B7CA26B51951p-97", 53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDNU],
+      [[514],                            ["-0x8DBE5135A8CA9p-96",  53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDND],
+      [["0x100FFFFFFFFFFFp-43", 53, 16], ["-0x11B7CA26B51951p-97", 53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDZN],
+      [["0x100FFFFFFFFFFFp-43", 53, 16], ["-0x11B7CA26B51951p-97", 53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDZZ],
+      [["0x100FFFFFFFFFFFp-43", 53, 16], ["-0x11B7CA26B51951p-97", 53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDZU],
+      [["0x100FFFFFFFFFFFp-43", 53, 16], ["-0x8DBE5135A8CA9p-96",  53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDZD],
+      [[514],                            ["-0x11B7CA26B51951p-97", 53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDUN],
+      [[514],                            ["-0x11B7CA26B51951p-97", 53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDUZ],
+      [[514],                            ["-0x11B7CA26B51951p-97", 53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDUU],
+      [[514],                            ["-0x8DBE5135A8CA9p-96",  53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDUD],
+      [["0x100FFFFFFFFFFFp-43", 53, 16], ["-0x11B7CA26B51951p-97", 53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDDN],
+      [["0x100FFFFFFFFFFFp-43", 53, 16], ["-0x11B7CA26B51951p-97", 53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDDZ],
+      [["0x100FFFFFFFFFFFp-43", 53, 16], ["-0x11B7CA26B51951p-97", 53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDDU],
+      [["0x100FFFFFFFFFFFp-43", 53, 16], ["-0x8DBE5135A8CA9p-96",  53, 16], ["0x3243F6A8885A3p-49", 53, 16], ["-0x1BBDD1808C59A3p-50", 53, 16], MPC::MPC_RNDDD]
+    ]
+    data.each do |expected_real, expected_imag, input_real, input_imag, rounding_mode|
+      actual = MPC.new([GMP::F.new(*input_real), GMP::F.new(*input_imag)]).sin(rounding_mode)
+      actual.real.should eq GMP::F.new(*expected_real)
+      actual.imag.should eq GMP::F.new(*expected_imag)
+    end
+  end
+end

data/spec/sinh_spec.rb

@@ -0,0 +1,65 @@
+require File.join(File.dirname(__FILE__), 'spec_helper')
+
+# All tests adapted from MPC 1.0.1's tests/sinh.dat
+describe MPC, '#sinh' do
+  it 'should calculate the hyperbolic sine of a pure real argument' do
+    data = [
+      ["-0x12cd9fc44eb98p-48",   MPC.new([GMP::F(-1, 7), GMP::F(0, 7)])],
+      [ "0x12cd9fc44eb98p-48",   MPC.new([GMP::F( 1, 7), GMP::F(0, 7)])]
+    ]
+    data.each do |expected, input|
+      actual = input.sinh(MPC::MPC_RNDNN, 50, 50)
+      actual.real.should eq GMP::F.new(expected, 50, 16)
+      actual.imag.should eq GMP::F(0)
+    end
+  end
+
+  it 'should calculate the hyperbolic sine of a pure real argument, using hash arguments' do
+    data = [
+      ["-0x12cd9fc44eb98p-48",   MPC.new([GMP::F(-1, 7), GMP::F(0, 7)])],
+      [ "0x12cd9fc44eb98p-48",   MPC.new([GMP::F( 1, 7), GMP::F(0, 7)])]
+    ]
+    data.each do |expected, input|
+      actual = input.sinh(:rounding_mode => MPC::MPC_RNDNN, :precision => 50)
+      actual.real.should eq GMP::F.new(expected, 50, 16)
+      actual.imag.should eq GMP::F(0)
+    end
+  end
+
+  it 'should calculate the hyperbolic sine of a pure imaginary argument' do
+    data = [
+      ["-0xd76aa47848678p-52",   MPC.new([GMP::F(0, 7), GMP::F(-1, 7)])],
+      [ "0xd76aa47848678p-52",   MPC.new([GMP::F(0, 7), GMP::F( 1, 7)])]
+    ]
+    data.each do |expected, input|
+      actual = input.sinh(MPC::MPC_RNDNN, 50, 50)
+      actual.real.should eq GMP::F(0)
+      actual.imag.should eq GMP::F.new(expected, 50, 16)
+    end
+  end
+
+  it 'should calculate the hyperbolic sine of an argument with +1 and -1' do
+    data = [
+      ["-0xa28cfec023fc8p-52", "-0x14c67b74f6cc5p-48",  MPC.new([GMP::F(-1, 7), GMP::F(-1, 7)])],
+      ["-0xa28cfec023fc8p-52",  "0x14c67b74f6cc5p-48",  MPC.new([GMP::F(-1, 7), GMP::F( 1, 7)])],
+      [ "0xa28cfec023fc8p-52", "-0x14c67b74f6cc5p-48",  MPC.new([GMP::F( 1, 7), GMP::F(-1, 7)])],
+      [ "0xa28cfec023fc8p-52",  "0x14c67b74f6cc5p-48",  MPC.new([GMP::F( 1, 7), GMP::F( 1, 7)])],
+    ]
+    data.each do |expected_real, expected_imag, input|
+      actual = input.sinh(MPC::MPC_RNDNN, 50, 50)
+      actual.real.should eq GMP::F.new(expected_real, 50, 16)
+      actual.imag.should eq GMP::F.new(expected_imag, 50, 16)
+    end
+  end
+
+  it 'should calculate the hyperbolic sine of ieee-754 double precision numbers' do
+    data = [
+      [["0xF48D4FDF29C53p-105", 53, 16], [2], ["0x15124271980435p-52", 53, 16], ["0x3243F6A8885A3p-49", 53, 16], MPC::MPC_RNDNN],
+    ]
+    data.each do |expected_real, expected_imag, input_real, input_imag, rounding_mode|
+      actual = MPC.new([GMP::F.new(*input_real), GMP::F.new(*input_imag)]).sinh(rounding_mode)
+      actual.real.should eq GMP::F.new(*expected_real)
+      actual.imag.should eq GMP::F.new(*expected_imag)
+    end
+  end
+end

data/spec/spec_helper.rb

@@ -0,0 +1 @@
+require File.join(File.dirname(__FILE__), '..', 'lib', 'mpc')

data/spec/sqr_spec.rb

@@ -0,0 +1,96 @@
+require File.join(File.dirname(__FILE__), 'spec_helper')
+
+# All tests adapted from MPC 1.0.1's tests/sqr.dat
+describe MPC, '#sqr' do
+  it 'should calculate the square of a pure real argument' do
+    data = [
+      [["0x12345676543230p+52", 53, 16], [ 0,  2], MPC.new([GMP::F( "0x1111111000000f", 53, 16), GMP::F(0, 17)]), MPC::MPC_RNDNN],
+      [["0x1234567654322fp+52", 53, 16], [ 0,  3], MPC.new([GMP::F("-0x1111111000000f", 54, 16), GMP::F(0, 16)]), MPC::MPC_RNDZN],
+      [["0x12345676543230p+52", 53, 16], [ 0,  4], MPC.new([GMP::F( "0x1111111000000f", 55, 16), GMP::F(0, 15)]), MPC::MPC_RNDUN],
+      [["0x1234567654322fp+52", 53, 16], [ 0,  5], MPC.new([GMP::F("-0x1111111000000f", 56, 16), GMP::F(0, 14)]), MPC::MPC_RNDDN],
+      [["0x1234567654322fp+52", 53, 16], [ 0,  6], MPC.new([GMP::F( "0x1111111000000f", 57, 16), GMP::F(0, 13)]), MPC::MPC_RNDZZ],
+      [["0x12345676543230p+52", 53, 16], [ 0,  7], MPC.new([GMP::F("-0x1111111000000f", 58, 16), GMP::F(0, 12)]), MPC::MPC_RNDUZ],
+      [["0x1234567654322fp+52", 53, 16], [ 0,  8], MPC.new([GMP::F( "0x1111111000000f", 59, 16), GMP::F(0, 11)]), MPC::MPC_RNDDZ],
+      [["0x12345676543230p+52", 53, 16], [ 0,  9], MPC.new([GMP::F("-0x1111111000000f", 60, 16), GMP::F(0, 10)]), MPC::MPC_RNDNZ],
+      [["0x12345676543230p+52", 53, 16], [ 0, 10], MPC.new([GMP::F( "0x1111111000000f", 61, 16), GMP::F(0,  9)]), MPC::MPC_RNDUU],
+      [["0x1234567654322fp+52", 53, 16], [ 0, 11], MPC.new([GMP::F("-0x1111111000000f", 62, 16), GMP::F(0,  8)]), MPC::MPC_RNDDU],
+      [["0x12345676543230p+52", 53, 16], [ 0, 12], MPC.new([GMP::F( "0x1111111000000f", 63, 16), GMP::F(0,  7)]), MPC::MPC_RNDNU],
+      [["0x1234567654322fp+52", 53, 16], [ 0, 13], MPC.new([GMP::F("-0x1111111000000f", 64, 16), GMP::F(0,  6)]), MPC::MPC_RNDZU],
+      [["0x1234567654322fp+52", 53, 16], [ 0, 14], MPC.new([GMP::F( "0x1111111000000f", 65, 16), GMP::F(0,  5)]), MPC::MPC_RNDDD],
+      [["0x12345676543230p+52", 53, 16], [ 0, 15], MPC.new([GMP::F("-0x1111111000000f", 66, 16), GMP::F(0,  4)]), MPC::MPC_RNDND],
+      [["0x1234567654322fp+52", 53, 16], [ 0, 16], MPC.new([GMP::F( "0x1111111000000f", 67, 16), GMP::F(0,  3)]), MPC::MPC_RNDZD],
+      [["0x12345676543230p+52", 53, 16], [ 0, 17], MPC.new([GMP::F("-0x1111111000000f", 68, 16), GMP::F(0,  2)]), MPC::MPC_RNDUD],
+    ]
+    data.each do |expected_real, expected_imag, input, rounding_mode|
+      actual = input.sqr(rounding_mode)
+      actual.real.should eq GMP::F.new(*expected_real)
+      actual.imag.should eq GMP::F.new(*expected_imag)
+    end
+  end
+
+  it 'should calculate the square of a pure imaginary argument' do
+    data = [
+      [["-0xE1000002000000p+56", 53, 16], [ 0, 53], MPC.new([GMP::F(0, 53), GMP::F( "0xf0000001111111", 53, 16)]), MPC::MPC_RNDNN],
+      [["-0xe1000001fffff8p+56", 53, 16], [ 0, 52], MPC.new([GMP::F(0, 51), GMP::F( "0xf0000001111111", 54, 16)]), MPC::MPC_RNDZN],
+      [["-0xe1000001fffff8p+56", 53, 16], [ 0, 51], MPC.new([GMP::F(0, 49), GMP::F("-0xf0000001111111", 55, 16)]), MPC::MPC_RNDUN],
+      [["-0xE1000002000000p+56", 53, 16], [ 0, 50], MPC.new([GMP::F(0, 47), GMP::F("-0xf0000001111111", 56, 16)]), MPC::MPC_RNDDN],
+      [["-0xe1000001fffff8p+56", 53, 16], [ 0, 49], MPC.new([GMP::F(0, 45), GMP::F( "0xf0000001111111", 57, 16)]), MPC::MPC_RNDZZ],
+      [["-0xe1000001fffff8p+56", 53, 16], [ 0, 48], MPC.new([GMP::F(0, 43), GMP::F( "0xf0000001111111", 58, 16)]), MPC::MPC_RNDUZ],
+      [["-0xE1000002000000p+56", 53, 16], [ 0, 47], MPC.new([GMP::F(0, 41), GMP::F("-0xf0000001111111", 59, 16)]), MPC::MPC_RNDDZ],
+      [["-0xE1000002000000p+56", 53, 16], [ 0, 46], MPC.new([GMP::F(0, 39), GMP::F("-0xf0000001111111", 60, 16)]), MPC::MPC_RNDNZ],
+      [["-0xe1000001fffff8p+56", 53, 16], [ 0, 45], MPC.new([GMP::F(0, 37), GMP::F( "0xf0000001111111", 61, 16)]), MPC::MPC_RNDUU],
+      [["-0xE1000002000000p+56", 53, 16], [ 0, 44], MPC.new([GMP::F(0, 35), GMP::F( "0xf0000001111111", 62, 16)]), MPC::MPC_RNDDU],
+      [["-0xE1000002000000p+56", 53, 16], [ 0, 43], MPC.new([GMP::F(0, 33), GMP::F("-0xf0000001111111", 63, 16)]), MPC::MPC_RNDNU],
+      [["-0xe1000001fffff8p+56", 53, 16], [ 0, 42], MPC.new([GMP::F(0, 31), GMP::F("-0xf0000001111111", 64, 16)]), MPC::MPC_RNDZU],
+      [["-0xE1000002000000p+56", 53, 16], [ 0, 41], MPC.new([GMP::F(0, 29), GMP::F("-0xf0000001111111", 65, 16)]), MPC::MPC_RNDDD],
+      [["-0xE1000002000000p+56", 53, 16], [ 0, 40], MPC.new([GMP::F(0, 27), GMP::F("-0xf0000001111111", 66, 16)]), MPC::MPC_RNDND],
+      [["-0xe1000001fffff8p+56", 53, 16], [ 0, 39], MPC.new([GMP::F(0, 25), GMP::F("-0xf0000001111111", 67, 16)]), MPC::MPC_RNDZD],
+      [["-0xe1000001fffff8p+56", 53, 16], [ 0, 38], MPC.new([GMP::F(0, 23), GMP::F("-0xf0000001111111", 68, 16)]), MPC::MPC_RNDUD],
+    ]
+    data.each do |expected_real, expected_imag, input, rounding_mode|
+      actual = input.sqr(rounding_mode)
+      actual.real.should eq GMP::F.new(*expected_real)
+      actual.imag.should eq GMP::F.new(*expected_imag)
+    end
+  end
+
+  it 'should calculate the square of a double precision number' do
+    data = [
+      [[ "0x10000000020000p+04",   53, 16], [ "0x10000000effff",        53, 16], [ "0x400008000180fp-22",   53, 16], [ "0x7ffff0077efcbp-32",   53, 16], MPC::MPC_RNDNN],
+      [[ "0x3ffffffffffffd",       53, 16], [ "0x7ffffffffffff4p+52",   53, 16], [ "0x1fffffffffffff",      53, 16], [ "0x1ffffffffffffe",      53, 16], MPC::MPC_RNDZN],
+      [[ "0x1c16e5d4c4d5e7p-45",   53, 16], ["-0x7ffffff800007p-47",    53, 16], [ "0xf",                   53, 16], ["-0x1111111000000fp-53",  53, 16], MPC::MPC_RNDUN],
+      [[ "0xfdbac097c8dc50p+2096", 53, 16], [ "0x7f6e5d4c3b2a2p+1036",  53, 16], [ "0xfedcba9876543p+1024", 53, 16], [ "0x10000000000001p-42",  53, 16], MPC::MPC_RNDDN],
+      [["-0x10000000020000p+04",   53, 16], [ "0x10000000efffefp-04",   53, 16], [ "0x7ffff0077efcbp-32",   53, 16], [ "0x400008000180fp-22",   53, 16], MPC::MPC_RNDZZ],
+      [[ "0x3ffffffffffffe",       53, 16], ["-0x7ffffffffffff4p+52",   53, 16], [ "0x1fffffffffffff",      53, 16], ["-0x1ffffffffffffe",      53, 16], MPC::MPC_RNDUZ],
+      [[ "0xe0b72ea626af3p-44",    53, 16], [ "0x7ffffff800007p-47",    53, 16], [ "0xf",                   53, 16], [ "0x1111111000000fp-53",  53, 16], MPC::MPC_RNDDZ],
+      [["-0xfdbac097c8dc58p+2096", 53, 16], [ "0x7f6e5d4c3b2a1cp+1032", 53, 16], ["-0x10000000000001p-42",  53, 16], ["-0xfedcba9876543p+1024", 53, 16], MPC::MPC_RNDNZ],
+      [[ "0x10000000020001p+04",   53, 16], ["-0x10000000efffefp-04",   53, 16], [ "0x400008000180fp-22",   53, 16], ["-0x7ffff0077efcbp-32",   53, 16], MPC::MPC_RNDUU],
+      [["-0x3ffffffffffffe",       53, 16], ["-0x7ffffffffffff4p+52",   53, 16], ["-0x1ffffffffffffe",      53, 16], [ "0x1fffffffffffff",      53, 16], MPC::MPC_RNDDU],
+      [["-0x1C16E5D4C4D5E7p-45",   53, 16], [ "0x1ffffffe00001dp-49",   53, 16], ["-0x1111111000000fp-53",  53, 16], ["-0xf",                   53, 16], MPC::MPC_RNDNU],
+      [["-0xfdbac097c8dc50p+2096", 53, 16], ["-0x7f6e5d4c3b2a1cp+1032", 53, 16], [ "0x10000000000001p-42",  53, 16], ["-0xfedcba9876543p+1024", 53, 16], MPC::MPC_RNDZU],
+      [["-0x10000000020001p+04",   53, 16], ["-0x10000000effff",        53, 16], ["-0x7ffff0077efcbp-32",   53, 16], [ "0x400008000180fp-22",   53, 16], MPC::MPC_RNDDD],
+      [[ "0x3ffffffffffffd",       53, 16], ["-0x7ffffffffffff8p+52",   53, 16], ["-0x1fffffffffffff",      53, 16], [ "0x1ffffffffffffe",      53, 16], MPC::MPC_RNDND],
+      [["-0xE0B72EA626AF3p-44",    53, 16], ["-0x1FFFFFFE00001Dp-49",   53, 16], [ "0x1111111000000fp-53",  53, 16], ["-0xf",                   53, 16], MPC::MPC_RNDZD],
+      [[ "0xfdbac097c8dc58p+2096", 53, 16], ["-0x7f6e5d4c3b2a2p+1036",  53, 16], ["-0xfedcba9876543p+1024", 53, 16], ["0x10000000000001p-42",   53, 16], MPC::MPC_RNDUD]
+    ]
+    data.each do |expected_real, expected_imag, input_real, input_imag, rounding_mode|
+      actual = MPC.new([GMP::F.new(*input_real), GMP::F(*input_imag)]).sqr(rounding_mode)
+      actual.real.should eq GMP::F.new(*expected_real)
+      actual.imag.should eq GMP::F.new(*expected_imag)
+    end
+  end
+
+  it 'should calculate the square: "improve test coverage"' do
+    actual = MPC.new([GMP::F(17592186044416, 30), GMP::F(536870913, 30)], 30).sqr
+    actual.real.should eq GMP::F.new(309485009533114692573069312, 30)
+    actual.imag.should eq GMP::F.new(18889465966662952943616, 30)
+  end
+
+  it 'should calculate the square: "improve test coverage"; a case where x+y or x-y are 0' do
+    actual = MPC.new([GMP::F(1, 4), GMP::F(1, 4)], 4).sqr
+    actual.real.should eq GMP::F.new(0, 4)
+    actual.imag.should eq GMP::F.new(2, 4)
+  end
+
+  ### There are still more tests in sqr.data, each a result of some bug or other thing
+end