Method Solve

kratos/utilities/math_utils.cpp:28–41

template<class TDataType>
void MathUtils<TDataType>::Solve(
    MatrixType A,
    VectorType& rX,
    const VectorType& rB
    )
{
    const SizeType size1 = A.size1();
    rX = rB;
    typedef permutation_matrix<SizeType> pmatrix;
    pmatrix pm(size1);
    int singular = lu_factorize(A,pm);
    KRATOS_DEBUG_ERROR_IF(singular == 1) << "Matrix is singular: " << A << std::endl;
    lu_substitute(A, pm, rX);
}

/// Explicit instantation
template class MathUtils<double>;