gnu_mpc 0.8.0

data/manual_template.latex

@@ -0,0 +1,177 @@
+\documentclass[$if(fontsize)$$fontsize$,$endif$$if(lang)$$lang$,$endif$$if(papersize)$$papersize$,$endif$]{$documentclass$}
+\usepackage[T1]{fontenc}
+\usepackage{lmodern}
+\usepackage{amssymb,amsmath}
+\usepackage{booktabs}
+\usepackage{ifxetex,ifluatex}
+\usepackage{fixltx2e} % provides \textsubscript
+% use microtype if available
+\IfFileExists{microtype.sty}{\usepackage{microtype}}{}
+\ifnum 0\ifxetex 1\fi\ifluatex 1\fi=0 % if pdftex
+  \usepackage[utf8]{inputenc}
+$if(euro)$
+  \usepackage{eurosym}
+$endif$
+\else % if luatex or xelatex
+  \usepackage{fontspec}
+  \ifxetex
+    \usepackage{xltxtra,xunicode}
+  \fi
+  \defaultfontfeatures{Mapping=tex-text,Scale=MatchLowercase}
+  \newcommand{\euro}{€}
+$if(mainfont)$
+    \setmainfont{$mainfont$}
+$endif$
+$if(sansfont)$
+    \setsansfont{$sansfont$}
+$endif$
+$if(monofont)$
+    \setmonofont{$monofont$}
+$endif$
+$if(mathfont)$
+    \setmathfont{$mathfont$}
+$endif$
+\fi
+$if(geometry)$
+\usepackage[$for(geometry)$$geometry$$sep$,$endfor$]{geometry}
+$else$
+\usepackage[top=1in, bottom=0.8in, left=1.2in, right=1.2in]{geometry}
+$endif$
+$if(natbib)$
+\usepackage{natbib}
+\bibliographystyle{plainnat}
+$endif$
+$if(biblatex)$
+\usepackage{biblatex}
+$if(biblio-files)$
+\bibliography{$biblio-files$}
+$endif$
+$endif$
+$if(listings)$
+\usepackage{listings}
+$endif$
+$if(lhs)$
+\lstnewenvironment{code}{\lstset{language=Haskell,basicstyle=\small\ttfamily}}{}
+$endif$
+$if(highlighting-macros)$
+$highlighting-macros$
+$endif$
+$if(verbatim-in-note)$
+\usepackage{fancyvrb}
+$endif$
+$if(tables)$
+\usepackage{longtable}
+$endif$
+$if(graphics)$
+\usepackage{graphicx}
+% We will generate all images so they have a width \maxwidth. This means
+% that they will get their normal width if they fit onto the page, but
+% are scaled down if they would overflow the margins.
+\makeatletter
+\def\maxwidth{\ifdim\Gin@nat@width>\linewidth\linewidth
+\else\Gin@nat@width\fi}
+\makeatother
+\let\Oldincludegraphics\includegraphics
+\renewcommand{\includegraphics}[1]{\Oldincludegraphics[width=\maxwidth]{#1}}
+$endif$
+\ifxetex
+  \usepackage[setpagesize=false, % page size defined by xetex
+              unicode=false, % unicode breaks when used with xetex
+              xetex]{hyperref}
+\else
+  \usepackage[unicode=true]{hyperref}
+\fi
+\hypersetup{breaklinks=true,
+            bookmarks=true,
+            pdfauthor={$author-meta$},
+            pdftitle={$title-meta$},
+            colorlinks=true,
+            urlcolor=$if(urlcolor)$$urlcolor$$else$blue$endif$,
+            linkcolor=$if(linkcolor)$$linkcolor$$else$magenta$endif$,
+            pdfborder={0 0 0}}
+$if(links-as-notes)$
+% Make links footnotes instead of hotlinks:
+\renewcommand{\href}[2]{#2\footnote{\url{#1}}}
+$endif$
+$if(strikeout)$
+\usepackage[normalem]{ulem}
+% avoid problems with \sout in headers with hyperref:
+\pdfstringdefDisableCommands{\renewcommand{\sout}{}}
+$endif$
+\setlength{\parindent}{0pt}
+\setlength{\parskip}{6pt plus 2pt minus 1pt}
+\setlength{\emergencystretch}{3em}  % prevent overfull lines
+$if(numbersections)$
+$else$
+\setcounter{secnumdepth}{0}
+$endif$
+$if(verbatim-in-note)$
+\VerbatimFootnotes % allows verbatim text in footnotes
+$endif$
+$if(lang)$
+\ifxetex
+  \usepackage{polyglossia}
+  \setmainlanguage{$mainlang$}
+\else
+  \usepackage[$lang$]{babel}
+\fi
+$endif$
+$for(header-includes)$
+$header-includes$
+$endfor$
+
+\newlength{\methwidth}
+\setlength{\methwidth}{0.8in}
+\newlength{\defnwidth}
+\setlength{\defnwidth}{\textwidth}
+\addtolength{\defnwidth}{-2.0in}
+\newlength{\returnwidth}
+\setlength{\returnwidth}{0.8in}
+\def\qquad{\quad\quad}
+\def\qqqquad{\quad\quad\quad\quad}
+
+$if(title)$
+\title{$title$}
+$endif$
+\author{$for(author)$$author$$sep$ \and $endfor$}
+\date{$date$}
+
+\begin{document}
+$if(title)$
+\maketitle
+$endif$
+
+$for(include-before)$
+$include-before$
+
+$endfor$
+$if(toc)$
+{
+\hypersetup{linkcolor=black}
+\tableofcontents
+}
+$endif$
+$body$
+
+$if(natbib)$
+$if(biblio-files)$
+$if(biblio-title)$
+$if(book-class)$
+\renewcommand\bibname{$biblio-title$}
+$else$
+\renewcommand\refname{$biblio-title$}
+$endif$
+$endif$
+\bibliography{$biblio-files$}
+
+$endif$
+$endif$
+$if(biblatex)$
+\printbibliography$if(biblio-title)$[title=$biblio-title$]$endif$
+
+$endif$
+$for(include-after)$
+$include-after$
+
+$endfor$
+\end{document}

data/spec/acos_spec.rb

@@ -0,0 +1,66 @@
+require File.join(File.dirname(__FILE__), 'spec_helper')
+
+# All tests adapted from MPC 1.0.1's tests/acos.dat
+describe MPC, '#acos' do
+  it 'should calculate the inverse cosine of a pure real argument' do
+    data = [
+      [[ "0x1921FB54442D18p-51", 53, 16], [ "0x13D2B7539DBA4Cp-51", 53, 16], MPC.new([GMP::F(-6  ), -GMP::F(0)])],
+      [[ "0x1921FB54442D18p-51", 53, 16], ["-0x13D2B7539DBA4Cp-51", 53, 16], MPC.new([GMP::F(-6  ),  GMP::F(0)])],
+      [[ "0x1921FB54442D18p-51", 53, 16], [ "0x15124271980435p-52", 53, 16], MPC.new([GMP::F(-2  ), -GMP::F(0)])],
+      [[ "0x1921FB54442D18p-51", 53, 16], ["-0x15124271980435p-52", 53, 16], MPC.new([GMP::F(-2  ),  GMP::F(0)])],
+      [[ "0x1921FB54442D18p-51", 53, 16], [0],                               MPC.new([GMP::F(-1  ), -GMP::F(0)])],
+      [[ "0x1921FB54442D18p-51", 53, 16], [0],                               MPC.new([GMP::F(-1  ),  GMP::F(0)])],
+      [[ "0x10C152382D7366p-51", 53, 16], [0],                               MPC.new([GMP::F(-0.5), -GMP::F(0)])],
+      [[ "0x10C152382D7366p-51", 53, 16], [0],                               MPC.new([GMP::F(-0.5),  GMP::F(0)])],
+      [[ "0x10C152382D7366p-52", 53, 16], [0],                               MPC.new([GMP::F( 0.5),  GMP::F(0)])],
+      [[ "0x10C152382D7366p-52", 53, 16], [0],                               MPC.new([GMP::F( 0.5), -GMP::F(0)])],
+      [[0],                               [0],                               MPC.new([GMP::F( 1  ),  GMP::F(0)])],
+      [[0],                               [0],                               MPC.new([GMP::F( 1  ), -GMP::F(0)])],
+      [[0],                               ["-0x15124271980435p-52", 53, 16], MPC.new([GMP::F( 2  ),  GMP::F(0)])],
+      [[0],                               [ "0x15124271980435p-52", 53, 16], MPC.new([GMP::F( 2  ), -GMP::F(0)])],
+      [[0],                               ["-0x13D2B7539DBA4Cp-51", 53, 16], MPC.new([GMP::F( 6  ),  GMP::F(0)])],
+      [[0],                               [ "0x13D2B7539DBA4Cp-51", 53, 16], MPC.new([GMP::F( 6  ), -GMP::F(0)])]
+    ]
+    data.each do |expected_real, expected_imag, input|
+      actual = input.acos
+      actual.real.should eq GMP::F.new(*expected_real)
+      actual.imag.should eq GMP::F.new(*expected_imag)
+    end
+  end
+
+  it 'should calculate the inverse cosine of a pure imaginary argument' do
+    data = [
+      [["0x1921FB54442D18p-52", 53, 16], [ "0x1D185B507EDC0Ep-52", 53, 16], MPC.new([-GMP::F(0), GMP::F(-3   )])],
+      [["0x1921FB54442D18p-52", 53, 16], [ "0x1D185B507EDC0Ep-52", 53, 16], MPC.new([ GMP::F(0), GMP::F(-3   )])],
+      [["0x1921FB54442D18p-52", 53, 16], [ "0x1FACFB2399E637p-55", 53, 16], MPC.new([-GMP::F(0), GMP::F(-0.25)])],
+      [["0x1921FB54442D18p-52", 53, 16], [ "0x1FACFB2399E637p-55", 53, 16], MPC.new([ GMP::F(0), GMP::F(-0.25)])],
+      [["0x1921FB54442D18p-52", 53, 16], ["-0x1FACFB2399E637p-55", 53, 16], MPC.new([-GMP::F(0), GMP::F( 0.25)])],
+      [["0x1921FB54442D18p-52", 53, 16], ["-0x1FACFB2399E637p-55", 53, 16], MPC.new([ GMP::F(0), GMP::F( 0.25)])],
+      [["0x1921FB54442D18p-52", 53, 16], ["-0x1D185B507EDC0Ep-52", 53, 16], MPC.new([-GMP::F(0), GMP::F( 3   )])],
+      [["0x1921FB54442D18p-52", 53, 16], ["-0x1D185B507EDC0Ep-52", 53, 16], MPC.new([ GMP::F(0), GMP::F( 3   )])],
+    ]
+    data.each do |expected_real, expected_imag, input|
+      actual = input.acos
+      actual.real.should eq GMP::F.new(*expected_real)
+      actual.imag.should eq GMP::F.new(*expected_imag)
+    end
+  end
+
+  it 'should calculate the inverse cosine with other precisions' do
+    actual = MPC.new([GMP::F.new(2), GMP::F.new(1)]).acos(MPC::MPC_RNDNZ, 2)
+    actual.real.should eq GMP::F.new( 0.5, 2)
+    actual.imag.should eq GMP::F.new(-1,   2)
+  end
+
+  it 'should calculate the inverse cosine with other precisions' do
+    actual = MPC.new([GMP::F.new(8.5), GMP::F.new(-71)]).acos(MPC::MPC_RNDNU, 9)
+    actual.real.should eq GMP::F.new("0x5Dp-6", 9, 16)
+    actual.imag.should eq GMP::F.new("0x9Fp-5", 9, 16)
+  end
+
+  it 'should calculate the inverse cosine with other precisions' do
+    actual = MPC.new([GMP::F.new("0x3243F6A8885A3p-48", 53, 16), GMP::F.new("0x162E42FEFA39EFp-53", 53, 16)]).acos
+    actual.real.should eq GMP::F.new("0x74C141310E695p-53", 53, 16)
+    actual.imag.should eq GMP::F.new("-0x1D6D2CFA9F3F11p-52", 53, 16)
+  end
+end

data/spec/add_args_spec.rb

@@ -0,0 +1,51 @@
+require File.join(File.dirname(__FILE__), 'spec_helper')
+
+describe MPC, '#add' do
+  it 'should calculate the sum of an MPC and a Fixnum' do
+    op1 = MPC.new([7, 16])
+    actual = op1.add(0)
+    actual.real.should eq 7
+    actual.imag.should eq 16
+
+    actual = op1.add(6)
+    actual.real.should eq 13
+    actual.imag.should eq 16
+
+    actual = op1.add(-1000)
+    actual.real.should eq -993
+    actual.imag.should eq 16
+
+    op1 = MPC.new([GMP::F(3.1), GMP::F(16)])
+    actual = op1.add(3)
+    actual.real.should eq GMP::F(6.1)
+    actual.imag.should eq 16
+  end
+
+  it 'should calculate the sum of an MPC and a GMP::Z' do
+    op1 = MPC.new([7, 16])
+    actual = op1.add(GMP::Z(0))
+    actual.real.should eq 7
+    actual.imag.should eq 16
+
+    actual = op1.add(GMP::Z(-1000))
+    actual.real.should eq -993
+    actual.imag.should eq 16
+
+    op1 = MPC.new([GMP::F(3.1), GMP::F(16)])
+    actual = op1.add(GMP::Z(3))
+    actual.real.should eq GMP::F(6.1)
+    actual.imag.should eq 16
+  end
+
+  it 'should calculate the sum of an MPC and a Bignum' do
+    op1 = MPC.new([7, 16])
+    actual = op1.add(2**65)  # 36893488147419103232
+    actual.real.should eq 36893488147419103239
+    actual.imag.should eq 16
+
+    op1 = MPC.new([GMP::F(3.1), GMP::F(16)])
+    actual = op1.add(2**65)  # 36893488147419103232
+    actual.real.should eq GMP::F("36893488147419103235.1")
+    actual.imag.should eq 16
+  end
+end

data/spec/add_fr_spec.rb

@@ -0,0 +1,54 @@
+require File.join(File.dirname(__FILE__), 'spec_helper')
+
+# All tests adapted from MPC 1.0.1's tests/add.dat
+describe MPC, '#add' do
+  it 'should calculate the sum of two real MPCs' do
+    op1 = MPC.new([1, 0])
+    op2_106 = GMP::F("0x10000000000001p-106", 53, 16)
+    actual = op1.add(op2_106, MPC::MPC_RNDNN)
+    actual.real.should eq GMP::F("0x10000000000000p-52", 53, 16)
+    actual.imag.should eq GMP::F(0)
+
+    op2_105 = GMP::F("0x10000000000001p-105", 53, 16)
+    actual = op1.add(op2_105, MPC::MPC_RNDNN)
+    actual.real.should eq GMP::F("0x10000000000001p-52", 53, 16)
+    actual.imag.should eq GMP::F(0)
+
+    op2_104 = GMP::F("0x10000000000001p-104", 53, 16)
+    actual = op1.add(op2_104, MPC::MPC_RNDNN)
+    actual.real.should eq GMP::F("0x10000000000001p-52", 53, 16)
+    actual.imag.should eq GMP::F(0)
+
+    actual = op1.add(op2_105, MPC::MPC_RNDZZ)
+    actual.real.should eq GMP::F("0x10000000000000p-52", 53, 16)
+    actual.imag.should eq GMP::F(0)
+
+    actual = op1.add(op2_105, MPC::MPC_RNDUU)
+    actual.real.should eq GMP::F("0x10000000000001p-52", 53, 16)
+    actual.imag.should eq GMP::F(0)
+
+    actual = op1.add(op2_105, MPC::MPC_RNDDD)
+    actual.real.should eq GMP::F("0x10000000000000p-52", 53, 16)
+    actual.imag.should eq GMP::F(0)
+  end
+
+  it 'should calculate the sum of an imaginary argument with a complex one.' do
+    op1 = MPC.new([0, 1])
+    op2 = GMP::F(1)
+    actual = op1.add(op2, MPC::MPC_RNDNN)
+    actual.real.should eq GMP::F(1)
+    actual.imag.should eq GMP::F(1)
+
+    actual = op1.add(op2, MPC::MPC_RNDZZ)
+    actual.real.should eq GMP::F(1)
+    actual.imag.should eq GMP::F(1)
+
+    actual = op1.add(op2, MPC::MPC_RNDUU)
+    actual.real.should eq GMP::F(1)
+    actual.imag.should eq GMP::F(1)
+
+    actual = op1.add(op2, MPC::MPC_RNDDD)
+    actual.real.should eq GMP::F(1)
+    actual.imag.should eq GMP::F(1)
+  end
+end

data/spec/add_spec.rb

@@ -0,0 +1,32 @@
+require File.join(File.dirname(__FILE__), 'spec_helper')
+
+# All tests adapted from MPC 1.0.1's tests/add.dat
+describe MPC, '#add' do
+  it 'should calculate the sum of two real MPCs' do
+    op1 = MPC.new([1, 0])
+    op2 = MPC.new([GMP::F("0x10000000000001p-105", 53, 16), GMP::F(0)])
+    actual = op1.add(op2, MPC::MPC_RNDNN)
+    actual.real.should eq GMP::F("0x10000000000001p-52", 53, 16)
+    actual.imag.should eq GMP::F(0)
+
+    actual = op1.add(op2, MPC::MPC_RNDZZ)
+    actual.real.should eq GMP::F("0x10000000000000p-52", 53, 16)
+    actual.imag.should eq GMP::F(0)
+
+    actual = op1.add(op2, MPC::MPC_RNDUU)
+    actual.real.should eq GMP::F("0x10000000000001p-52", 53, 16)
+    actual.imag.should eq GMP::F(0)
+
+    actual = op1.add(op2, MPC::MPC_RNDDD)
+    actual.real.should eq GMP::F("0x10000000000000p-52", 53, 16)
+    actual.imag.should eq GMP::F(0)
+  end
+
+  it 'should calculate the sum of two real MPCs' do
+    op2 = MPC.new([GMP::F(0), GMP::F("0x10000000000001p-105", 53, 16)])
+    op1 = MPC.new([0, 1])
+    actual = op1.add(op2, MPC::MPC_RNDNN)
+    actual.real.should eq GMP::F(0)
+    actual.imag.should eq GMP::F("0x10000000000001p-52", 53, 16)
+  end
+end

data/spec/asin_spec.rb

@@ -0,0 +1,66 @@
+require File.join(File.dirname(__FILE__), 'spec_helper')
+
+# All tests adapted from MPC 1.0.1's tests/asin.dat
+describe MPC, '#sin' do
+  it 'should calculate the inverse sine of a pure real argument' do
+    data = [
+      [["-0x1921FB54442D18p-52", 53, 16], ["-0x1ECC2CAEC5160Ap-53", 53, 16], MPC.new([GMP::F(-1.5), -GMP::F(0)])],
+      [["-0x1921FB54442D18p-52", 53, 16], [ "0x1ECC2CAEC5160Ap-53", 53, 16], MPC.new([GMP::F(-1.5),  GMP::F(0)])],
+      [["-0x1921FB54442D18p-52", 53, 16], [0],                               MPC.new([GMP::F(-1  ), -GMP::F(0)])],
+      [["-0x1921FB54442D18p-52", 53, 16], [0],                               MPC.new([GMP::F(-1  ),  GMP::F(0)])],
+      [["-0x10C152382D7366p-53", 53, 16], [0],                               MPC.new([GMP::F(-0.5), -GMP::F(0)])],
+      [["-0x10C152382D7366p-53", 53, 16], [0],                               MPC.new([GMP::F(-0.5),  GMP::F(0)])],
+      [[ "0x10C152382D7366p-53", 53, 16], [0],                               MPC.new([GMP::F( 0.5), -GMP::F(0)])],
+      [[ "0x10C152382D7366p-53", 53, 16], [0],                               MPC.new([GMP::F( 0.5),  GMP::F(0)])],
+      [[ "0x1921FB54442D18p-52", 53, 16], [0],                               MPC.new([GMP::F( 1  ), -GMP::F(0)])],
+      [[ "0x1921FB54442D18p-52", 53, 16], [0],                               MPC.new([GMP::F( 1  ),  GMP::F(0)])],
+      [[ "0x1921FB54442D18p-52", 53, 16], ["-0x1ECC2CAEC5160Ap-53", 53, 16], MPC.new([GMP::F( 1.5), -GMP::F(0)])],
+      [[ "0x1921FB54442D18p-52", 53, 16], [ "0x1ECC2CAEC5160Ap-53", 53, 16], MPC.new([GMP::F( 1.5),  GMP::F(0)])]
+    ]
+    data.each do |expected_real, expected_imag, input|
+      actual = input.asin
+      actual.real.should eq GMP::F.new(*expected_real)
+      actual.imag.should eq GMP::F.new(*expected_imag)
+    end
+  end
+
+  it 'should calculate the inverse sine of a pure imaginary argument' do
+    data = [
+      [[-GMP::F(0)], ["-0x131DC0090B63D8p-52", 53, 16], MPC.new([-GMP::F(0), GMP::F(-1.5)])],
+      [[ GMP::F(0)], ["-0x131DC0090B63D8p-52", 53, 16], MPC.new([ GMP::F(0), GMP::F(-1.5)])],
+      [[-GMP::F(0)], ["-0x1C34366179D427p-53", 53, 16], MPC.new([-GMP::F(0), GMP::F(-1  )])],
+      [[ GMP::F(0)], ["-0x1C34366179D427p-53", 53, 16], MPC.new([ GMP::F(0), GMP::F(-1  )])],
+      [[-GMP::F(0)], ["-0x1ECC2CAEC5160Ap-54", 53, 16], MPC.new([-GMP::F(0), GMP::F(-0.5)])],
+      [[ GMP::F(0)], ["-0x1ECC2CAEC5160Ap-54", 53, 16], MPC.new([ GMP::F(0), GMP::F(-0.5)])],
+      [[-GMP::F(0)], [ "0x1ECC2CAEC5160Ap-54", 53, 16], MPC.new([-GMP::F(0), GMP::F( 0.5)])],
+      [[ GMP::F(0)], [ "0x1ECC2CAEC5160Ap-54", 53, 16], MPC.new([ GMP::F(0), GMP::F( 0.5)])],
+      [[-GMP::F(0)], [ "0x1C34366179D427p-53", 53, 16], MPC.new([-GMP::F(0), GMP::F( 1  )])],
+      [[ GMP::F(0)], [ "0x1C34366179D427p-53", 53, 16], MPC.new([ GMP::F(0), GMP::F( 1  )])],
+      [[-GMP::F(0)], [ "0x131DC0090B63D8p-52", 53, 16], MPC.new([-GMP::F(0), GMP::F( 1.5)])],
+      [[ GMP::F(0)], [ "0x131DC0090B63D8p-52", 53, 16], MPC.new([ GMP::F(0), GMP::F( 1.5)])]
+    ]
+    data.each do |expected_real, expected_imag, input|
+      actual = input.asin
+      actual.real.should eq GMP::F.new(*expected_real)
+      actual.imag.should eq GMP::F.new(*expected_imag)
+    end
+  end
+
+  it 'should calculate the inverse sine of ieee-754 double precision numbers' do
+    actual = MPC.new([GMP::F.new(17), GMP::F.new(42)]).asin(MPC::MPC_RNDNN)
+    actual.real.should eq GMP::F.new("0x189BF9EC7FCD5Bp-54", 53, 16)
+    actual.imag.should eq GMP::F.new("0x1206ECFA94614Bp-50", 53, 16)
+  end
+
+  it 'should calculate the inverse sine with other precisions' do
+    actual = MPC.new([GMP::F.new(96), GMP::F.new("0x1p-8", 2, 16)]).asin(MPC::MPC_RNDNN, 2)
+    actual.real.should eq GMP::F.new(1.5, 2)
+    actual.imag.should eq GMP::F.new(6,   2)
+  end
+
+  it 'should calculate the inverse sine with other precisions' do
+    actual = MPC.new([GMP::F.new(96), GMP::F.new("0x1p-8", 2, 16)]).asin(MPC::MPC_RNDNN, 8)
+    actual.real.should eq GMP::F.new("0xC9p-7", 8, 16)
+    actual.imag.should eq GMP::F.new("0x15p-2", 8, 16)
+  end
+end