Python qrの例

コード例 #1

    def test_qr(self):
        M = 100
        N = 10
        Y = expr.randn(M, N, tile_hint=(M / 4, N)).force()

        Q1, R1 = qr(Y)
        Q1 = Q1.glom()
        Q2, R2 = linalg.qr(Y.glom(), mode="economic")

        assert np.allclose(np.absolute(Q1), np.absolute(Q2))
        assert np.allclose(np.absolute(R1), np.absolute(R2))

コード例 #2

ファイル: test_qr.py プロジェクト: MaggieQi/spartan

 def test_qr(self):
   M = 100
   N = 10
   Y = expr.randn(M, N, tile_hint=(M/4, N)).force()
   
   Q1, R1 = qr(Y)
   Q1 = Q1.glom()
   Q2, R2 = linalg.qr(Y.glom(), mode="economic")
   
   assert np.allclose(np.absolute(Q1), np.absolute(Q2))
   assert np.allclose(np.absolute(R1), np.absolute(R2))

コード例 #3

def benchmark_qr(ctx, timer):
    M = 1280
    N = 1280
    Y = np.random.randn(M, N)
    Y = expr.from_numpy(Y)
    #Y = expr.randn(M, N)

    t1 = datetime.now()
    Q, R = qr(Y)
    t2 = datetime.now()
    cost_time = millis(t1, t2)

    print "total cost time:%s ms" % (cost_time)

コード例 #4

ファイル: test_qr.py プロジェクト: MaggieQi/spartan

def benchmark_qr(ctx, timer):
  M = 1280
  N = 1280
  Y = np.random.randn(M, N)
  Y = expr.from_numpy(Y)
  #Y = expr.randn(M, N)

  t1 = datetime.now()
  Q, R = qr(Y)
  t2 = datetime.now()
  cost_time = millis(t1, t2)
    
  print "total cost time:%s ms" % (cost_time)