adaptive_precision_utilities.hpp

Go to the documentation of this file.

//This file is part of Bertini 2.
//
//include/bertini2/trackers/adaptive_precision_utilities.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.
//
//include/bertini2/trackers/adaptive_precision_utilities.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 include/bertini2/trackers/adaptive_precision_utilities.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/trackers/config.hpp"
 
 
namespace bertini{ namespace tracking { namespace adaptive{

inline
void SetPrecision(SampCont<complex_mp> & samples, unsigned prec)
{
    for (auto& s : samples)
        for (unsigned ii=0; ii<s.size(); ii++)
            s(ii).precision(prec);
}

inline
void SetPrecision(TimeCont<complex_mp> & times, unsigned prec)
{
    for (auto& t : times)
        t.precision(prec);
}

inline
unsigned MaxPrecision(SampCont<complex_mp> const& samples)
{
    unsigned max_precision = 0;
    for (const auto& s : samples)
        if(Precision(s(0)) > max_precision)
            max_precision = Precision(s(0));
    return max_precision;
}

inline
unsigned MaxPrecision(TimeCont<complex_mp> const& times)
{
    unsigned max_precision = 0;
    for (const auto& t : times)
        if(Precision(t) > max_precision)
            max_precision = Precision(t);
    return max_precision;
}
 

inline
unsigned EnsureAtUniformPrecision(TimeCont<complex_dbl> & /*times*/, SampCont<complex_dbl> & /*samples*/)
{
    return DoublePrecision();
}
 

inline
unsigned EnsureAtUniformPrecision(TimeCont<complex_mp> & times, SampCont<complex_mp> & samples)
{
    auto def_prec = ThreadPrecision();
    if (std::any_of(begin(times),end(times),[=](auto const& p){return Precision(p)!=def_prec;})
        ||
        std::any_of(begin(samples),end(samples),[=](auto const& p){return Precision(p)!=def_prec;}))
    {
        auto max_precision = max(MaxPrecision(samples), MaxPrecision(times));
 
        SetThreadPrecision(max_precision);
        SetPrecision(times, max_precision);
        SetPrecision(samples, max_precision);
        return max_precision;
    }
    return def_prec;
 
}
 

inline
unsigned EnsureAtUniformPrecision(TimeCont<complex_mp> & times, SampCont<complex_mp> & samples, SampCont<complex_mp> & derivatives)
{
    auto def_prec = ThreadPrecision();
    if (std::any_of(begin(samples),end(samples),[=](auto const& p){return Precision(p)!=def_prec;})
        ||
        std::any_of(begin(derivatives),end(derivatives),[=](auto const& p){return Precision(p)!=def_prec;}))
    {
        auto max_precision = max(MaxPrecision(samples),MaxPrecision(times),MaxPrecision(derivatives));
 
        SetThreadPrecision(max_precision);
 
        SetPrecision(times, max_precision);
        SetPrecision(samples, max_precision);
        SetPrecision(derivatives, max_precision);
        return max_precision;
    }
    return def_prec;
}
 
 
 
}}} // re: namespaces