ruby-lapack 1.3 → 1.4

data/Rakefile

@@ -40,7 +40,7 @@ DLLIB = "ext/#{NAME}.so"
 task :default => DLLIB
 
 desc "building extensions"
-file DLLIB do
+file DLLIB => "ext/Makefile" do
   system("cd ext; make")
 end
 so_file = File.join("lib", target_prefix, "#{NAME}.so")
@@ -48,6 +48,9 @@ file so_file => DLLIB do
   mkdir File.dirname(so_file)
   cp DLLIB, so_file
 end
+file "ext/Makefile" do
+  system("cd ext; ruby extconf.rb")
+end
 
 desc "install files to system"
 task :install => [:install_so, :install_rb]
@@ -77,7 +80,7 @@ TEST_FILES = FileList["tests/**/*.rb"]
 
 spec = Gem::Specification.new do |s|
   s.name = "ruby-lapack"
-  s.version = "1.3"
+  s.version = "1.4"
   s.summary = "A Ruby wrapper of Lapack"
   s.description = <<EOL
 Ruby-LAPACK is a Ruby wrapper of Lapack, which is a linear algebra package (http://www.netlib.org/lapack/).

data/dev/make_csrc.rb

@@ -704,7 +704,11 @@ EOF
 
   code << <<"EOF"
 void
-init_lapack_#{sub_name}(VALUE mLapack){
+init_lapack_#{sub_name}(VALUE mLapack, VALUE sH, VALUE sU, VALUE zero){
+  sHelp = sH;
+  sUsage = sU;
+  #{RBPREFIX}ZERO = zero;
+
   rb_define_module_function(mLapack, \"#{sub_name}\", #{RBPREFIX}#{sub_name}, -1);
 }
 EOF
@@ -790,8 +794,8 @@ extern integer ilatrans_(char* trans);
 extern integer ilaenv_(integer* ispec, char* name, char* opts, integer* n1, integer* n2, integer* n3, integer* n4);
 
 
-VALUE sHelp, sUsage;
-VALUE #{RBPREFIX}ZERO;
+static VALUE sHelp, sUsage;
+static VALUE #{RBPREFIX}ZERO;
 
 EOF
   }
@@ -805,7 +809,7 @@ EOF
 EOF
 
     sub_names.each{|sname|
-      file.print "extern void init_lapack_#{sname}(VALUE mLapack);\n"
+      file.print "extern void init_lapack_#{sname}(VALUE mLapack, VALUE sHelp, VALUE sUsage, VALUE #{RBPREFIX}ZERO);\n"
     }
 
     file.print <<"EOF"
@@ -826,7 +830,7 @@ void Init_lapack(){
 
 EOF
     sub_names.each{|sname|
-      file.print "  init_lapack_#{sname}(mLapack);\n"
+      file.print "  init_lapack_#{sname}(mLapack, sHelp, sUsage, #{RBPREFIX}ZERO);\n"
     }
     file.print "}\n"
   }

data/doc/dsy.html

@@ -451,7 +451,7 @@ FORTRAN MANUAL
     <H2>dsyevr</H2>
     <PRE>
 USAGE:
-  m, w, z, isuppz, work, iwork, info, a = NumRu::Lapack.dsyevr( jobz, range, uplo, a, vl, vu, il, iu, abstol, liwork, [:lwork => lwork, :usage => usage, :help => help])
+  m, w, z, isuppz, work, iwork, info, a = NumRu::Lapack.dsyevr( jobz, range, uplo, a, vl, vu, il, iu, abstol, [:lwork => lwork, :liwork => liwork, :usage => usage, :help => help])
 
 
 FORTRAN MANUAL

data/doc/others.html

@@ -6,158 +6,256 @@
     <A NAME="top"></A>
     <H1>other routines</H1>
     <UL>
-      <LI><A HREF="#ilaslc">ilaslc</A></LI>
-      <LI><A HREF="#ilatrans">ilatrans</A></LI>
-      <LI><A HREF="#ilazlc">ilazlc</A></LI>
-      <LI><A HREF="#iladlr">iladlr</A></LI>
-      <LI><A HREF="#ilaver">ilaver</A></LI>
-      <LI><A HREF="#iladiag">iladiag</A></LI>
+      <LI><A HREF="#iparmq">iparmq</A></LI>
       <LI><A HREF="#lsamen">lsamen</A></LI>
       <LI><A HREF="#ilazlr">ilazlr</A></LI>
-      <LI><A HREF="#iparmq">iparmq</A></LI>
-      <LI><A HREF="#xerbla">xerbla</A></LI>
+      <LI><A HREF="#iladiag">iladiag</A></LI>
+      <LI><A HREF="#ilaslc">ilaslc</A></LI>
+      <LI><A HREF="#ilaver">ilaver</A></LI>
+      <LI><A HREF="#ilaslr">ilaslr</A></LI>
+      <LI><A HREF="#ilauplo">ilauplo</A></LI>
       <LI><A HREF="#xerbla_array">xerbla_array</A></LI>
+      <LI><A HREF="#iladlr">iladlr</A></LI>
+      <LI><A HREF="#ilazlc">ilazlc</A></LI>
+      <LI><A HREF="#icmax1">icmax1</A></LI>
+      <LI><A HREF="#ilaclr">ilaclr</A></LI>
+      <LI><A HREF="#iladlc">iladlc</A></LI>
+      <LI><A HREF="#ilatrans">ilatrans</A></LI>
+      <LI><A HREF="#ilaenv">ilaenv</A></LI>
       <LI><A HREF="#ilaclc">ilaclc</A></LI>
+      <LI><A HREF="#xerbla">xerbla</A></LI>
+      <LI><A HREF="#ieeeck">ieeeck</A></LI>
       <LI><A HREF="#izmax1">izmax1</A></LI>
       <LI><A HREF="#ilaprec">ilaprec</A></LI>
-      <LI><A HREF="#ilaenv">ilaenv</A></LI>
-      <LI><A HREF="#ieeeck">ieeeck</A></LI>
-      <LI><A HREF="#iladlc">iladlc</A></LI>
-      <LI><A HREF="#ilauplo">ilauplo</A></LI>
-      <LI><A HREF="#icmax1">icmax1</A></LI>
-      <LI><A HREF="#ilaslr">ilaslr</A></LI>
-      <LI><A HREF="#ilaclr">ilaclr</A></LI>
     </UL>
 
-    <A NAME="ilaslc"></A>
-    <H2>ilaslc</H2>
+    <A NAME="iparmq"></A>
+    <H2>iparmq</H2>
     <PRE>
 USAGE:
-  __out__ = NumRu::Lapack.ilaslc( m, a, [:usage => usage, :help => help])
+  __out__ = NumRu::Lapack.iparmq( ispec, name, opts, n, ilo, ihi, lwork, [:usage => usage, :help => help])
 
 
 FORTRAN MANUAL
-      INTEGER FUNCTION ILASLC( M, N, A, LDA )
+      INTEGER FUNCTION IPARMQ( ISPEC, NAME, OPTS, N, ILO, IHI, LWORK )
 
 *  Purpose
 *  =======
 *
-*  ILASLC scans A for its last non-zero column.
+*       This program sets problem and machine dependent parameters
+*       useful for xHSEQR and its subroutines. It is called whenever 
+*       ILAENV is called with 12 <= ISPEC <= 16
 *
 
 *  Arguments
 *  =========
 *
-*  M       (input) INTEGER
-*          The number of rows of the matrix A.
+*       ISPEC  (input) integer scalar
+*              ISPEC specifies which tunable parameter IPARMQ should
+*              return.
 *
-*  N       (input) INTEGER
-*          The number of columns of the matrix A.
+*              ISPEC=12: (INMIN)  Matrices of order nmin or less
+*                        are sent directly to xLAHQR, the implicit
+*                        double shift QR algorithm.  NMIN must be
+*                        at least 11.
 *
-*  A       (input) REAL array, dimension (LDA,N)
-*          The m by n matrix A.
+*              ISPEC=13: (INWIN)  Size of the deflation window.
+*                        This is best set greater than or equal to
+*                        the number of simultaneous shifts NS.
+*                        Larger matrices benefit from larger deflation
+*                        windows.
 *
-*  LDA     (input) INTEGER
-*          The leading dimension of the array A. LDA >= max(1,M).
+*              ISPEC=14: (INIBL) Determines when to stop nibbling and
+*                        invest in an (expensive) multi-shift QR sweep.
+*                        If the aggressive early deflation subroutine
+*                        finds LD converged eigenvalues from an order
+*                        NW deflation window and LD.GT.(NW*NIBBLE)/100,
+*                        then the next QR sweep is skipped and early
+*                        deflation is applied immediately to the
+*                        remaining active diagonal block.  Setting
+*                        IPARMQ(ISPEC=14) = 0 causes TTQRE to skip a
+*                        multi-shift QR sweep whenever early deflation
+*                        finds a converged eigenvalue.  Setting
+*                        IPARMQ(ISPEC=14) greater than or equal to 100
+*                        prevents TTQRE from skipping a multi-shift
+*                        QR sweep.
 *
-
-*  =====================================================================
+*              ISPEC=15: (NSHFTS) The number of simultaneous shifts in
+*                        a multi-shift QR iteration.
 *
-
-
-    </PRE>
-    <A HREF="#top">go to the page top</A>
-
-    <A NAME="ilatrans"></A>
-    <H2>ilatrans</H2>
-    <PRE>
-USAGE:
-  __out__ = NumRu::Lapack.ilatrans( trans, [:usage => usage, :help => help])
-
-
-FORTRAN MANUAL
-      INTEGER FUNCTION ILATRANS( TRANS )
-
-*  Purpose
-*  =======
+*              ISPEC=16: (IACC22) IPARMQ is set to 0, 1 or 2 with the
+*                        following meanings.
+*                        0:  During the multi-shift QR sweep,
+*                            xLAQR5 does not accumulate reflections and
+*                            does not use matrix-matrix multiply to
+*                            update the far-from-diagonal matrix
+*                            entries.
+*                        1:  During the multi-shift QR sweep,
+*                            xLAQR5 and/or xLAQRaccumulates reflections and uses
+*                            matrix-matrix multiply to update the
+*                            far-from-diagonal matrix entries.
+*                        2:  During the multi-shift QR sweep.
+*                            xLAQR5 accumulates reflections and takes
+*                            advantage of 2-by-2 block structure during
+*                            matrix-matrix multiplies.
+*                        (If xTRMM is slower than xGEMM, then
+*                        IPARMQ(ISPEC=16)=1 may be more efficient than
+*                        IPARMQ(ISPEC=16)=2 despite the greater level of
+*                        arithmetic work implied by the latter choice.)
 *
-*  This subroutine translates from a character string specifying a
-*  transposition operation to the relevant BLAST-specified integer
-*  constant.
+*       NAME    (input) character string
+*               Name of the calling subroutine
 *
-*  ILATRANS returns an INTEGER.  If ILATRANS < 0, then the input is not
-*  a character indicating a transposition operator.  Otherwise ILATRANS
-*  returns the constant value corresponding to TRANS.
+*       OPTS    (input) character string
+*               This is a concatenation of the string arguments to
+*               TTQRE.
+*
+*       N       (input) integer scalar
+*               N is the order of the Hessenberg matrix H.
+*
+*       ILO     (input) INTEGER
+*       IHI     (input) INTEGER
+*               It is assumed that H is already upper triangular
+*               in rows and columns 1:ILO-1 and IHI+1:N.
+*
+*       LWORK   (input) integer scalar
+*               The amount of workspace available.
 *
 
-*  Arguments
-*  =========
-*  TRANS   (input) CHARACTER*1
-*          Specifies the form of the system of equations:
-*          = 'N':  No transpose
-*          = 'T':  Transpose
-*          = 'C':  Conjugate transpose
-
-*  =====================================================================
+*  Further Details
+*  ===============
+*
+*       Little is known about how best to choose these parameters.
+*       It is possible to use different values of the parameters
+*       for each of CHSEQR, DHSEQR, SHSEQR and ZHSEQR.
+*
+*       It is probably best to choose different parameters for
+*       different matrices and different parameters at different
+*       times during the iteration, but this has not been
+*       implemented --- yet.
+*
+*
+*       The best choices of most of the parameters depend
+*       in an ill-understood way on the relative execution
+*       rate of xLAQR3 and xLAQR5 and on the nature of each
+*       particular eigenvalue problem.  Experiment may be the
+*       only practical way to determine which choices are most
+*       effective.
+*
+*       Following is a list of default values supplied by IPARMQ.
+*       These defaults may be adjusted in order to attain better
+*       performance in any particular computational environment.
+*
+*       IPARMQ(ISPEC=12) The xLAHQR vs xLAQR0 crossover point.
+*                        Default: 75. (Must be at least 11.)
+*
+*       IPARMQ(ISPEC=13) Recommended deflation window size.
+*                        This depends on ILO, IHI and NS, the
+*                        number of simultaneous shifts returned
+*                        by IPARMQ(ISPEC=15).  The default for
+*                        (IHI-ILO+1).LE.500 is NS.  The default
+*                        for (IHI-ILO+1).GT.500 is 3*NS/2.
 *
+*       IPARMQ(ISPEC=14) Nibble crossover point.  Default: 14.
+*
+*       IPARMQ(ISPEC=15) Number of simultaneous shifts, NS.
+*                        a multi-shift QR iteration.
+*
+*                        If IHI-ILO+1 is ...
+*
+*                        greater than      ...but less    ... the
+*                        or equal to ...      than        default is
+*
+*                                0               30       NS =   2+
+*                               30               60       NS =   4+
+*                               60              150       NS =  10
+*                              150              590       NS =  **
+*                              590             3000       NS =  64
+*                             3000             6000       NS = 128
+*                             6000             infinity   NS = 256
+*
+*                    (+)  By default matrices of this order are
+*                         passed to the implicit double shift routine
+*                         xLAHQR.  See IPARMQ(ISPEC=12) above.   These
+*                         values of NS are used only in case of a rare
+*                         xLAHQR failure.
+*
+*                    (**) The asterisks (**) indicate an ad-hoc
+*                         function increasing from 10 to 64.
+*
+*       IPARMQ(ISPEC=16) Select structured matrix multiply.
+*                        (See ISPEC=16 above for details.)
+*                        Default: 3.
+*
+*     ================================================================
 
 
     </PRE>
     <A HREF="#top">go to the page top</A>
 
-    <A NAME="ilazlc"></A>
-    <H2>ilazlc</H2>
+    <A NAME="lsamen"></A>
+    <H2>lsamen</H2>
     <PRE>
 USAGE:
-  __out__ = NumRu::Lapack.ilazlc( m, a, [:usage => usage, :help => help])
+  __out__ = NumRu::Lapack.lsamen( n, ca, cb, [:usage => usage, :help => help])
 
 
 FORTRAN MANUAL
-      INTEGER FUNCTION ILAZLC( M, N, A, LDA )
+      LOGICAL          FUNCTION LSAMEN( N, CA, CB )
 
 *  Purpose
 *  =======
 *
-*  ILAZLC scans A for its last non-zero column.
+*  LSAMEN  tests if the first N letters of CA are the same as the
+*  first N letters of CB, regardless of case.
+*  LSAMEN returns .TRUE. if CA and CB are equivalent except for case
+*  and .FALSE. otherwise.  LSAMEN also returns .FALSE. if LEN( CA )
+*  or LEN( CB ) is less than N.
 *
 
 *  Arguments
 *  =========
 *
-*  M       (input) INTEGER
-*          The number of rows of the matrix A.
-*
 *  N       (input) INTEGER
-*          The number of columns of the matrix A.
-*
-*  A       (input) COMPLEX*16 array, dimension (LDA,N)
-*          The m by n matrix A.
+*          The number of characters in CA and CB to be compared.
 *
-*  LDA     (input) INTEGER
-*          The leading dimension of the array A. LDA >= max(1,M).
+*  CA      (input) CHARACTER*(*)
+*  CB      (input) CHARACTER*(*)
+*          CA and CB specify two character strings of length at least N.
+*          Only the first N characters of each string will be accessed.
 *
 
-*  =====================================================================
+* =====================================================================
 *
+*     .. Local Scalars ..
+      INTEGER            I
+*     ..
+*     .. External Functions ..
+      LOGICAL            LSAME
+      EXTERNAL           LSAME
+*     ..
+*     .. Intrinsic Functions ..
+      INTRINSIC          LEN
+*     ..
 
 
     </PRE>
     <A HREF="#top">go to the page top</A>
 
-    <A NAME="iladlr"></A>
-    <H2>iladlr</H2>
+    <A NAME="ilazlr"></A>
+    <H2>ilazlr</H2>
     <PRE>
 USAGE:
-  __out__ = NumRu::Lapack.iladlr( m, a, [:usage => usage, :help => help])
+  __out__ = NumRu::Lapack.ilazlr( m, a, [:usage => usage, :help => help])
 
 
 FORTRAN MANUAL
-      INTEGER FUNCTION ILADLR( M, N, A, LDA )
+      INTEGER FUNCTION ILAZLR( M, N, A, LDA )
 
 *  Purpose
 *  =======
 *
-*  ILADLR scans A for its last non-zero row.
+*  ILAZLR scans A for its last non-zero row.
 *
 
 *  Arguments
@@ -169,7 +267,7 @@ FORTRAN MANUAL
 *  N       (input) INTEGER
 *          The number of columns of the matrix A.
 *
-*  A       (input) DOUBLE PRECISION array, dimension (LDA,N)
+*  A       (input) COMPLEX*16 array, dimension (LDA,N)
 *          The m by n matrix A.
 *
 *  LDA     (input) INTEGER
@@ -180,47 +278,6 @@ FORTRAN MANUAL
 *
 
 
-    </PRE>
-    <A HREF="#top">go to the page top</A>
-
-    <A NAME="ilaver"></A>
-    <H2>ilaver</H2>
-    <PRE>
-USAGE:
-  vers_major, vers_minor, vers_patch = NumRu::Lapack.ilaver( , [:usage => usage, :help => help])
-
-
-FORTRAN MANUAL
-      SUBROUTINE ILAVER( VERS_MAJOR, VERS_MINOR, VERS_PATCH )
-
-*  Purpose
-*  =======
-*
-*  This subroutine return the Lapack version.
-*
-
-*  Arguments
-*  =========
-*  VERS_MAJOR   (output) INTEGER
-*      return the lapack major version
-*  VERS_MINOR   (output) INTEGER
-*      return the lapack minor version from the major version
-*  VERS_PATCH   (output) INTEGER
-*      return the lapack patch version from the minor version
-
-*  =====================================================================
-*
-      INTEGER VERS_MAJOR, VERS_MINOR, VERS_PATCH
-*  =====================================================================
-      VERS_MAJOR = 3
-      VERS_MINOR = 3
-      VERS_PATCH = 0
-*  =====================================================================
-*
-      RETURN
-      END
-
-
     </PRE>
     <A HREF="#top">go to the page top</A>
 
@@ -259,69 +316,20 @@ FORTRAN MANUAL
     </PRE>
     <A HREF="#top">go to the page top</A>
 
-    <A NAME="lsamen"></A>
-    <H2>lsamen</H2>
-    <PRE>
-USAGE:
-  __out__ = NumRu::Lapack.lsamen( n, ca, cb, [:usage => usage, :help => help])
-
-
-FORTRAN MANUAL
-      LOGICAL          FUNCTION LSAMEN( N, CA, CB )
-
-*  Purpose
-*  =======
-*
-*  LSAMEN  tests if the first N letters of CA are the same as the
-*  first N letters of CB, regardless of case.
-*  LSAMEN returns .TRUE. if CA and CB are equivalent except for case
-*  and .FALSE. otherwise.  LSAMEN also returns .FALSE. if LEN( CA )
-*  or LEN( CB ) is less than N.
-*
-
-*  Arguments
-*  =========
-*
-*  N       (input) INTEGER
-*          The number of characters in CA and CB to be compared.
-*
-*  CA      (input) CHARACTER*(*)
-*  CB      (input) CHARACTER*(*)
-*          CA and CB specify two character strings of length at least N.
-*          Only the first N characters of each string will be accessed.
-*
-
-* =====================================================================
-*
-*     .. Local Scalars ..
-      INTEGER            I
-*     ..
-*     .. External Functions ..
-      LOGICAL            LSAME
-      EXTERNAL           LSAME
-*     ..
-*     .. Intrinsic Functions ..
-      INTRINSIC          LEN
-*     ..
-
-
-    </PRE>
-    <A HREF="#top">go to the page top</A>
-
-    <A NAME="ilazlr"></A>
-    <H2>ilazlr</H2>
+    <A NAME="ilaslc"></A>
+    <H2>ilaslc</H2>
     <PRE>
 USAGE:
-  __out__ = NumRu::Lapack.ilazlr( m, a, [:usage => usage, :help => help])
+  __out__ = NumRu::Lapack.ilaslc( m, a, [:usage => usage, :help => help])
 
 
 FORTRAN MANUAL
-      INTEGER FUNCTION ILAZLR( M, N, A, LDA )
+      INTEGER FUNCTION ILASLC( M, N, A, LDA )
 
 *  Purpose
 *  =======
 *
-*  ILAZLR scans A for its last non-zero row.
+*  ILASLC scans A for its last non-zero column.
 *
 
 *  Arguments
@@ -333,7 +341,7 @@ FORTRAN MANUAL
 *  N       (input) INTEGER
 *          The number of columns of the matrix A.
 *
-*  A       (input) COMPLEX*16 array, dimension (LDA,N)
+*  A       (input) REAL array, dimension (LDA,N)
 *          The m by n matrix A.
 *
 *  LDA     (input) INTEGER
@@ -347,207 +355,116 @@ FORTRAN MANUAL
     </PRE>
     <A HREF="#top">go to the page top</A>
 
-    <A NAME="iparmq"></A>
-    <H2>iparmq</H2>
+    <A NAME="ilaver"></A>
+    <H2>ilaver</H2>
     <PRE>
 USAGE:
-  __out__ = NumRu::Lapack.iparmq( ispec, name, opts, n, ilo, ihi, lwork, [:usage => usage, :help => help])
+  vers_major, vers_minor, vers_patch = NumRu::Lapack.ilaver( , [:usage => usage, :help => help])
 
 
 FORTRAN MANUAL
-      INTEGER FUNCTION IPARMQ( ISPEC, NAME, OPTS, N, ILO, IHI, LWORK )
+      SUBROUTINE ILAVER( VERS_MAJOR, VERS_MINOR, VERS_PATCH )
 
 *  Purpose
 *  =======
 *
-*       This program sets problem and machine dependent parameters
-*       useful for xHSEQR and its subroutines. It is called whenever 
-*       ILAENV is called with 12 <= ISPEC <= 16
+*  This subroutine return the Lapack version.
 *
 
 *  Arguments
 *  =========
+*  VERS_MAJOR   (output) INTEGER
+*      return the lapack major version
+*  VERS_MINOR   (output) INTEGER
+*      return the lapack minor version from the major version
+*  VERS_PATCH   (output) INTEGER
+*      return the lapack patch version from the minor version
+
+*  =====================================================================
 *
-*       ISPEC  (input) integer scalar
-*              ISPEC specifies which tunable parameter IPARMQ should
-*              return.
-*
-*              ISPEC=12: (INMIN)  Matrices of order nmin or less
-*                        are sent directly to xLAHQR, the implicit
-*                        double shift QR algorithm.  NMIN must be
-*                        at least 11.
-*
-*              ISPEC=13: (INWIN)  Size of the deflation window.
-*                        This is best set greater than or equal to
-*                        the number of simultaneous shifts NS.
-*                        Larger matrices benefit from larger deflation
-*                        windows.
-*
-*              ISPEC=14: (INIBL) Determines when to stop nibbling and
-*                        invest in an (expensive) multi-shift QR sweep.
-*                        If the aggressive early deflation subroutine
-*                        finds LD converged eigenvalues from an order
-*                        NW deflation window and LD.GT.(NW*NIBBLE)/100,
-*                        then the next QR sweep is skipped and early
-*                        deflation is applied immediately to the
-*                        remaining active diagonal block.  Setting
-*                        IPARMQ(ISPEC=14) = 0 causes TTQRE to skip a
-*                        multi-shift QR sweep whenever early deflation
-*                        finds a converged eigenvalue.  Setting
-*                        IPARMQ(ISPEC=14) greater than or equal to 100
-*                        prevents TTQRE from skipping a multi-shift
-*                        QR sweep.
-*
-*              ISPEC=15: (NSHFTS) The number of simultaneous shifts in
-*                        a multi-shift QR iteration.
-*
-*              ISPEC=16: (IACC22) IPARMQ is set to 0, 1 or 2 with the
-*                        following meanings.
-*                        0:  During the multi-shift QR sweep,
-*                            xLAQR5 does not accumulate reflections and
-*                            does not use matrix-matrix multiply to
-*                            update the far-from-diagonal matrix
-*                            entries.
-*                        1:  During the multi-shift QR sweep,
-*                            xLAQR5 and/or xLAQRaccumulates reflections and uses
-*                            matrix-matrix multiply to update the
-*                            far-from-diagonal matrix entries.
-*                        2:  During the multi-shift QR sweep.
-*                            xLAQR5 accumulates reflections and takes
-*                            advantage of 2-by-2 block structure during
-*                            matrix-matrix multiplies.
-*                        (If xTRMM is slower than xGEMM, then
-*                        IPARMQ(ISPEC=16)=1 may be more efficient than
-*                        IPARMQ(ISPEC=16)=2 despite the greater level of
-*                        arithmetic work implied by the latter choice.)
-*
-*       NAME    (input) character string
-*               Name of the calling subroutine
-*
-*       OPTS    (input) character string
-*               This is a concatenation of the string arguments to
-*               TTQRE.
-*
-*       N       (input) integer scalar
-*               N is the order of the Hessenberg matrix H.
-*
-*       ILO     (input) INTEGER
-*       IHI     (input) INTEGER
-*               It is assumed that H is already upper triangular
-*               in rows and columns 1:ILO-1 and IHI+1:N.
-*
-*       LWORK   (input) integer scalar
-*               The amount of workspace available.
-*
-
-*  Further Details
-*  ===============
-*
-*       Little is known about how best to choose these parameters.
-*       It is possible to use different values of the parameters
-*       for each of CHSEQR, DHSEQR, SHSEQR and ZHSEQR.
-*
-*       It is probably best to choose different parameters for
-*       different matrices and different parameters at different
-*       times during the iteration, but this has not been
-*       implemented --- yet.
-*
-*
-*       The best choices of most of the parameters depend
-*       in an ill-understood way on the relative execution
-*       rate of xLAQR3 and xLAQR5 and on the nature of each
-*       particular eigenvalue problem.  Experiment may be the
-*       only practical way to determine which choices are most
-*       effective.
-*
-*       Following is a list of default values supplied by IPARMQ.
-*       These defaults may be adjusted in order to attain better
-*       performance in any particular computational environment.
-*
-*       IPARMQ(ISPEC=12) The xLAHQR vs xLAQR0 crossover point.
-*                        Default: 75. (Must be at least 11.)
-*
-*       IPARMQ(ISPEC=13) Recommended deflation window size.
-*                        This depends on ILO, IHI and NS, the
-*                        number of simultaneous shifts returned
-*                        by IPARMQ(ISPEC=15).  The default for
-*                        (IHI-ILO+1).LE.500 is NS.  The default
-*                        for (IHI-ILO+1).GT.500 is 3*NS/2.
+      INTEGER VERS_MAJOR, VERS_MINOR, VERS_PATCH
+*  =====================================================================
+      VERS_MAJOR = 3
+      VERS_MINOR = 3
+      VERS_PATCH = 0
+*  =====================================================================
 *
-*       IPARMQ(ISPEC=14) Nibble crossover point.  Default: 14.
+      RETURN
+      END
+
+
+    </PRE>
+    <A HREF="#top">go to the page top</A>
+
+    <A NAME="ilaslr"></A>
+    <H2>ilaslr</H2>
+    <PRE>
+USAGE:
+  __out__ = NumRu::Lapack.ilaslr( m, a, [:usage => usage, :help => help])
+
+
+FORTRAN MANUAL
+      INTEGER FUNCTION ILASLR( M, N, A, LDA )
+
+*  Purpose
+*  =======
 *
-*       IPARMQ(ISPEC=15) Number of simultaneous shifts, NS.
-*                        a multi-shift QR iteration.
+*  ILASLR scans A for its last non-zero row.
 *
-*                        If IHI-ILO+1 is ...
+
+*  Arguments
+*  =========
 *
-*                        greater than      ...but less    ... the
-*                        or equal to ...      than        default is
+*  M       (input) INTEGER
+*          The number of rows of the matrix A.
 *
-*                                0               30       NS =   2+
-*                               30               60       NS =   4+
-*                               60              150       NS =  10
-*                              150              590       NS =  **
-*                              590             3000       NS =  64
-*                             3000             6000       NS = 128
-*                             6000             infinity   NS = 256
+*  N       (input) INTEGER
+*          The number of columns of the matrix A.
 *
-*                    (+)  By default matrices of this order are
-*                         passed to the implicit double shift routine
-*                         xLAHQR.  See IPARMQ(ISPEC=12) above.   These
-*                         values of NS are used only in case of a rare
-*                         xLAHQR failure.
+*  A       (input) REAL array, dimension (LDA,N)
+*          The m by n matrix A.
 *
-*                    (**) The asterisks (**) indicate an ad-hoc
-*                         function increasing from 10 to 64.
+*  LDA     (input) INTEGER
+*          The leading dimension of the array A. LDA >= max(1,M).
 *
-*       IPARMQ(ISPEC=16) Select structured matrix multiply.
-*                        (See ISPEC=16 above for details.)
-*                        Default: 3.
+
+*  =====================================================================
 *
-*     ================================================================
 
 
     </PRE>
     <A HREF="#top">go to the page top</A>
 
-    <A NAME="xerbla"></A>
-    <H2>xerbla</H2>
+    <A NAME="ilauplo"></A>
+    <H2>ilauplo</H2>
     <PRE>
 USAGE:
-   = NumRu::Lapack.xerbla( srname, info, [:usage => usage, :help => help])
+  __out__ = NumRu::Lapack.ilauplo( uplo, [:usage => usage, :help => help])
 
 
 FORTRAN MANUAL
-      SUBROUTINE XERBLA( SRNAME, INFO )
+      INTEGER FUNCTION ILAUPLO( UPLO )
 
 *  Purpose
 *  =======
 *
-*  XERBLA  is an error handler for the LAPACK routines.
-*  It is called by an LAPACK routine if an input parameter has an
-*  invalid value.  A message is printed and execution stops.
+*  This subroutine translated from a character string specifying a
+*  upper- or lower-triangular matrix to the relevant BLAST-specified
+*  integer constant.
 *
-*  Installers may consider modifying the STOP statement in order to
-*  call system-specific exception-handling facilities.
+*  ILAUPLO returns an INTEGER.  If ILAUPLO < 0, then the input is not
+*  a character indicating an upper- or lower-triangular matrix.
+*  Otherwise ILAUPLO returns the constant value corresponding to UPLO.
 *
 
 *  Arguments
 *  =========
-*
-*  SRNAME  (input) CHARACTER*(*)
-*          The name of the routine which called XERBLA.
-*
-*  INFO    (input) INTEGER
-*          The position of the invalid parameter in the parameter list
-*          of the calling routine.
-*
+*  UPLO    (input) CHARACTER
+*          = 'U':  A is upper triangular;
+*          = 'L':  A is lower triangular.
 
-* =====================================================================
+*  =====================================================================
 *
-*     .. Intrinsic Functions ..
-      INTRINSIC          LEN_TRIM
-*     ..
 
 
     </PRE>
@@ -620,20 +537,20 @@ FORTRAN MANUAL
     </PRE>
     <A HREF="#top">go to the page top</A>
 
-    <A NAME="ilaclc"></A>
-    <H2>ilaclc</H2>
+    <A NAME="iladlr"></A>
+    <H2>iladlr</H2>
     <PRE>
 USAGE:
-  __out__ = NumRu::Lapack.ilaclc( m, a, [:usage => usage, :help => help])
+  __out__ = NumRu::Lapack.iladlr( m, a, [:usage => usage, :help => help])
 
 
 FORTRAN MANUAL
-      INTEGER FUNCTION ILACLC( M, N, A, LDA )
+      INTEGER FUNCTION ILADLR( M, N, A, LDA )
 
 *  Purpose
 *  =======
 *
-*  ILACLC scans A for its last non-zero column.
+*  ILADLR scans A for its last non-zero row.
 *
 
 *  Arguments
@@ -645,7 +562,7 @@ FORTRAN MANUAL
 *  N       (input) INTEGER
 *          The number of columns of the matrix A.
 *
-*  A       (input) COMPLEX array, dimension (LDA,N)
+*  A       (input) DOUBLE PRECISION array, dimension (LDA,N)
 *          The m by n matrix A.
 *
 *  LDA     (input) INTEGER
@@ -659,26 +576,65 @@ FORTRAN MANUAL
     </PRE>
     <A HREF="#top">go to the page top</A>
 
-    <A NAME="izmax1"></A>
-    <H2>izmax1</H2>
+    <A NAME="ilazlc"></A>
+    <H2>ilazlc</H2>
     <PRE>
 USAGE:
-  __out__ = NumRu::Lapack.izmax1( cx, incx, [:usage => usage, :help => help])
+  __out__ = NumRu::Lapack.ilazlc( m, a, [:usage => usage, :help => help])
 
 
 FORTRAN MANUAL
-      INTEGER          FUNCTION IZMAX1( N, CX, INCX )
+      INTEGER FUNCTION ILAZLC( M, N, A, LDA )
 
 *  Purpose
 *  =======
 *
-*  IZMAX1 finds the index of the element whose real part has maximum
+*  ILAZLC scans A for its last non-zero column.
+*
+
+*  Arguments
+*  =========
+*
+*  M       (input) INTEGER
+*          The number of rows of the matrix A.
+*
+*  N       (input) INTEGER
+*          The number of columns of the matrix A.
+*
+*  A       (input) COMPLEX*16 array, dimension (LDA,N)
+*          The m by n matrix A.
+*
+*  LDA     (input) INTEGER
+*          The leading dimension of the array A. LDA >= max(1,M).
+*
+
+*  =====================================================================
+*
+
+
+    </PRE>
+    <A HREF="#top">go to the page top</A>
+
+    <A NAME="icmax1"></A>
+    <H2>icmax1</H2>
+    <PRE>
+USAGE:
+  __out__ = NumRu::Lapack.icmax1( cx, incx, [:usage => usage, :help => help])
+
+
+FORTRAN MANUAL
+      INTEGER          FUNCTION ICMAX1( N, CX, INCX )
+
+*  Purpose
+*  =======
+*
+*  ICMAX1 finds the index of the element whose real part has maximum
 *  absolute value.
 *
-*  Based on IZAMAX from Level 1 BLAS.
+*  Based on ICAMAX from Level 1 BLAS.
 *  The change is to use the 'genuine' absolute value.
 *
-*  Contributed by Nick Higham for use with ZLACON.
+*  Contributed by Nick Higham for use with CLACON.
 *
 
 *  Arguments
@@ -687,7 +643,7 @@ FORTRAN MANUAL
 *  N       (input) INTEGER
 *          The number of elements in the vector CX.
 *
-*  CX      (input) COMPLEX*16 array, dimension (N)
+*  CX      (input) COMPLEX array, dimension (N)
 *          The vector whose elements will be summed.
 *
 *  INCX    (input) INTEGER
@@ -698,55 +654,132 @@ FORTRAN MANUAL
 *
 *     .. Local Scalars ..
       INTEGER            I, IX
-      DOUBLE PRECISION   SMAX
-      COMPLEX*16         ZDUM
+      REAL               SMAX
+      COMPLEX            ZDUM
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS
 *     ..
 *     .. Statement Functions ..
-      DOUBLE PRECISION   CABS1
+      REAL               CABS1
 *     ..
 *     .. Statement Function definitions ..
 *
-*     NEXT LINE IS THE ONLY MODIFICATION.
-      CABS1( ZDUM ) = ABS( ZDUM )
-*     ..
+*     NEXT LINE IS THE ONLY MODIFICATION.
+      CABS1( ZDUM ) = ABS( ZDUM )
+*     ..
+
+
+    </PRE>
+    <A HREF="#top">go to the page top</A>
+
+    <A NAME="ilaclr"></A>
+    <H2>ilaclr</H2>
+    <PRE>
+USAGE:
+  __out__ = NumRu::Lapack.ilaclr( m, a, [:usage => usage, :help => help])
+
+
+FORTRAN MANUAL
+      INTEGER FUNCTION ILACLR( M, N, A, LDA )
+
+*  Purpose
+*  =======
+*
+*  ILACLR scans A for its last non-zero row.
+*
+
+*  Arguments
+*  =========
+*
+*  M       (input) INTEGER
+*          The number of rows of the matrix A.
+*
+*  N       (input) INTEGER
+*          The number of columns of the matrix A.
+*
+*  A       (input) COMPLEX          array, dimension (LDA,N)
+*          The m by n matrix A.
+*
+*  LDA     (input) INTEGER
+*          The leading dimension of the array A. LDA >= max(1,M).
+*
+
+*  =====================================================================
+*
+
+
+    </PRE>
+    <A HREF="#top">go to the page top</A>
+
+    <A NAME="iladlc"></A>
+    <H2>iladlc</H2>
+    <PRE>
+USAGE:
+  __out__ = NumRu::Lapack.iladlc( m, a, [:usage => usage, :help => help])
+
+
+FORTRAN MANUAL
+      INTEGER FUNCTION ILADLC( M, N, A, LDA )
+
+*  Purpose
+*  =======
+*
+*  ILADLC scans A for its last non-zero column.
+*
+
+*  Arguments
+*  =========
+*
+*  M       (input) INTEGER
+*          The number of rows of the matrix A.
+*
+*  N       (input) INTEGER
+*          The number of columns of the matrix A.
+*
+*  A       (input) DOUBLE PRECISION array, dimension (LDA,N)
+*          The m by n matrix A.
+*
+*  LDA     (input) INTEGER
+*          The leading dimension of the array A. LDA >= max(1,M).
+*
+
+*  =====================================================================
+*
 
 
     </PRE>
     <A HREF="#top">go to the page top</A>
 
-    <A NAME="ilaprec"></A>
-    <H2>ilaprec</H2>
+    <A NAME="ilatrans"></A>
+    <H2>ilatrans</H2>
     <PRE>
 USAGE:
-  __out__ = NumRu::Lapack.ilaprec( prec, [:usage => usage, :help => help])
+  __out__ = NumRu::Lapack.ilatrans( trans, [:usage => usage, :help => help])
 
 
 FORTRAN MANUAL
-      INTEGER FUNCTION ILAPREC( PREC )
+      INTEGER FUNCTION ILATRANS( TRANS )
 
 *  Purpose
 *  =======
 *
-*  This subroutine translated from a character string specifying an
-*  intermediate precision to the relevant BLAST-specified integer
+*  This subroutine translates from a character string specifying a
+*  transposition operation to the relevant BLAST-specified integer
 *  constant.
 *
-*  ILAPREC returns an INTEGER.  If ILAPREC < 0, then the input is not a
-*  character indicating a supported intermediate precision.  Otherwise
-*  ILAPREC returns the constant value corresponding to PREC.
+*  ILATRANS returns an INTEGER.  If ILATRANS < 0, then the input is not
+*  a character indicating a transposition operator.  Otherwise ILATRANS
+*  returns the constant value corresponding to TRANS.
 *
 
 *  Arguments
 *  =========
-*  PREC    (input) CHARACTER
+*  TRANS   (input) CHARACTER*1
 *          Specifies the form of the system of equations:
-*          = 'S':  Single
-*          = 'D':  Double
-*          = 'I':  Indigenous
-*          = 'X', 'E':  Extra
+*          = 'N':  No transpose
+*          = 'T':  Transpose
+*          = 'C':  Conjugate transpose
 
 *  =====================================================================
 *
@@ -877,149 +910,156 @@ FORTRAN MANUAL
     </PRE>
     <A HREF="#top">go to the page top</A>
 
-    <A NAME="ieeeck"></A>
-    <H2>ieeeck</H2>
+    <A NAME="ilaclc"></A>
+    <H2>ilaclc</H2>
     <PRE>
 USAGE:
-  __out__ = NumRu::Lapack.ieeeck( ispec, zero, one, [:usage => usage, :help => help])
+  __out__ = NumRu::Lapack.ilaclc( m, a, [:usage => usage, :help => help])
 
 
 FORTRAN MANUAL
-      INTEGER          FUNCTION IEEECK( ISPEC, ZERO, ONE )
+      INTEGER FUNCTION ILACLC( M, N, A, LDA )
 
 *  Purpose
 *  =======
 *
-*  IEEECK is called from the ILAENV to verify that Infinity and
-*  possibly NaN arithmetic is safe (i.e. will not trap).
+*  ILACLC scans A for its last non-zero column.
 *
 
 *  Arguments
 *  =========
 *
-*  ISPEC   (input) INTEGER
-*          Specifies whether to test just for inifinity arithmetic
-*          or whether to test for infinity and NaN arithmetic.
-*          = 0: Verify infinity arithmetic only.
-*          = 1: Verify infinity and NaN arithmetic.
+*  M       (input) INTEGER
+*          The number of rows of the matrix A.
 *
-*  ZERO    (input) REAL
-*          Must contain the value 0.0
-*          This is passed to prevent the compiler from optimizing
-*          away this code.
+*  N       (input) INTEGER
+*          The number of columns of the matrix A.
 *
-*  ONE     (input) REAL
-*          Must contain the value 1.0
-*          This is passed to prevent the compiler from optimizing
-*          away this code.
+*  A       (input) COMPLEX array, dimension (LDA,N)
+*          The m by n matrix A.
 *
-*  RETURN VALUE:  INTEGER
-*          = 0:  Arithmetic failed to produce the correct answers
-*          = 1:  Arithmetic produced the correct answers
+*  LDA     (input) INTEGER
+*          The leading dimension of the array A. LDA >= max(1,M).
+*
+
+*  =====================================================================
 *
-*     .. Local Scalars ..
-      REAL               NAN1, NAN2, NAN3, NAN4, NAN5, NAN6, NEGINF,
-     $                   NEGZRO, NEWZRO, POSINF
-*     ..
 
 
     </PRE>
     <A HREF="#top">go to the page top</A>
 
-    <A NAME="iladlc"></A>
-    <H2>iladlc</H2>
+    <A NAME="xerbla"></A>
+    <H2>xerbla</H2>
     <PRE>
 USAGE:
-  __out__ = NumRu::Lapack.iladlc( m, a, [:usage => usage, :help => help])
+   = NumRu::Lapack.xerbla( srname, info, [:usage => usage, :help => help])
 
 
 FORTRAN MANUAL
-      INTEGER FUNCTION ILADLC( M, N, A, LDA )
+      SUBROUTINE XERBLA( SRNAME, INFO )
 
 *  Purpose
 *  =======
 *
-*  ILADLC scans A for its last non-zero column.
+*  XERBLA  is an error handler for the LAPACK routines.
+*  It is called by an LAPACK routine if an input parameter has an
+*  invalid value.  A message is printed and execution stops.
+*
+*  Installers may consider modifying the STOP statement in order to
+*  call system-specific exception-handling facilities.
 *
 
 *  Arguments
 *  =========
 *
-*  M       (input) INTEGER
-*          The number of rows of the matrix A.
-*
-*  N       (input) INTEGER
-*          The number of columns of the matrix A.
-*
-*  A       (input) DOUBLE PRECISION array, dimension (LDA,N)
-*          The m by n matrix A.
+*  SRNAME  (input) CHARACTER*(*)
+*          The name of the routine which called XERBLA.
 *
-*  LDA     (input) INTEGER
-*          The leading dimension of the array A. LDA >= max(1,M).
+*  INFO    (input) INTEGER
+*          The position of the invalid parameter in the parameter list
+*          of the calling routine.
 *
 
-*  =====================================================================
+* =====================================================================
 *
+*     .. Intrinsic Functions ..
+      INTRINSIC          LEN_TRIM
+*     ..
 
 
     </PRE>
     <A HREF="#top">go to the page top</A>
 
-    <A NAME="ilauplo"></A>
-    <H2>ilauplo</H2>
+    <A NAME="ieeeck"></A>
+    <H2>ieeeck</H2>
     <PRE>
 USAGE:
-  __out__ = NumRu::Lapack.ilauplo( uplo, [:usage => usage, :help => help])
+  __out__ = NumRu::Lapack.ieeeck( ispec, zero, one, [:usage => usage, :help => help])
 
 
 FORTRAN MANUAL
-      INTEGER FUNCTION ILAUPLO( UPLO )
+      INTEGER          FUNCTION IEEECK( ISPEC, ZERO, ONE )
 
 *  Purpose
 *  =======
 *
-*  This subroutine translated from a character string specifying a
-*  upper- or lower-triangular matrix to the relevant BLAST-specified
-*  integer constant.
-*
-*  ILAUPLO returns an INTEGER.  If ILAUPLO < 0, then the input is not
-*  a character indicating an upper- or lower-triangular matrix.
-*  Otherwise ILAUPLO returns the constant value corresponding to UPLO.
+*  IEEECK is called from the ILAENV to verify that Infinity and
+*  possibly NaN arithmetic is safe (i.e. will not trap).
 *
 
 *  Arguments
 *  =========
-*  UPLO    (input) CHARACTER
-*          = 'U':  A is upper triangular;
-*          = 'L':  A is lower triangular.
-
-*  =====================================================================
 *
+*  ISPEC   (input) INTEGER
+*          Specifies whether to test just for inifinity arithmetic
+*          or whether to test for infinity and NaN arithmetic.
+*          = 0: Verify infinity arithmetic only.
+*          = 1: Verify infinity and NaN arithmetic.
+*
+*  ZERO    (input) REAL
+*          Must contain the value 0.0
+*          This is passed to prevent the compiler from optimizing
+*          away this code.
+*
+*  ONE     (input) REAL
+*          Must contain the value 1.0
+*          This is passed to prevent the compiler from optimizing
+*          away this code.
+*
+*  RETURN VALUE:  INTEGER
+*          = 0:  Arithmetic failed to produce the correct answers
+*          = 1:  Arithmetic produced the correct answers
+*
+*     .. Local Scalars ..
+      REAL               NAN1, NAN2, NAN3, NAN4, NAN5, NAN6, NEGINF,
+     $                   NEGZRO, NEWZRO, POSINF
+*     ..
 
 
     </PRE>
     <A HREF="#top">go to the page top</A>
 
-    <A NAME="icmax1"></A>
-    <H2>icmax1</H2>
+    <A NAME="izmax1"></A>
+    <H2>izmax1</H2>
     <PRE>
 USAGE:
-  __out__ = NumRu::Lapack.icmax1( cx, incx, [:usage => usage, :help => help])
+  __out__ = NumRu::Lapack.izmax1( cx, incx, [:usage => usage, :help => help])
 
 
 FORTRAN MANUAL
-      INTEGER          FUNCTION ICMAX1( N, CX, INCX )
+      INTEGER          FUNCTION IZMAX1( N, CX, INCX )
 
 *  Purpose
 *  =======
 *
-*  ICMAX1 finds the index of the element whose real part has maximum
+*  IZMAX1 finds the index of the element whose real part has maximum
 *  absolute value.
 *
-*  Based on ICAMAX from Level 1 BLAS.
+*  Based on IZAMAX from Level 1 BLAS.
 *  The change is to use the 'genuine' absolute value.
 *
-*  Contributed by Nick Higham for use with CLACON.
+*  Contributed by Nick Higham for use with ZLACON.
 *
 
 *  Arguments
@@ -1028,7 +1068,7 @@ FORTRAN MANUAL
 *  N       (input) INTEGER
 *          The number of elements in the vector CX.
 *
-*  CX      (input) COMPLEX array, dimension (N)
+*  CX      (input) COMPLEX*16 array, dimension (N)
 *          The vector whose elements will be summed.
 *
 *  INCX    (input) INTEGER
@@ -1039,14 +1079,14 @@ FORTRAN MANUAL
 *
 *     .. Local Scalars ..
       INTEGER            I, IX
-      REAL               SMAX
-      COMPLEX            ZDUM
+      DOUBLE PRECISION   SMAX
+      COMPLEX*16         ZDUM
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS
 *     ..
 *     .. Statement Functions ..
-      REAL               CABS1
+      DOUBLE PRECISION   CABS1
 *     ..
 *     .. Statement Function definitions ..
 *
@@ -1058,76 +1098,36 @@ FORTRAN MANUAL
     </PRE>
     <A HREF="#top">go to the page top</A>
 
-    <A NAME="ilaslr"></A>
-    <H2>ilaslr</H2>
+    <A NAME="ilaprec"></A>
+    <H2>ilaprec</H2>
     <PRE>
 USAGE:
-  __out__ = NumRu::Lapack.ilaslr( m, a, [:usage => usage, :help => help])
+  __out__ = NumRu::Lapack.ilaprec( prec, [:usage => usage, :help => help])
 
 
 FORTRAN MANUAL
-      INTEGER FUNCTION ILASLR( M, N, A, LDA )
+      INTEGER FUNCTION ILAPREC( PREC )
 
 *  Purpose
 *  =======
 *
-*  ILASLR scans A for its last non-zero row.
-*
-
-*  Arguments
-*  =========
-*
-*  M       (input) INTEGER
-*          The number of rows of the matrix A.
-*
-*  N       (input) INTEGER
-*          The number of columns of the matrix A.
-*
-*  A       (input) REAL array, dimension (LDA,N)
-*          The m by n matrix A.
-*
-*  LDA     (input) INTEGER
-*          The leading dimension of the array A. LDA >= max(1,M).
-*
-
-*  =====================================================================
-*
-
-
-    </PRE>
-    <A HREF="#top">go to the page top</A>
-
-    <A NAME="ilaclr"></A>
-    <H2>ilaclr</H2>
-    <PRE>
-USAGE:
-  __out__ = NumRu::Lapack.ilaclr( m, a, [:usage => usage, :help => help])
-
-
-FORTRAN MANUAL
-      INTEGER FUNCTION ILACLR( M, N, A, LDA )
-
-*  Purpose
-*  =======
+*  This subroutine translated from a character string specifying an
+*  intermediate precision to the relevant BLAST-specified integer
+*  constant.
 *
-*  ILACLR scans A for its last non-zero row.
+*  ILAPREC returns an INTEGER.  If ILAPREC < 0, then the input is not a
+*  character indicating a supported intermediate precision.  Otherwise
+*  ILAPREC returns the constant value corresponding to PREC.
 *
 
 *  Arguments
 *  =========
-*
-*  M       (input) INTEGER
-*          The number of rows of the matrix A.
-*
-*  N       (input) INTEGER
-*          The number of columns of the matrix A.
-*
-*  A       (input) COMPLEX          array, dimension (LDA,N)
-*          The m by n matrix A.
-*
-*  LDA     (input) INTEGER
-*          The leading dimension of the array A. LDA >= max(1,M).
-*
+*  PREC    (input) CHARACTER
+*          Specifies the form of the system of equations:
+*          = 'S':  Single
+*          = 'D':  Double
+*          = 'I':  Indigenous
+*          = 'X', 'E':  Extra
 
 *  =====================================================================
 *