;ELC
;;; Compiled
;;; in Emacs version 28.2
;;; with all optimizations.



(byte-code "\300\301!\210\300\302!\207" [require calc-ext calc-macs] 2)
#@12 

(fn ARG)
(defalias 'calc-mdet #[257 "\300\301\302\"`\"\207" [calc-do make-closure #[0 "\301\302\303\300#\207" [V0 calc-unary-op "mdet" calcFunc-det] 4]] 5 (#$ . 156) "P"])
#@12 

(fn ARG)
(defalias 'calc-mtrace #[257 "\300\301\302\"`\"\207" [calc-do make-closure #[0 "\301\302\303\300#\207" [V0 calc-unary-op "mtr" calcFunc-tr] 4]] 5 (#$ . 337) "P"])
#@12 

(fn ARG)
(defalias 'calc-mlud #[257 "\300\301\302\"`\"\207" [calc-do make-closure #[0 "\301\302\303\300#\207" [V0 calc-unary-op "mlud" calcFunc-lud] 4]] 5 (#$ . 518) "P"])
#@10 

(fn A)
(defalias 'math-row-matrix #[257 "\211\211\242\300=\262
\203\301
!\204\300
D\207\207" [vec math-matrixp] 4 (#$ . 699)])
#@10 

(fn A)
(defalias 'math-col-matrix #[257 "\211\211\242\300=\262
\203\301
!\204\300\302\303A\"B\207\207" [vec math-matrixp mapcar #[257 "\300
D\207" [vec] 3 "\n\n(fn X)"]] 5 (#$ . 838)])
#@12 

(fn A B)
(defalias 'math-mul-mats #[514 "\300
A@G\300\211\211\211\211A\211\262\n\203i\262\300\262S\211\262\301V\203^@A\262A\262\302@@8\"\262
A\262
A\211\262\203V\303
\302@@8\"\"\262
\2028\211B\262\202\304BB\262\202\n\304\237B\207" [nil 0 math-mul math-add vec] 15 (#$ . 1036)])
#@12 

(fn A B)
(defalias 'math-mul-mat-vec #[514 "\300\301\302\303\"A\"B\207" [vec mapcar make-closure #[257 "\301
\300\"\207" [V0 math-dot-product] 4 "\n\n(fn ROW)"]] 7 (#$ . 1361)])
#@12 

(fn MAT)
(defalias 'calcFunc-tr #[257 "\300
!\203\301\302GS\211A@A@$\207\303
\304\"\207" [math-square-matrixp math-matrix-trace-step 2 math-reject-arg square-matrixp] 6 (#$ . 1549)])
#@23 

(fn N SIZE MAT SUM)
(defalias 'math-matrix-trace-step #[1028 "X\203\300T\301	\211	88\"$\207\207" [math-matrix-trace-step math-add] 13 (#$ . 1744)])
#@10 

(fn M)
(defalias 'math-matrix-inv-raw #[257 "\211GS\211\300X\203h
\301!\302
!?\205g
\303\304U\203\304\202e
\305U\203I\306\211\305\21188\307\305A@8!E\306\307\3058A@!A@A@EE\202e
\300U\205e
\306\211\310\311\300\211	88\305\211\n88\"\311\300\305\n88\305\30088\"\"\310\311\300	A@8\305\30088\"\311\300\21188\305A@8\"\"\310\311\300\30588\305A@8\"\311\300A@8\305\211
88\"\"F\306\310\311\300\305\n88\300\n8A@\"\311\300\21188\3058A@\"\"\310\311\300\21188\nA@A@\"\311\300A@8\300\f8A@\"\"\310\311\300A@8\305\f8A@\"\311\300\305
88\fA@A@\"\"F\306\310\311\305\30088\3058A@\"\311\305\211\f88\300\f8A@\"\"\310\311\305A@8\300\f8A@\"\311\305\300
88\fA@A@\"\"\310\311\305\211
88\fA@A@\"\311\305
A@8\3058A@\"\"FF\"\207\312!\211\205v
\313
\314\304\"\"\207" [3 math-det-raw math-zerop math-div 1 2 vec math-neg math-sub math-mul math-matrix-lud math-lud-solve calcFunc-idn] 16 (#$ . 1910)])
#@10 

(fn M)
(defalias 'calcFunc-det #[257 "\301
!\203\302\303\\\304!)!\207\211\242\305=\203-\306
A@!\204)\307
A@\310\"\203-\211A@\207\311
\312\"\207" [calc-internal-prec math-square-matrixp math-normalize 2 math-det-raw calcFunc-idn math-zerop math-equal-int 1 math-reject-arg square-matrixp] 4 (#$ . 2855)])
#@10 

(fn M)
(defalias 'math-det-raw #[257 "\211GS\211\301U\203
A@A@\207\211\302U\2032\303\304A@A@\302\21188\"\304\302A@8\3028A@\"\"\207\211\305U\203\305\303\211\211\306\211\304A@A@\304\302\21188\305\211\f88\"\"\304\302	A@8\304\305\302\f88\305\f8A@\"\"\"\304\305A@8\304\302\n8A@\302\305\f88\"\"\"\304\305A@8\304\302\211\n88\305\n8A@\"\"\"\304A@A@\304\305\302	88\302\305\n88\"\"\"\304\302A@8\304\3028A@\305\211	88\"\"\"\207\307!\211\203\327\211@\310\3028\")\207\311\207" [math-det-lu 1 2 math-sub math-mul 3 math-add math-matrix-lud math-det-step 0] 14 (#$ . 3175)])
#@15 

(fn N PROD)
(defalias 'math-det-step #[514 "
\301V\203\302S\303\21188\"\"\207\207" [math-det-lu 0 math-det-step math-mul] 9 (#$ . 3781)])
(defvar math-lud-cache nil)
#@10 

(fn M)
(defalias 'math-matrix-lud #[257 "\303
\"	\nD
\203
A@
\232\203
AA\207\3042\305!0

B\203-
\241\210\2023
BB
\207" [math-lud-cache calc-internal-prec calc-prefer-frac assoc singular math-do-matrix-lud] 7 (#$ . 3961)])
#@98 Determine a useful value for checking the size of potential pivots
in LUD decomposition.

(fn A)
(defalias 'math-lud-pivot-check #[257 "\211\242\300=\203(\301
A@!\203&\301\3028!\203&\303
A@\3028\"\304=\203&\304\207\305\207\306
!\207" [mod math-integerp 2 math-gcd 1 0 math-abs-approx] 5 (#$ . 4206)])
#@10 

(fn M)
(defalias 'math-do-matrix-lud #[257 "\301
!\211GS\302\303\302\211\211\211\303\302	X\203[
\303\262\304\262\262W\203~\305=\2039\306\307\310\311\n\f#\"\21088\262\303\262W\203k\312\31388\n\n88\"\"\262T\262\202F8\233\240\210T\262\202 	X\203\372\305=\203\227\306\307\310\311\n\f#\"\21088\262\303\262W\203\311\312\31388\n\n88\"\"\262T\262\202\2448\233\240\210\314!\315!\204\352
\316\"\317U\266\202\203\361\211\262\262\210T\262\202~V\203
\320\n#\262\n
[\262
B\262
\21188\315
!\203*
\321\322\302\"\266\202S
\262	T\211\262\n\nX\203R
	\f8\233\323	88\"\240\210\202.
\210T\262\202	
\237E\207" [calc-display-working-message math-copy-matrix nil 1 0 lots math-do-working "LUD step" format "%d/%d" math-sub math-mul math-lud-pivot-check math-zerop math-compare -1 math-swap-rows throw singular math-div] 18 (#$ . 4519)])
#@16 

(fn M R1 R2)
(defalias 'math-swap-rows #[771 "

U\204&
S\233\211AS\233\211A\211A\241\210\241\210
A\241\210
\241\266\207" [] 10 (#$ . 5489)])
#@29 

(fn LUD B &optional NEED)
(defalias 'math-lud-solve #[770 "\203	
\301!\211GS
A@GS@\302\303\211\211\211\211\211	X\203
\304=\203,\305\306\"\210\302\262\303\262
A@\262\nX\203\253
@\262
A\262\f88\262
\f8\233
88\240\210\204q\307
!\204\231\262\202\231\262W\203\231\310
\311	\f88\n	88\"\"\262
T\262\202t\f8\233
\240\210T\262\2028S\211\262\302Y\203\376\f88\262
\262T\211\262\nX\203\347\310
\311	\f88\n	88\"\"\262
\202\302\f8\233\312\211\f88\"\240\210\202\253T\262\202\n\207\211\205
\313
\314\"\207" [calc-display-working-message math-copy-matrix 1 nil lots math-do-working "LUD solver step" math-zerop math-sub math-mul math-div math-reject-arg "*Singular matrix"] 21 (#$ . 5650)])
#@10 

(fn M)
(defalias 'calcFunc-lud #[257 "\301
!\203\240\302\303\\\304!\211\205\224\305
@!\305@!@GS\306\307\"\310\307\211X\203i\307\262

W\203C\2118\233\311\240\210
T\262\202-\211
8\233\307\240\210
T\211\262X\203b\2118\233\311\240\210\202L\211T\262
\202$\211S\211\262\307Y\203\213\211SA@8\211U\204\207\312#\262\210\202i\313F\266\206\262
)!\206\244\314
\315\"\207\314
\316\"\207" [calc-internal-prec math-square-matrixp math-normalize 2 math-matrix-lud math-copy-matrix calcFunc-idn 1 nil 0 math-swap-rows vec math-reject-arg "*Singular matrix" square-matrixp] 14 (#$ . 6436)])
(provide 'calc-mtx)