https://mediawiki.zeropage.org/index.php?action=history&feed=atom&title=CubicSpline%2F1002%2Ftest_lu.py 2026-05-15T20:58:32Z Revision history for this page on the wiki MediaWiki 1.39.8 https://mediawiki.zeropage.org/index.php?title=CubicSpline/1002/test_lu.py&diff=31027&oldid=prev 2021-02-07T05:23:03Z

<p></p> <p><b>New page</b></p><div> import unittest<br /> from Numeric import *<br /> from Matrix import *<br /> from LuDecomposition import *<br /> <br /> class TestLuDecomposition(unittest.TestCase):<br /> def setUp(self):<br /> self.a = [[1.0,2.0,3.0],[4.0,5.0,6.0],[7.0,8.0,9.0]]<br /> self.matrixA = Matrix(array(self.a))<br /> <br /> def testLwrite(self):<br /> ld = LuDecomposition(self.a)<br /> col = 0 <br /> ld._writeL(col)<br /> self.assertEquals(ld._getL(), [[1,0,0],[4,0,0],[7,0,0]])<br /> <br /> def testUwrite(self):<br /> ld = LuDecomposition(self.a)<br /> col = 0<br /> ld._writeL(col)<br /> row = 0<br /> ld._writeU(row)<br /> self.assertEquals(ld._getU(), [[1,2.0,3.0],[0,1,0], [0,0,1]])<br /> <br /> def testLuDecMore(self):<br /> ld = LuDecomposition(self.a)<br /> ld._writeL(0)<br /> ld._writeU(0)<br /> self.assertEquals(ld._getL(), [[1,0,0],[4,0,0],[7,0,0]])<br /> self.assertEquals(ld._getU(), [[1,2,3],[0,1,0], [0,0,1]])<br /> ld._writeL(1)<br /> ld._writeU(1)<br /> self.assertEquals(ld._getL(), [[1,0,0],[4,-3,0],[7,-6,0]])<br /> self.assertEquals(ld._getU(), [[1,2,3],[0,1,2], [0,0,1]])<br /> ld._writeL(2)<br /> ld._writeU(2)<br /> self.assertEquals(ld._getL(), [[1,0,0],[4,-3,0],[7,-6,0]])<br /> self.assertEquals(ld._getU(), [[1,2,3],[0,1,2], [0,0,1]])<br /> <br /> def testLuDecomposition(self):<br /> l, u = LuDecomposition(self.a).perform()<br /> matrixL = Matrix(array(l))<br /> matrixU = Matrix(array(u))<br /> actual = matrixL * matrixU<br /> self.assertEquals(str(actual), str(self.matrixA))<br /> <br /> <br /> if __name__==&quot;__main__&quot;:<br /> unittest.main()</div>

imported>Unknown