CompoundConvergenceCriterion.hpp

Go to the documentation of this file.

// This file is part of GQCG-GQCP.
//
// Copyright (C) 2017-2020  the GQCG developers
//
// GQCG-GQCP is free software: you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// GQCG-GQCP 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 Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with GQCG-GQCP.  If not, see <http://www.gnu.org/licenses/>.
 
#pragma once
 
 
#include "Mathematical/Algorithm/ConvergenceCriterion.hpp"
 
#include <boost/format.hpp>
 
#include <memory>
#include <vector>
 
 
namespace GQCP {
 
 
template <typename _Environment>
class CompoundConvergenceCriterion:
    public ConvergenceCriterion<_Environment> {
 
public:
    using Environment = _Environment;
    using Self = CompoundConvergenceCriterion<Environment>;
 
 
private:
    std::vector<std::shared_ptr<ConvergenceCriterion<Environment>>> criteria;
 
 
public:
    /*
     *  CONSTRUCTORS
     */
 
    template <typename C1, typename C2>
    CompoundConvergenceCriterion(const C1& criterion1, const C2& criterion2) :
        criteria {std::make_shared<C1>(criterion1), std::make_shared<C2>(criterion2)} {}
 
 
    std::string description() const override {
 
        std::string description_string = (boost::format("A compound convergence criterion step consisting of %s combined criteria:\n") % this->count()).str();
 
        for (size_t i = 0; i < this->count(); i++) {
            const auto& criterion = this->criteria[i];
            description_string += (boost::format("\t%s. %s\n") % std::to_string(i + 1) % criterion->description()).str();
        }
        return description_string;
    }
 
 
    bool isFulfilled(Environment& environment) override {
 
        // Check if every criterion is fulfilled.
        for (const auto& criterion : this->criteria) {
            if (!criterion->isFulfilled(environment)) {
                return false;
            }
        }
        return true;
    }
 
 
    /*
     *  PUBLIC METHODS
     */
 
    template <typename C>
    Self add(const C& criterion) {
        this->criteria.push_back(std::make_shared<C>(criterion));
        return *this;
    }
 
 
    size_t count() const { return this->criteria.size(); }
};
 
 
}  // namespace GQCP