Python augmented_matrix Examples

Example #1

def solve(a, b):
    print "Coefficients:\n", a
    print "Values:\n", b
    p = mo.augmented_matrix(a, b)
    print 'Augmented:\n', p
    p = ge.row_scaling(p)
    print 'Row scaled:\n', p
    # p = partial_pivoting(p)
    # print "Partially pivoted:\n", p
    p = transform_to_identity(p)
    print 'Transformed to identity form:\n', p
    q = ge.backward_substitution(p)
    print "On backward substitution:\n", q

Example #2

File: gauss_jordan_elimination.py Project: archphy/numerical-methods-python

def solve(a,b):
	print "Coefficients:\n", a
	print "Values:\n", b
	p = mo.augmented_matrix(a,b)
	print 'Augmented:\n', p
	p = ge.row_scaling(p)
	print 'Row scaled:\n', p
	# p = partial_pivoting(p)
	# print "Partially pivoted:\n", p
	p = transform_to_identity(p)
	print 'Transformed to identity form:\n', p
	q = ge.backward_substitution(p)
	print "On backward substitution:\n", q

Example #3

File: thomas-tridiagonal.py Project: archphy/numerical-methods-python

def solve(a,b):
	if check_tridiagonal(a):
		(l,u) = thomas_decomposition(a)
	print 'Lower triangular matrix:\n', l
	print 'Upper triangular matrix:\n', u
	y = dd.forward_substitution(l,b)
	print 'On forward substitution:\n', y
	s = ge.backward_substitution(mo.augmented_matrix(u,y))
	print 'On backward substitution:\n', s

# test
# A = array([[3,-1,0,0],[2,-3,2,0],[0,1,2,5],[0,0,1,-1]], dtype='float')
# B = mo.transpose(array([1,2,3,4]))
# solve(A,B)

Example #4

File: dolittle_decomposition.py Project: archphy/numerical-methods-python

def solve(a,b):
	(l,u) = ge.transform_to_row_echleon(a, True)
	print 'Lower triangular matrix:\n', l
	print 'Upper triangular matrix:\n', u
	y = forward_substitution(l,b)
	print 'On forward substitution:\n', y
	s = ge.backward_substitution(mo.augmented_matrix(u,y))
	print 'On backward substitution:\n', s

# test
# A = array([[4,1],[1,3]], dtype='float')
# B = mo.transpose(array([1,2]))
# solve(A,B)
# A = array([[2,3,-1],[4,4,-3],[-2,1,-1]])
# B = mo.transpose(array([5,3,-3]))
# solve(A,B)

Example #5

def solve(a, b):
    (l, u) = ge.transform_to_row_echleon(a, True)
    print 'Lower triangular matrix:\n', l
    print 'Upper triangular matrix:\n', u
    y = forward_substitution(l, b)
    print 'On forward substitution:\n', y
    s = ge.backward_substitution(mo.augmented_matrix(u, y))
    print 'On backward substitution:\n', s


# test
# A = array([[4,1],[1,3]], dtype='float')
# B = mo.transpose(array([1,2]))
# solve(A,B)
# A = array([[2,3,-1],[4,4,-3],[-2,1,-1]])
# B = mo.transpose(array([5,3,-3]))
# solve(A,B)

Example #6

File: gauss_elimination.py Project: archphy/numerical-methods-python

def solve(a,b):
	print "Coefficients:\n", a
	print "Values:\n", b
	p = mo.augmented_matrix(a,b)
	print 'Augmented:\n', p
	p = row_scaling(p)
	print 'Row scaled:\n', p
	# p = partial_pivoting(p)
	# print "Partially pivoted:\n", p
	p = transform_to_row_echleon(p)
	print 'Transformed to Row-echleon form:\n', p
	q = backward_substitution(p)
	print "On backward substitution:\n", q

# test
# A = array([[2,3,-1],[4,4,-3],[-2,1,-1]])
# B = mo.transpose(array([5,3,-3]))
# A = array([[4,1],[1,3]])
# B = mo.transpose(array([1,2]))
# solve(A,B)

Example #7

File: gauss_elimination.py Project: archphy/numerical-methods-python

def solve(a, b):
    print "Coefficients:\n", a
    print "Values:\n", b
    p = mo.augmented_matrix(a, b)
    print 'Augmented:\n', p
    p = row_scaling(p)
    print 'Row scaled:\n', p
    # p = partial_pivoting(p)
    # print "Partially pivoted:\n", p
    p = transform_to_row_echleon(p)
    print 'Transformed to Row-echleon form:\n', p
    q = backward_substitution(p)
    print "On backward substitution:\n", q


# test
# A = array([[2,3,-1],[4,4,-3],[-2,1,-1]])
# B = mo.transpose(array([5,3,-3]))
# A = array([[4,1],[1,3]])
# B = mo.transpose(array([1,2]))
# solve(A,B)