;ELC ;;; compiled by amu@mit.edu on Wed Dec 18 01:19:04 1996 ;;; from file /afs/sipb.mit.edu/contrib/emacs/packages/calc-2.02e/calc-mat.el ;;; emacs version 19.30.1. ;;; bytecomp version FSF 2.10 ;;; optimization is on. ;;; this file uses opcodes which do not exist in Emacs 18. (if (and (boundp 'emacs-version) (or (and (boundp 'epoch::version) epoch::version) (string-lessp emacs-version "19.29"))) (error "`calc-mat.el' was compiled for Emacs 19.29 or later")) (byte-code "!!" [require calc-ext calc-macs] 2) (defalias 'calc-Need-calc-mat #[nil "" [nil] 1]) (defalias 'calc-mdet #[(arg) "\"" [calc-do #[nil " #" [calc-unary-op "mdet" calcFunc-det arg] 4] 1247] 3 nil "P"]) (defalias 'calc-mtrace #[(arg) "\"" [calc-do #[nil " #" [calc-unary-op "mtr" calcFunc-tr arg] 4] 1359] 3 nil "P"]) (defalias 'calc-mlud #[(arg) "\"" [calc-do #[nil " #" [calc-unary-op "mlud" calcFunc-lud arg] 4] 1471] 3 nil "P"]) (defalias 'math-row-matrix #[(a) ":@=!D" [a vec math-matrixp] 2]) (defalias 'math-col-matrix #[(a) ":@=!A\"B" [a vec math-matrixp mapcar #[(x) " D" [vec x] 2]] 4]) (defalias 'math-mul-mats #[(a b) " A@G A k SV_ @A A\f@ @8\"\fA AV\n\f@ @8\"\"=\nB#BB.B" [nil b accum bp ap col row cols mat a 0 math-mul math-add vec] 7]) (defalias 'math-mul-mat-vec #[(a b) " A\"B" [vec mapcar #[(row) " \n\"" [math-dot-product row b] 3] a] 4]) (defalias 'calcFunc-tr #[(mat) " ! GS A@A@$ \"" [math-square-matrixp mat math-matrix-trace-step 2 math-reject-arg square-matrixp] 5]) (defalias 'math-matrix-trace-step #[(n size mat sum) " XT  88\"$ " [n size math-matrix-trace-step mat math-add sum] 9]) (defalias 'math-matrix-inv-raw #[(m) "GSXE!\f!?A U ǂ? UGɉȉ88A@8!E8A@!A@A@EE? U?ɉ‰88ȉ88\"8888\"\"A@888\"‰88A@8\"\"88A@8\"A@8ȉ88\"\"F888A@\"‰888A@\"\"‰88A@A@\"A@88A@\"\"A@88A@\"88A@A@\"\"F888A@\"ȉ888A@\"\"A@88A@\"88A@A@\"\"ȉ88A@A@\"A@88A@\"\"FF\f\")W!V \"\"))" [m n 3 math-det-raw det math-zerop math-div 1 2 vec math-neg math-sub math-mul math-matrix-lud lud math-lud-solve calcFunc-idn] 14]) (defalias 'calcFunc-det #[(m) " ! \\ !)! =- A@!) A@\"- A@ \"" [math-square-matrixp m math-normalize calc-internal-prec 2 math-det-raw calcFunc-idn math-zerop math-equal-int 1 math-reject-arg square-matrixp] 3]) (defalias 'math-det-raw #[(m) "GSUA@A@ U5A@A@É88\"A@88A@\"\" UĉljA@A@É88Ɖ88\"\"A@8888A@\"\"\"A@88A@88\"\"\"A@8É888A@\"\"\"A@A@8888\"\"\"A@88A@Ɖ88\"\"\"!  @\n  8\")))" [m n 1 2 math-sub math-mul 3 math-add math-matrix-lud lud lu math-det-step 0] 13]) (defalias 'math-det-step #[(n prod) "VS\f 88\"\"\f" [n 0 math-det-step math-mul prod lu] 7]) (defalias 'math-matrix-lud #[(m) " \n\" \fDA@ AAAȍ  B\n7\n> \nB\nB **" [assoc m math-lud-cache calc-internal-prec calc-prefer-frac context old singular (math-do-matrix-lud m) lud entry] 4]) (byte-code "! B‡" [boundp math-lud-cache nil current-load-list] 2) (defalias 'math-do-matrix-lud #[(m) " !GS \n \f  XV  W=G #\" \n88\n Wv\n \n88\n88\"\"\nTS \n8\n T0 X= #\" \n88\nW\n \n88\n88\"\"\nT \n8\n\n! \"U  ) T V \n #\f[\f  B \n88!+\"L T XL \n8 \n88\".)T!\n \fE.\n" [math-copy-matrix m lu n nil i 1 j k imax sum big d index 0 calc-display-working-message lots math-do-working "LUD step" format "%d/%d" math-sub math-mul math-abs-approx dum math-compare -1 math-swap-rows pivot math-zerop throw singular math-div] 8]) (defalias 'math-swap-rows #[(m r1 r2) " U0S\nA S\nAA \f \fA\f-\n" [r1 r2 m r1prev row1 r2prev row2 r2next] 3]) (defalias 'math-lud-solve #[(lud b &optional need) "'\n!GS A@GS@\n \f  X#=?\"\n A@\n\fX@\fA\f 88\f 8\n 88 !\n   \nW \n88 88\"\" T \n 8\nT\nL\nS\nY\n 88\n  T \fX \n88 88\"\"\n 8\n88\"T- . 1\"" [lud math-copy-matrix b x n m lu 1 col nil i j ip ii index sum calc-display-working-message lots math-do-working "LUD solver step" math-zerop math-sub math-mul math-div need math-reject-arg "*Singular matrix"] 8]) (defalias 'calcFunc-lud #[(m) " ! \\ !@!@! @GS\n\n\" \nXW[8ѠTA8̠T\nX 8ѠeT7SYSA@8U # )  F.*! \" \"" [math-square-matrixp m math-normalize calc-internal-prec 2 math-matrix-lud lud math-copy-matrix lmat umat n calcFunc-idn 1 perm nil i j 0 pos math-swap-rows vec math-reject-arg "*Singular matrix" square-matrixp] 6])