narumij https://blog.hatena.ne.jp/narumij/ ちょっとしたメモなど https://narumij.hatenadiary.org/ ニューメリカルレシピinCを参考に書いてみました。 ublasにピボット無しのLU分解のメソッドがあることに気付いたので、 また逆行列やってみるときはそこら辺かな。※何度か修正してます。 (in-package #:cl) (defpackage #:lu (:use :cl) (:export #:m-inv)) (in-package #:lu) (defmethod matrix-at ((m array) i j) (aref m i j)) (defmethod (setf matrix-at) (value (m array) i j) (setf (aref m i j) val… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fnarumij.hatenadiary.org%2Fentry%2F20100408%2F1270715432" title="逆行列再び - ちょっとしたメモなど" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> Hatena Blog https://hatena.blog 2010-04-08 17:30:32 rich https://narumij.hatenadiary.org/entry/20100408/1270715432 1.0 100%