linalg 1.0.2

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 90308570ba2bd52acca9489497ecf195ac1dccda
+  data.tar.gz: 98273ca4701588085a03034fc04fd286289770c0
+SHA512:
+  metadata.gz: 86139429d64a55b96f1d3c509797bc1b189e3e75e4c30629dce989b7284d5f706bbea878fcbe99bd90c96df88f9a7e2db825c52167ae145b9db699b78f544298
+  data.tar.gz: 78239441788657d0ea82dcc46ac891ac5fbbb5e78c59c1ab9f8c311fb8d4cdce182dbd9213b966003ac5282f5f7c927c6d62cfd017d0c6de9b1dd45e76e8e7c6

data/LICENSE

@@ -0,0 +1,24 @@
+
+"linalg -- Ruby Linear Algebra Package" is copyrighted free software
+by James M. Lawrence.  Distributed under the MIT license.
+
+Copyright (c) 2004-2008 James M. Lawrence
+
+Permission is hereby granted, free of charge, to any person obtaining
+a copy of this software and associated documentation files (the
+"Software"), to deal in the Software without restriction, including
+without limitation the rights to use, copy, modify, merge, publish,
+distribute, sublicense, and/or sell copies of the Software, and to
+permit persons to whom the Software is furnished to do so, subject to
+the following conditions:
+
+The above copyright notice and this permission notice shall be
+included in all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
+NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
+LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
+WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

data/ext/g2c_typedefs.h

@@ -0,0 +1,5 @@
+typedef int __g77_integer;
+typedef unsigned int __g77_uinteger;
+typedef long int __g77_longint;
+typedef unsigned long int __g77_ulongint;
+

data/ext/lapack/extconf.rb

@@ -0,0 +1,309 @@
+#
+# Copyright (c) 2004-2008 by James M. Lawrence
+#
+# See LICENSE
+#
+
+require 'mkmf'
+require 'enumerator'
+
+module Enumerable
+   def map_with_index
+      arr = []
+      each_with_index { |e, i|
+         arr << yield(e, i)
+      }
+      arr
+   end
+end
+
+class Integer
+   def enum
+      enum_for(:times)
+   end
+end
+
+module GenerateHeader
+   def self.write(dir, outname, version = nil)
+      outname = File.expand_path(outname)
+      Dir.chdir(dir) {
+         Dir.glob('*.P') { |f| File.unlink(f) }
+         if system("f2c -R -P -\!c *.f")
+            File.open(outname, "w") { |outfile|
+               macro = File.basename(outname.upcase.gsub(".","_"))
+               outfile.puts "#ifndef #{macro}"
+               outfile.puts "#define #{macro}"
+               outfile.puts "/* version #{version} */ " if version
+               outfile.puts "/* generated #{Time.now} */"
+               outfile.puts
+               Dir.glob('*.P') { |pname|
+                  File.open(pname) { |pfile|
+                     pfile.each { |line|
+                        next unless line =~ %r!^extern!
+
+                        # I don't understand why f2c inserts these ftnlen params
+                        while line.gsub!(%r!, ftnlen [a-z]+_len\)!, ")") ; end
+
+                        outfile.write line
+                     }
+                  }
+               }
+
+               outfile.puts
+               outfile.puts "#endif"
+            }
+         end
+      }
+   end
+end
+
+module ScanHeader
+   extend self
+
+   def ctype_canonical_form(ctype)
+      ctype.gsub(%r!\*!, " * ").gsub(%r!\s+!, " ").strip
+   end
+
+   def parse_prototype(s, argnames = true)
+      func, args = s.strip.gsub(%r!\s+!, " ").split(%r![\(\)]!)
+      ret, name =
+         func.
+         strip.
+         match(%r!^(\w[*\s\w]*[*\s])(\w+)$!).
+         captures
+      desc = {}
+      desc[:name] = name
+      desc[:return] = ctype_canonical_form(ret)
+      desc[:args] = []
+      unless args == "" or args.strip == "void"
+         args.split(",").each { |arg|
+            arg.sub!(%r!\w+\s*$!, "") if argnames
+            desc[:args].push(ctype_canonical_form(arg))
+         }
+      end
+      desc
+   end
+
+   def scan_header_io(io, ppdefs = nil)
+      s = io.read
+      s.gsub!(%r!^\s*\#.*$!, "")
+      s.gsub!(%r!//.*$!, "")
+      s.gsub!(%r!/\*.*?\*/!m, "")
+      if ppdefs
+         s.gsub!(%r!\w+!) { |word|
+            if ppdefs[word]
+               ppdefs[word]
+            else
+               word
+            end
+         }
+      end
+      s.scan(%r!\w[*\w\s]*?\w+\s*\(.*?\)\s*\;!m) { |func|
+         yield parse_prototype(func)
+      }
+   end
+
+   def scan_header(header, ppdefs = nil, &block)
+      File.open(header) { |file|
+         scan_header_io(file, ppdefs, &block)
+      }
+   end
+end
+
+module FunctionDB
+   def self.read
+      db = {}
+
+      cppdefs = {
+         "extern" => "",
+
+         # these are in f2c.h
+         "C_f" => "void",
+         "Z_f" => "void",
+      }
+
+      ["include/LAPACK.h", "include/BLAS.h"].each { |filename|
+         print "reading #{filename}..."
+         ScanHeader.scan_header(filename, cppdefs) { |desc|
+            db[desc[:name]] = desc
+         }
+         puts "ok"
+      }
+
+      db
+   end
+end
+
+module Lapack_c
+   extend self
+
+   def write(db)
+      {
+
+         :s => /^s/,
+         :d => /^d/,
+         :c => /^c/,
+         :z => /^z/,
+         :x => /^[^sdcz]/,
+
+      }.each_pair { |letter, regex|
+         part = {}
+         db.keys.grep(regex).each { |key|
+            part[key] = db.delete(key)
+         }
+
+         write_file(part,
+                    "rb_lapack_#{letter}",
+                    "define_lapack_#{letter}")
+      }
+   end
+
+   def write_file(db, filename, funcname)
+      File.open("#{filename}.c", "w+") { |file|
+         file.puts %{#include "rb_lapack.h"\n\n}
+         db.keys.sort.each { |func|
+            Lapack_c.write_function(db[func], file)
+         }
+
+         file.puts "void #{funcname}()\n{\n"
+
+         db.keys.sort.each { |func|
+         file.puts <<EOS
+    rb_define_singleton_method(rb_cLapack,
+                               "#{func.sub(%r!_$!, "")}",
+                               #{Lapack_c.rbname(func)},
+                               #{-1}) ;
+EOS
+         }
+
+         file.puts "\n}"
+      }
+   end
+
+   def rbname(name)
+      "rb_" + name.sub(%r!_$!, "")
+   end
+
+   def write_function(desc, out)
+     rbname = rbname(desc[:name])
+     args = ['int argc', 'VALUE *argv', 'VALUE self']
+     conv = desc[:args].map { |type|
+       "to_" + case type
+               when %r!^(\w+)$!
+                 $1
+               when %r!^(\w+)\s*\*$!
+                 "#{$1}_ptr"
+               else
+                 raise "unknown arg type in #{desc[:name]}: '#{type}'"
+               end.gsub(" ","_")
+     }
+     assign = desc[:args].map_with_index { |type, i|
+       "#{type} p#{i} = (#{type})NUM2ULONG(rb_funcall(v#{i}, #{conv[i]}, 0))"
+      }
+      ftncall = desc[:name] + "(" + desc[:args].size.enum.inject("") { |acc, i|
+         acc + "p#{i}, "
+      }.chop.chop + ")"
+      ret = if desc[:return] =~ %r!^void$!i
+               "#{ftncall} ; return Qnil ;"
+            else
+               "return " +
+                  case desc[:return]
+                  when %r!(^real$|^doublereal$)!
+                     "rb_float_new((double)"
+                  when %r!(^int|^logical$)!
+                     "LONG2NUM((long)"
+                  else
+                     raise "unknown type return type in #{desc[:name]}: #{desc[:return]}"
+                  end +
+                  "#{ftncall})"
+            end
+     define_values = []
+     (0..desc[:args].size).map do |j|
+       define_values << "VALUE v#{j};"
+     end
+
+      out.puts <<EOS
+VALUE #{rbname}(#{args.join(",")})
+{
+    #{define_values.join("\n    ") }
+    rb_scan_args(argc, argv, "#{desc[:args].size}*", #{(0..desc[:args].size).map{ |j| "&v#{j}" }.join(', ') });
+    #{assign.join(" ;\n    ")} ;
+
+    #{ret} ;
+}
+
+EOS
+   end
+end
+
+module PullCode
+   def self.pull(sofile, *reps)
+      `nm -D #{sofile}`.each { |line|
+         line =~ %r!U (\w+)!
+         next unless $1
+         func = $1
+         reps.each { |d|
+            d = d.sub %r!/$!, ""
+            [
+               "#{d}/#{func}.c",
+               d + "/" + func.sub(%r!_$!, "") + ".c",
+            ].each { |fname|
+               if FileTest.exist? fname
+                  FileUtils.cp fname, "."
+                  puts "-->copied #{fname}"
+                  break
+               end
+            }
+         }
+      }
+   end
+end
+
+module Main
+   extend self
+
+   def genheader(*args)
+      GenerateHeader.write(*args)
+   end
+
+   def pull(*args)
+      PullCode.pull(*args)
+   end
+
+   def write_c
+      db = FunctionDB.read
+      Lapack_c.write(db)
+   end
+
+   def config
+      unless have_header("f2c.h") and
+            have_library("blas") and
+            have_library("lapack")
+         puts "A full LAPACK installation was not found."
+         exit(-1)
+      end
+
+      $distcleanfiles = [
+         "rb_lapack_s.c",
+         "rb_lapack_d.c",
+         "rb_lapack_c.c",
+         "rb_lapack_z.c",
+         "rb_lapack_x.c",
+      ]
+   end
+
+   def create
+      create_makefile("lapack")
+   end
+
+   STDOUT.sync = true
+   action = ARGV.shift
+
+   if action
+      send(action, *ARGV)
+   else
+      config
+      write_c
+      create
+   end
+end

data/ext/lapack/include/BLAS.h

@@ -0,0 +1,154 @@
+/*
+ * Copyright (c) 2004-2008 by James M. Lawrence
+ *
+ * See LICENSE
+ *
+ */
+#ifndef BLAS_H
+#define BLAS_H
+/* version 20031025 */ 
+/* generated Fri Jun 04 05:58:15 EDT 2004-2008 */
+
+extern int caxpy_(integer *n, complex *ca, complex *cx, integer *incx, complex *cy, integer *incy);
+extern int ccopy_(integer *n, complex *cx, integer *incx, complex *cy, integer *incy);
+extern C_f cdotc_(complex * ret_val, integer *n, complex *cx, integer *incx, complex *cy, integer *incy);
+extern C_f cdotu_(complex * ret_val, integer *n, complex *cx, integer *incx, complex *cy, integer *incy);
+extern int cgbmv_(char *trans, integer *m, integer *n, integer *kl, integer *ku, complex *alpha, complex *a, integer *lda, complex *x, integer *incx, complex *beta, complex *y, integer *incy);
+extern int cgemm_(char *transa, char *transb, integer *m, integer *n, integer *k, complex *alpha, complex *a, integer *lda, complex *b, integer *ldb, complex *beta, complex *c__, integer *ldc);
+extern int cgemv_(char *trans, integer *m, integer *n, complex *alpha, complex *a, integer *lda, complex *x, integer *incx, complex *beta, complex *y, integer *incy);
+extern int cgerc_(integer *m, integer *n, complex *alpha, complex *x, integer *incx, complex *y, integer *incy, complex *a, integer *lda);
+extern int cgeru_(integer *m, integer *n, complex *alpha, complex *x, integer *incx, complex *y, integer *incy, complex *a, integer *lda);
+extern int chbmv_(char *uplo, integer *n, integer *k, complex *alpha, complex *a, integer *lda, complex *x, integer *incx, complex *beta, complex *y, integer *incy);
+extern int chemm_(char *side, char *uplo, integer *m, integer *n, complex *alpha, complex *a, integer *lda, complex *b, integer *ldb, complex *beta, complex *c__, integer *ldc);
+extern int chemv_(char *uplo, integer *n, complex *alpha, complex *a, integer *lda, complex *x, integer *incx, complex *beta, complex *y, integer *incy);
+extern int cher_(char *uplo, integer *n, real *alpha, complex *x, integer *incx, complex *a, integer *lda);
+extern int cher2_(char *uplo, integer *n, complex *alpha, complex *x, integer *incx, complex *y, integer *incy, complex *a, integer *lda);
+extern int cher2k_(char *uplo, char *trans, integer *n, integer *k, complex *alpha, complex *a, integer *lda, complex *b, integer *ldb, real *beta, complex *c__, integer *ldc);
+extern int cherk_(char *uplo, char *trans, integer *n, integer *k, real *alpha, complex *a, integer *lda, real *beta, complex *c__, integer *ldc);
+extern int chpmv_(char *uplo, integer *n, complex *alpha, complex *ap, complex *x, integer *incx, complex *beta, complex *y, integer *incy);
+extern int chpr_(char *uplo, integer *n, real *alpha, complex *x, integer *incx, complex *ap);
+extern int chpr2_(char *uplo, integer *n, complex *alpha, complex *x, integer *incx, complex *y, integer *incy, complex *ap);
+extern int crotg_(complex *ca, complex *cb, real *c__, complex *s);
+extern int cscal_(integer *n, complex *ca, complex *cx, integer *incx);
+extern int csscal_(integer *n, real *sa, complex *cx, integer *incx);
+extern int cswap_(integer *n, complex *cx, integer *incx, complex *cy, integer *incy);
+extern int csymm_(char *side, char *uplo, integer *m, integer *n, complex *alpha, complex *a, integer *lda, complex *b, integer *ldb, complex *beta, complex *c__, integer *ldc);
+extern int csyr2k_(char *uplo, char *trans, integer *n, integer *k, complex *alpha, complex *a, integer *lda, complex *b, integer *ldb, complex *beta, complex *c__, integer *ldc);
+extern int csyrk_(char *uplo, char *trans, integer *n, integer *k, complex *alpha, complex *a, integer *lda, complex *beta, complex *c__, integer *ldc);
+extern int ctbmv_(char *uplo, char *trans, char *diag, integer *n, integer *k, complex *a, integer *lda, complex *x, integer *incx);
+extern int ctbsv_(char *uplo, char *trans, char *diag, integer *n, integer *k, complex *a, integer *lda, complex *x, integer *incx);
+extern int ctpmv_(char *uplo, char *trans, char *diag, integer *n, complex *ap, complex *x, integer *incx);
+extern int ctpsv_(char *uplo, char *trans, char *diag, integer *n, complex *ap, complex *x, integer *incx);
+extern int ctrmm_(char *side, char *uplo, char *transa, char *diag, integer *m, integer *n, complex *alpha, complex *a, integer *lda, complex *b, integer *ldb);
+extern int ctrmv_(char *uplo, char *trans, char *diag, integer *n, complex *a, integer *lda, complex *x, integer *incx);
+extern int ctrsm_(char *side, char *uplo, char *transa, char *diag, integer *m, integer *n, complex *alpha, complex *a, integer *lda, complex *b, integer *ldb);
+extern int ctrsv_(char *uplo, char *trans, char *diag, integer *n, complex *a, integer *lda, complex *x, integer *incx);
+extern doublereal dasum_(integer *n, doublereal *dx, integer *incx);
+extern int daxpy_(integer *n, doublereal *da, doublereal *dx, integer *incx, doublereal *dy, integer *incy);
+extern doublereal dcabs1_(doublecomplex *z__);
+extern int dcopy_(integer *n, doublereal *dx, integer *incx, doublereal *dy, integer *incy);
+extern doublereal ddot_(integer *n, doublereal *dx, integer *incx, doublereal *dy, integer *incy);
+extern int dgbmv_(char *trans, integer *m, integer *n, integer *kl, integer *ku, doublereal *alpha, doublereal *a, integer *lda, doublereal *x, integer *incx, doublereal *beta, doublereal *y, integer *incy);
+extern int dgemm_(char *transa, char *transb, integer *m, integer *n, integer *k, doublereal *alpha, doublereal *a, integer *lda, doublereal *b, integer *ldb, doublereal *beta, doublereal *c__, integer *ldc);
+extern int dgemv_(char *trans, integer *m, integer *n, doublereal *alpha, doublereal *a, integer *lda, doublereal *x, integer *incx, doublereal *beta, doublereal *y, integer *incy);
+extern int dger_(integer *m, integer *n, doublereal *alpha, doublereal *x, integer *incx, doublereal *y, integer *incy, doublereal *a, integer *lda);
+extern doublereal dnrm2_(integer *n, doublereal *x, integer *incx);
+extern int drot_(integer *n, doublereal *dx, integer *incx, doublereal *dy, integer *incy, doublereal *c__, doublereal *s);
+extern int drotg_(doublereal *da, doublereal *db, doublereal *c__, doublereal *s);
+extern int dsbmv_(char *uplo, integer *n, integer *k, doublereal *alpha, doublereal *a, integer *lda, doublereal *x, integer *incx, doublereal *beta, doublereal *y, integer *incy);
+extern int dscal_(integer *n, doublereal *da, doublereal *dx, integer *incx);
+extern int dspmv_(char *uplo, integer *n, doublereal *alpha, doublereal *ap, doublereal *x, integer *incx, doublereal *beta, doublereal *y, integer *incy);
+extern int dspr_(char *uplo, integer *n, doublereal *alpha, doublereal *x, integer *incx, doublereal *ap);
+extern int dspr2_(char *uplo, integer *n, doublereal *alpha, doublereal *x, integer *incx, doublereal *y, integer *incy, doublereal *ap);
+extern int dswap_(integer *n, doublereal *dx, integer *incx, doublereal *dy, integer *incy);
+extern int dsymm_(char *side, char *uplo, integer *m, integer *n, doublereal *alpha, doublereal *a, integer *lda, doublereal *b, integer *ldb, doublereal *beta, doublereal *c__, integer *ldc);
+extern int dsymv_(char *uplo, integer *n, doublereal *alpha, doublereal *a, integer *lda, doublereal *x, integer *incx, doublereal *beta, doublereal *y, integer *incy);
+extern int dsyr_(char *uplo, integer *n, doublereal *alpha, doublereal *x, integer *incx, doublereal *a, integer *lda);
+extern int dsyr2_(char *uplo, integer *n, doublereal *alpha, doublereal *x, integer *incx, doublereal *y, integer *incy, doublereal *a, integer *lda);
+extern int dsyr2k_(char *uplo, char *trans, integer *n, integer *k, doublereal *alpha, doublereal *a, integer *lda, doublereal *b, integer *ldb, doublereal *beta, doublereal *c__, integer *ldc);
+extern int dsyrk_(char *uplo, char *trans, integer *n, integer *k, doublereal *alpha, doublereal *a, integer *lda, doublereal *beta, doublereal *c__, integer *ldc);
+extern int dtbmv_(char *uplo, char *trans, char *diag, integer *n, integer *k, doublereal *a, integer *lda, doublereal *x, integer *incx);
+extern int dtbsv_(char *uplo, char *trans, char *diag, integer *n, integer *k, doublereal *a, integer *lda, doublereal *x, integer *incx);
+extern int dtpmv_(char *uplo, char *trans, char *diag, integer *n, doublereal *ap, doublereal *x, integer *incx);
+extern int dtpsv_(char *uplo, char *trans, char *diag, integer *n, doublereal *ap, doublereal *x, integer *incx);
+extern int dtrmm_(char *side, char *uplo, char *transa, char *diag, integer *m, integer *n, doublereal *alpha, doublereal *a, integer *lda, doublereal *b, integer *ldb);
+extern int dtrmv_(char *uplo, char *trans, char *diag, integer *n, doublereal *a, integer *lda, doublereal *x, integer *incx);
+extern int dtrsm_(char *side, char *uplo, char *transa, char *diag, integer *m, integer *n, doublereal *alpha, doublereal *a, integer *lda, doublereal *b, integer *ldb);
+extern int dtrsv_(char *uplo, char *trans, char *diag, integer *n, doublereal *a, integer *lda, doublereal *x, integer *incx);
+extern doublereal dzasum_(integer *n, doublecomplex *zx, integer *incx);
+extern doublereal dznrm2_(integer *n, doublecomplex *x, integer *incx);
+extern integer icamax_(integer *n, complex *cx, integer *incx);
+extern integer idamax_(integer *n, doublereal *dx, integer *incx);
+extern integer isamax_(integer *n, real *sx, integer *incx);
+extern integer izamax_(integer *n, doublecomplex *zx, integer *incx);
+extern logical lsame_(char *ca, char *cb);
+extern real sasum_(integer *n, real *sx, integer *incx);
+extern int saxpy_(integer *n, real *sa, real *sx, integer *incx, real *sy, integer *incy);
+extern real scasum_(integer *n, complex *cx, integer *incx);
+extern real scnrm2_(integer *n, complex *x, integer *incx);
+extern int scopy_(integer *n, real *sx, integer *incx, real *sy, integer *incy);
+extern real sdot_(integer *n, real *sx, integer *incx, real *sy, integer *incy);
+extern int sgbmv_(char *trans, integer *m, integer *n, integer *kl, integer *ku, real *alpha, real *a, integer *lda, real *x, integer *incx, real *beta, real *y, integer *incy);
+extern int sgemm_(char *transa, char *transb, integer *m, integer *n, integer *k, real *alpha, real *a, integer *lda, real *b, integer *ldb, real *beta, real *c__, integer *ldc);
+extern int sgemv_(char *trans, integer *m, integer *n, real *alpha, real *a, integer *lda, real *x, integer *incx, real *beta, real *y, integer *incy);
+extern int sger_(integer *m, integer *n, real *alpha, real *x, integer *incx, real *y, integer *incy, real *a, integer *lda);
+extern real snrm2_(integer *n, real *x, integer *incx);
+extern int srot_(integer *n, real *sx, integer *incx, real *sy, integer *incy, real *c__, real *s);
+extern int srotg_(real *sa, real *sb, real *c__, real *s);
+extern int ssbmv_(char *uplo, integer *n, integer *k, real *alpha, real *a, integer *lda, real *x, integer *incx, real *beta, real *y, integer *incy);
+extern int sscal_(integer *n, real *sa, real *sx, integer *incx);
+extern int sspmv_(char *uplo, integer *n, real *alpha, real *ap, real *x, integer *incx, real *beta, real *y, integer *incy);
+extern int sspr_(char *uplo, integer *n, real *alpha, real *x, integer *incx, real *ap);
+extern int sspr2_(char *uplo, integer *n, real *alpha, real *x, integer *incx, real *y, integer *incy, real *ap);
+extern int sswap_(integer *n, real *sx, integer *incx, real *sy, integer *incy);
+extern int ssymm_(char *side, char *uplo, integer *m, integer *n, real *alpha, real *a, integer *lda, real *b, integer *ldb, real *beta, real *c__, integer *ldc);
+extern int ssymv_(char *uplo, integer *n, real *alpha, real *a, integer *lda, real *x, integer *incx, real *beta, real *y, integer *incy);
+extern int ssyr_(char *uplo, integer *n, real *alpha, real *x, integer *incx, real *a, integer *lda);
+extern int ssyr2_(char *uplo, integer *n, real *alpha, real *x, integer *incx, real *y, integer *incy, real *a, integer *lda);
+extern int ssyr2k_(char *uplo, char *trans, integer *n, integer *k, real *alpha, real *a, integer *lda, real *b, integer *ldb, real *beta, real *c__, integer *ldc);
+extern int ssyrk_(char *uplo, char *trans, integer *n, integer *k, real *alpha, real *a, integer *lda, real *beta, real *c__, integer *ldc);
+extern int stbmv_(char *uplo, char *trans, char *diag, integer *n, integer *k, real *a, integer *lda, real *x, integer *incx);
+extern int stbsv_(char *uplo, char *trans, char *diag, integer *n, integer *k, real *a, integer *lda, real *x, integer *incx);
+extern int stpmv_(char *uplo, char *trans, char *diag, integer *n, real *ap, real *x, integer *incx);
+extern int stpsv_(char *uplo, char *trans, char *diag, integer *n, real *ap, real *x, integer *incx);
+extern int strmm_(char *side, char *uplo, char *transa, char *diag, integer *m, integer *n, real *alpha, real *a, integer *lda, real *b, integer *ldb);
+extern int strmv_(char *uplo, char *trans, char *diag, integer *n, real *a, integer *lda, real *x, integer *incx);
+extern int strsm_(char *side, char *uplo, char *transa, char *diag, integer *m, integer *n, real *alpha, real *a, integer *lda, real *b, integer *ldb);
+extern int strsv_(char *uplo, char *trans, char *diag, integer *n, real *a, integer *lda, real *x, integer *incx);
+extern int xerbla_(char *srname, integer *info);
+extern int zaxpy_(integer *n, doublecomplex *za, doublecomplex *zx, integer *incx, doublecomplex *zy, integer *incy);
+extern int zcopy_(integer *n, doublecomplex *zx, integer *incx, doublecomplex *zy, integer *incy);
+extern Z_f zdotc_(doublecomplex * ret_val, integer *n, doublecomplex *zx, integer *incx, doublecomplex *zy, integer *incy);
+extern Z_f zdotu_(doublecomplex * ret_val, integer *n, doublecomplex *zx, integer *incx, doublecomplex *zy, integer *incy);
+extern int zdscal_(integer *n, doublereal *da, doublecomplex *zx, integer *incx);
+extern int zgbmv_(char *trans, integer *m, integer *n, integer *kl, integer *ku, doublecomplex *alpha, doublecomplex *a, integer *lda, doublecomplex *x, integer *incx, doublecomplex *beta, doublecomplex *y, integer *incy);
+extern int zgemm_(char *transa, char *transb, integer *m, integer *n, integer *k, doublecomplex *alpha, doublecomplex *a, integer *lda, doublecomplex *b, integer *ldb, doublecomplex *beta, doublecomplex *c__, integer *ldc);
+extern int zgemv_(char *trans, integer *m, integer *n, doublecomplex *alpha, doublecomplex *a, integer *lda, doublecomplex *x, integer *incx, doublecomplex *beta, doublecomplex *y, integer *incy);
+extern int zgerc_(integer *m, integer *n, doublecomplex *alpha, doublecomplex *x, integer *incx, doublecomplex *y, integer *incy, doublecomplex *a, integer *lda);
+extern int zgeru_(integer *m, integer *n, doublecomplex *alpha, doublecomplex *x, integer *incx, doublecomplex *y, integer *incy, doublecomplex *a, integer *lda);
+extern int zhbmv_(char *uplo, integer *n, integer *k, doublecomplex *alpha, doublecomplex *a, integer *lda, doublecomplex *x, integer *incx, doublecomplex *beta, doublecomplex *y, integer *incy);
+extern int zhemm_(char *side, char *uplo, integer *m, integer *n, doublecomplex *alpha, doublecomplex *a, integer *lda, doublecomplex *b, integer *ldb, doublecomplex *beta, doublecomplex *c__, integer *ldc);
+extern int zhemv_(char *uplo, integer *n, doublecomplex *alpha, doublecomplex *a, integer *lda, doublecomplex *x, integer *incx, doublecomplex *beta, doublecomplex *y, integer *incy);
+extern int zher_(char *uplo, integer *n, doublereal *alpha, doublecomplex *x, integer *incx, doublecomplex *a, integer *lda);
+extern int zher2_(char *uplo, integer *n, doublecomplex *alpha, doublecomplex *x, integer *incx, doublecomplex *y, integer *incy, doublecomplex *a, integer *lda);
+extern int zher2k_(char *uplo, char *trans, integer *n, integer *k, doublecomplex *alpha, doublecomplex *a, integer *lda, doublecomplex *b, integer *ldb, doublereal *beta, doublecomplex *c__, integer *ldc);
+extern int zherk_(char *uplo, char *trans, integer *n, integer *k, doublereal *alpha, doublecomplex *a, integer *lda, doublereal *beta, doublecomplex *c__, integer *ldc);
+extern int zhpmv_(char *uplo, integer *n, doublecomplex *alpha, doublecomplex *ap, doublecomplex *x, integer *incx, doublecomplex *beta, doublecomplex *y, integer *incy);
+extern int zhpr_(char *uplo, integer *n, doublereal *alpha, doublecomplex *x, integer *incx, doublecomplex *ap);
+extern int zhpr2_(char *uplo, integer *n, doublecomplex *alpha, doublecomplex *x, integer *incx, doublecomplex *y, integer *incy, doublecomplex *ap);
+extern int zrotg_(doublecomplex *ca, doublecomplex *cb, doublereal *c__, doublecomplex *s);
+extern int zscal_(integer *n, doublecomplex *za, doublecomplex *zx, integer *incx);
+extern int zswap_(integer *n, doublecomplex *zx, integer *incx, doublecomplex *zy, integer *incy);
+extern int zsymm_(char *side, char *uplo, integer *m, integer *n, doublecomplex *alpha, doublecomplex *a, integer *lda, doublecomplex *b, integer *ldb, doublecomplex *beta, doublecomplex *c__, integer *ldc);
+extern int zsyr2k_(char *uplo, char *trans, integer *n, integer *k, doublecomplex *alpha, doublecomplex *a, integer *lda, doublecomplex *b, integer *ldb, doublecomplex *beta, doublecomplex *c__, integer *ldc);
+extern int zsyrk_(char *uplo, char *trans, integer *n, integer *k, doublecomplex *alpha, doublecomplex *a, integer *lda, doublecomplex *beta, doublecomplex *c__, integer *ldc);
+extern int ztbmv_(char *uplo, char *trans, char *diag, integer *n, integer *k, doublecomplex *a, integer *lda, doublecomplex *x, integer *incx);
+extern int ztbsv_(char *uplo, char *trans, char *diag, integer *n, integer *k, doublecomplex *a, integer *lda, doublecomplex *x, integer *incx);
+extern int ztpmv_(char *uplo, char *trans, char *diag, integer *n, doublecomplex *ap, doublecomplex *x, integer *incx);
+extern int ztpsv_(char *uplo, char *trans, char *diag, integer *n, doublecomplex *ap, doublecomplex *x, integer *incx);
+extern int ztrmm_(char *side, char *uplo, char *transa, char *diag, integer *m, integer *n, doublecomplex *alpha, doublecomplex *a, integer *lda, doublecomplex *b, integer *ldb);
+extern int ztrmv_(char *uplo, char *trans, char *diag, integer *n, doublecomplex *a, integer *lda, doublecomplex *x, integer *incx);
+extern int ztrsm_(char *side, char *uplo, char *transa, char *diag, integer *m, integer *n, doublecomplex *alpha, doublecomplex *a, integer *lda, doublecomplex *b, integer *ldb);
+extern int ztrsv_(char *uplo, char *trans, char *diag, integer *n, doublecomplex *a, integer *lda, doublecomplex *x, integer *incx);
+
+#endif