b2static_linear_p_refinement_solver.H

//------------------------------------------------------------------------
// b2static_linear_p_refinement_solver.H --
//
//
// written by Mathias Doreille
//            Neda Ebrahimi Pour <neda.ebrahimipour@dlr.de>
//            Thomas Blome <thomas.blome@dlr.de>
//
// (c) 2004-2012 SMR Engineering & Development SA
//               2502 Bienne, Switzerland
//
// (c) 2023 Deutsches Zentrum für Luft- und Raumfahrt (DLR) e.V.
//          Linder Höhe, 51147 Köln
//
// All Rights Reserved.  Proprietary source code.  The contents of
// this file may not be disclosed to third parties, copied or
// duplicated in any form, in whole or in part, without the prior
// written permission of SMR.
//------------------------------------------------------------------------
 
#ifndef B2STATIC_LINEAR_P_REFINEMENT_SOLVER_H_
#define B2STATIC_LINEAR_P_REFINEMENT_SOLVER_H_
 
#include "b2static_linear_solver.H"
 
namespace b2000::solver {
 
template <typename T, typename MATRIX_FORMAT>
class StaticLinearPRefinementSolver : public StaticLinearSolver<T, MATRIX_FORMAT> {
public:
    using type_t = ObjectTypeComplete<
          StaticLinearPRefinementSolver, typename StaticLinearSolver<T, MATRIX_FORMAT>::type_t>;
    using StaticLinearSolver<T, MATRIX_FORMAT>::case_iterator_;
 
    void solve() {
        while (solve_iteration()) { ; }
    }
 
    bool solve_iteration() {
        if (case_iterator_.get() == nullptr) {
            case_iterator_ =
                  StaticLinearSolver<T, MATRIX_FORMAT>::model_->get_case_list().get_case_iterator();
            this->case_ = case_iterator_->next();
        }
        if (this->case_ == nullptr) {
            // All the case are solved.
            return false;
        }
        StaticLinearSolver<T, MATRIX_FORMAT>::solve_on_case(*(this->case_));
        if (!StaticLinearSolver<T, MATRIX_FORMAT>::model_->get_solution().is_need_refinement()) {
            // Go to the next case
            this->case_ = case_iterator_->next();
        }
        return this->case_ != nullptr;
    }
 
    static type_t type;
};
 
template <typename T, typename MATRIX_FORMAT>
typename StaticLinearPRefinementSolver<T, MATRIX_FORMAT>::type_t
      StaticLinearPRefinementSolver<T, MATRIX_FORMAT>::type(
            "StaticLinearPRefinementSolver", type_name<T, MATRIX_FORMAT>(),
            StringList("LINEAR_P_REFINEMENT"), module, &Solver::type);
 
}  // namespace b2000::solver
 
#endif  // B2STATIC_LINEAR_P_REFINEMENT_SOLVER_H_