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
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
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)
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)
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)
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)
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)