240
edits
Interval package: Difference between revisions
Jump to navigation
Jump to search
→Matrix operations: added dot and matrix multiplication with less accuracy and fixed accuracy statements
(→Matrix operations: Again, did some correctness and accuracy changes to the code) |
(→Matrix operations: added dot and matrix multiplication with less accuracy and fixed accuracy statements) |
||
| Line 236: | Line 236: | ||
Above mentioned operations can also be applied element-wise to interval vectors and matrices. Many operations use [http://www.gnu.org/software/octave/doc/interpreter/Vectorization-and-Faster-Code-Execution.html#Vectorization-and-Faster-Code-Execution vectorization techniques]. | Above mentioned operations can also be applied element-wise to interval vectors and matrices. Many operations use [http://www.gnu.org/software/octave/doc/interpreter/Vectorization-and-Faster-Code-Execution.html#Vectorization-and-Faster-Code-Execution vectorization techniques]. | ||
In addition, there are matrix operations on interval matrices. These operations comprise: | In addition, there are matrix operations on interval matrices. These operations comprise: tight dot product, fast dot product, tight matrix multiplication, fast matrix multiplication, tight vector sums, matrix inversion, matrix powers, and solving linear systems. As a result of missing hardware / low-level library support and missing optimizations, these operations are quite slow compared to familiar operations in floating-point arithmetic. | ||
octave:1> A = infsup ([1, 2, 3; 4, 0, 0; 0, 0, 1]); A (2, 3) = "[0, 6]" | octave:1> A = infsup ([1, 2, 3; 4, 0, 0; 0, 0, 1]); A (2, 3) = "[0, 6]" | ||
A = 3×3 interval matrix | A = 3×3 interval matrix | ||