Line data    Source code

       1              : // Copyright (C) 2023 Helmut Schmidt University
       2              : // This file is part of the Mamico project. For conditions of distribution
       3              : // and use, please see the copyright notice in Mamico's main folder
       4              : 
       5              : #include "coupling/solvers/LBCouetteSolver.h"
       6              : 
       7              : using State = coupling::interface::PintableMacroSolverState;
       8              : 
       9           16 : std::unique_ptr<State> coupling::solvers::LBCouetteSolverState::operator+(const State& rhs) {
      10           16 :   const LBCouetteSolverState* other = dynamic_cast<const LBCouetteSolverState*>(&rhs);
      11              : 
      12              : #if (COUPLING_MD_ERROR == COUPLING_MD_YES)
      13           16 :   if (other == nullptr) {
      14            0 :     std::cout << "ERROR LBCouetteSolverState operator+ type mismatch" << std::endl;
      15            0 :     exit(EXIT_FAILURE);
      16              :   }
      17           16 :   if (other->_pdf.size() != _pdf.size()) {
      18            0 :     std::cout << "ERROR LBCouetteSolverState operator+ size mismatch" << std::endl;
      19            0 :     exit(EXIT_FAILURE);
      20              :   }
      21              : #endif
      22              : 
      23           16 :   std::unique_ptr<LBCouetteSolverState> res = std::make_unique<LBCouetteSolverState>(*this);
      24           64 :   for (std::vector<double>::size_type i = 0; i < _pdf.size(); ++i)
      25           48 :     res->_pdf[i] += other->_pdf[i];
      26              : 
      27           16 :   return res;
      28           16 : }
      29              : 
      30            8 : std::unique_ptr<State> coupling::solvers::LBCouetteSolverState::operator-(const State& rhs) {
      31            8 :   const LBCouetteSolverState* other = dynamic_cast<const LBCouetteSolverState*>(&rhs);
      32              : 
      33              : #if (COUPLING_MD_ERROR == COUPLING_MD_YES)
      34            8 :   if (other == nullptr) {
      35            0 :     std::cout << "ERROR LBCouetteSolverState operator- type mismatch" << std::endl;
      36            0 :     exit(EXIT_FAILURE);
      37              :   }
      38            8 :   if (other->_pdf.size() != _pdf.size()) {
      39            0 :     std::cout << "ERROR LBCouetteSolverState operator- size mismatch" << std::endl;
      40            0 :     exit(EXIT_FAILURE);
      41              :   }
      42              : #endif
      43              : 
      44            8 :   std::unique_ptr<LBCouetteSolverState> res = std::make_unique<LBCouetteSolverState>(*this);
      45           32 :   for (std::vector<double>::size_type i = 0; i < _pdf.size(); ++i)
      46           24 :     res->_pdf[i] -= other->_pdf[i];
      47              : 
      48            8 :   return res;
      49            8 : }