start_base.hpp

Go to the documentation of this file.

//This file is part of Bertini 2.
//
//start_base.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.
//
//start_base.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 start_base.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/system.hpp"
 
#include <functional>
#include <memory>
 
 
namespace bertini
{
    namespace start_system{
 
        // forward declare so the StartSystemFactory alias / MakeStartFactory below can be defined
        // here (they only need the StartSystem base, declared just below).
        class StartSystem;

        template<typename SystemType>
        using StartSystemFactory =
            std::function<std::shared_ptr<StartSystem>(SystemType const&)>;

        template<typename StartType, typename SystemType = bertini::System>
        StartSystemFactory<SystemType> MakeStartFactory()
        {
            return [](SystemType const& target) -> std::shared_ptr<StartSystem> {
                return std::make_shared<StartType>(target);
            };
        }

        class StartSystem : public System
        {
 
        public:
 

            virtual unsigned long long NumStartPoints() const = 0;

            template<typename T>
            Vec<T> StartPoint(unsigned long long index) const
            {
                return GenerateStartPoint(T(),index);
            }
 
            virtual ~StartSystem() = default;
 
        private:
            virtual Vec<complex_dbl> GenerateStartPoint(complex_dbl,unsigned long long index) const = 0;
            virtual Vec<complex_mp> GenerateStartPoint(complex_mp,unsigned long long index) const = 0;
 
            friend class boost::serialization::access;

            template <typename Archive>
            void serialize(Archive& ar, const unsigned /*version*/) {
                ar & boost::serialization::base_object<System>(*this);
            }
 
        };
 
    }
}