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 : }