picard_matrix_solver.h

// This file is part of Hermes2D
//
// Copyright (c) 2009 hp-FEM group at the University of Nevada, Reno (UNR).
// Email: hpfem-group@unr.edu, home page: http://www.hpfem.org/.
//
// Hermes2D 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 2 of the License,
// or (at your option) any later version.
//
// Hermes2D 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 Hermes2D; if not, write to the Free Software
// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
#ifndef __HERMES_COMMON_PICARD_MATRIX_SOLVER_H_
#define __HERMES_COMMON_PICARD_MATRIX_SOLVER_H_

#include "solvers/nonlinear_matrix_solver.h"

namespace Hermes
{
  namespace Solvers
  {
    template<typename Scalar>
    class HERMES_API PicardMatrixSolver : public NonlinearMatrixSolver < Scalar >
    {
    public:
      PicardMatrixSolver();
      virtual ~PicardMatrixSolver() {};

      void use_overloaded_damping_factor_condition(bool onOff);

#pragma region anderson-public
      void use_Anderson_acceleration(bool to_set);

      void set_num_last_vector_used(int num);

      void set_anderson_beta(double beta);
#pragma endregion

    protected:
      virtual double update_solution_return_change_norm(Scalar* linear_system_solution);

      virtual void init_solving(Scalar*& coeff_vec);

      virtual void deinit_solving();

      virtual void solve_linear_system();

      virtual double calculate_residual_norm();

      virtual bool damping_factor_condition();

      void init_picard();

      inline std::string getClassName() const { return "PicardMatrixSolver"; }

      bool damping_factor_condition_overloaded;
#pragma region anderson-private
      // Anderson.
      int num_last_vectors_used;
      bool anderson_is_on;
      double anderson_beta;
      Scalar** previous_vectors;
      Scalar* anderson_coeffs;

      void init_anderson();
      void deinit_anderson();

      void handle_previous_vectors();
      void calculate_anderson_coeffs();

      Scalar* previous_Anderson_sln_vector;

      int vec_in_memory;
#pragma endregion
    };
  }
}
#endif