[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[no subject]



In System 78.34, ZMail 38.5, microcode 836, 60.45hz, on Lisp Machine Twenty:

I may be using math:decompose solve etc wrong but:

base
10.

ibase
10.

(setq a (make-array '(2 2)))
#<ART-Q-2.-2. 23376503>

(setq b (make-array 2))
#<ART-Q-2. 23376544>

(math:fill-2d-array a '((102 -52) (200 -2)))
NIL

(aset 5777 b 0)
5777.

(aset 19901 b 1)
19901.

(multiple-value (lu ps) (math:decompose a))
#<ART-Q-2.-2. 23376723>

(setq x (math:solve lu ps b))
#<ART-Q-2. 23376767>

(setq wrong-b (math:multiply-matrices a x))
#<ART-Q-2. 23377237>

(aref wrong-b 0)
5777.0s0

(aref wrong-b 1)
39564.5s0	;<<<<<<<this clearly should be 19901.