ConstantDamper.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
 
 
namespace GQCP {
 
 
class ConstantDamper {
private:
    double alpha;  // the damping factor
 
 
public:
    /*
     *  CONSTRUCTORS
     */
 
    ConstantDamper(const double alpha) :
        alpha {alpha} {
 
        if ((this->alpha > 1.0) || (this->alpha < 0)) {
            throw std::invalid_argument("ConstantDamper::ConstantDamper(const double): The given damping factor must be between 0.0 and 1.0.");
        }
    }
 
 
    /*
     *  PUBLIC METHODS
     */
 
    template <typename Subject>
    Subject accelerate(const Subject& last, const Subject& next_to_last) {
 
        Subject accelerated_subject = last * this->alpha + next_to_last * (1 - this->alpha);
        return accelerated_subject;
    }
};
 
 
}  // namespace GQCP