DirectTransferStrategy.h

// Copyright (C) 2015 Technische Universitaet Muenchen
// This file is part of the Mamico project. For conditions of distribution
// and use, please see the copyright notice in Mamico's main folder, or at
// www5.in.tum.de/mamico
#ifndef _MOLECULARDYNAMICS_COUPLING_TRANSFERSTRATEGIES_DIRECTTRANSFERSTRATEGY_H_
#define _MOLECULARDYNAMICS_COUPLING_TRANSFERSTRATEGIES_DIRECTTRANSFERSTRATEGY_H_
 
#include "coupling/cell-mappings/ComputeMassMapping.h"
#include "coupling/cell-mappings/ComputeMomentumMapping.h"
#include "coupling/transferstrategies/TransferStrategy.h"
#include "tarch/la/Vector.h"
 
namespace coupling {
namespace transferstrategies {
template <class LinkedCell, unsigned int dim> class DirectTransferStrategy;
}
} // namespace coupling

template <class LinkedCell, unsigned int dim>
class coupling::transferstrategies::DirectTransferStrategy : public coupling::transferstrategies::TransferStrategy<LinkedCell, dim> {
public:
  DirectTransferStrategy(coupling::interface::MDSolverInterface<LinkedCell, dim>* const mdSolverInterface)
      : coupling::transferstrategies::TransferStrategy<LinkedCell, dim>(mdSolverInterface), _massMapping(mdSolverInterface),
        _momentumMapping(mdSolverInterface) {}

  virtual ~DirectTransferStrategy() {}

  void processInnerCouplingCellBeforeReceivingMacroscopicSolverData(coupling::datastructures::CouplingCellWithLinkedCells<LinkedCell, dim>& cell,
                                                                    I02 index) override {
    // reset quantities
    const tarch::la::Vector<dim, double> zero(0.0);
    cell.setMicroscopicMass(0.0);
    cell.setMicroscopicMomentum(zero);
  }

  void processOuterCouplingCellBeforeReceivingMacroscopicSolverData(coupling::datastructures::CouplingCellWithLinkedCells<LinkedCell, dim>& cell,
                                                                    I02 index) override {
    // reset quantities
    const tarch::la::Vector<dim, double> zero(0.0);
    cell.setMicroscopicMass(0.0);
    cell.setMicroscopicMomentum(zero);
  }

  void processInnerCouplingCellBeforeSendingMDSolverData(coupling::datastructures::CouplingCellWithLinkedCells<LinkedCell, dim>& cell, I02 index) override {
    cell.iterateConstCells(_massMapping);
    cell.iterateConstCells(_momentumMapping);
    cell.setMacroscopicMass(_massMapping.getMass());
    cell.setMacroscopicMomentum(_momentumMapping.getMomentum());
  }
 
private:
  coupling::cellmappings::ComputeMassMapping<LinkedCell, dim> _massMapping;
  coupling::cellmappings::ComputeMomentumMapping<LinkedCell, dim> _momentumMapping;
};
#endif // _MOLECULARDYNAMICS_COUPLING_TRANSFERSTRATEGIES_DIRECTTRANSFERSTRATEGY_H_