mhom.hpp

Go to the documentation of this file.

//This file is part of Bertini 2.
//
//mhom.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.
//
//mhom.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 mhom.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 MHomogeneous : public StartSystem
        {
        public:
            MHomogeneous() = default;
            virtual ~MHomogeneous() = default;

            MHomogeneous(System const& s);

            void CreateDegreeMatrix(System const& s);

            void GenerateValidPartitions(System const& s);

            int ChooseColumnInRow(System const& s,Vec<int>& variable_group_counter, int row, int column);

            unsigned long long NumStartPoints() const override;

            unsigned long long NumStartPointsForPartition(Vec<int> partition) const;

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

            Mat<int> degree_matrix_; // stores degrees of all functions in all homogeneous variable groups.

            std::deque< Vec<int> > valid_partitions_;
 
 
        private:

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

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

            template<typename T>
            void GenerateStartPointT(Vec<T>& start_point, unsigned long long index) const;
 
 
 
 
            std::vector< VariableGroup > var_groups_;
            Mat<Mat<complex_mp>> linear_coeffs_;
            std::vector< std::vector<size_t> > variable_cols_; 
            size_t num_hom_groups_ = 0; 
 
            mutable Vec<complex_mp> temp_v_mp_;
 
            friend class boost::serialization::access;
 
            template <typename Archive>
            void serialize(Archive& ar, const unsigned /*version*/)
            {
                // Serialize ALL persistent state.  This object is broadcast to MPI workers via
                // ZeroDimSolver::DistributeSystems (a polymorphic shared_ptr<StartSystem> archive);
                // dropping any member leaves a worker with a hollow start system whose
                // NumStartPoints() returns 0, sizing the per-path metadata to nothing and crashing
                // the parallel solve.  Order base-first so Boost object tracking re-links the
                // var_groups_ shared_ptr<Variable> entries to the same Variable nodes serialized in
                // the base System (variable_groups_/hom_variable_groups_), preserving node identity.
                ar & boost::serialization::base_object<StartSystem>(*this);
                ar & degree_matrix_;
                ar & valid_partitions_;
                ar & var_groups_;
                ar & linear_coeffs_;   // random factor coefficients: PRIMARY data, not recomputable
                ar & variable_cols_;
                ar & num_hom_groups_;
            }
 
        };
    }//end start_system namespace
}//end bertini namespace