total_degree_binomial.hpp

Go to the documentation of this file.

//This file is part of Bertini 2.
//
//total_degree_binomial.hpp is free software: you can redistribute it and/or modify
//it under the terms of the GNU General Public License as published by
//the Free Software Foundation, either version 3 of the License, or
//(at your option) any later version.
//
//total_degree_binomial.hpp is distributed in the hope that it will be useful,
//but WITHOUT ANY WARRANTY; without even the implied warranty of
//MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//GNU General Public License for more details.
//
//You should have received a copy of the GNU General Public License
//along with total_degree_binomial.hpp.  If not, see <http://www.gnu.org/licenses/>.
//
// Copyright(C) Bertini2 Development Team
//
// See <http://www.gnu.org/licenses/> for a copy of the license,
// as well as COPYING.  Bertini2 is provided with permitted
// additional terms in the b2/licenses/ directory.

 
#pragma once
 
#include "bertini2/system/start_base.hpp"
#include "bertini2/system/start/utility.hpp"
 
namespace bertini
{
    namespace start_system{
 

        class TotalDegreeBinomial : public StartSystem
        {
        public:
            TotalDegreeBinomial() = default;
            virtual ~TotalDegreeBinomial() = default;

            TotalDegreeBinomial(System const& s);
 

            template<typename NumT>
            NumT RandomValue(size_t index) const
            {
                // A direct read of the literal constant --- no node-level evaluation.  The Complex
                // node stores a complex_mp at its (max) creation precision; convert to NumT.
                auto const& v = random_values_[index]->GetValue();
                if constexpr (std::is_same<NumT, complex_dbl>::value)
                    return complex_dbl(double(v.real()), double(v.imag()));
                else
                    return NumT(v);
            }
 

            std::vector<std::shared_ptr<node::Complex> > const& RandomValues()
            {
                return random_values_;
            }
 

            unsigned long long NumStartPoints() const override;

            TotalDegreeBinomial& operator*=(Nd const& n);
 
            TotalDegreeBinomial& operator+=(System const& sys) = delete;

            void SanityChecks(System const& s);
 
        private:

            void CopyDegrees(System const& s);

            void SeedRandomValues(int num_functions);

            void GenerateFunctions();
 

            Vec<complex_dbl> GenerateStartPoint(complex_dbl,unsigned long long index) const override;

            Vec<complex_mp> GenerateStartPoint(complex_mp,unsigned long long index) const override;
 
            std::vector<std::shared_ptr<node::Complex> > random_values_; 
            std::vector<unsigned long long> degrees_; 
 
 
            friend class boost::serialization::access;
 
            template <typename Archive>
            void serialize(Archive& ar, const unsigned /*version*/) {
                ar & boost::serialization::base_object<StartSystem>(*this);
                ar & random_values_;
                ar & degrees_;
            }
 
        };
    }
}