示例#1
0
 def test_lu(self):
     tm,tmr,tma=self.tm,self.tmr,self.tma
     print 'checking for LU decomposition.'
     t0=time.time()
     lusys=lin.trilu(tm)
     L,U=lusys.L,lusys.U
     t1=time.time()
     L1,U=L.tocsr(),U.tocsc()
     lures=L1.dot(U)
     diff=abs(lures-tmr).toarray()
     print 'err -> %s, Elapse -> %s'%(diff.sum(),t1-t0)
     assert_allclose(diff,0,atol=1e-10)
示例#2
0
 def test_lu(self):
     tm, tmr, tma = self.tm, self.tmr, self.tma
     print 'checking for LU decomposition.'
     t0 = time.time()
     lusys = lin.trilu(tm)
     L, U = lusys.L, lusys.U
     t1 = time.time()
     L1, U = L.tocsr(), U.tocsc()
     lures = L1.dot(U)
     diff = abs(lures - tmr).toarray()
     print 'err -> %s, Elapse -> %s' % (diff.sum(), t1 - t0)
     assert_allclose(diff, 0, atol=1e-10)
示例#3
0
 def test_ul(self):
     print 'checking for UL decomposition.'
     tm,tmr,tma=self.tm,self.tmr,self.tma
     lusys=lin.trilu(tm)
     t0=time.time()
     ulsys=lin.triul(tm)
     U,L=lusys.L,lusys.U
     t1=time.time()
     U2,L=U.tocsr(),L.tocsc()
     ulres=U2.dot(L)
     diff=abs(ulres-tmr).toarray()
     print 'err -> %s, Elapse -> %s'%(diff.sum(),t1-t0)
     assert_allclose(diff,0,atol=1e-10)
示例#4
0
 def test_ul(self):
     print 'checking for UL decomposition.'
     tm, tmr, tma = self.tm, self.tmr, self.tma
     lusys = lin.trilu(tm)
     t0 = time.time()
     ulsys = lin.triul(tm)
     U, L = lusys.L, lusys.U
     t1 = time.time()
     U2, L = U.tocsr(), L.tocsc()
     ulres = U2.dot(L)
     diff = abs(ulres - tmr).toarray()
     print 'err -> %s, Elapse -> %s' % (diff.sum(), t1 - t0)
     assert_allclose(diff, 0, atol=1e-10)