Line data    Source code

       1             : // Copyright (C) 2015 Technische Universitaet Muenchen
       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, or at
       4             : // www5.in.tum.de/mamico
       5             : 
       6             : template <class LinkedCell, unsigned int dim>
       7           0 : void coupling::transferstrategies::TransferStrategy4SchwarzCoupling<LinkedCell, dim>::beginProcessInnerCouplingCellsBeforeReceivingMacroscopicSolverData() {
       8             :   // reset sample counter
       9           0 :   _sampleCounter = 0;
      10           0 : }
      11             : 
      12             : template <class LinkedCell, unsigned int dim>
      13           0 : void coupling::transferstrategies::TransferStrategy4SchwarzCoupling<LinkedCell, dim>::processInnerCouplingCellAfterReceivingMacroscopicSolverData(
      14             :     coupling::datastructures::CouplingCellWithLinkedCells<LinkedCell, dim>& cell, I02 index) {
      15             :   // convert momentum to velocity values for cont->MD transfer
      16           0 :   if (cell.getMicroscopicMass() == 0.0) {
      17           0 :     cell.setMicroscopicMomentum(tarch::la::Vector<dim, double>(0.0));
      18             :   } else {
      19           0 :     cell.setMicroscopicMomentum((1.0 / cell.getMicroscopicMass()) * cell.getMicroscopicMomentum());
      20             :   }
      21             :   // currently: no mass transfer anywhere
      22           0 :   cell.setMicroscopicMass(0.0);
      23             :   // reset macroscopic solver values (from averaging) and store old mass
      24           0 :   cell.setMacroscopicMass(0.0);
      25           0 :   cell.setMacroscopicMomentum(tarch::la::Vector<dim, double>(0.0));
      26           0 : }
      27             : 
      28             : template <class LinkedCell, unsigned int dim>
      29           0 : void coupling::transferstrategies::TransferStrategy4SchwarzCoupling<LinkedCell, dim>::processOuterCouplingCellAfterReceivingMacroscopicSolverData(
      30             :     coupling::datastructures::CouplingCellWithLinkedCells<LinkedCell, dim>& cell, I02 index) {
      31           0 :   if (cell.getMicroscopicMass() == 0.0) {
      32           0 :     cell.setMicroscopicMomentum(tarch::la::Vector<dim, double>(0.0));
      33             :   } else {
      34           0 :     cell.setMicroscopicMomentum((1.0 / cell.getMicroscopicMass()) * cell.getMicroscopicMomentum());
      35             :   }
      36           0 :   cell.setMicroscopicMass(0.0);
      37           0 :   cell.setMacroscopicMass(0.0);
      38           0 :   cell.setMacroscopicMomentum(tarch::la::Vector<dim, double>(0.0));
      39           0 : }
      40             : 
      41             : template <class LinkedCell, unsigned int dim>
      42           0 : void coupling::transferstrategies::TransferStrategy4SchwarzCoupling<LinkedCell, dim>::processInnerCouplingCellBeforeSendingMDSolverData(
      43             :     coupling::datastructures::CouplingCellWithLinkedCells<LinkedCell, dim>& cell, I02 index) {
      44             : #if (COUPLING_MD_DEBUG == COUPLING_MD_YES)
      45             :   std::cout << "processInnerCouplingCellBeforeSendingMDSolverData(): Data "
      46             :                "before averaging from cell "
      47             :             << index << "Mass: " << cell.getMacroscopicMass() << " , momentum: " << cell.getMacroscopicMomentum() << std::endl;
      48             : #endif
      49             :   // average sampled data from MD
      50           0 :   if (_sampleCounter != 0) {
      51           0 :     const double mass = cell.getMacroscopicMass() / ((double)_sampleCounter);
      52           0 :     const tarch::la::Vector<dim, double> momentum = (1.0 / ((double)_sampleCounter)) * cell.getMacroscopicMomentum();
      53           0 :     cell.setMacroscopicMass(mass);
      54           0 :     cell.setMacroscopicMomentum(momentum);
      55             :   } else {
      56           0 :     cell.setMacroscopicMass(0.0);
      57           0 :     cell.setMacroscopicMomentum(tarch::la::Vector<dim, double>(0.0));
      58             :   }
      59             : #if (COUPLING_MD_DEBUG == COUPLING_MD_YES)
      60             :   std::cout << "Samplecounter=" << _sampleCounter << std::endl;
      61             :   std::cout << "processInnerCouplingCellBeforeSendingMDSolverData(): Data "
      62             :                "from cell "
      63             :             << index << "Mass: " << cell.getMacroscopicMass() << " , momentum: " << cell.getMacroscopicMomentum() << std::endl;
      64             : #endif
      65           0 : }
      66             : 
      67             : template <class LinkedCell, unsigned int dim>
      68           0 : void coupling::transferstrategies::TransferStrategy4SchwarzCoupling<LinkedCell, dim>::processOuterCouplingCellBeforeSendingMDSolverData(
      69             :     coupling::datastructures::CouplingCellWithLinkedCells<LinkedCell, dim>& cell, I02 index) {
      70           0 :   cell.setMacroscopicMass(0.0);
      71           0 :   cell.setMacroscopicMomentum(tarch::la::Vector<dim, double>(0.0));
      72           0 : }
      73             : 
      74             : template <class LinkedCell, unsigned int dim>
      75           0 : void coupling::transferstrategies::TransferStrategy4SchwarzCoupling<LinkedCell, dim>::beginProcessInnerCouplingCellsAfterMDTimestep() {
      76           0 :   _timestepCounter++;
      77           0 :   if (sample()) {
      78           0 :     _sampleCounter++;
      79             :   }
      80           0 : }
      81             : 
      82             : template <class LinkedCell, unsigned int dim>
      83           0 : void coupling::transferstrategies::TransferStrategy4SchwarzCoupling<LinkedCell, dim>::processInnerCouplingCellAfterMDTimestep(
      84             :     coupling::datastructures::CouplingCellWithLinkedCells<LinkedCell, dim>& cell, I02 index) {
      85             :   // averaging: evaluate current state
      86           0 :   if (sample()) {
      87           0 :     cell.iterateConstCells(_massMapping);
      88           0 :     const double mass = _massMapping.getMass();
      89           0 :     cell.iterateConstCells(_momentumMapping);
      90           0 :     const tarch::la::Vector<dim, double> momentum = _momentumMapping.getMomentum();
      91           0 :     cell.addMacroscopicMass(mass);
      92           0 :     cell.addMacroscopicMomentum(momentum);
      93             :   }
      94           0 : }
      95             : 
      96           0 : template <class LinkedCell, unsigned int dim> bool coupling::transferstrategies::TransferStrategy4SchwarzCoupling<LinkedCell, dim>::sample() const {
      97             :   // sample strategy: for a hard-coded time-limit, sample only over 20% of
      98             :   // coupling cycle
      99           0 :   const unsigned int timestepInThisCouplingCycle = _timestepCounter % _numberMDSteps;
     100           0 :   return ((timestepInThisCouplingCycle > 0.8 * _numberMDSteps - 1) && ((_timestepCounter - 1) % _sampleEveryTimestep == 0));
     101             : }