Esempi in Python per setdiag

Esempio n. 1

File: test_linalg.py Progetto: yuyunli2/projects

def test_setdiag():
    """
    Testing setdiag.
    """

    A = np.random.rand(30, 20)
    x = np.random.rand(20)
    linalg.setdiag(A, x)
    assert np.sum(np.abs(x - np.diag(A))) < 1e-10
    A = np.random.rand(30, 20).T
    x = np.random.rand(20)
    linalg.setdiag(A, x)
    assert np.sum(np.abs(x - np.diag(A))) < 1e-10

Esempio n. 2

File: test_linalg.py Progetto: goelhardik/projects

def test_setdiag():
    """
    Testing setdiag.
    """

    A = np.random.rand(30,20)
    x = np.random.rand(20)
    linalg.setdiag(A,x)
    assert np.sum(np.abs(x-np.diag(A))) < 1e-10
    A = np.random.rand(30,20).T
    x = np.random.rand(20)
    linalg.setdiag(A,x)
    assert np.sum(np.abs(x-np.diag(A))) < 1e-10

Esempio n. 3

print "NumPy vs mathutils.linalg diff.:", np.sum(np.abs(C - np.dot(A, B)))

print "Testing getdiag"
A = np.random.rand(30, 20)
x = np.zeros((20))
linalg.getdiag(A, x)
print "NumPy vs mathutils.linalg diff.:", np.sum(np.abs(x - np.diag(A)))
A = np.random.rand(30, 20).T
x = np.zeros((20))
linalg.getdiag(A, x)
print "NumPy vs mathutils.linalg diff.:", np.sum(np.abs(x - np.diag(A)))

print "Testing setdiag"
A = np.random.rand(30, 20)
x = np.random.rand(20)
linalg.setdiag(A, x)
print "NumPy vs mathutils.linalg diff.:", np.sum(np.abs(x - np.diag(A)))
A = np.random.rand(30, 20).T
x = np.random.rand(20)
linalg.setdiag(A, x)
print "NumPy vs mathutils.linalg diff.:", np.sum(np.abs(x - np.diag(A)))

print "Testing multiple_row_accumulate"
X = np.random.rand(30, 20)
Y = np.zeros((10, 20))
rows = np.hstack([
    np.arange(10, dtype='i'),
    np.arange(10, dtype='i'),
    np.arange(10, dtype='i'),
    np.arange(10, dtype='i'),
    np.arange(10, dtype='i')

Esempio n. 4

File: mathutils.py Progetto: MultiPath/TMBP

print "NumPy vs mathutils.linalg diff.:",np.sum(np.abs(C - np.dot(A,B)))

print "Testing getdiag"
A = np.random.rand(30,20)
x = np.zeros((20))
linalg.getdiag(A,x)
print "NumPy vs mathutils.linalg diff.:",np.sum(np.abs(x-np.diag(A)))
A = np.random.rand(30,20).T
x = np.zeros((20))
linalg.getdiag(A,x)
print "NumPy vs mathutils.linalg diff.:",np.sum(np.abs(x-np.diag(A)))

print "Testing setdiag"
A = np.random.rand(30,20)
x = np.random.rand(20)
linalg.setdiag(A,x)
print "NumPy vs mathutils.linalg diff.:",np.sum(np.abs(x-np.diag(A)))
A = np.random.rand(30,20).T
x = np.random.rand(20)
linalg.setdiag(A,x)
print "NumPy vs mathutils.linalg diff.:",np.sum(np.abs(x-np.diag(A)))

print "Testing multiple_row_accumulate"
X = np.random.rand(30,20)
Y = np.zeros((10,20))
rows = np.hstack([np.arange(10,dtype='i'),np.arange(10,dtype='i'),np.arange(10,dtype='i'),np.arange(10,dtype='i'),np.arange(10,dtype='i')])
np.random.shuffle(rows)
rows = rows[:30]
linalg.multiple_row_accumulate(X,rows,Y)
Y2 = np.zeros((10,20))
for i,r in enumerate(rows):