solve

AnsSolver.solve(b, x=None)

Solve a linear system.

Parameters:

bAnsVec

AnsMath vector.

xAnsVec, optional

AnsMath vector to place the solution into.

Returns:

AnsVec

Solution vector, which is identical to the x parameter if supplied.

Examples

>>> k = mm.stiff(fname='PRSMEMB.full')
>>> s = mm.factorize(k)
>>> b = mm.get_vec(fname='PRSMEMB.full', mat_id="RHS")
>>> x = s.solve(b)
>>> x
AnsMath vector size 20000