Post Reply 
41CL Quiz: Determinant of 30x30 anti-Identity matrix
05-11-2018, 07:34 PM (This post was last modified: 05-11-2018 07:39 PM by Ángel Martin.)
Post: #3
RE: 41CL Quiz: Determinant of 30x30 anti-Identity matrix
(05-11-2018 06:06 PM)Valentin Albillo Wrote:  
(05-11-2018 09:10 AM)Ángel Martin Wrote:  This quiz is just for you, proud owners of a 41CL: -

Sorry, Ángel, you know I couldn't resist ...

Quote:Calculate the determinant of a 30x30 matrix with all elements equal to "1" except the diagonal, which has all zeros. I call this type "anti-Identity", but it's totally made up for you purists out there so don't nail me down just yet.

The exact result is Det = -29

My solution returns -28.999999977 in 11 minutes, give or take...

I know, I know, this is just for the 41CL but nevertheless here's the result for the HP-71B under J-F Garnier's Emu71. From the command line:

>DESTROY ALL @ OPTION BASE 1 @ DIM A(30,30),C(30,30)
>MAT A=CON @ MAT C=IDN @ MAT A=A-C @ DET(A)

      -29.0000000001

The timing is 0.15 seconds, give or take a hundredth ...

Have a nice weekend and best regards.
V.

Not a problem Valentín - I'm very glad you chimed in.
Ah the beauty of Saturn MCODE, I'm sooo jealous of such sheer power ;-)

The original goal of this exercise was to show that it's possible to deal with 900 element matrices using the expanded registers on the CL board. Obviously the bottleneck here is the FOCAL program that calculates the determinant. I'm toying with the idea of modifying the MCODE-based functions in the Advantage to also support the Y-registers, but it's a tall order...

Saludos y buen finde.
ÁM.

PS. I'm curious to know: does EMU-71 run faster than the real machine (being a PC application), or it also emulates the Saturn CPU speed?

"To live or die by your own sword one must first learn to wield it aptly."
Find all posts by this user
Quote this message in a reply
Post Reply 


Messages In This Thread
RE: 41CL Quiz: Determinant of 30x30 anti-Identity matrix - Ángel Martin - 05-11-2018 07:34 PM



User(s) browsing this thread: 2 Guest(s)