polynomial_utilities.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Jordi Cotela
//
 
#pragma once
 
// System includes
#include <array>
#include <vector>
 
// External includes
 
// Project includes
#include "includes/kratos_export_api.h"
 
namespace Kratos::PolynomialUtilities
{
    using PolynomialType = std::vector<double>;
    using IntervalType = std::array<double, 2>;
 
    std::size_t KRATOS_API(KRATOS_CORE) Degree(const PolynomialType& rPolynomial);
 
    double KRATOS_API(KRATOS_CORE) Evaluate(const PolynomialType& rPolynomial, double x);
 
    PolynomialType KRATOS_API(KRATOS_CORE) Differentiate(const PolynomialType& rPolynomial);
 
    PolynomialType KRATOS_API(KRATOS_CORE) Multiply(const PolynomialType& rA, const PolynomialType& rB);
 
    void KRATOS_API(KRATOS_CORE) Divide(
        PolynomialType& rQuotient,
        PolynomialType& rRemainder,
        const PolynomialType& rA,
        const PolynomialType& rB
        );
 
    void KRATOS_API(KRATOS_CORE) IsolateRoots(
        std::vector<IntervalType>& rRootIntervals,
        const PolynomialType& rPolynomial,
        const IntervalType& rRange
        );
 
    double KRATOS_API(KRATOS_CORE) FindRoot(
        const PolynomialType& rPolynomial,
        const IntervalType& rRange
        );
 
}