Talk:Gauss-Jordan matrix inversion: Difference between revisions
Talk:Gauss-Jordan matrix inversion (view source)
Revision as of 19:58, 3 November 2018
, 5 years ago→what's wrong (with my REXX program)?: thanks. program fixed for the specal case
Thundergnat (talk | contribs) (Seems to work with the Perl 6 implementation.) |
Walterpachl (talk | contribs) (→what's wrong (with my REXX program)?: thanks. program fixed for the specal case) |
||
Line 145:
::::I ran it through the Perl 6 implementation and it works and is nominally correct. (If I invert the inverted matrix, I end up with the original.) I added it as another test case. See the last two tests under [[Gauss-Jordan_matrix_inversion#Perl_6]] --[[User:Thundergnat|Thundergnat]] ([[User talk:Thundergnat|talk]]) 18:29, 3 November 2018 (UTC)
:::::Thank you, that helped. I fixed my program for the case that a.2.2=0.
My inverse matrix is the same as Perl 6's
<pre>
show 6 The inverse matrix
-4525/6238 2529/6238 -233/3119 1481/3119 -639/6238
1033/6238 -1075/6238 342/3119 -447/3119 871/6238
1299/6238 -289/6238 -204/3119 -390/3119 739/6238
782/3119 -222/3119 237/3119 -556/3119 -177/3119
-474/3119 -17/3119 -24/3119 688/3119 -140/3119
</pre>
<br>
Generalization of the REXX algorithm for ANY pivot=0 must be worked out (tbd) --Walter Pachl 19:58, 3 November 2018 (UTC)
| |||