PolynomialCoefficient.cpp

Go to the documentation of this file.

namespace interpolation {
 
    std::ostream& operator<<(std::ostream& out, const PolynomialCoefficient& coeff) {
        std::vector<int> in_var = coeff.InVariables();
        for (size_t i = 0; i < coeff.InVariables().size(); ++i)
            out << in_var[i] << " ";
        out << "| " << coeff.OutVariable();
        out << " | " << coeff.Coefficient();
        return out;
    }
 
    void PolynomialCoefficient::SpaceTransform(
            std::vector<int> space_in, std::vector<int> space_out) {
        std::map<int, int> mapping;  // probably optimise this
        for (unsigned int i = 0; i < space_out.size(); i++)
            for (unsigned int j = 0; j < space_in.size(); j++)
                if (space_out[i] == space_in[j])
                    mapping[j] = i;  // mapping[space_in_index] returns space_out_index
 
        std::vector<int> in_variables(_inVarByVec.size());
        for (unsigned int con = 0; con < in_variables.size(); con++) {
            if (mapping.find(in_variables[con]) != mapping.end())
                in_variables[con] = mapping[in_variables[con]];
            else
                throw(GeneralOpalException(
                        "PolynomialVector::PolynomialCoefficient::SpaceTransform",
                        "Input variable not found in space transform"));
        }
 
        if (mapping.find(_outVar) != mapping.end())
            _outVar = mapping[_outVar];
        else
            throw(GeneralOpalException(
                    "PolynomialVector::PolynomialCoefficient::SpaceTransform",
                    "Output variable not found in space transform"));
        _inVarByVec = in_variables;
    }
}  // namespace interpolation