constant_relaxation_convergence_accelerator.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Ruben Zorrilla
//
//
 
#if !defined(KRATOS_CONSTANT_RELAXATION_CONVERGENCE_ACCELERATOR)
#define  KRATOS_CONSTANT_RELAXATION_CONVERGENCE_ACCELERATOR
 
/* System includes */
 
 
/* External includes */
 
 
/* Project includes */
#include "solving_strategies/convergence_accelerators/convergence_accelerator.h"
 
 
namespace Kratos
{
 
 
 
 
 
 
template<class TSparseSpace, class TDenseSpace>
class ConstantRelaxationConvergenceAccelerator: public ConvergenceAccelerator<TSparseSpace, TDenseSpace>
{
public:
 
 
    KRATOS_CLASS_POINTER_DEFINITION( ConstantRelaxationConvergenceAccelerator );
 
    typedef ConvergenceAccelerator<TSparseSpace, TDenseSpace>              BaseType;
 
    typedef typename BaseType::Pointer                              BaseTypePointer;
 
    typedef typename BaseType::VectorType                                VectorType;
    typedef typename BaseType::VectorPointerType                  VectorPointerType;
 
 
    static Parameters GetDefaultParameters()
    {
        Parameters default_settings(R"(
        {
            "solver_type": "constant_relaxation",
            "w": 0.5,
            "interface_block_newton" : false
        }
        )");
 
        return default_settings;
    }
 
    ConstantRelaxationConvergenceAccelerator(Parameters rConvAcceleratorParameters)
    : BaseType(),
      mOmega([] (Parameters x) -> double {x.ValidateAndAssignDefaults(GetDefaultParameters()); return x["w"].GetDouble();} (rConvAcceleratorParameters))
    {
    }
 
    ConstantRelaxationConvergenceAccelerator(const double rOmega = 0.5)
    : BaseType(),
      mOmega(rOmega)
    {
    }
 
    ConstantRelaxationConvergenceAccelerator(const ConstantRelaxationConvergenceAccelerator &rOther) = delete;
 
    virtual ~ConstantRelaxationConvergenceAccelerator() = default;
 
 
 
    void UpdateSolution(const VectorType &rResidualVector,
                        VectorType &rIterationGuess) override
    {
        KRATOS_TRY;
 
        TSparseSpace::UnaliasedAdd(rIterationGuess, mOmega, rResidualVector);
 
        KRATOS_CATCH("");
    }
 
 
 
 
 
private:
 
 
    const double mOmega;
 
 
 
 
 
 
 
}; /* Class ConstantRelaxationConvergenceAccelerator */
 
 
}  /* namespace Kratos.*/
 
#endif /* KRATOS_CONSTANT_RELAXATION_CONVERGENCE_ACCELERATOR defined */