87ed8700e2
* first fix to start branch PR * explicit conversion from sizet to int * Array.h and DiagonalMatrix.h * many sizet to int conversions * removed some unused variables and added expl overrides * removed many unused parameters * more typing * even more * last of the easy changes
123 lines
3.6 KiB
C++
123 lines
3.6 KiB
C++
/***************************************************************************
|
|
* Copyright (c) 2023 Ondsel, Inc. *
|
|
* *
|
|
* This file is part of OndselSolver. *
|
|
* *
|
|
* See LICENSE file for details about copyright. *
|
|
***************************************************************************/
|
|
|
|
#include <algorithm>
|
|
|
|
#include "Polynomial.h"
|
|
#include "Constant.h"
|
|
#include "Sum.h"
|
|
|
|
using namespace MbD;
|
|
|
|
MbD::Polynomial::Polynomial(Symsptr var, std::shared_ptr<std::vector<double>> coefficients)
|
|
{
|
|
assert(!coefficients->empty());
|
|
xx = var;
|
|
std::transform(coefficients->begin(), coefficients->end(), coeffs->begin(),
|
|
[&](auto& coeff) { return sptrConstant(coeff); }
|
|
);
|
|
}
|
|
MbD::Polynomial::Polynomial(Symsptr var, std::shared_ptr<std::vector<Symsptr>> coefficients)
|
|
{
|
|
assert(!coefficients->empty());
|
|
xx = var;
|
|
coeffs->insert(coeffs->end(), coefficients->begin(), coefficients->end());
|
|
}
|
|
|
|
Symsptr MbD::Polynomial::expandUntil(Symsptr, std::shared_ptr<std::unordered_set<Symsptr>> set)
|
|
{
|
|
auto newCoeffs = std::make_shared<std::vector<Symsptr>>();
|
|
for (int i = 0; i < (int)coeffs->size(); i++)
|
|
{
|
|
auto coeff = coeffs->at(i);
|
|
auto newCoeff = coeff->expandUntil(coeff, set);
|
|
newCoeffs->push_back(newCoeff);
|
|
}
|
|
return std::make_shared<Polynomial>(xx, newCoeffs);
|
|
}
|
|
|
|
Symsptr MbD::Polynomial::simplifyUntil(Symsptr, std::shared_ptr<std::unordered_set<Symsptr>> set)
|
|
{
|
|
auto newCoeffs = std::make_shared<std::vector<Symsptr>>();
|
|
for (int i = 0; i < (int)coeffs->size(); i++)
|
|
{
|
|
auto coeff = coeffs->at(i);
|
|
auto newCoeff = coeff->simplifyUntil(coeff, set);
|
|
newCoeffs->push_back(newCoeff);
|
|
}
|
|
return std::make_shared<Polynomial>(xx, newCoeffs);
|
|
}
|
|
|
|
Symsptr MbD::Polynomial::differentiateWRTx()
|
|
{
|
|
//Differentiate powers
|
|
if (coeffs->size() == 1) return sptrConstant(0.0);
|
|
auto newCoeffs = std::make_shared<std::vector<Symsptr>>();
|
|
for (int i = 1; i < (int)coeffs->size(); i++)
|
|
{
|
|
auto newCoeff = i * coeffs->at(i)->getValue();
|
|
newCoeffs->push_back(sptrConstant(newCoeff));
|
|
}
|
|
auto poly1 = std::make_shared<Polynomial>(xx, newCoeffs);
|
|
//Differentiate coeffs
|
|
auto coeffDerivs = std::make_shared<std::vector<Symsptr>>();
|
|
std::transform(coeffs->begin(),
|
|
coeffs->end(),
|
|
std::back_inserter(*coeffDerivs),
|
|
[&](auto& coeff) { return coeff->differentiateWRT(xx); }
|
|
);
|
|
auto poly2 = std::make_shared<Polynomial>(xx, coeffDerivs);
|
|
return std::make_shared<Sum>(poly1, poly2);
|
|
}
|
|
|
|
Symsptr MbD::Polynomial::integrateWRT(Symsptr var)
|
|
{
|
|
assert(xx == var);
|
|
auto newCoeffs = std::make_shared<std::vector<Symsptr>>();
|
|
newCoeffs->push_back(sptrConstant(0.0));
|
|
for (int i = 0; i < (int)coeffs->size(); i++)
|
|
{
|
|
auto newCoeff = coeffs->at(i)->getValue() / (i + 1);
|
|
newCoeffs->push_back(sptrConstant(newCoeff));
|
|
}
|
|
return std::make_shared<Polynomial>(var, newCoeffs);
|
|
}
|
|
|
|
double MbD::Polynomial::getValue()
|
|
{
|
|
auto xvalue = xx->getValue();
|
|
auto xpower = 1.0;
|
|
auto answer = 0.0;
|
|
for (int i = 0; i < (int)coeffs->size(); i++)
|
|
{
|
|
answer += coeffs->at(i)->getValue() * xpower;
|
|
xpower *= xvalue;
|
|
}
|
|
return answer;
|
|
}
|
|
|
|
void MbD::Polynomial::integrationConstant(double integConstant)
|
|
{
|
|
auto coeff0 = coeffs->at(0);
|
|
coeff0->setValue(coeff0->getValue() + integConstant);
|
|
}
|
|
|
|
std::ostream& MbD::Polynomial::printOn(std::ostream& s) const
|
|
{
|
|
s << "Polynomial(";
|
|
s << *xx << ", ";
|
|
s << "coeffs{";
|
|
s << *coeffs->at(0);
|
|
for (int i = 1; i < (int)coeffs->size(); i++)
|
|
{
|
|
s << ", " << *coeffs->at(i);
|
|
}
|
|
s << "})";
|
|
return s;
|
|
}
|