nmatrix 0.0.1 → 0.0.2
Sign up to get free protection for your applications and to get access to all the features.
- data/.gitignore +27 -0
- data/.rspec +2 -0
- data/Gemfile +3 -5
- data/Guardfile +6 -0
- data/History.txt +33 -0
- data/Manifest.txt +41 -38
- data/README.rdoc +88 -11
- data/Rakefile +35 -53
- data/ext/nmatrix/data/complex.h +372 -0
- data/ext/nmatrix/data/data.cpp +275 -0
- data/ext/nmatrix/data/data.h +707 -0
- data/ext/nmatrix/data/rational.h +421 -0
- data/ext/nmatrix/data/ruby_object.h +446 -0
- data/ext/nmatrix/extconf.rb +101 -51
- data/ext/nmatrix/new_extconf.rb +56 -0
- data/ext/nmatrix/nmatrix.cpp +1609 -0
- data/ext/nmatrix/nmatrix.h +265 -849
- data/ext/nmatrix/ruby_constants.cpp +134 -0
- data/ext/nmatrix/ruby_constants.h +103 -0
- data/ext/nmatrix/storage/common.cpp +70 -0
- data/ext/nmatrix/storage/common.h +170 -0
- data/ext/nmatrix/storage/dense.cpp +665 -0
- data/ext/nmatrix/storage/dense.h +116 -0
- data/ext/nmatrix/storage/list.cpp +1088 -0
- data/ext/nmatrix/storage/list.h +129 -0
- data/ext/nmatrix/storage/storage.cpp +658 -0
- data/ext/nmatrix/storage/storage.h +99 -0
- data/ext/nmatrix/storage/yale.cpp +1601 -0
- data/ext/nmatrix/storage/yale.h +208 -0
- data/ext/nmatrix/ttable_helper.rb +126 -0
- data/ext/nmatrix/{yale/smmp1_header.template.c → types.h} +36 -9
- data/ext/nmatrix/util/io.cpp +295 -0
- data/ext/nmatrix/util/io.h +117 -0
- data/ext/nmatrix/util/lapack.h +1175 -0
- data/ext/nmatrix/util/math.cpp +557 -0
- data/ext/nmatrix/util/math.h +1363 -0
- data/ext/nmatrix/util/sl_list.cpp +475 -0
- data/ext/nmatrix/util/sl_list.h +255 -0
- data/ext/nmatrix/util/util.h +78 -0
- data/lib/nmatrix/blas.rb +70 -0
- data/lib/nmatrix/io/mat5_reader.rb +567 -0
- data/lib/nmatrix/io/mat_reader.rb +162 -0
- data/lib/{string.rb → nmatrix/monkeys.rb} +49 -2
- data/lib/nmatrix/nmatrix.rb +199 -0
- data/lib/nmatrix/nvector.rb +103 -0
- data/lib/nmatrix/version.rb +27 -0
- data/lib/nmatrix.rb +22 -230
- data/nmatrix.gemspec +59 -0
- data/scripts/mac-brew-gcc.sh +47 -0
- data/spec/4x4_sparse.mat +0 -0
- data/spec/4x5_dense.mat +0 -0
- data/spec/blas_spec.rb +47 -0
- data/spec/elementwise_spec.rb +164 -0
- data/spec/io_spec.rb +60 -0
- data/spec/lapack_spec.rb +52 -0
- data/spec/math_spec.rb +96 -0
- data/spec/nmatrix_spec.rb +93 -89
- data/spec/nmatrix_yale_spec.rb +52 -36
- data/spec/nvector_spec.rb +1 -1
- data/spec/slice_spec.rb +257 -0
- data/spec/spec_helper.rb +51 -0
- data/spec/utm5940.mtx +83844 -0
- metadata +113 -71
- data/.autotest +0 -23
- data/.gemtest +0 -0
- data/ext/nmatrix/cblas.c +0 -150
- data/ext/nmatrix/dense/blas_header.template.c +0 -52
- data/ext/nmatrix/dense/elementwise.template.c +0 -107
- data/ext/nmatrix/dense/gemm.template.c +0 -159
- data/ext/nmatrix/dense/gemv.template.c +0 -130
- data/ext/nmatrix/dense/rationalmath.template.c +0 -68
- data/ext/nmatrix/dense.c +0 -307
- data/ext/nmatrix/depend +0 -18
- data/ext/nmatrix/generator/syntax_tree.rb +0 -481
- data/ext/nmatrix/generator.rb +0 -594
- data/ext/nmatrix/list.c +0 -774
- data/ext/nmatrix/nmatrix.c +0 -1977
- data/ext/nmatrix/rational.c +0 -98
- data/ext/nmatrix/yale/complexmath.template.c +0 -71
- data/ext/nmatrix/yale/elementwise.template.c +0 -46
- data/ext/nmatrix/yale/elementwise_op.template.c +0 -73
- data/ext/nmatrix/yale/numbmm.template.c +0 -94
- data/ext/nmatrix/yale/smmp1.template.c +0 -21
- data/ext/nmatrix/yale/smmp2.template.c +0 -43
- data/ext/nmatrix/yale/smmp2_header.template.c +0 -46
- data/ext/nmatrix/yale/sort_columns.template.c +0 -56
- data/ext/nmatrix/yale/symbmm.template.c +0 -54
- data/ext/nmatrix/yale/transp.template.c +0 -68
- data/ext/nmatrix/yale.c +0 -726
- data/lib/array.rb +0 -67
- data/spec/syntax_tree_spec.rb +0 -46
@@ -0,0 +1,129 @@
|
|
1
|
+
/////////////////////////////////////////////////////////////////////
|
2
|
+
// = NMatrix
|
3
|
+
//
|
4
|
+
// A linear algebra library for scientific computation in Ruby.
|
5
|
+
// NMatrix is part of SciRuby.
|
6
|
+
//
|
7
|
+
// NMatrix was originally inspired by and derived from NArray, by
|
8
|
+
// Masahiro Tanaka: http://narray.rubyforge.org
|
9
|
+
//
|
10
|
+
// == Copyright Information
|
11
|
+
//
|
12
|
+
// SciRuby is Copyright (c) 2010 - 2012, Ruby Science Foundation
|
13
|
+
// NMatrix is Copyright (c) 2012, Ruby Science Foundation
|
14
|
+
//
|
15
|
+
// Please see LICENSE.txt for additional copyright notices.
|
16
|
+
//
|
17
|
+
// == Contributing
|
18
|
+
//
|
19
|
+
// By contributing source code to SciRuby, you agree to be bound by
|
20
|
+
// our Contributor Agreement:
|
21
|
+
//
|
22
|
+
// * https://github.com/SciRuby/sciruby/wiki/Contributor-Agreement
|
23
|
+
//
|
24
|
+
// == list.h
|
25
|
+
//
|
26
|
+
// List-of-lists n-dimensional matrix storage. Uses singly-linked
|
27
|
+
// lists.
|
28
|
+
|
29
|
+
#ifndef LIST_H
|
30
|
+
#define LIST_H
|
31
|
+
|
32
|
+
/*
|
33
|
+
* Standard Includes
|
34
|
+
*/
|
35
|
+
|
36
|
+
#include <stdlib.h>
|
37
|
+
|
38
|
+
/*
|
39
|
+
* Project Includes
|
40
|
+
*/
|
41
|
+
|
42
|
+
#include "types.h"
|
43
|
+
|
44
|
+
#include "data/data.h"
|
45
|
+
|
46
|
+
#include "common.h"
|
47
|
+
|
48
|
+
#include "util/sl_list.h"
|
49
|
+
|
50
|
+
#include "nmatrix.h"
|
51
|
+
|
52
|
+
/*
|
53
|
+
* Macros
|
54
|
+
*/
|
55
|
+
|
56
|
+
/*
|
57
|
+
* Types
|
58
|
+
*/
|
59
|
+
|
60
|
+
|
61
|
+
|
62
|
+
/*
|
63
|
+
* Data
|
64
|
+
*/
|
65
|
+
|
66
|
+
extern "C" {
|
67
|
+
|
68
|
+
/*
|
69
|
+
* Functions
|
70
|
+
*/
|
71
|
+
|
72
|
+
////////////////
|
73
|
+
// Lifecycle //
|
74
|
+
///////////////
|
75
|
+
|
76
|
+
LIST_STORAGE* nm_list_storage_create(dtype_t dtype, size_t* shape, size_t dim, void* init_val);
|
77
|
+
void nm_list_storage_delete(STORAGE* s);
|
78
|
+
void nm_list_storage_delete_ref(STORAGE* s);
|
79
|
+
void nm_list_storage_mark(void*);
|
80
|
+
|
81
|
+
///////////////
|
82
|
+
// Accessors //
|
83
|
+
///////////////
|
84
|
+
|
85
|
+
void* nm_list_storage_ref(STORAGE* s, SLICE* slice);
|
86
|
+
void* nm_list_storage_get(STORAGE* s, SLICE* slice);
|
87
|
+
void* nm_list_storage_insert(STORAGE* s, SLICE* slice, void* val);
|
88
|
+
void* nm_list_storage_remove(STORAGE* s, SLICE* slice);
|
89
|
+
|
90
|
+
///////////
|
91
|
+
// Tests //
|
92
|
+
///////////
|
93
|
+
|
94
|
+
bool nm_list_storage_eqeq(const STORAGE* left, const STORAGE* right);
|
95
|
+
|
96
|
+
//////////
|
97
|
+
// Math //
|
98
|
+
//////////
|
99
|
+
|
100
|
+
STORAGE* nm_list_storage_ew_op(nm::ewop_t op, const STORAGE* left, const STORAGE* right);
|
101
|
+
STORAGE* nm_list_storage_matrix_multiply(const STORAGE_PAIR& casted_storage, size_t* resulting_shape, bool vector);
|
102
|
+
|
103
|
+
|
104
|
+
/////////////
|
105
|
+
// Utility //
|
106
|
+
/////////////
|
107
|
+
|
108
|
+
size_t nm_list_storage_count_elements_r(const LIST* l, size_t recursions);
|
109
|
+
size_t nm_list_storage_count_nd_elements(const LIST_STORAGE* s);
|
110
|
+
|
111
|
+
/*
|
112
|
+
* Count non-zero elements. See also count_list_storage_nd_elements.
|
113
|
+
*/
|
114
|
+
inline size_t nm_list_storage_count_elements(const LIST_STORAGE* s) {
|
115
|
+
return nm_list_storage_count_elements_r(s->rows, s->dim - 1);
|
116
|
+
}
|
117
|
+
|
118
|
+
/////////////////////////
|
119
|
+
// Copying and Casting //
|
120
|
+
/////////////////////////
|
121
|
+
|
122
|
+
LIST_STORAGE* nm_list_storage_copy(const LIST_STORAGE* rhs);
|
123
|
+
STORAGE* nm_list_storage_copy_transposed(const STORAGE* rhs_base);
|
124
|
+
STORAGE* nm_list_storage_cast_copy(const STORAGE* rhs, dtype_t new_dtype);
|
125
|
+
VALUE nm_list_storage_to_hash(const LIST_STORAGE* s, const dtype_t dtype);
|
126
|
+
|
127
|
+
} // end of extern "C" block
|
128
|
+
|
129
|
+
#endif // LIST_H
|