numo-gsl-alt 0.1.3

data/ext/numo/gsl/spmatrix/erbpp_spmatrix.rb

@@ -0,0 +1,101 @@
+require_relative "../gen/erbpp_gsl"
+require_relative "parse_spmatrix"
+
+class_list = [
+ ["SpMatrix","spmatrix",[]],
+ ["SpBlas","spblas",[]],
+ ["IterSolve","splinalg_itersolve",[]],
+]
+list = {}
+ErbppGsl.read_func_pattern(*class_list.map{|name,base|
+  [/gsl_#{base}_/, list[name]=[]]
+})
+const_list = ErbppGsl.read_const
+
+DefLib.new do
+  set line_number: true
+  set erb_dir: %w[tmpl ../gen/tmpl]
+  set erb_suffix: ".c"
+  set ns_var: "mG"
+
+  ErbPP.new(self,"cast_1d_contiguous")
+  ErbPP.new(self,"cast_2d_contiguous")
+
+  cname = "SpMatrix"
+  cbase = cname.downcase
+  set file_name: "gsl_#{cname}.c"
+  set include_files: %w[gsl/gsl_spmatrix.h gsl/gsl_spblas.h gsl/gsl_splinalg.h]
+  set lib_name: cname.downcase
+
+  DefSpMatrix.new(self) do
+    name = cname
+    base = cbase
+    set name: base
+    set class_name: name
+    set class_var: "c"+name
+    set full_class_name: "Numo::GSL::"+name
+    set struct: "gsl_"+base
+
+    const_list.each do |a|
+      m = a[0].sub(/^GSL_SPMATRIX_/,"")
+      v = "INT2FIX(#{a[0]})"
+      def_const(m, v, desc:a[1]||"")
+    end
+
+    undef_alloc_func
+    list[name].each do |h|
+      check_func(h)
+    end
+  end
+
+  DefSpBlas.new(self) do
+    name = "SpBlas"
+    base = name.downcase
+    set name: base
+    set module_name: name
+    set module_var: "m"+name
+    set full_module_name: "Numo::GSL::"+name
+    set struct: "gsl_spmatrix"
+
+    set class_name: cname
+    set class_var: "c"+cname
+    set full_class_name: "Numo::GSL::"+cname
+    set data_type_var: cbase+"_data_type"
+
+    def_const "NOTRANS", "INT2FIX(CblasNoTrans)"
+    def_const "TRANS", "INT2FIX(CblasTrans)"
+
+    list[name].each do |h|
+      check_func(h)
+    end
+  end
+
+  def_module do
+    name = "SpLinalg"
+    base = name.downcase
+    set name: base
+    set module_name: name
+    set module_var: "m"+name
+    set full_module_name: "Numo::GSL::"+name
+
+    DefIterSolve.new(self) do
+      set ns_var: "mSpLinalg"
+      name = "IterSolve"
+      base = "splinalg_itersolve"
+      set name: base
+      set class_name: name
+      set class_var: "c"+name
+      set full_class_name: "Numo::GSL::SpLinalg::"+name
+      set struct: "gsl_"+base
+
+      set sm_struct: "gsl_"+cbase
+      set sm_data_type_var: cbase+"_data_type"
+
+      undef_alloc_func
+      list[name].each do |h|
+        check_func(h)
+      end
+    end
+  end
+
+end.run

data/ext/numo/gsl/spmatrix/extconf.rb

@@ -0,0 +1,9 @@
+require_relative '../extconf_gsl.rb'
+
+# source file to compile
+srcs = %w(
+ gsl_spmatrix
+)
+$objs = srcs.collect{|i| i+".o"}
+
+create_makefile('numo/gsl/spmatrix')

data/ext/numo/gsl/spmatrix/func_2.2.1.rb

@@ -0,0 +1,297 @@
+[{:func_name=>"gsl_spblas_dgemm",
+  :func_type=>"int",
+  :args=>
+   [["const double", "alpha"],
+    ["const gsl_spmatrix *", "A"],
+    ["const gsl_spmatrix *", "B"],
+    ["gsl_spmatrix *", "C"]],
+  :desc=>
+   "This function computes the sparse matrix-matrix product\n" +
+   "C = \\alpha A B. The matrices must be in compressed format."},
+ {:func_name=>"gsl_spblas_dgemv",
+  :func_type=>"int",
+  :args=>
+   [["const CBLAS_TRANSPOSE_t", "TransA"],
+    ["const double", "alpha"],
+    ["const gsl_spmatrix *", "A"],
+    ["const gsl_vector *", "x"],
+    ["const double", "beta"],
+    ["gsl_vector *", "y"]],
+  :desc=>
+   "This function computes the matrix-vector product and sum\n" +
+   "y \\leftarrow \\alpha op(A) x + \\beta y, where\n" +
+   "op(A) = A, A^T for TransA = CblasNoTrans,\n" +
+   "CblasTrans. In-place computations are not supported, so\n" +
+   "x and y must be distinct vectors.\n" +
+   "The matrix A may be in triplet or compressed format."},
+ {:func_name=>"gsl_splinalg_itersolve_alloc",
+  :func_type=>"gsl_splinalg_itersolve *",
+  :args=>
+   [["const gsl_splinalg_itersolve_type *", "T"],
+    ["const size_t", "n"],
+    ["const size_t", "m"]],
+  :desc=>
+   "This function allocates a workspace for the iterative solution of\n" +
+   "n-by-n sparse matrix systems. The iterative solver type\n" +
+   "is specified by T. The argument m specifies the size\n" +
+   "of the solution candidate subspace {\\cal K}_m. The dimension\n" +
+   "m may be set to 0 in which case a reasonable default value is used."},
+ {:func_name=>"gsl_splinalg_itersolve_free",
+  :func_type=>"void",
+  :args=>[["gsl_splinalg_itersolve *", "w"]],
+  :desc=>"This function frees the memory associated with the workspace w."},
+ {:func_name=>"gsl_splinalg_itersolve_iterate",
+  :func_type=>"int",
+  :args=>
+   [["const gsl_spmatrix *", "A"],
+    ["const gsl_vector *", "b"],
+    ["const double", "tol"],
+    ["gsl_vector *", "x"],
+    ["gsl_splinalg_itersolve *", "w"]],
+  :desc=>
+   "This function performs one iteration of the iterative method for\n" +
+   "the sparse linear system specfied by the matrix A, right hand\n" +
+   "side vector b and solution vector x. On input, x\n" +
+   "must be set to an initial guess for the solution. On output,\n" +
+   "x is updated to give the current solution estimate. The\n" +
+   "parameter tol specifies the relative tolerance between the residual\n" +
+   "norm and norm of b in order to check for convergence.\n" +
+   "When the following condition is satisfied:\n" +
+   "|| A x - b || <= tol * || b ||\n" +
+   "the method has converged, the function returns GSL_SUCCESS and\n" +
+   "the final solution is provided in x. Otherwise, the function\n" +
+   "returns GSL_CONTINUE to signal that more iterations are\n" +
+   "required. Here, || \\cdot || represents the Euclidean norm.\n" +
+   "The input matrix A may be in triplet or compressed column\n" +
+   "format."},
+ {:func_name=>"gsl_splinalg_itersolve_name",
+  :func_type=>"const char *",
+  :args=>[["const gsl_splinalg_itersolve *", "w"]],
+  :desc=>"This function returns a string pointer to the name of the solver."},
+ {:func_name=>"gsl_splinalg_itersolve_normr",
+  :func_type=>"double",
+  :args=>[["const gsl_splinalg_itersolve *", "w"]],
+  :desc=>
+   "This function returns the current residual norm\n" +
+   "||r|| = ||A x - b||, which is updated after each call to\n" +
+   "gsl_splinalg_itersolve_iterate."},
+ {:func_name=>"gsl_spmatrix_add",
+  :func_type=>"int",
+  :args=>
+   [["gsl_spmatrix *", "c"],
+    ["const gsl_spmatrix *", "a"],
+    ["const gsl_spmatrix *", "b"]],
+  :desc=>
+   "This function computes the sum c = a + b. The three matrices must\n" +
+   "have the same dimensions and be stored in a compressed format."},
+ {:func_name=>"gsl_spmatrix_alloc",
+  :func_type=>"gsl_spmatrix *",
+  :args=>[["const size_t", "n1"], ["const size_t", "n2"]],
+  :desc=>
+   "This function allocates a sparse matrix of size n1-by-n2 and\n" +
+   "initializes it to all zeros. If the size of the matrix is not known at allocation\n" +
+   "time, both n1 and n2 may be set to 1, and they will automatically\n" +
+   "grow as elements are added to the matrix. This function sets the\n" +
+   "matrix to the triplet representation, which is the easiest for adding\n" +
+   "and accessing matrix elements. This function tries to make a reasonable guess\n" +
+   "for the number of non-zero elements (nzmax) which will be added to the matrix by\n" +
+   "assuming a sparse density of 10\\%. The function\n" +
+   "gsl_spmatrix_alloc_nzmax can be used if this number is known more\n" +
+   "accurately. The workspace is of size O(nzmax)."},
+ {:func_name=>"gsl_spmatrix_alloc_nzmax",
+  :func_type=>"gsl_spmatrix *",
+  :args=>
+   [["const size_t", "n1"],
+    ["const size_t", "n2"],
+    ["const size_t", "nzmax"],
+    ["const size_t", "sptype"]],
+  :desc=>
+   "This function allocates a sparse matrix of size n1-by-n2 and\n" +
+   "initializes it to all zeros. If the size of the matrix is not known at allocation\n" +
+   "time, both n1 and n2 may be set to 1, and they will automatically\n" +
+   "grow as elements are added to the matrix. The parameter nzmax specifies\n" +
+   "the maximum number of non-zero elements which will be added to the matrix.\n" +
+   "It does not need to be precisely known in advance, since storage space will \n" +
+   "automatically grow using gsl_spmatrix_realloc if nzmax is not\n" +
+   "large enough. Accurate knowledge of this parameter reduces the number of\n" +
+   "reallocation calls required. The parameter sptype specifies the\n" +
+   "storage format of the sparse matrix. Possible values are\n" +
+   "The allocated gsl_spmatrix structure is of size O(nzmax)."},
+ {:func_name=>"gsl_spmatrix_ccs",
+  :func_type=>"gsl_spmatrix *",
+  :args=>[["const gsl_spmatrix *", "T"]],
+  :desc=>
+   "This function creates a sparse matrix in compressed column format\n" +
+   "from the input sparse matrix T which must be in triplet format.\n" +
+   "A pointer to a newly allocated matrix is returned. The calling function\n" +
+   "should free the newly allocated matrix when it is no longer needed."},
+ {:func_name=>"gsl_spmatrix_crs",
+  :func_type=>"gsl_spmatrix *",
+  :args=>[["const gsl_spmatrix *", "T"]],
+  :desc=>
+   "This function creates a sparse matrix in compressed row format\n" +
+   "from the input sparse matrix T which must be in triplet format.\n" +
+   "A pointer to a newly allocated matrix is returned. The calling function\n" +
+   "should free the newly allocated matrix when it is no longer needed."},
+ {:func_name=>"gsl_spmatrix_d2sp",
+  :func_type=>"int",
+  :args=>[["gsl_spmatrix *", "S"], ["const gsl_matrix *", "A"]],
+  :desc=>
+   "This function converts the dense matrix A into sparse triplet format\n" +
+   "and stores the result in S."},
+ {:func_name=>"gsl_spmatrix_equal",
+  :func_type=>"int",
+  :args=>[["const gsl_spmatrix *", "a"], ["const gsl_spmatrix *", "b"]],
+  :desc=>
+   "This function returns 1 if the matrices a and b are equal (by comparison of\n" +
+   "element values) and 0 otherwise. The matrices a and b must be in the same\n" +
+   "sparse storage format for comparison."},
+ {:func_name=>"gsl_spmatrix_fprintf",
+  :func_type=>"int",
+  :args=>
+   [["FILE *", "stream"],
+    ["const gsl_spmatrix *", "m"],
+    ["const char *", "format"]],
+  :desc=>
+   "This function writes the elements of the matrix m line-by-line to\n" +
+   "the stream stream using the format specifier format, which\n" +
+   "should be one of the %g, %e or %f formats for\n" +
+   "floating point numbers.  The function returns 0 for success and\n" +
+   "GSL_EFAILED if there was a problem writing to the file. The\n" +
+   "input matrix m may be in any storage format, and the output file\n" +
+   "will be written in MatrixMarket format."},
+ {:func_name=>"gsl_spmatrix_fread",
+  :func_type=>"int",
+  :args=>[["FILE *", "stream"], ["gsl_spmatrix *", "m"]],
+  :desc=>
+   "This function reads into the matrix m from the open stream\n" +
+   "stream in binary format.  The matrix m must be preallocated\n" +
+   "with the correct storage format, dimensions and have a sufficiently large nzmax\n" +
+   "in order to read in all matrix elements, otherwise GSL_EBADLEN\n" +
+   "is returned. The return value is 0 for success and\n" +
+   "GSL_EFAILED if there was a problem reading from the file.  The\n" +
+   "data is assumed to have been written in the native binary format on the\n" +
+   "same architecture."},
+ {:func_name=>"gsl_spmatrix_free",
+  :func_type=>"void",
+  :args=>[["gsl_spmatrix *", "m"]],
+  :desc=>
+   "This function frees the memory associated with the sparse matrix m."},
+ {:func_name=>"gsl_spmatrix_fscanf",
+  :func_type=>"gsl_spmatrix *",
+  :args=>[["FILE *", "stream"]],
+  :desc=>
+   "This function reads sparse matrix data in the MatrixMarket format\n" +
+   "from the stream stream and stores it in a newly allocated matrix\n" +
+   "which is returned in triplet format.  The function returns 0 for success and\n" +
+   "GSL_EFAILED if there was a problem reading from the file. The\n" +
+   "user should free the returned matrix when it is no longer needed."},
+ {:func_name=>"gsl_spmatrix_fwrite",
+  :func_type=>"int",
+  :args=>[["FILE *", "stream"], ["const gsl_spmatrix *", "m"]],
+  :desc=>
+   "This function writes the elements of the matrix m to the stream\n" +
+   "stream in binary format.  The return value is 0 for success and\n" +
+   "GSL_EFAILED if there was a problem writing to the file.  Since the\n" +
+   "data is written in the native binary format it may not be portable\n" +
+   "between different architectures."},
+ {:func_name=>"gsl_spmatrix_get",
+  :func_type=>"double",
+  :args=>
+   [["const gsl_spmatrix *", "m"],
+    ["const size_t", "i"],
+    ["const size_t", "j"]],
+  :desc=>
+   "This function returns element (i,j) of the matrix m.\n" +
+   "The matrix may be in triplet or compressed format."},
+ {:func_name=>"gsl_spmatrix_memcpy",
+  :func_type=>"int",
+  :args=>[["gsl_spmatrix *", "dest"], ["const gsl_spmatrix *", "src"]],
+  :desc=>
+   "This function copies the elements of the sparse matrix src into\n" +
+   "dest. The two matrices must have the same dimensions and be in the\n" +
+   "same storage format."},
+ {:func_name=>"gsl_spmatrix_minmax",
+  :func_type=>"int",
+  :args=>
+   [["const gsl_spmatrix *", "m"],
+    ["double *", "min_out"],
+    ["double *", "max_out"]],
+  :desc=>
+   "This function returns the minimum and maximum elements of the matrix\n" +
+   "m, storing them in min_out and max_out, and searching\n" +
+   "only the non-zero values."},
+ {:func_name=>"gsl_spmatrix_nnz",
+  :func_type=>"size_t",
+  :args=>[["const gsl_spmatrix *", "m"]],
+  :desc=>"This function returns the number of non-zero elements in m."},
+ {:func_name=>"gsl_spmatrix_ptr",
+  :func_type=>"double *",
+  :args=>
+   [["gsl_spmatrix *", "m"], ["const size_t", "i"], ["const size_t", "j"]],
+  :desc=>
+   "This function returns a pointer to the (i,j) element of the matrix m.\n" +
+   "If the (i,j) element is not explicitly stored in the matrix,\n" +
+   "a null pointer is returned."},
+ {:func_name=>"gsl_spmatrix_realloc",
+  :func_type=>"int",
+  :args=>[["const size_t", "nzmax"], ["gsl_spmatrix *", "m"]],
+  :desc=>
+   "This function reallocates the storage space for m to accomodate\n" +
+   "nzmax non-zero elements. It is typically called internally by\n" +
+   "gsl_spmatrix_set if the user wants to add more elements to the\n" +
+   "sparse matrix than the previously specified nzmax."},
+ {:func_name=>"gsl_spmatrix_scale",
+  :func_type=>"int",
+  :args=>[["gsl_spmatrix *", "m"], ["const double", "x"]],
+  :desc=>
+   "This function scales all elements of the matrix m by the constant\n" +
+   "factor x. The result m(i,j) \\leftarrow x m(i,j) is stored in m."},
+ {:func_name=>"gsl_spmatrix_set",
+  :func_type=>"int",
+  :args=>
+   [["gsl_spmatrix *", "m"],
+    ["const size_t", "i"],
+    ["const size_t", "j"],
+    ["const double", "x"]],
+  :desc=>
+   "This function sets element (i,j) of the matrix m to\n" +
+   "the value x. The matrix must be in triplet representation."},
+ {:func_name=>"gsl_spmatrix_set_zero",
+  :func_type=>"int",
+  :args=>[["gsl_spmatrix *", "m"]],
+  :desc=>
+   "This function sets (or resets) all the elements of the matrix m to zero."},
+ {:func_name=>"gsl_spmatrix_sp2d",
+  :func_type=>"int",
+  :args=>[["gsl_matrix *", "A"], ["const gsl_spmatrix *", "S"]],
+  :desc=>
+   "This function converts the sparse matrix S into a dense matrix and\n" +
+   "stores the result in A. S must be in triplet format."},
+ {:func_name=>"gsl_spmatrix_transpose",
+  :func_type=>"int",
+  :args=>[["gsl_spmatrix *", "m"]],
+  :desc=>
+   "This function replaces the matrix m by its transpose,\n" +
+   "preserving the storage format of the input matrix. Currently,\n" +
+   "only triplet matrix inputs are supported."},
+ {:func_name=>"gsl_spmatrix_transpose2",
+  :func_type=>"int",
+  :args=>[["gsl_spmatrix *", "m"]],
+  :desc=>
+   "This function replaces the matrix m by its transpose, but\n" +
+   "changes the storage format for compressed matrix inputs. Since\n" +
+   "compressed column storage is the transpose of compressed row storage,\n" +
+   "this function simply converts a CCS matrix to CRS and vice versa.\n" +
+   "This is the most efficient way to transpose a compressed storage\n" +
+   "matrix, but the user should note that the storage format of their\n" +
+   "compressed matrix will change on output. For triplet matrices,\n" +
+   "the output matrix is also in triplet storage."},
+ {:func_name=>"gsl_spmatrix_transpose_memcpy",
+  :func_type=>"int",
+  :args=>[["gsl_spmatrix *", "dest"], ["const gsl_spmatrix *", "src"]],
+  :desc=>
+   "This function copies the transpose of the sparse matrix src into\n" +
+   "dest. The dimensions of dest must match the transpose of the\n" +
+   "matrix src. Also, both matrices must use the same sparse storage\n" +
+   "format."}]