linmeric 0.1.0 → 0.2.0

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (109) hide show
  1. checksums.yaml +4 -4
  2. data/Gemfile +6 -0
  3. data/LICENSE.txt +21 -0
  4. data/README.md +51 -0
  5. data/Rakefile +10 -0
  6. data/bin/help/Help.rb +184 -0
  7. data/bin/help/Help_inst.rb +244 -0
  8. data/bin/linguide +23 -0
  9. data/bin/linmeric +36 -70
  10. data/doc/Archive.html +352 -0
  11. data/doc/Calculator/Evaluator.html +477 -0
  12. data/doc/Calculator/Lexer.html +186 -0
  13. data/doc/Calculator/Token.html +257 -0
  14. data/doc/Calculator.html +181 -0
  15. data/doc/Dim.html +257 -0
  16. data/doc/Filename.html +246 -0
  17. data/doc/Fixnum.html +253 -0
  18. data/doc/Float.html +253 -0
  19. data/doc/Function.html +784 -0
  20. data/doc/InputError.html +102 -0
  21. data/doc/Instructions_en.txt +6 -7
  22. data/doc/Instructions_it.txt +6 -9
  23. data/doc/Integrators.html +436 -0
  24. data/doc/LU.html +435 -0
  25. data/doc/Lexer.html +261 -0
  26. data/doc/Linmeric.html +109 -0
  27. data/doc/Listener.html +605 -0
  28. data/doc/Matrix.html +1719 -0
  29. data/doc/MyArgError.html +102 -0
  30. data/doc/NilClass.html +202 -0
  31. data/doc/Numeric.html +251 -0
  32. data/doc/Parser.html +389 -0
  33. data/doc/PrintError.html +622 -0
  34. data/doc/README_md.html +142 -0
  35. data/doc/Scp.html +530 -0
  36. data/doc/Sizer.html +652 -0
  37. data/doc/String.html +540 -0
  38. data/doc/Token.html +341 -0
  39. data/doc/Tool.html +394 -0
  40. data/doc/created.rid +18 -0
  41. data/doc/css/fonts.css +167 -0
  42. data/doc/css/rdoc.css +590 -0
  43. data/doc/fonts/Lato-Light.ttf +0 -0
  44. data/doc/fonts/Lato-LightItalic.ttf +0 -0
  45. data/doc/fonts/Lato-Regular.ttf +0 -0
  46. data/doc/fonts/Lato-RegularItalic.ttf +0 -0
  47. data/doc/fonts/SourceCodePro-Bold.ttf +0 -0
  48. data/doc/fonts/SourceCodePro-Regular.ttf +0 -0
  49. data/doc/images/add.png +0 -0
  50. data/doc/images/arrow_up.png +0 -0
  51. data/doc/images/brick.png +0 -0
  52. data/doc/images/brick_link.png +0 -0
  53. data/doc/images/bug.png +0 -0
  54. data/doc/images/bullet_black.png +0 -0
  55. data/doc/images/bullet_toggle_minus.png +0 -0
  56. data/doc/images/bullet_toggle_plus.png +0 -0
  57. data/doc/images/date.png +0 -0
  58. data/doc/images/delete.png +0 -0
  59. data/doc/images/find.png +0 -0
  60. data/doc/images/loadingAnimation.gif +0 -0
  61. data/doc/images/macFFBgHack.png +0 -0
  62. data/doc/images/package.png +0 -0
  63. data/doc/images/page_green.png +0 -0
  64. data/doc/images/page_white_text.png +0 -0
  65. data/doc/images/page_white_width.png +0 -0
  66. data/doc/images/plugin.png +0 -0
  67. data/doc/images/ruby.png +0 -0
  68. data/doc/images/tag_blue.png +0 -0
  69. data/doc/images/tag_green.png +0 -0
  70. data/doc/images/transparent.png +0 -0
  71. data/doc/images/wrench.png +0 -0
  72. data/doc/images/wrench_orange.png +0 -0
  73. data/doc/images/zoom.png +0 -0
  74. data/doc/index.html +190 -0
  75. data/doc/js/darkfish.js +161 -0
  76. data/doc/js/jquery.js +9404 -0
  77. data/doc/js/navigation.js +142 -0
  78. data/doc/js/navigation.js.gz +0 -0
  79. data/doc/js/search.js +109 -0
  80. data/doc/js/search_index.js +1 -0
  81. data/doc/js/search_index.js.gz +0 -0
  82. data/doc/js/searcher.js +228 -0
  83. data/doc/js/searcher.js.gz +0 -0
  84. data/doc/table_of_contents.html +834 -0
  85. data/lib/linmeric/Archive.rb +22 -1
  86. data/lib/linmeric/Calculator.rb +252 -0
  87. data/lib/linmeric/CnGal_Matrix_class.rb +130 -49
  88. data/lib/linmeric/CnGal_new_classes.rb +139 -37
  89. data/lib/linmeric/CnGal_tools.rb +23 -40
  90. data/lib/linmeric/Error_print.rb +35 -1
  91. data/lib/linmeric/Function_class.rb +70 -11
  92. data/lib/linmeric/Integrators.rb +81 -35
  93. data/lib/linmeric/LU.rb +26 -5
  94. data/lib/linmeric/Lexer.rb +19 -1
  95. data/lib/linmeric/Listener.rb +25 -5
  96. data/lib/linmeric/Parser.rb +23 -1
  97. data/lib/linmeric/Scopify.rb +52 -31
  98. data/lib/linmeric/Sizer.rb +43 -10
  99. data/lib/linmeric/Token.rb +23 -4
  100. data/lib/linmeric/version.rb +3 -0
  101. data/lib/linmeric.rb +10 -8
  102. data/lib/linmeric_bin.rb +12 -0
  103. metadata +126 -22
  104. data/doc/Instructions_en.html +0 -231
  105. data/doc/Instructions_it.html +0 -231
  106. data/doc/README_en.html +0 -177
  107. data/doc/README_en.txt +0 -30
  108. data/doc/README_it.html +0 -187
  109. data/doc/README_it.txt +0 -32
data/doc/LU.html ADDED
@@ -0,0 +1,435 @@
1
+ <!DOCTYPE html>
2
+
3
+ <html>
4
+ <head>
5
+ <meta charset="UTF-8">
6
+
7
+ <title>module LU - RDoc Documentation</title>
8
+
9
+ <script type="text/javascript">
10
+ var rdoc_rel_prefix = "./";
11
+ </script>
12
+
13
+ <script src="./js/jquery.js"></script>
14
+ <script src="./js/darkfish.js"></script>
15
+
16
+ <link href="./css/fonts.css" rel="stylesheet">
17
+ <link href="./css/rdoc.css" rel="stylesheet">
18
+
19
+
20
+
21
+ <body id="top" role="document" class="module">
22
+ <nav role="navigation">
23
+ <div id="project-navigation">
24
+ <div id="home-section" role="region" title="Quick navigation" class="nav-section">
25
+ <h2>
26
+ <a href="./index.html" rel="home">Home</a>
27
+ </h2>
28
+
29
+ <div id="table-of-contents-navigation">
30
+ <a href="./table_of_contents.html#pages">Pages</a>
31
+ <a href="./table_of_contents.html#classes">Classes</a>
32
+ <a href="./table_of_contents.html#methods">Methods</a>
33
+ </div>
34
+ </div>
35
+
36
+ <div id="search-section" role="search" class="project-section initially-hidden">
37
+ <form action="#" method="get" accept-charset="utf-8">
38
+ <div id="search-field-wrapper">
39
+ <input id="search-field" role="combobox" aria-label="Search"
40
+ aria-autocomplete="list" aria-controls="search-results"
41
+ type="text" name="search" placeholder="Search" spellcheck="false"
42
+ title="Type to search, Up and Down to navigate, Enter to load">
43
+ </div>
44
+
45
+ <ul id="search-results" aria-label="Search Results"
46
+ aria-busy="false" aria-expanded="false"
47
+ aria-atomic="false" class="initially-hidden"></ul>
48
+ </form>
49
+ </div>
50
+
51
+ </div>
52
+
53
+
54
+
55
+ <div id="class-metadata">
56
+
57
+
58
+
59
+
60
+ <!-- Method Quickref -->
61
+ <div id="method-list-section" class="nav-section">
62
+ <h3>Methods</h3>
63
+
64
+ <ul class="link-list" role="directory">
65
+
66
+ <li ><a href="#method-c-L">::L</a>
67
+
68
+ <li ><a href="#method-c-U">::U</a>
69
+
70
+ <li ><a href="#method-c-factorize">::factorize</a>
71
+
72
+ <li ><a href="#method-c-reset">::reset</a>
73
+
74
+ <li ><a href="#method-c-solve">::solve</a>
75
+
76
+ <li ><a href="#method-c-swap">::swap</a>
77
+
78
+ </ul>
79
+ </div>
80
+
81
+ </div>
82
+ </nav>
83
+
84
+ <main role="main" aria-labelledby="module-LU">
85
+ <h1 id="module-LU" class="module">
86
+ module LU
87
+ </h1>
88
+
89
+ <section class="description">
90
+
91
+ <p>This module provides some method to perform <a href="LU.html">LU</a>
92
+ factorization on squared matrices.</p>
93
+ <dl class="rdoc-list note-list"><dt>Author
94
+ <dd>
95
+ <p>Massimiliano Dal Mas (<a
96
+ href="mailto:max.codeware@gmail.com">max.codeware@gmail.com</a>)</p>
97
+ </dd><dt>License
98
+ <dd>
99
+ <p>Distributed under MIT license</p>
100
+ </dd></dl>
101
+
102
+ </section>
103
+
104
+
105
+
106
+
107
+ <section id="5Buntitled-5D" class="documentation-section">
108
+
109
+
110
+
111
+
112
+
113
+
114
+
115
+
116
+
117
+ <section id="public-class-5Buntitled-5D-method-details" class="method-section">
118
+ <header>
119
+ <h3>Public Class Methods</h3>
120
+ </header>
121
+
122
+
123
+ <div id="method-c-L" class="method-detail ">
124
+
125
+ <div class="method-heading">
126
+ <span class="method-name">L</span><span
127
+ class="method-args">()</span>
128
+
129
+ <span class="method-click-advice">click to toggle source</span>
130
+
131
+ </div>
132
+
133
+
134
+ <div class="method-description">
135
+
136
+ <p>Returns the `L` matrix</p>
137
+
138
+
139
+
140
+
141
+ <div class="method-source-code" id="L-source">
142
+ <pre><span class="ruby-comment"># File lib/linmeric/LU.rb, line 27</span>
143
+ <span class="ruby-keyword">def</span> <span class="ruby-keyword">self</span>.<span class="ruby-constant">L</span>()
144
+ <span class="ruby-keyword">return</span> <span class="ruby-ivar">@L</span>
145
+ <span class="ruby-keyword">end</span></pre>
146
+ </div>
147
+
148
+ </div>
149
+
150
+
151
+
152
+
153
+ </div>
154
+
155
+
156
+ <div id="method-c-U" class="method-detail ">
157
+
158
+ <div class="method-heading">
159
+ <span class="method-name">U</span><span
160
+ class="method-args">()</span>
161
+
162
+ <span class="method-click-advice">click to toggle source</span>
163
+
164
+ </div>
165
+
166
+
167
+ <div class="method-description">
168
+
169
+ <p>Returns the `U` matrix</p>
170
+
171
+
172
+
173
+
174
+ <div class="method-source-code" id="U-source">
175
+ <pre><span class="ruby-comment"># File lib/linmeric/LU.rb, line 32</span>
176
+ <span class="ruby-keyword">def</span> <span class="ruby-keyword">self</span>.<span class="ruby-constant">U</span>()
177
+ <span class="ruby-keyword">return</span> <span class="ruby-ivar">@U</span>
178
+ <span class="ruby-keyword">end</span></pre>
179
+ </div>
180
+
181
+ </div>
182
+
183
+
184
+
185
+
186
+ </div>
187
+
188
+
189
+ <div id="method-c-factorize" class="method-detail ">
190
+
191
+ <div class="method-heading">
192
+ <span class="method-name">factorize</span><span
193
+ class="method-args">(mx,sol)</span>
194
+
195
+ <span class="method-click-advice">click to toggle source</span>
196
+
197
+ </div>
198
+
199
+
200
+ <div class="method-description">
201
+
202
+ <p>Performs <a href="LU.html">LU</a> factorization calculating <a
203
+ href="LU.html#method-c-L">::L</a> and <a href="LU.html#method-c-U">::U</a>
204
+ matrices</p>
205
+ <ul><li>
206
+ <p>*<strong>argument</strong>*: squared matrix to factorize</p>
207
+ </li><li>
208
+ <p>*<strong>argument</strong>*: n x 1 matrix of known values of the linear
209
+ system</p>
210
+ </li><li>
211
+ <p>*<strong>returns</strong>*: same as <a
212
+ href="LU.html#method-c-solve">::solve</a></p>
213
+ </li></ul>
214
+
215
+
216
+
217
+
218
+ <div class="method-source-code" id="factorize-source">
219
+ <pre><span class="ruby-comment"># File lib/linmeric/LU.rb, line 47</span>
220
+ <span class="ruby-keyword">def</span> <span class="ruby-keyword">self</span>.<span class="ruby-identifier">factorize</span>(<span class="ruby-identifier">mx</span>,<span class="ruby-identifier">sol</span>)
221
+ [<span class="ruby-identifier">mx</span>,<span class="ruby-identifier">sol</span>].<span class="ruby-identifier">each</span> <span class="ruby-keyword">do</span> <span class="ruby-operator">|</span><span class="ruby-identifier">vec</span><span class="ruby-operator">|</span>
222
+ <span class="ruby-keyword">return</span> <span class="ruby-keyword">nil</span> <span class="ruby-keyword">unless</span> <span class="ruby-identifier">vec</span>.<span class="ruby-identifier">is_a?</span> <span class="ruby-constant">Matrix</span>
223
+ <span class="ruby-keyword">end</span>
224
+ <span class="ruby-identifier">sol</span> = <span class="ruby-identifier">sol</span>.<span class="ruby-identifier">tr</span> <span class="ruby-keyword">if</span> <span class="ruby-identifier">sol</span>.<span class="ruby-identifier">getCls</span> <span class="ruby-operator">&gt;</span> <span class="ruby-value">1</span>
225
+ <span class="ruby-keyword">return</span> <span class="ruby-keyword">nil</span> <span class="ruby-keyword">unless</span> <span class="ruby-identifier">sol</span>.<span class="ruby-identifier">getCls</span> <span class="ruby-operator">==</span> <span class="ruby-value">1</span>
226
+ <span class="ruby-keyword">return</span> <span class="ruby-keyword">nil</span> <span class="ruby-keyword">unless</span> <span class="ruby-identifier">mx</span>.<span class="ruby-identifier">is_squared?</span>
227
+ <span class="ruby-keyword">return</span> <span class="ruby-keyword">nil</span> <span class="ruby-keyword">unless</span> <span class="ruby-identifier">mx</span>.<span class="ruby-identifier">getRws</span> <span class="ruby-operator">==</span> <span class="ruby-identifier">sol</span>.<span class="ruby-identifier">getRws</span>
228
+ <span class="ruby-identifier">rows</span> = <span class="ruby-identifier">mx</span>.<span class="ruby-identifier">getRws</span>
229
+ <span class="ruby-keyword">for</span> <span class="ruby-identifier">k</span> <span class="ruby-keyword">in</span> <span class="ruby-value">0</span><span class="ruby-operator">...</span><span class="ruby-identifier">rows</span> <span class="ruby-keyword">do</span>
230
+ <span class="ruby-identifier">column</span> = <span class="ruby-identifier">mx</span>[<span class="ruby-identifier">k</span><span class="ruby-operator">...</span><span class="ruby-identifier">mx</span>.<span class="ruby-identifier">getRws</span>,<span class="ruby-identifier">k</span>].<span class="ruby-identifier">export</span>.<span class="ruby-identifier">map!</span> { <span class="ruby-operator">|</span><span class="ruby-identifier">val</span><span class="ruby-operator">|</span> <span class="ruby-identifier">val</span>.<span class="ruby-identifier">abs</span>}
231
+ <span class="ruby-identifier">max_index</span> = <span class="ruby-identifier">column</span>.<span class="ruby-identifier">index</span>(<span class="ruby-identifier">column</span>.<span class="ruby-identifier">max</span>)
232
+ <span class="ruby-identifier">mx</span> = <span class="ruby-keyword">self</span>.<span class="ruby-identifier">swap</span>(<span class="ruby-identifier">mx</span>,<span class="ruby-identifier">k</span>,<span class="ruby-identifier">k</span> <span class="ruby-operator">+</span> <span class="ruby-identifier">max_index</span>) <span class="ruby-keyword">unless</span> <span class="ruby-identifier">k</span> <span class="ruby-operator">==</span> <span class="ruby-identifier">max_index</span> <span class="ruby-operator">+</span> <span class="ruby-identifier">k</span>
233
+ <span class="ruby-identifier">sol</span> = <span class="ruby-keyword">self</span>.<span class="ruby-identifier">swap</span>(<span class="ruby-identifier">sol</span>,<span class="ruby-value">0</span>,<span class="ruby-identifier">k</span> <span class="ruby-operator">+</span> <span class="ruby-identifier">max_index</span>) <span class="ruby-keyword">unless</span> <span class="ruby-identifier">k</span> <span class="ruby-operator">==</span> <span class="ruby-identifier">max_index</span> <span class="ruby-operator">+</span> <span class="ruby-identifier">k</span>
234
+ <span class="ruby-keyword">for</span> <span class="ruby-identifier">i</span> <span class="ruby-keyword">in</span> (<span class="ruby-identifier">k</span><span class="ruby-operator">+</span><span class="ruby-value">1</span>)<span class="ruby-operator">...</span><span class="ruby-identifier">mx</span>.<span class="ruby-identifier">getRws</span> <span class="ruby-keyword">do</span>
235
+ <span class="ruby-identifier">alpha</span> = (<span class="ruby-identifier">mx</span>[<span class="ruby-identifier">k</span>,<span class="ruby-identifier">k</span>] <span class="ruby-operator">!=</span> <span class="ruby-value">0</span>) <span class="ruby-operator">?</span> (<span class="ruby-identifier">mx</span>[<span class="ruby-identifier">i</span>,<span class="ruby-identifier">k</span>] <span class="ruby-operator">/</span> <span class="ruby-identifier">mx</span>[<span class="ruby-identifier">k</span>,<span class="ruby-identifier">k</span>].<span class="ruby-identifier">to_f</span>) <span class="ruby-operator">:</span> <span class="ruby-value">0</span>
236
+ <span class="ruby-identifier">mx</span>[<span class="ruby-identifier">i</span>,<span class="ruby-identifier">k</span>] = <span class="ruby-identifier">alpha</span>
237
+ <span class="ruby-keyword">for</span> <span class="ruby-identifier">j</span> <span class="ruby-keyword">in</span> (<span class="ruby-identifier">k</span><span class="ruby-operator">+</span><span class="ruby-value">1</span>)<span class="ruby-operator">...</span><span class="ruby-identifier">mx</span>.<span class="ruby-identifier">getCls</span> <span class="ruby-keyword">do</span>
238
+ <span class="ruby-identifier">mx</span>[<span class="ruby-identifier">i</span>,<span class="ruby-identifier">j</span>] <span class="ruby-operator">-=</span> <span class="ruby-identifier">alpha</span> <span class="ruby-operator">*</span> <span class="ruby-identifier">mx</span>[<span class="ruby-identifier">k</span>,<span class="ruby-identifier">j</span>]
239
+ <span class="ruby-keyword">end</span>
240
+ <span class="ruby-keyword">end</span>
241
+ <span class="ruby-keyword">end</span>
242
+ <span class="ruby-ivar">@L</span> = <span class="ruby-constant">Matrix</span>.<span class="ruby-identifier">identity</span>(<span class="ruby-identifier">mx</span>.<span class="ruby-identifier">getRws</span>) <span class="ruby-operator">+</span> <span class="ruby-constant">Matrix</span>.<span class="ruby-identifier">new</span>(<span class="ruby-identifier">mx</span>.<span class="ruby-identifier">getRws</span>,<span class="ruby-identifier">mx</span>.<span class="ruby-identifier">getCls</span>){ <span class="ruby-operator">|</span><span class="ruby-identifier">i</span>,<span class="ruby-identifier">j</span><span class="ruby-operator">|</span> (<span class="ruby-identifier">i</span> <span class="ruby-operator">&gt;</span> <span class="ruby-identifier">j</span>) <span class="ruby-operator">?</span> <span class="ruby-identifier">mx</span>[<span class="ruby-identifier">i</span>,<span class="ruby-identifier">j</span>] <span class="ruby-operator">:</span> <span class="ruby-value">0</span>}
243
+ <span class="ruby-ivar">@U</span> = <span class="ruby-constant">Matrix</span>.<span class="ruby-identifier">new</span>(<span class="ruby-identifier">mx</span>.<span class="ruby-identifier">getRws</span>,<span class="ruby-identifier">mx</span>.<span class="ruby-identifier">getCls</span>){ <span class="ruby-operator">|</span><span class="ruby-identifier">i</span>,<span class="ruby-identifier">j</span><span class="ruby-operator">|</span> (<span class="ruby-identifier">i</span> <span class="ruby-operator">&lt;=</span> <span class="ruby-identifier">j</span>) <span class="ruby-operator">?</span> <span class="ruby-identifier">mx</span>[<span class="ruby-identifier">i</span>,<span class="ruby-identifier">j</span>] <span class="ruby-operator">:</span> <span class="ruby-value">0</span>}
244
+ <span class="ruby-keyword">return</span> <span class="ruby-identifier">solve</span>(<span class="ruby-identifier">sol</span>)
245
+ <span class="ruby-keyword">end</span>
246
+
247
+ <span class="ruby-comment"># Finds the solutions of the linear system</span>
248
+ <span class="ruby-comment">#</span>
249
+ <span class="ruby-comment"># * **argument**: n x 1 matrix of known values of the linear system</span>
250
+ <span class="ruby-comment"># * **returns**: n x 1 matrix with the solutions of the system</span>
251
+ <span class="ruby-keyword">def</span> <span class="ruby-keyword">self</span>.<span class="ruby-identifier">solve</span>(<span class="ruby-identifier">sol</span>)
252
+ <span class="ruby-identifier">z</span> = <span class="ruby-constant">Matrix</span>.<span class="ruby-identifier">new</span>(<span class="ruby-identifier">sol</span>.<span class="ruby-identifier">getRws</span>,<span class="ruby-value">1</span>) {<span class="ruby-value">0</span>}
253
+ <span class="ruby-identifier">x</span> = <span class="ruby-constant">Matrix</span>.<span class="ruby-identifier">new</span>(<span class="ruby-identifier">sol</span>.<span class="ruby-identifier">getRws</span>,<span class="ruby-value">1</span>) {<span class="ruby-value">0</span>}
254
+ <span class="ruby-keyword">for</span> <span class="ruby-identifier">i</span> <span class="ruby-keyword">in</span> <span class="ruby-value">0</span><span class="ruby-operator">...</span><span class="ruby-identifier">sol</span>.<span class="ruby-identifier">getRws</span> <span class="ruby-keyword">do</span>
255
+ <span class="ruby-identifier">z</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>] = <span class="ruby-identifier">sol</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>]
256
+ <span class="ruby-keyword">for</span> <span class="ruby-identifier">j</span> <span class="ruby-keyword">in</span> <span class="ruby-value">0</span><span class="ruby-operator">...</span><span class="ruby-identifier">i</span> <span class="ruby-keyword">do</span>
257
+ <span class="ruby-identifier">z</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>] <span class="ruby-operator">-=</span> <span class="ruby-ivar">@L</span>[<span class="ruby-identifier">i</span>,<span class="ruby-identifier">j</span>] <span class="ruby-operator">*</span> (<span class="ruby-identifier">z</span>[<span class="ruby-identifier">j</span>,<span class="ruby-value">0</span>] <span class="ruby-operator">||</span> <span class="ruby-value">1</span>)
258
+ <span class="ruby-keyword">end</span>
259
+ <span class="ruby-keyword">end</span>
260
+
261
+ (<span class="ruby-identifier">sol</span>.<span class="ruby-identifier">getRws</span> <span class="ruby-operator">-</span> <span class="ruby-value">1</span>).<span class="ruby-identifier">downto</span>(<span class="ruby-value">0</span>) <span class="ruby-keyword">do</span> <span class="ruby-operator">|</span><span class="ruby-identifier">i</span><span class="ruby-operator">|</span>
262
+ <span class="ruby-identifier">x</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>] = <span class="ruby-identifier">z</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>]
263
+ <span class="ruby-keyword">for</span> <span class="ruby-identifier">j</span> <span class="ruby-keyword">in</span> <span class="ruby-identifier">i</span><span class="ruby-operator">...</span><span class="ruby-identifier">sol</span>.<span class="ruby-identifier">getRws</span>
264
+ <span class="ruby-identifier">x</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>] <span class="ruby-operator">-=</span> (<span class="ruby-ivar">@U</span>[<span class="ruby-identifier">i</span>,<span class="ruby-identifier">j</span><span class="ruby-operator">+</span><span class="ruby-value">1</span>] <span class="ruby-operator">||</span> <span class="ruby-value">0</span>) <span class="ruby-operator">*</span> (<span class="ruby-identifier">x</span>[<span class="ruby-identifier">j</span><span class="ruby-operator">+</span><span class="ruby-value">1</span>,<span class="ruby-value">0</span>] <span class="ruby-operator">||</span> <span class="ruby-value">0</span>)
265
+ <span class="ruby-keyword">end</span>
266
+ <span class="ruby-identifier">x</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>] = <span class="ruby-identifier">x</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>] <span class="ruby-operator">/</span> <span class="ruby-ivar">@U</span>[<span class="ruby-identifier">i</span>,<span class="ruby-identifier">i</span>].<span class="ruby-identifier">to_f</span>
267
+ <span class="ruby-keyword">end</span>
268
+ <span class="ruby-keyword">return</span> <span class="ruby-identifier">x</span>
269
+ <span class="ruby-keyword">end</span>
270
+
271
+
272
+ <span class="ruby-keyword">end</span>
273
+
274
+
275
+
276
+
277
+
278
+
279
+ </pre>
280
+ </div>
281
+
282
+ </div>
283
+
284
+
285
+
286
+
287
+ </div>
288
+
289
+
290
+ <div id="method-c-reset" class="method-detail ">
291
+
292
+ <div class="method-heading">
293
+ <span class="method-name">reset</span><span
294
+ class="method-args">()</span>
295
+
296
+ <span class="method-click-advice">click to toggle source</span>
297
+
298
+ </div>
299
+
300
+
301
+ <div class="method-description">
302
+
303
+ <p>Sets to nil @L and @U variables</p>
304
+
305
+
306
+
307
+
308
+ <div class="method-source-code" id="reset-source">
309
+ <pre><span class="ruby-comment"># File lib/linmeric/LU.rb, line 37</span>
310
+ <span class="ruby-keyword">def</span> <span class="ruby-keyword">self</span>.<span class="ruby-identifier">reset</span>
311
+ <span class="ruby-ivar">@L</span> = <span class="ruby-keyword">nil</span>
312
+ <span class="ruby-ivar">@U</span> = <span class="ruby-keyword">nil</span>
313
+ <span class="ruby-keyword">end</span></pre>
314
+ </div>
315
+
316
+ </div>
317
+
318
+
319
+
320
+
321
+ </div>
322
+
323
+
324
+ <div id="method-c-solve" class="method-detail ">
325
+
326
+ <div class="method-heading">
327
+ <span class="method-name">solve</span><span
328
+ class="method-args">(sol)</span>
329
+
330
+ <span class="method-click-advice">click to toggle source</span>
331
+
332
+ </div>
333
+
334
+
335
+ <div class="method-description">
336
+
337
+ <p>Finds the solutions of the linear system</p>
338
+ <ul><li>
339
+ <p>*<strong>argument</strong>*: n x 1 matrix of known values of the linear
340
+ system</p>
341
+ </li><li>
342
+ <p>*<strong>returns</strong>*: n x 1 matrix with the solutions of the system</p>
343
+ </li></ul>
344
+
345
+
346
+
347
+
348
+ <div class="method-source-code" id="solve-source">
349
+ <pre><span class="ruby-comment"># File lib/linmeric/LU.rb, line 78</span>
350
+ <span class="ruby-keyword">def</span> <span class="ruby-keyword">self</span>.<span class="ruby-identifier">solve</span>(<span class="ruby-identifier">sol</span>)
351
+ <span class="ruby-identifier">z</span> = <span class="ruby-constant">Matrix</span>.<span class="ruby-identifier">new</span>(<span class="ruby-identifier">sol</span>.<span class="ruby-identifier">getRws</span>,<span class="ruby-value">1</span>) {<span class="ruby-value">0</span>}
352
+ <span class="ruby-identifier">x</span> = <span class="ruby-constant">Matrix</span>.<span class="ruby-identifier">new</span>(<span class="ruby-identifier">sol</span>.<span class="ruby-identifier">getRws</span>,<span class="ruby-value">1</span>) {<span class="ruby-value">0</span>}
353
+ <span class="ruby-keyword">for</span> <span class="ruby-identifier">i</span> <span class="ruby-keyword">in</span> <span class="ruby-value">0</span><span class="ruby-operator">...</span><span class="ruby-identifier">sol</span>.<span class="ruby-identifier">getRws</span> <span class="ruby-keyword">do</span>
354
+ <span class="ruby-identifier">z</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>] = <span class="ruby-identifier">sol</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>]
355
+ <span class="ruby-keyword">for</span> <span class="ruby-identifier">j</span> <span class="ruby-keyword">in</span> <span class="ruby-value">0</span><span class="ruby-operator">...</span><span class="ruby-identifier">i</span> <span class="ruby-keyword">do</span>
356
+ <span class="ruby-identifier">z</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>] <span class="ruby-operator">-=</span> <span class="ruby-ivar">@L</span>[<span class="ruby-identifier">i</span>,<span class="ruby-identifier">j</span>] <span class="ruby-operator">*</span> (<span class="ruby-identifier">z</span>[<span class="ruby-identifier">j</span>,<span class="ruby-value">0</span>] <span class="ruby-operator">||</span> <span class="ruby-value">1</span>)
357
+ <span class="ruby-keyword">end</span>
358
+ <span class="ruby-keyword">end</span>
359
+
360
+ (<span class="ruby-identifier">sol</span>.<span class="ruby-identifier">getRws</span> <span class="ruby-operator">-</span> <span class="ruby-value">1</span>).<span class="ruby-identifier">downto</span>(<span class="ruby-value">0</span>) <span class="ruby-keyword">do</span> <span class="ruby-operator">|</span><span class="ruby-identifier">i</span><span class="ruby-operator">|</span>
361
+ <span class="ruby-identifier">x</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>] = <span class="ruby-identifier">z</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>]
362
+ <span class="ruby-keyword">for</span> <span class="ruby-identifier">j</span> <span class="ruby-keyword">in</span> <span class="ruby-identifier">i</span><span class="ruby-operator">...</span><span class="ruby-identifier">sol</span>.<span class="ruby-identifier">getRws</span>
363
+ <span class="ruby-identifier">x</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>] <span class="ruby-operator">-=</span> (<span class="ruby-ivar">@U</span>[<span class="ruby-identifier">i</span>,<span class="ruby-identifier">j</span><span class="ruby-operator">+</span><span class="ruby-value">1</span>] <span class="ruby-operator">||</span> <span class="ruby-value">0</span>) <span class="ruby-operator">*</span> (<span class="ruby-identifier">x</span>[<span class="ruby-identifier">j</span><span class="ruby-operator">+</span><span class="ruby-value">1</span>,<span class="ruby-value">0</span>] <span class="ruby-operator">||</span> <span class="ruby-value">0</span>)
364
+ <span class="ruby-keyword">end</span>
365
+ <span class="ruby-identifier">x</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>] = <span class="ruby-identifier">x</span>[<span class="ruby-identifier">i</span>,<span class="ruby-value">0</span>] <span class="ruby-operator">/</span> <span class="ruby-ivar">@U</span>[<span class="ruby-identifier">i</span>,<span class="ruby-identifier">i</span>].<span class="ruby-identifier">to_f</span>
366
+ <span class="ruby-keyword">end</span>
367
+ <span class="ruby-keyword">return</span> <span class="ruby-identifier">x</span>
368
+ <span class="ruby-keyword">end</span></pre>
369
+ </div>
370
+
371
+ </div>
372
+
373
+
374
+
375
+
376
+ </div>
377
+
378
+
379
+ <div id="method-c-swap" class="method-detail ">
380
+
381
+ <div class="method-heading">
382
+ <span class="method-name">swap</span><span
383
+ class="method-args">(mx,r1,r2)</span>
384
+
385
+ <span class="method-click-advice">click to toggle source</span>
386
+
387
+ </div>
388
+
389
+
390
+ <div class="method-description">
391
+
392
+ <p>Swaps two rows of a matrix (pivoting)</p>
393
+ <ul><li>
394
+ <p>*<strong>argument</strong>*: matrix the rows must be swapped on</p>
395
+ </li><li>
396
+ <p>*<strong>argument</strong>*: first row</p>
397
+ </li><li>
398
+ <p>*<strong>argument</strong>*: second row</p>
399
+ </li><li>
400
+ <p>*<strong>returns</strong>*: new matrix with swapped rows</p>
401
+ </li></ul>
402
+
403
+
404
+
405
+
406
+ <div class="method-source-code" id="swap-source">
407
+ <pre><span class="ruby-comment"># File lib/linmeric/LU.rb, line 19</span>
408
+ <span class="ruby-keyword">def</span> <span class="ruby-keyword">self</span>.<span class="ruby-identifier">swap</span>(<span class="ruby-identifier">mx</span>,<span class="ruby-identifier">r1</span>,<span class="ruby-identifier">r2</span>)
409
+ <span class="ruby-keyword">for</span> <span class="ruby-identifier">i</span> <span class="ruby-keyword">in</span> <span class="ruby-value">0</span><span class="ruby-operator">...</span><span class="ruby-identifier">mx</span>.<span class="ruby-identifier">getCls</span> <span class="ruby-keyword">do</span>
410
+ <span class="ruby-identifier">mx</span>[<span class="ruby-identifier">r1</span>,<span class="ruby-identifier">i</span>],<span class="ruby-identifier">mx</span>[<span class="ruby-identifier">r2</span>,<span class="ruby-identifier">i</span>] = <span class="ruby-identifier">mx</span>[<span class="ruby-identifier">r2</span>,<span class="ruby-identifier">i</span>],<span class="ruby-identifier">mx</span>[<span class="ruby-identifier">r1</span>,<span class="ruby-identifier">i</span>]
411
+ <span class="ruby-keyword">end</span>
412
+ <span class="ruby-keyword">return</span> <span class="ruby-identifier">mx</span>
413
+ <span class="ruby-keyword">end</span></pre>
414
+ </div>
415
+
416
+ </div>
417
+
418
+
419
+
420
+
421
+ </div>
422
+
423
+
424
+ </section>
425
+
426
+ </section>
427
+ </main>
428
+
429
+
430
+ <footer id="validator-badges" role="contentinfo">
431
+ <p><a href="http://validator.w3.org/check/referer">Validate</a>
432
+ <p>Generated by <a href="http://docs.seattlerb.org/rdoc/">RDoc</a> 4.2.0.
433
+ <p>Based on <a href="http://deveiate.org/projects/Darkfish-RDoc/">Darkfish</a> by <a href="http://deveiate.org">Michael Granger</a>.
434
+ </footer>
435
+