gsl 1.15.3 → 1.16.0.6

data/rdoc/linalg_complex.rdoc

@@ -83,6 +83,6 @@
 # * GSL::Linalg::Complex::LU::LUMatrix#sgndet(signum)
 #
 #
-# {back}[link:files/rdoc/linalg_rdoc.html]
+# {back}[link:rdoc/linalg_rdoc.html]
 #
 #

data/rdoc/math.rdoc

@@ -1,18 +1,18 @@
 #
 # = Mathematical Functions
 # Contents:
-# 1. {Mathematical Constants}[link:files/rdoc/math_rdoc.html#1]
-#    1. {Infinities and Not-a-number}[link:files/rdoc/math_rdoc.html#2]
-#    1. {Constants}[link:files/rdoc/math_rdoc.html#2.1]
-# 1. {Module functions}[link:files/rdoc/math_rdoc.html#2.2]
-# 1. {Elementary Functions}[link:files/rdoc/math_rdoc.html#3]
-# 1. {Small Integer Powers}[link:files/rdoc/math_rdoc.html#4]
-# 1. {Testing the Sign of Numbers}[link:files/rdoc/math_rdoc.html#5]
-# 1. {Testing for Odd and Even Numbers}[link:files/rdoc/math_rdoc.html#6]
-# 1. {Maximum and Minimum functions}[link:files/rdoc/math_rdoc.html#7]
-# 1. {Approximate Comparison of Floating Point Numbers}[link:files/rdoc/math_rdoc.html#8]
-#
-# == {}[link:index.html"name="1] Mathematical Constants
+# 1. {Mathematical Constants}[link:rdoc/math_rdoc.html#label-Mathematical+Constants]
+#    1. {Infinities and Not-a-number}[link:rdoc/math_rdoc.html#label-Infinities+and+Not-a-number]
+#    1. {Constants}[link:rdoc/math_rdoc.html#label-Constants]
+# 1. {Module functions}[link:rdoc/math_rdoc.html#label-Module+functions]
+# 1. {Elementary Functions}[link:rdoc/math_rdoc.html#label-Elementary+Functions]
+# 1. {Small Integer Powers}[link:rdoc/math_rdoc.html#label-Small+Integer+Powers]
+# 1. {Testing the Sign of Numbers}[link:rdoc/math_rdoc.html#label-Testing+the+Sign+of+Numbers]
+# 1. {Testing for Odd and Even Numbers}[link:rdoc/math_rdoc.html#label-Testing+for+Odd+and+Even+Numbers]
+# 1. {Maximum and Minimum functions}[link:rdoc/math_rdoc.html#label-Maximum+and+Minimum+functions]
+# 1. {Approximate Comparison of Floating Point Numbers}[link:rdoc/math_rdoc.html#label-Approximate+Comparison+of+Floating+Point+Numbers]
+#
+# == Mathematical Constants
 # ---
 # * GSL::M_E
 #
@@ -82,18 +82,18 @@
 #
 #   Euler's constant
 #
-# == {}[link:index.html"name="2] Infinities and Not-a-number
+# == Infinities and Not-a-number
 #
-# === {}[link:index.html"name="2.1] Constants
+# === Constants
 # ---
 # * GSL::POSINF
 #
-#   The IEEE representation of positive infinity, 
+#   The IEEE representation of positive infinity,
 #   computed from the expression +1.0/0.0.
 # ---
 # * GSL::NEGINF
 #
-#   The IEEE representation of negative infinity, 
+#   The IEEE representation of negative infinity,
 #   computed from the expression -1.0/0.0.
 # ---
 # * GSL::NAN
@@ -101,7 +101,7 @@
 #   The IEEE representation of the Not-a-Number symbol,
 #   computed from the ratio 0.0/0.0.
 #
-# === {}[link:index.html"name="2.2] Module functions
+# === Module functions
 # ---
 # * GSL::isnan(x)
 #
@@ -113,59 +113,59 @@
 # ---
 # * GSL::isinf(x)
 #
-#   This returns +1 if <tt>x</tt>  is positive infinity, 
+#   This returns +1 if <tt>x</tt>  is positive infinity,
 #   -1 if <tt>x</tt>  is negative infinity and 0 otherwise.
 #   NOTE: In Darwin9.5.0-gcc4.0.1, this method returns 1 for -inf.
 # ---
 # * GSL::isinf?(x)
 #
-#   This returns <tt>true</tt> if <tt>x</tt> is positive or negative infinity, 
+#   This returns <tt>true</tt> if <tt>x</tt> is positive or negative infinity,
 #   and <tt>false</tt> otherwise.
 # ---
 # * GSL::finite(x)
 #
-#   This returns 1 if <tt>x</tt> is a real number, 
+#   This returns 1 if <tt>x</tt> is a real number,
 #   and 0 if it is infinite or not-a-number.
 # ---
 # * GSL::finite?(x)
 #
-#   This returns <tt>true</tt> if <tt>x</tt> is a real number, 
+#   This returns <tt>true</tt> if <tt>x</tt> is a real number,
 #   and <tt>false</tt> if it is infinite or not-a-number.
 #
-# == {}[link:index.html"name="3] Elementary Functions
+# == Elementary Functions
 # ---
 # * GSL::log1p(x)
 #
-#   This method computes the value of log(1+x) 
-#   in a way that is accurate for small <tt>x</tt>. It provides an alternative 
+#   This method computes the value of log(1+x)
+#   in a way that is accurate for small <tt>x</tt>. It provides an alternative
 #   to the BSD math function log1p(x).
 # ---
 # * GSL::expm1(x)
 #
-#   This method computes the value of exp(x)-1 
-#   in a way that is accurate for small <tt>x</tt>. It provides an alternative 
+#   This method computes the value of exp(x)-1
+#   in a way that is accurate for small <tt>x</tt>. It provides an alternative
 #   to the BSD math function expm1(x).
 # ---
 # * GSL::hypot(x, y)
 #
-#   This method computes the value of sqrt{x^2 + y^2} in a way that 
+#   This method computes the value of sqrt{x^2 + y^2} in a way that
 #   avoids overflow.
 # ---
-# * GSL::hypot3(x, y, z) 
+# * GSL::hypot3(x, y, z)
 #
-#   Computes the value of sqrt{x^2 + y^2 + z^2} in a way that avoids overflow. 
+#   Computes the value of sqrt{x^2 + y^2 + z^2} in a way that avoids overflow.
 # ---
 # * GSL::acosh(x)
 #
-#   This method computes the value of arccosh(x). 
+#   This method computes the value of arccosh(x).
 # ---
 # * GSL::asinh(x)
 #
-#   This method computes the value of arcsinh(x). 
+#   This method computes the value of arcsinh(x).
 # ---
 # * GSL::atanh(x)
 #
-#   This method computes the value of arctanh(x). 
+#   This method computes the value of arctanh(x).
 #
 #   These methods above can take argument <tt>x</tt> of
 #   Integer, Float, Array, Vector or Matrix.
@@ -173,22 +173,22 @@
 # ---
 # * GSL::ldexp(x)
 #
-#   This method computes the value of x * 2^e. 
+#   This method computes the value of x * 2^e.
 # ---
 # * GSL::frexp(x)
 #
-#   This method splits the number <tt>x</tt> into its normalized fraction 
-#   f and exponent e, such that x = f * 2^e and 0.5 <= f < 1. 
-#   The method returns f and the exponent e as an array, [f, e]. 
-#   If <tt>x</tt> is zero, both f and e are set to zero. 
+#   This method splits the number <tt>x</tt> into its normalized fraction
+#   f and exponent e, such that x = f * 2^e and 0.5 <= f < 1.
+#   The method returns f and the exponent e as an array, [f, e].
+#   If <tt>x</tt> is zero, both f and e are set to zero.
 #
-# == {}[link:index.html"name="4] Small Integer Powers
+# == Small Integer Powers
 # ---
 # * GSL::pow_int(x, n)
 #
-#   This routine computes the power <tt>x^n</tt> for integer <tt>n</tt>. 
-#   The power is computed efficiently -- for example, x^8 is computed as 
-#   ((x^2)^2)^2, requiring only 3 multiplications. 
+#   This routine computes the power <tt>x^n</tt> for integer <tt>n</tt>.
+#   The power is computed efficiently -- for example, x^8 is computed as
+#   ((x^2)^2)^2, requiring only 3 multiplications.
 #
 # ---
 # * GSL::pow_2(x)
@@ -200,25 +200,25 @@
 # * GSL::pow_8(x)
 # * GSL::pow_9(x)
 #
-#   These methods can be used to compute small integer powers x^2, x^3, etc. 
+#   These methods can be used to compute small integer powers x^2, x^3, etc.
 #   efficiently.
 #
-# == {}[link:index.html"name="5] Testing the Sign of Numbers
+# == Testing the Sign of Numbers
 # ---
 # * GSL::SIGN(x)
 # * GSL::sign(x)
 #
-#   Return the sign of <tt>x</tt>. 
-#   It is defined as ((x) >= 0 ? 1 : -1). 
-#   Note that with this definition the sign of zero is positive 
+#   Return the sign of <tt>x</tt>.
+#   It is defined as ((x) >= 0 ? 1 : -1).
+#   Note that with this definition the sign of zero is positive
 #   (regardless of its IEEE sign bit).
 #
-# == {}[link:index.html"name="6] Testing for Odd and Even Numbers
+# == Testing for Odd and Even Numbers
 # ---
 # * GSL::is_odd(n)
 # * GSL::IS_ODD(n)
 #
-#   Evaluate to 1 if <tt>n</tt> is odd and 0 if <tt>n</tt> is even. 
+#   Evaluate to 1 if <tt>n</tt> is odd and 0 if <tt>n</tt> is even.
 #   The argument <tt>n</tt> must be of Fixnum type.
 # ---
 # * GSL::is_odd?(n)
@@ -229,7 +229,7 @@
 # * GSL::is_even(n)
 # * GSL::IS_EVEN(n)
 #
-#   Evaluate to 1 if <tt>n</tt> is even and 0 if <tt>n</tt> is odd. 
+#   Evaluate to 1 if <tt>n</tt> is even and 0 if <tt>n</tt> is odd.
 #   The argument <tt>n</tt> must be of Fixnum type.
 # ---
 # * GSL::is_even?(n)
@@ -237,25 +237,25 @@
 #
 #   Return <tt>true</tt> if <tt>n</tt> is even and <tt>false</tt> if odd.
 #
-# == {}[link:index.html"name="7] Maximum and Minimum functions
+# == Maximum and Minimum functions
 # ---
 # * GSL::max(a, b)
 # * GSL::MAX(a, b)
 # * GSL::min(a, b)
 # * GSL::MIN(a, b)
 #
-#   
-# == {}[link:index.html"name="8] Approximate Comparison of Floating Point Numbers
+#
+# == Approximate Comparison of Floating Point Numbers
 # ---
 # * GSL::fcmp(a, b, epsilon = 1e-10)
 #
-#   This method determines whether <tt>x</tt> and <tt>y</tt> are approximately equal to a 
+#   This method determines whether <tt>x</tt> and <tt>y</tt> are approximately equal to a
 #   relative accuracy <tt>epsilon</tt>.
 # ---
 # * GSL::equal?(a, b, epsilon = 1e-10)
 #
 #
-# == {}[link:index.html"name="9] Module Constants
+# == Module Constants
 # ---
 # * GSL::VERSION
 #
@@ -267,10 +267,10 @@
 #
 #   Ruby/GSL version
 #
-# {prev}[link:files/rdoc/ehandling_rdoc.html]
-# {next}[link:files/rdoc/complex_rdoc.html]
+# {prev}[link:rdoc/ehandling_rdoc.html]
+# {next}[link:rdoc/complex_rdoc.html]
 #
-# {Reference index}[link:files/rdoc/ref_rdoc.html]
-# {top}[link:files/rdoc/index_rdoc.html]
+# {Reference index}[link:rdoc/ref_rdoc.html]
+# {top}[link:index.html]
 #
 #

data/rdoc/matrix.rdoc

@@ -1,29 +1,29 @@
 #
 # = Matrices
 # Contents:
-# 1. {Class methods}[link:files/rdoc/matrix_rdoc.html#1]
-# 1. {Instance methods}[link:files/rdoc/matrix_rdoc.html#2]
-#    1. {Accessing matrix elements}[link:files/rdoc/matrix_rdoc.html#2.1]
-#    1. {Initializing matrix elements}[link:files/rdoc/matrix_rdoc.html#2.2]
-#    1. {IO}[link:files/rdoc/matrix_rdoc.html#2.3]
-#    1. {Matrix views}[link:files/rdoc/matrix_rdoc.html#2.4]
-#    1. {Creating row and column views}[link:files/rdoc/matrix_rdoc.html#2.5]
-#    1. {Iterators}[link:files/rdoc/matrix_rdoc.html#2.6]
-#    1. {Copying matrices}[link:files/rdoc/matrix_rdoc.html#2.7]
-#    1. {Copying rows and columns}[link:files/rdoc/matrix_rdoc.html#2.8]
-#    1. {Exchanging rows and columns}[link:files/rdoc/matrix_rdoc.html#2.9]
-#    1. {Matrix operations}[link:files/rdoc/matrix_rdoc.html#2.10]
-#    1. {Finding maximum and minimum elements of matrices}[link:files/rdoc/matrix_rdoc.html#2.11]
-#    1. {Matrix properties}[link:files/rdoc/matrix_rdoc.html#2.12]
-# 1. {NArray}[link:files/rdoc/matrix_rdoc.html#3]
-# 1. {Special matrices}[link:files/rdoc/matrix_rdoc.html#4]
-#
-# == {}[link:index.html"name="1] Class methods
+# 1. {Class methods}[link:rdoc/matrix_rdoc.html#label-Class+methods]
+# 1. {Instance methods}[link:rdoc/matrix_rdoc.html#label-Instance+Methods]
+#    1. {Accessing matrix elements}[link:rdoc/matrix_rdoc.html#label-Accessing+matrix+elements]
+#    1. {Initializing matrix elements}[link:rdoc/matrix_rdoc.html#label-Initializing+matrix+elements]
+#    1. {IO}[link:rdoc/matrix_rdoc.html#label-IO]
+#    1. {Matrix views}[link:rdoc/matrix_rdoc.html#label-Matrix+views]
+#    1. {Creating row and column views}[link:rdoc/matrix_rdoc.html#label-Creating+row+and+column+views]
+#    1. {Iterators}[link:rdoc/matrix_rdoc.html#label-Iterators]
+#    1. {Copying matrices}[link:rdoc/matrix_rdoc.html#label-Copying+matrices]
+#    1. {Copying rows and columns}[link:rdoc/matrix_rdoc.html#label-Copying+rows+and+columns]
+#    1. {Exchanging rows and columns}[link:rdoc/matrix_rdoc.html#label-Exchanging+rows+and+columns]
+#    1. {Matrix operations}[link:rdoc/matrix_rdoc.html#label-Matrix+operations]
+#    1. {Finding maximum and minimum elements of matrices}[link:rdoc/matrix_rdoc.html#label-Finding+maximum+and+minimum+elements+of+matrices]
+#    1. {Matrix properties}[link:rdoc/matrix_rdoc.html#label-Matrix+properties]
+# 1. {NArray}[link:rdoc/matrix_rdoc.html#label-NArray]
+# 1. {Special matrices}[link:rdoc/matrix_rdoc.html#label-Special+matrices]
+#
+# == Class methods
 #
 # ---
 # * GSL::Matrix.alloc(n)
 # * GSL::Matrix.alloc(size1, size2)
-# * GSL::Matrix.alloc(array) 
+# * GSL::Matrix.alloc(array)
 # * GSL::Matrix.alloc(arrays)
 # * GSL::Matrix.alloc( ... )
 # * GSL::Matrix[ ... ]
@@ -32,21 +32,21 @@
 #
 # 1. From arrays
 #        >> m = GSL::Matrix[[1, 2, 3], [4, 5, 6], [7, 8, 9]]
-#        => GSL::Matrix 
-#        [ 1.000e+00 2.000e+00 3.000e+00 
-#          4.000e+00 5.000e+00 6.000e+00 
+#        => GSL::Matrix
+#        [ 1.000e+00 2.000e+00 3.000e+00
+#          4.000e+00 5.000e+00 6.000e+00
 #          7.000e+00 8.000e+00 9.000e+00 ]
-#     
+#
 # 1. With an array and rows&cols,
 #           m = GSL::Matrix.alloc([1, 2, 3, 4, 5, 6, 7, 8, 9], 3, 3)
 #
 # 1. With Range objects,
 #           >> m = GSL::Matrix.alloc(1..3, 4..6, 7..9)
-#           [ 1.000e+00 2.000e+00 3.000e+00 
-#             4.000e+00 5.000e+00 6.000e+00 
+#           [ 1.000e+00 2.000e+00 3.000e+00
+#             4.000e+00 5.000e+00 6.000e+00
 #             7.000e+00 8.000e+00 9.000e+00 ]
 #           >> m2 = GSL::Matrix[1..6, 2, 3]
-#           [ 1.000e+00 2.000e+00 3.000e+00 
+#           [ 1.000e+00 2.000e+00 3.000e+00
 #             4.000e+00 5.000e+00 6.000e+00 ]
 #
 # ---
@@ -56,14 +56,14 @@
 #   Examples:
 #     >> m = GSL::Matrix::Int.eye(3)
 #     => GSL::Matrix::Int
-#     [ 1 0 0 
-#       0 1 0 
+#     [ 1 0 0
+#       0 1 0
 #       0 0 1 ]
 #     >> m = GSL::Matrix::Int.eye(2, 4)
 #     => GSL::Matrix::Int
-#     [ 1 0 0 0 
+#     [ 1 0 0 0
 #       0 1 0 0 ]
-#  
+#
 # ---
 # * GSL::Matrix.identity(n)
 # * GSL::Matrix.scalar(n)
@@ -71,7 +71,7 @@
 # * GSL::Matrix.I(n)
 #
 #   Create diagonal matrices of dimensions n*n, of values 1.0.
-#   
+#
 # ---
 # * GSL::Matrix.diagonal(a, b, c, ...)
 # * GSL::Matrix.diagonal(Ary)
@@ -83,14 +83,14 @@
 #   Example:
 #     >> GSL::Matrix::Int.diagonal(1..4)
 #     => GSL::Matrix::Int
-#     [ 1 0 0 0 
-#       0 2 0 0 
-#       0 0 3 0 
+#     [ 1 0 0 0
+#       0 2 0 0
+#       0 0 3 0
 #       0 0 0 4 ]
 #     >> GSL::Matrix::Int.diagonal(2, 5, 3)
 #     => GSL::Matrix::Int
-#     [ 2 0 0 
-#       0 5 0 
+#     [ 2 0 0
+#       0 5 0
 #       0 0 3 ]
 #
 # ---
@@ -111,26 +111,26 @@
 #   Example:
 #
 #     >> m = GSL::Matrix::Int.indgen(3, 5)
-#     => GSL::Matrix::Int 
-#     [  0  1  2  3  4 
-#        5  6  7  8  9 
+#     => GSL::Matrix::Int
+#     [  0  1  2  3  4
+#        5  6  7  8  9
 #       10 11 12 13 14 ]
 #     >> m = GSL::Matrix::Int.indgen(3, 5, 2)
-#     => GSL::Matrix::Int 
-#     [  2  3  4  5  6 
-#        7  8  9 10 11 
+#     => GSL::Matrix::Int
+#     [  2  3  4  5  6
+#        7  8  9 10 11
 #       12 13 14 15 16 ]
 #     >> m = GSL::Matrix.indgen(2, 3, 4.5, 6.7)
 #     => GSL::Matrix
-#     [  4.500e+00  1.120e+01  1.790e+01 
+#     [  4.500e+00  1.120e+01  1.790e+01
 #        2.460e+01  3.130e+01  3.800e+01 ]
 #
-# === {}[link:index.html"name="1.1] NOTE:
+# === NOTE:
 # Matrix dimensions are limited within the range of Fixnum.
 # For 32-bit CPU, the maximum of matrix dimension is 2^30 ~ 1e9.
 #
-# == {}[link:index.html"name="2] Instance Methods
-# === {}[link:index.html"name="2.1] Accessing matrix elements 
+# == Instance Methods
+# === Accessing matrix elements
 #
 # ---
 # * GSL::Matrix#size1
@@ -156,7 +156,7 @@
 #
 # ---
 # * GSL::Matrix#set(args, val)
-# * GSL::Matrix#[args]=val
+# * \GSL::Matrix#[args]=val
 #
 #   If <tt>args</tt> is empty and <tt>val</tt> is an Array (i.e. called with just a
 #   single Array argument), the Array's elements are taken as row contents.
@@ -173,7 +173,7 @@
 #   as indexes.
 #
 #   If <tt>args</tt> are two <tt>Fixnums</tt>, <tt>i</tt> and <tt>j</tt>, this method
-#   sets the <tt>(i,j)</tt>-th element of the matrix <tt>self</tt> to <tt>val</tt>. 
+#   sets the <tt>(i,j)</tt>-th element of the matrix <tt>self</tt> to <tt>val</tt>.
 #
 #   If <tt>args</tt> is a single <tt>Fixnum</tt>, <tt>i</tt>, this method sets the
 #   element at row <tt>i</tt>/<tt>size2</tt>, column <tt>i</tt>%<tt>size2</tt> to
@@ -194,18 +194,18 @@
 #   NOTE: GSL does not provide a matrix copy function that properly copies data
 #   across overlapping memory regions, so watch out if assigning to part of a
 #   Matrix from another part of itself (see <tt>#set</tt> example of
-#   {GSL::Vector}[link:files/rdoc/vector_rdoc.html]).
+#   {GSL::Vector}[link:rdoc/vector_rdoc.html]).
 #
 # ---
 # * GSL::Matrix#get(args)
-# * GSL::Matrix#[args]
+# * \GSL::Matrix#[args]
 #
 #   If <tt>args</tt> are two <tt>Fixnums</tt>, <tt>i</tt> and <tt>j</tt>, this method
-#   returns the <tt>(i,j)</tt>-th element of the matrix <tt>self</tt>. 
+#   returns the <tt>(i,j)</tt>-th element of the matrix <tt>self</tt>.
 #
 #   If <tt>args</tt> is a single <tt>Fixnum</tt>, <tt>i</tt>, this method returns the
 #   element at row <tt>i</tt>/<tt>size2</tt>, column <tt>i</tt>%<tt>size2</tt>.
-#   
+#
 #   All other forms of <tt>args</tt> are treated as with <tt>Matrix#submatrix</tt>
 #   and a View object is returned.
 #
@@ -226,8 +226,8 @@
 #   Examples:
 #     >> m = GSL::Matrix[1..9, 3, 3]
 #     => GSL::Matrix
-#     [ 1.000e+00 2.000e+00 3.000e+00 
-#       4.000e+00 5.000e+00 6.000e+00 
+#     [ 1.000e+00 2.000e+00 3.000e+00
+#       4.000e+00 5.000e+00 6.000e+00
 #       7.000e+00 8.000e+00 9.000e+00 ]
 #     >> m[1, 2]
 #     => 6.0
@@ -235,23 +235,23 @@
 #     => 123
 #     >> m
 #     => GSL::Matrix
-#     [ 1.000e+00 2.000e+00 3.000e+00 
-#       4.000e+00 5.000e+00 1.230e+02 
+#     [ 1.000e+00 2.000e+00 3.000e+00
+#       4.000e+00 5.000e+00 1.230e+02
 #       7.000e+00 8.000e+00 9.000e+00 ]
 #     >> m[1]
 #     => 2.0
 #     >> m.set([3, 5, 2], [4, 5, 3], [7, 1, 5])
 #     => GSL::Matrix
-#     [ 3.000e+00 5.000e+00 2.000e+00 
-#       4.000e+00 5.000e+00 3.000e+00 
+#     [ 3.000e+00 5.000e+00 2.000e+00
+#       4.000e+00 5.000e+00 3.000e+00
 #       7.000e+00 1.000e+00 5.000e+00 ]
 #     >> m[1][1]   # old/unsupported form
 #     NoMethodError: undefined method `[]' for 2.0:Float
 #             from (irb):8
 #     >> m = GSL::Matrix::Int[1..9, 3, 3]
 #     => GSL::Matrix::Int
-#     [ 1 2 3 
-#       4 5 6 
+#     [ 1 2 3
+#       4 5 6
 #       7 8 9 ]
 #     >> m[1]      # m[0,1]
 #     => 2
@@ -274,7 +274,7 @@
 #     >> GSL::Matrix.eye(3).to_a
 #     => [[1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0]]
 #
-# === {}[link:index.html"name="2.2] Initializing matrix elements 
+# === Initializing matrix elements
 # ---
 # * GSL::Matrix#set_all(x)
 #
@@ -283,16 +283,16 @@
 # ---
 # * GSL::Matrix#set_zero
 #
-#   This method sets all the elements of the matrix to zero. 
+#   This method sets all the elements of the matrix to zero.
 #
 # ---
 # * GSL::Matrix#set_identity
 #
-#   This method sets the elements of the matrix to the corresponding 
-#   elements of the identity matrix, i.e. a unit diagonal with all off-diagonal 
-#   elements zero. This applies to both square and rectangular matrices. 
+#   This method sets the elements of the matrix to the corresponding
+#   elements of the identity matrix, i.e. a unit diagonal with all off-diagonal
+#   elements zero. This applies to both square and rectangular matrices.
 #
-# === {}[link:index.html"name="2.3] IO
+# === IO
 # ---
 # * GSL::Matrix#fwrite(io)
 # * GSL::Matrix#fwrite(filename)
@@ -304,7 +304,7 @@
 # * GSL::Matrix#fscanf(filename)
 #
 #
-# === {}[link:index.html"name="2.4] Matrix views
+# === Matrix views
 # The <tt>GSL::Matrix::View</tt> class is defined to be used as "references" to
 # matrices. The <tt>Matrix::View</tt> class is a subclass of <tt>Matrix</tt>, and an
 # instance of the <tt>View</tt> class created by slicing a <tt>Matrix</tt> object can
@@ -394,97 +394,97 @@
 #
 #   Ex:
 #      >> v = Vector[1..9]
-#      => GSL::Vector 
+#      => GSL::Vector
 #      [ 1.000e+00 2.000e+00 3.000e+00 4.000e+00 5.000e+00 6.000e+00 7.000e+00 8.000e+00 9.000e+00 ]
 #      >> m = v.matrix_view(3, 3)
-#      => GSL::Matrix::View 
-#      [ 1.000e+00 2.000e+00 3.000e+00 
-#        4.000e+00 5.000e+00 6.000e+00 
+#      => GSL::Matrix::View
+#      [ 1.000e+00 2.000e+00 3.000e+00
+#        4.000e+00 5.000e+00 6.000e+00
 #        7.000e+00 8.000e+00 9.000e+00 ]
 #      >> m[1][1] = 99.99
 #      => 99.99
 #      >> v
-#      => GSL::Vector 
+#      => GSL::Vector
 #      [ 1.000e+00 2.000e+00 3.000e+00 4.000e+00 9.999e+01 6.000e+00 7.000e+00 8.000e+00 9.000e+00 ]
-#      >> 
+#      >>
 #
 #
-# === {}[link:index.html"name="2.5] Creating row and column views
+# === Creating row and column views
 #
 # ---
 # * GSL::Matrix#row(i)
 #
 #   These methods return <tt>i</tt>-th row of the matrix as a <tt>Vector::View</tt>
 #   object. Any modifications to the <tt>Vectror::View</tt> object returned by this method
-#   propagate to the original matrix. 
-#   
+#   propagate to the original matrix.
+#
 # ---
 # * GSL::Matrix#column(i)
-# * GSL::Matrix#col(i) 
+# * GSL::Matrix#col(i)
 #
 #   These methods return a vector view of the <tt>j</tt>-th column of the matrix.
 #
 # ---
 # * GSL::Matrix#subrow(i, offset, n)
 #
-#   Returns a vector view of the <tt>i</tt>-th row of the matrix <tt>self</tt> 
-#   beginning at <tt>offset</tt> elements past the first column 
+#   Returns a vector view of the <tt>i</tt>-th row of the matrix <tt>self</tt>
+#   beginning at <tt>offset</tt> elements past the first column
 #   and containing <tt>n</tt> elements. (>= GSL-1.10)
 #
 # ---
 # * GSL::Matrix#subcolumn(j, offset, n)
 #
-#   Returns a vector view of the <tt>j</tt>-th column of the matrix <tt>self</tt> 
-#   beginning at <tt>offset</tt> elements past the first row 
+#   Returns a vector view of the <tt>j</tt>-th column of the matrix <tt>self</tt>
+#   beginning at <tt>offset</tt> elements past the first row
 #   and containing <tt>n</tt> elements. (>= GSL-1.10)
 #
 # ---
 # * GSL::Matrix#diag
-# * GSL::Matrix#diagonal 
+# * GSL::Matrix#diagonal
 #
 #   This method returns a <tt>Vector::View</tt> of the diagonal of the matrix.
-#   The matrix is not required to be square. For a rectangular matrix the 
+#   The matrix is not required to be square. For a rectangular matrix the
 #   length of the diagonal is the same as the smaller dimension of the matrix.
 #
 #
 #   Ex:
 #     >> m = GSL::Matrix[1..9, 3, 3]
-#     => GSL::Matrix 
-#     [ 1.000e+00 2.000e+00 3.000e+00 
-#       4.000e+00 5.000e+00 6.000e+00 
+#     => GSL::Matrix
+#     [ 1.000e+00 2.000e+00 3.000e+00
+#       4.000e+00 5.000e+00 6.000e+00
 #       7.000e+00 8.000e+00 9.000e+00 ]
 #     >> m.row(1)
-#     => GSL::Vector::View 
+#     => GSL::Vector::View
 #     [ 4.000e+00 5.000e+00 6.000e+00 ]
 #     >> m.col(2)
-#     => GSL::Vector::Col::View 
-#     [ 3.000e+00 
-#       6.000e+00 
+#     => GSL::Vector::Col::View
+#     [ 3.000e+00
+#       6.000e+00
 #       9.000e+00 ]
 #     >> m.col(2)[2] = 123
 #     => 123
 #     >> m
-#     => GSL::Matrix 
-#     [ 1.000e+00 2.000e+00 3.000e+00 
-#       4.000e+00 5.000e+00 6.000e+00 
+#     => GSL::Matrix
+#     [ 1.000e+00 2.000e+00 3.000e+00
+#       4.000e+00 5.000e+00 6.000e+00
 #       7.000e+00 8.000e+00 1.230e+02 ]
 #     >> m.diagonal
-#     => GSL::Vector::View: 
+#     => GSL::Vector::View:
 #     [ 1.000e+00 5.000e+00 1.230e+02 ]
 #
 # ---
 # * GSL::Matrix#subdiagonal(k)
 #
-#   Returns a vector view view of the <tt>k</tt>-th subdiagonal 
-#   of the matrix <tt>self</tt>. 
-#   The matrix is not required to be square. The diagonal of the matrix 
+#   Returns a vector view view of the <tt>k</tt>-th subdiagonal
+#   of the matrix <tt>self</tt>.
+#   The matrix is not required to be square. The diagonal of the matrix
 #   corresponds to k = 0.
 #
 # ---
 # * GSL::Matrix#superdiagonal(k)
 #
-#   Returns a vector view of the <tt>k</tt>-th superdiagonal of the matrix <tt>self</tt>. 
-#   The matrix is not required to be square. The diagonal of the matrix 
+#   Returns a vector view of the <tt>k</tt>-th superdiagonal of the matrix <tt>self</tt>.
+#   The matrix is not required to be square. The diagonal of the matrix
 #   corresponds to k = 0.
 #
 # ---
@@ -493,22 +493,22 @@
 #   Creates a <tt>GSL::Vector</tt> object "flattening" the rows of the matrix <tt>self</tt>.
 #
 #       >> m = GSL::Matrix[1..6, 2, 3]
-#       => GSL::Matrix 
-#       [ 1.000e+00 2.000e+00 3.000e+00 
+#       => GSL::Matrix
+#       [ 1.000e+00 2.000e+00 3.000e+00
 #         4.000e+00 5.000e+00 6.000e+00 ]
 #       >> m.to_v
-#       => GSL::Vector 
+#       => GSL::Vector
 #       [ 1.000e+00 2.000e+00 3.000e+00 4.000e+00 5.000e+00 6.000e+00 ]
 #
-# === {}[link:index.html"name="2.6] Iterators
+# === Iterators
 # ---
 # * GSL::Matrix#each_row
 #
-#   Iterator for each of rows in the matrix <tt>self</tt>. 
+#   Iterator for each of rows in the matrix <tt>self</tt>.
 # ---
 # * GSL::Matrix#each_col
 #
-#   Iterator for each of columns in the matrix <tt>self</tt>. 
+#   Iterator for each of columns in the matrix <tt>self</tt>.
 #
 # ---
 # * GSL::Matrix#collect { |item| .. }
@@ -517,7 +517,7 @@
 # * GSL::Matrix#map! { |item| .. }
 #
 #
-# === {}[link:index.html"name="2.7] Copying matrices
+# === Copying matrices
 # ---
 # * GSL::Matrix#clone
 # * GSL::Matrix#duplicate
@@ -529,35 +529,35 @@
 # * GSL::Matrix.swap(dest, src)
 #
 #
-# === {}[link:index.html"name="2.8] Copying rows and columns
+# === Copying rows and columns
 #
 # ---
 # * GSL::Matrix#get_row(i)
 #
-#   This method returns a new vector (not a view) which contains the elements 
+#   This method returns a new vector (not a view) which contains the elements
 #   of the <tt>i</tt>-th row of the matrix <tt>self</tt>.
 #
 # ---
 # * GSL::Matrix#get_col(j)
 #
-#   This method returns a new vector (not a view) which contains the elements of the <tt>j</tt>-th 
+#   This method returns a new vector (not a view) which contains the elements of the <tt>j</tt>-th
 #   column of the matrix <tt>self</tt>.
 #
 # ---
 # * GSL::Matrix#set_row(i, v)
 #
-#   This method copies the elements of the vector <tt>v</tt> into the <tt>i</tt>-th 
-#   row of the matrix. 
-#   The length of the vector must be the same as the length of the row. 
+#   This method copies the elements of the vector <tt>v</tt> into the <tt>i</tt>-th
+#   row of the matrix.
+#   The length of the vector must be the same as the length of the row.
 #
 # ---
 # * GSL::Matrix#set_col(j, v)
 #
-#   This method copies the elements of the vector <tt>v</tt> into the <tt>j</tt>-th 
-#   column of the matrix. The length of the vector must be the same as the length 
-#   of the column. 
+#   This method copies the elements of the vector <tt>v</tt> into the <tt>j</tt>-th
+#   column of the matrix. The length of the vector must be the same as the length
+#   of the column.
 #
-# === {}[link:index.html"name="2.9] Exchanging rows and columns
+# === Exchanging rows and columns
 # ---
 # * GSL::Matrix#swap_rows!(i, j)
 #
@@ -570,7 +570,7 @@
 # ---
 # * GSL::Matrix#swap_columns!(i, j)
 #
-#   This method exchanges the <tt>i</tt>-th and <tt>j</tt>-th columns of the matrix <tt>in-place</tt>. 
+#   This method exchanges the <tt>i</tt>-th and <tt>j</tt>-th columns of the matrix <tt>in-place</tt>.
 # ---
 # * GSL::Matrix#swap_columns(i, j)
 #
@@ -592,9 +592,9 @@
 # ---
 # * GSL::Matrix#transpose!
 #
-#   This method replaces the matrix by its transpose by copying the 
+#   This method replaces the matrix by its transpose by copying the
 #   elements of the matrix <tt>in-place</tt>. The matrix must be square for this
-#   operation to be possible. 
+#   operation to be possible.
 #
 # ---
 # * GSL::Matrix#reverse_rows
@@ -602,14 +602,14 @@
 #
 #   Example:
 #       >> m = GSL::Matrix::Int[1..9, 3, 3]
-#       => GSL::Matrix::Int 
-#       [ 1 2 3 
-#         4 5 6 
+#       => GSL::Matrix::Int
+#       [ 1 2 3
+#         4 5 6
 #         7 8 9 ]
 #       >> m.reverse_rows
-#       => GSL::Matrix::Int 
-#       [ 7 8 9 
-#         4 5 6 
+#       => GSL::Matrix::Int
+#       [ 7 8 9
+#         4 5 6
 #         1 2 3 ]
 #
 # ---
@@ -619,74 +619,74 @@
 #   Example:
 #       >> m = GSL::Matrix::Int[1..9, 3, 3]
 #       => GSL::Matrix::Int
-#       [ 1 2 3 
-#         4 5 6 
+#       [ 1 2 3
+#         4 5 6
 #         7 8 9 ]
 #       >> m.reverse_rows.reverse_columns
-#       => GSL::Matrix::Int 
-#       [ 9 8 7 
-#         6 5 4 
+#       => GSL::Matrix::Int
+#       [ 9 8 7
+#         6 5 4
 #         3 2 1 ]
 #
 # ---
 # * GSL::Matrix#rot90(n = 1)
 #
-#   Return a copy of <tt>self</tt> with the elements rotated 
-#   counterclockwise in 90-degree increments. The argument <tt>n</tt> is 
-#   optional, and specifies how many 90-degree rotations are to be applied 
-#   (the default value is 1). 
+#   Return a copy of <tt>self</tt> with the elements rotated
+#   counterclockwise in 90-degree increments. The argument <tt>n</tt> is
+#   optional, and specifies how many 90-degree rotations are to be applied
+#   (the default value is 1).
 #   Negative values of <tt>n</tt> rotate the matrix in a clockwise direction.
 #
 #   Examples:
 #     >> m = GSL::Matrix::Int[1..6, 2, 3]
 #     => GSL::Matrix::Int
-#     [ 1 2 3 
+#     [ 1 2 3
 #       4 5 6 ]
 #     >> m.rot90
 #     => GSL::Matrix::Int
-#     [ 3 6 
-#       2 5 
+#     [ 3 6
+#       2 5
 #       1 4 ]
 #     >> m.rot90(2)
 #     => GSL::Matrix::Int
-#     [ 6 5 4 
+#     [ 6 5 4
 #       3 2 1 ]
 #     >> m.rot90(3)
 #     => GSL::Matrix::Int
-#     [ 4 1 
-#       5 2 
+#     [ 4 1
+#       5 2
 #       6 3 ]
 #     >> m.rot90(-1)
 #     => GSL::Matrix::Int
-#     [ 4 1 
-#       5 2 
+#     [ 4 1
+#       5 2
 #       6 3 ]
 #
 # ---
 # * GSL::Matrix#upper
 #
-#   This creates a matrix copying the upper half part of the matrix 
+#   This creates a matrix copying the upper half part of the matrix
 #   <tt>self</tt>, including the diagonal elements.
 # ---
 # * GSL::Matrix#lower
 #
-#   This creates a matrix copying the lower half part of the matrix 
+#   This creates a matrix copying the lower half part of the matrix
 #   <tt>self</tt>, including the diagonal elements.
 #
 #        >> m = GSL::Matrix[1..9, 3, 3]
 #        => GSL::Matrix
-#        [ 1.000e+00 2.000e+00 3.000e+00 
-#          4.000e+00 5.000e+00 6.000e+00 
+#        [ 1.000e+00 2.000e+00 3.000e+00
+#          4.000e+00 5.000e+00 6.000e+00
 #          7.000e+00 8.000e+00 9.000e+00 ]
 #        >> m.upper
-#        => GSL::Matrix 
-#        [ 1.000e+00 2.000e+00 3.000e+00 
-#          0.000e+00 5.000e+00 6.000e+00 
+#        => GSL::Matrix
+#        [ 1.000e+00 2.000e+00 3.000e+00
+#          0.000e+00 5.000e+00 6.000e+00
 #          0.000e+00 0.000e+00 9.000e+00 ]
 #        >> m.lower
-#        => GSL::Matrix 
-#        [ 1.000e+00 0.000e+00 0.000e+00 
-#          4.000e+00 5.000e+00 0.000e+00 
+#        => GSL::Matrix
+#        [ 1.000e+00 0.000e+00 0.000e+00
+#          4.000e+00 5.000e+00 0.000e+00
 #          7.000e+00 8.000e+00 9.000e+00 ]
 #
 # ---
@@ -699,49 +699,49 @@
 #     => true
 #     >> a = GSL::Matrix::Int[1..4, 2, 2]
 #     => GSL::Matrix::Int
-#     [ 1 2 
+#     [ 1 2
 #       3 4 ]
 #     >> b = GSL::Matrix::Int[5..10, 2, 3]
 #     => GSL::Matrix::Int
-#     [  5  6  7 
+#     [  5  6  7
 #        8  9 10 ]
 #     >> a.horzcat(b)
 #     => GSL::Matrix::Int
-#     [  1  2  5  6  7 
+#     [  1  2  5  6  7
 #        3  4  8  9 10 ]
 #
 # ---
 # * GSL::Matrix#vertcat(other)
 #
 #   Returns the vertical concatenation of <tt>self</tt> and <tt>other</tt>.
-#  
+#
 #   Ex:
 #     >> a = GSL::Matrix::Int[1..4, 2, 2]
 #     => GSL::Matrix::Int
-#     [ 1 2 
+#     [ 1 2
 #       3 4 ]
 #     >> b = GSL::Matrix::Int[5..10, 3, 2]
 #     => GSL::Matrix::Int
-#     [  5  6 
-#        7  8 
+#     [  5  6
+#        7  8
 #        9 10 ]
 #     >> a.vertcat(b)
 #     => GSL::Matrix::Int
-#     [  1  2 
-#        3  4 
-#        5  6 
-#        7  8 
+#     [  1  2
+#        3  4
+#        5  6
+#        7  8
 #        9 10 ]
 #
-# === {}[link:index.html"name="2.10] Matrix operations
+# === Matrix operations
 #
 # ---
 # * GSL::Matrix#add(b)
 # * GSL::Matrix#+(b)
 #
-#   This method adds the elements of matrix <tt>b</tt> 
-#   to the elements of the  matrix. 
-#   The two matrices must have the same dimensions. 
+#   This method adds the elements of matrix <tt>b</tt>
+#   to the elements of the  matrix.
+#   The two matrices must have the same dimensions.
 #
 #   If <tt>b</tt> is a scalar, these methods add it to all the elements
 #   of the matrix <tt>self</tt> (equivalent to the method <tt>add_constant</tt>).
@@ -750,29 +750,29 @@
 # * GSL::Matrix#sub(b)
 # * GSL::Matrix#-(b)
 #
-#   This method subtracts the elements of matrix <tt>b</tt> 
-#   from the elements of the  
-#   matrix. The two matrices must have the same dimensions. 
+#   This method subtracts the elements of matrix <tt>b</tt>
+#   from the elements of the
+#   matrix. The two matrices must have the same dimensions.
 #
 # ---
 # * GSL::Matrix#mul_elements(b)
 #
-#   This method multiplies the elements of the matrix by the elements of 
-#   matrix <tt>b</tt>. The two matrices must have the same dimensions. 
-#   If <tt>b</tt> is a scalar, the method <tt>scale</tt> (see below) 
+#   This method multiplies the elements of the matrix by the elements of
+#   matrix <tt>b</tt>. The two matrices must have the same dimensions.
+#   If <tt>b</tt> is a scalar, the method <tt>scale</tt> (see below)
 #   is called.
 #
 # ---
 # * GSL::Matrix#div_elements(b)
 #
 #
-#   This method divides the elements of the  matrix by the elements of 
-#   matrix <tt>b</tt>. The two matrices must have the same dimensions. 
+#   This method divides the elements of the  matrix by the elements of
+#   matrix <tt>b</tt>. The two matrices must have the same dimensions.
 #
 # ---
 # * GSL::Matrix#scale(x)
 #
-#   This method multiplies the elements of the  matrix by the constant 
+#   This method multiplies the elements of the  matrix by the constant
 #   factor <tt>x</tt>.
 #
 # ---
@@ -783,32 +783,32 @@
 # ---
 # * GSL::Matrix#*(b)
 #
-#   Matrix multiplication. 
+#   Matrix multiplication.
 #
 #   Ex:
 #
 #      >> a = GSL::Matrix[1..4, 2, 2]
-#      => GSL::Matrix 
-#      [ 1.000e+00 2.000e+00 
+#      => GSL::Matrix
+#      [ 1.000e+00 2.000e+00
 #        3.000e+00 4.000e+00 ]
 #      >> b = GSL::Matrix[5..8, 2, 2]
-#      => GSL::Matrix 
-#      [ 5.000e+00 6.000e+00 
+#      => GSL::Matrix
+#      [ 5.000e+00 6.000e+00
 #        7.000e+00 8.000e+00 ]
 #      >> a*b
-#      => GSL::Matrix 
-#      [ 1.900e+01 2.200e+01 
+#      => GSL::Matrix
+#      [ 1.900e+01 2.200e+01
 #        4.300e+01 5.000e+01 ]
 #      >> a*2
-#      => GSL::Matrix 
-#      [ 2.000e+00 4.000e+00 
+#      => GSL::Matrix
+#      [ 2.000e+00 4.000e+00
 #        6.000e+00 8.000e+00 ]
 #      >> c = Vector[1, 2]
-#      => GSL::Vector 
+#      => GSL::Vector
 #      [ 1.000e+00 2.000e+00 ]
 #      >> a*c.col
-#      => GSL::Vector::Col 
-#      [ 5.000e+00 
+#      => GSL::Vector::Col
+#      [ 5.000e+00
 #        1.100e+01 ]
 #
 # ---
@@ -819,26 +819,26 @@
 #   If a <tt>Vector::Col</tt> is given, this method solves the linear system
 #   by using LU decomposition.
 #
-#   Ex: 
+#   Ex:
 #       >> m = GSL::Matrix[1..4, 2, 2]
-#       => GSL::Matrix 
-#       [ 1.000e+00 2.000e+00   
+#       => GSL::Matrix
+#       [ 1.000e+00 2.000e+00
 #         3.000e+00 4.000e+00 ]
 #       >> m/3
-#       => GSL::Matrix 
+#       => GSL::Matrix
 #       [ 3.333e-01 6.667e-01         <--- 1/3, 2/3
 #         1.000e+00 1.333e+00 ]       <--- 3/3, 4/3
 #       >> b = Vector[5, 6].col
-#       => GSL::Vector::Col 
-#       [ 5.000e+00   
+#       => GSL::Vector::Col
+#       [ 5.000e+00
 #         6.000e+00 ]
 #       >> x = m/b          <--- Solve m (x,y) = b
-#       => GSL::Vector::Col 
+#       => GSL::Vector::Col
 #       [ -4.000e+00                  <--- x = -4
-#         4.500e+00 ]                 <--- y = 4.5 
+#         4.500e+00 ]                 <--- y = 4.5
 #       >> m*x
 #       => GSL::Vector::Col
-#       [  5.000e+00 
+#       [  5.000e+00
 #          6.000e+00 ]
 #
 # ---
@@ -846,7 +846,7 @@
 #
 #   Computes matrix power of <tt>b</tt>.
 #
-# === {}[link:index.html"name="2.11] Finding maximum and minimum elements of matrices
+# === Finding maximum and minimum elements of matrices
 #
 # ---
 # * GSL::Matrix#max
@@ -855,73 +855,73 @@
 #   These methods return the max/min value in the  matrix.
 #
 # ---
-# * GSL::Matrix#minmax 
+# * GSL::Matrix#minmax
 #
-#   This method returns a two elements array [min, max], 
+#   This method returns a two elements array [min, max],
 #   which contains the minimum
 #   and the maximum values in the  matrix.
 #
 # ---
-# * GSL::Matrix#max_index 
-# * GSL::Matrix#min_index 
+# * GSL::Matrix#max_index
+# * GSL::Matrix#min_index
 #
 #   These methods return the index of the max/min value in the matrix.
 #
 # ---
-# * GSL::Matrix#minmax_index 
+# * GSL::Matrix#minmax_index
 #
-#   This method returns a two elements array [imin, imax], 
+#   This method returns a two elements array [imin, imax],
 #   which contains the indices
 #   of the minimum and the maximum value in the matrix.
 #
-# === {}[link:index.html"name="2.12] Matrix properties
+# === Matrix properties
 # ---
 # * GSL::Matrix#isnull
 #
-#   This returns 1 if all the elements of the matrix <tt>self</tt> are zero, 
+#   This returns 1 if all the elements of the matrix <tt>self</tt> are zero,
 #   and 0 otherwise.
 #
 # ---
 # * GSL::Matrix#isnull?
 #
-#   This returns <tt>true</tt> if all the elements of the matrix <tt>self</tt> 
+#   This returns <tt>true</tt> if all the elements of the matrix <tt>self</tt>
 #   are zero, and <tt>false</tt> otherwise.
 #
 # ---
 # * GSL::Matrix#ispos
 # * GSL::Matrix#ispos?
 #
-#   (GSL-1.9 or later) Return 1 (true) if all the elements of the matrix <tt>self</tt> are strictly positive, and 0 (false) otherwise. 
+#   (GSL-1.9 or later) Return 1 (true) if all the elements of the matrix <tt>self</tt> are strictly positive, and 0 (false) otherwise.
 #
 # ---
 # * GSL::Matrix#isneg
 # * GSL::Matrix#isneg?
 #
-#   (GSL-1.9 or later) Return 1 (true) if all the elements of the matrix <tt>self</tt> are strictly negative, and 0 (false) otherwise. 
+#   (GSL-1.9 or later) Return 1 (true) if all the elements of the matrix <tt>self</tt> are strictly negative, and 0 (false) otherwise.
 #
 # ---
 # * GSL::Matrix#isnonneg
 # * GSL::Matrix#isnonneg?
 #
-#   (GSL-1.10 or later) Return 1 (true) if all the elements of the matrix <tt>self</tt> are non-negative , and 0 (false) otherwise. 
+#   (GSL-1.10 or later) Return 1 (true) if all the elements of the matrix <tt>self</tt> are non-negative , and 0 (false) otherwise.
 #
 # ---
 # * GSL::Matrix#any
 #
-#   Returns a Vector of ones and zeros with each element indicating 
-#   whether any of the elements of the corresponding column of the 
+#   Returns a Vector of ones and zeros with each element indicating
+#   whether any of the elements of the corresponding column of the
 #   matrix are nonzero.
 #
 # ---
 # * GSL::Matrix#all
 #
-#   Behaves like the method <tt>any</tt>, except that it returns 1 only if 
+#   Behaves like the method <tt>any</tt>, except that it returns 1 only if
 #   all the elements of the matrix.
 #
 # ---
 # * GSL:Matrix#trace
 #
-#   This returns trace of the matrix <tt>self</tt>, the sum of the diagonal 
+#   This returns trace of the matrix <tt>self</tt>, the sum of the diagonal
 #   elements.
 #
 # ---
@@ -944,13 +944,13 @@
 #   Example:
 #     >> m = GSL::Matrix::Int[-5..4, 3, 3]
 #     => GSL::Matrix::Int
-#     [ -5 -4 -3 
-#       -2 -1  0 
+#     [ -5 -4 -3
+#       -2 -1  0
 #        1  2  3 ]
 #     >> m.abs
 #     => GSL::Matrix::Int
-#     [ 5 4 3 
-#       2 1 0 
+#     [ 5 4 3
+#       2 1 0
 #       1 2 3 ]
 #
 # ---
@@ -958,15 +958,15 @@
 # * GSL::Matrix#==(other, eps = 1e-10)
 #
 #   Returns <tt>true</tt> if the matrices have same size and elements
-#   equal to absolute accurary <tt>eps</tt> for all the indices, 
+#   equal to absolute accurary <tt>eps</tt> for all the indices,
 #   and <tt>false</tt> otherwise.
 #
-# == {}[link:index.html"name="3] NArray
+# == NArray
 #
 # ---
 # * GSL::Matrix#to_na
 #
-#   The Matrix object <tt>self</tt> is converted into an <tt>NMatrix</tt> object. 
+#   The Matrix object <tt>self</tt> is converted into an <tt>NMatrix</tt> object.
 #   The matrix data are copied to newly allocated memory.
 #
 # ---
@@ -979,54 +979,54 @@
 # * NArray#to_gm_view
 # * NArray#to_gslm_view
 #
-#   A <tt>GSL::Matrix::View</tt> object is created from the NArray object <tt>na</tt>. 
-#   The data of <tt>na</tt> are 
-#   not copied, thus any modifications to the View object affect on the original 
-#   NArray object <tt>na</tt>. 
+#   A <tt>GSL::Matrix::View</tt> object is created from the NArray object <tt>na</tt>.
+#   The data of <tt>na</tt> are
+#   not copied, thus any modifications to the View object affect on the original
+#   NArray object <tt>na</tt>.
 #   The View object can be used as a reference to the NMatrix object.
 #
-# == {}[link:index.html"name="4] Special matrices
+# == Special matrices
 # ---
-# * GSL::Matrix.hirbert(n)
+# * GSL::Matrix.hilbert(n)
 #
 #   Returns the Hilbert matrix of order <tt>n</tt>. The <tt>ij</tt> element is
 #   defined as 1/(i+j+1).
 #
 # ---
-# * GSL::Matrix.invhirbert(n)
+# * GSL::Matrix.invhilbert(n)
 #
 #   Returns the inverse of a Hilbert matrix of order <tt>n</tt>.
 #
 #   Ex:
 #     >> m = GSL::Matrix.hilbert(4)
-#     => GSL::Matrix 
-#     [ 1.000e+00 5.000e-01 3.333e-01 2.500e-01 
-#       5.000e-01 3.333e-01 2.500e-01 2.000e-01 
-#       3.333e-01 2.500e-01 2.000e-01 1.667e-01 
+#     => GSL::Matrix
+#     [ 1.000e+00 5.000e-01 3.333e-01 2.500e-01
+#       5.000e-01 3.333e-01 2.500e-01 2.000e-01
+#       3.333e-01 2.500e-01 2.000e-01 1.667e-01
 #       2.500e-01 2.000e-01 1.667e-01 1.429e-01 ]
 #     >> invm = GSL::Matrix.invhilbert(4)
-#     => GSL::Matrix 
-#     [ 1.600e+01 -1.200e+02 2.400e+02 -1.400e+02 
-#       -1.200e+02 1.200e+03 -2.700e+03 1.680e+03 
-#       2.400e+02 -2.700e+03 6.480e+03 -4.200e+03 
+#     => GSL::Matrix
+#     [ 1.600e+01 -1.200e+02 2.400e+02 -1.400e+02
+#       -1.200e+02 1.200e+03 -2.700e+03 1.680e+03
+#       2.400e+02 -2.700e+03 6.480e+03 -4.200e+03
 #       -1.400e+02 1.680e+03 -4.200e+03 2.800e+03 ]
 #     >> invm2 = m.inv
-#     => GSL::Matrix 
-#     [ 1.600e+01 -1.200e+02 2.400e+02 -1.400e+02 
-#       -1.200e+02 1.200e+03 -2.700e+03 1.680e+03 
-#       2.400e+02 -2.700e+03 6.480e+03 -4.200e+03 
+#     => GSL::Matrix
+#     [ 1.600e+01 -1.200e+02 2.400e+02 -1.400e+02
+#       -1.200e+02 1.200e+03 -2.700e+03 1.680e+03
+#       2.400e+02 -2.700e+03 6.480e+03 -4.200e+03
 #       -1.400e+02 1.680e+03 -4.200e+03 2.800e+03 ]
 #     >> m*invm
-#     => GSL::Matrix 
-#     [ 1.000e+00 5.684e-14 -2.274e-13 1.137e-13 
-#       1.998e-15 1.000e+00 -4.663e-14 3.109e-14 
-#       3.664e-15 -7.239e-14 1.000e+00 -1.017e-13 
+#     => GSL::Matrix
+#     [ 1.000e+00 5.684e-14 -2.274e-13 1.137e-13
+#       1.998e-15 1.000e+00 -4.663e-14 3.109e-14
+#       3.664e-15 -7.239e-14 1.000e+00 -1.017e-13
 #       -2.442e-15 1.510e-14 -8.038e-14 1.000e+00 ]
 #     >> m*invm2
-#     => GSL::Matrix 
-#     [ 1.000e+00 0.000e+00 0.000e+00 0.000e+00 
-#       -1.554e-15 1.000e+00 -2.389e-14 8.349e-15 
-#       1.295e-15 3.405e-15 1.000e+00 -6.957e-15 
+#     => GSL::Matrix
+#     [ 1.000e+00 0.000e+00 0.000e+00 0.000e+00
+#       -1.554e-15 1.000e+00 -2.389e-14 8.349e-15
+#       1.295e-15 3.405e-15 1.000e+00 -6.957e-15
 #       1.110e-15 1.916e-14 1.707e-14 1.000e+00 ]
 #
 # ---
@@ -1035,23 +1035,23 @@
 #   Returns the Pascal matrix of order <tt>n</tt>, created from Pascal's triangle.
 #
 #     >> GSL::Matrix::Int.pascal(10)
-#     => GSL::Matrix::Int 
-#     [     1     1     1     1     1     1     1     1     1     1 
-#           1     2     3     4     5     6     7     8     9    10 
-#           1     3     6    10    15    21    28    36    45    55 
-#           1     4    10    20    35    56    84   120   165   220 
-#           1     5    15    35    70   126   210   330   495   715 
-#           1     6    21    56   126   252   462   792  1287  2002 
-#           1     7    28    84   210   462   924  1716  3003  5005 
-#           1     8    36   120   330   792  1716  3432  6435 11440 
-#           1     9    45   165   495  1287  3003  6435 12870 24310 
+#     => GSL::Matrix::Int
+#     [     1     1     1     1     1     1     1     1     1     1
+#           1     2     3     4     5     6     7     8     9    10
+#           1     3     6    10    15    21    28    36    45    55
+#           1     4    10    20    35    56    84   120   165   220
+#           1     5    15    35    70   126   210   330   495   715
+#           1     6    21    56   126   252   462   792  1287  2002
+#           1     7    28    84   210   462   924  1716  3003  5005
+#           1     8    36   120   330   792  1716  3432  6435 11440
+#           1     9    45   165   495  1287  3003  6435 12870 24310
 #           1    10    55   220   715  2002  5005 11440 24310 48620 ]
 #
 # ---
 # * GSL::Matrix.vandermonde(v)
 #
 #   Creates a Vendermonde matrix from a vector or an array <tt>v</tt>.
-#  
+#
 #      >> GSL::Matrix.vander([1, 2, 3, 4])
 #      => GSL::Matrix
 #      [ 1.000e+00 1.000e+00 1.000e+00 1.000e+00
@@ -1063,31 +1063,31 @@
 # * GSL::Matrix.toeplitz(v)
 #
 #   Creates a Toeplitz matrix from a vector or an array <tt>v</tt>.
-#  
+#
 #       >> GSL::Matrix::Int.toeplitz([1, 2, 3, 4, 5])
-#       => GSL::Matrix::Int 
-#       [ 1 2 3 4 5 
-#         2 1 2 3 4 
-#         3 2 1 2 3 
-#         4 3 2 1 2 
+#       => GSL::Matrix::Int
+#       [ 1 2 3 4 5
+#         2 1 2 3 4
+#         3 2 1 2 3
+#         4 3 2 1 2
 #         5 4 3 2 1 ]
 #
 # ---
 # * GSL::Matrix.circulant(v)
 #
 #   Creates a circulant matrix from a vector or an array <tt>v</tt>.
-#  
+#
 #       >> GSL::Matrix::Int.circulant([1, 2, 3, 4])
-#       => GSL::Matrix::Int 
-#       [ 4 1 2 3 
-#         3 4 1 2 
-#         2 3 4 1 
+#       => GSL::Matrix::Int
+#       [ 4 1 2 3
+#         3 4 1 2
+#         2 3 4 1
 #         1 2 3 4 ]
 #
-# {prev}[link:files/rdoc/vector_rdoc.html]
-# {next}[link:files/rdoc/perm_rdoc.html]
+# {prev}[link:rdoc/vector_rdoc.html]
+# {next}[link:rdoc/perm_rdoc.html]
+#
+# {Reference index}[link:rdoc/ref_rdoc.html]
+# {top}[link:index.html]
 #
-# {Reference index}[link:files/rdoc/ref_rdoc.html]
-# {top}[link:files/rdoc/index_rdoc.html]
-#   
 #