Line data    Source code

       1             : // This file is part of the Mamico project. For conditions of distribution
       2             : // and use, please see the copyright notice in Mamico's main folder, or at
       3             : // www5.in.tum.de/mamico
       4             : 
       5             : #pragma once
       6             : 
       7             : #define DEBUG_FILTER_JUNCTION_ASYM
       8             : 
       9             : #include "coupling/filtering/sequencing/FilterSequence.h"
      10             : 
      11             : // INCLUDE ALL RELEVANT JUNCTOR HEADERS HERE
      12             : #include "coupling/filtering/filters/WriteToFileJunctor.h"
      13             : 
      14             : /*
      15             :  * Filtering junction which does not use multiple identical input domains, but
      16             :  * rather uses one primary input partition to filter on while having access to a
      17             :  * secondary input data set. As no filtering is conducted on that secondary
      18             :  * input cell data set, it may be of any shape or size. The primary input
      19             :  * partition has to follow the restrictions to input data sets of
      20             :  * FilterSequences and FilterJunctions.
      21             :  *
      22             :  * @author Felix Maurer
      23             :  */
      24             : 
      25             : namespace coupling {
      26             : namespace filtering {
      27             : template <unsigned int dim> class AsymmetricalFilterJunction;
      28             : }
      29             : } // namespace coupling
      30             : 
      31             : template <unsigned int dim> class coupling::filtering::AsymmetricalFilterJunction : public coupling::filtering::FilterSequence<dim> {
      32             : public:
      33           0 :   AsymmetricalFilterJunction(const char* name,
      34             :                              const std::vector<coupling::datastructures::CouplingCell<dim>*> primaryInputCellVector,   // primary input of sequence.
      35             :                              const std::vector<coupling::datastructures::CouplingCell<dim>*> secondaryInputCellVector, // additional data, presented as
      36             :                                                                                                                        // macro cells as well
      37             : #if (COUPLING_MD_PARALLEL == COUPLING_MD_YES)
      38             :                              MPI_Comm comm,
      39             : #endif
      40             :                              std::array<bool, 7> filteredValues)
      41             :       : coupling::filtering::FilterSequence<dim>(name, primaryInputCellVector,
      42             : #if (COUPLING_MD_PARALLEL == COUPLING_MD_YES)
      43             :                                                  comm,
      44             : #endif
      45             :                                                  filteredValues),
      46           0 :         _inputCellVector_secondary(secondaryInputCellVector) {
      47             : #ifdef DEBUG_FILTER_JUNCTION_ASYM
      48           0 :     std::cout << PRINT_PREFIX() << "Begin initialization." << std::endl;
      49             : #endif
      50             : 
      51             :     // allocate and init secondary cell vectors
      52           0 :     for (auto cell : _inputCellVector_secondary) {
      53           0 :       _cellVector1_secondary.push_back(new coupling::datastructures::CouplingCell<dim>(*cell));
      54           0 :       _cellVector2_secondary.push_back(new coupling::datastructures::CouplingCell<dim>(*cell));
      55             :     }
      56             : #ifdef DEBUG_FILTER_JUNCTION_ASYM
      57           0 :     std::cout << PRINT_PREFIX() << "Initialized secondary cell vectors." << std::endl;
      58           0 :     std::cout << PRINT_PREFIX() << "First element of _cellVector1_secondary after init: " << _cellVector1_secondary[0] << std::endl;
      59           0 :     std::cout << PRINT_PREFIX() << "First element of _cellVector2_secondary after init: " << _cellVector2_secondary[0] << std::endl;
      60             : #endif
      61             : 
      62           0 :     coupling::filtering::FilterSequence<dim>::_isModifiable = false; // Dynamic filters are not yet supported. TODO
      63           0 :   }
      64             : 
      65           0 :   ~AsymmetricalFilterJunction() {
      66           0 :     for (auto secondarycell : _cellVector1_secondary)
      67           0 :       delete secondarycell;
      68           0 :     for (auto secondarycell : _cellVector2_secondary)
      69           0 :       delete secondarycell;
      70           0 :   }
      71             : 
      72             :   /*
      73             :    * This function is very similar to the interface's. Check
      74             :    * coupling::FilterSequence for more details.
      75             :    */
      76             :   int loadFiltersFromXML(tinyxml2::XMLElement* sequenceNode) override;
      77             : 
      78           0 :   void printFilters() override {
      79           0 :     std::cout << "Junctors in asymmetrical junction " << coupling::filtering::FilterSequence<dim>::_name << ": ";
      80           0 :     for (auto f : coupling::filtering::FilterSequence<dim>::_filters)
      81           0 :       std::cout << f->getType() << " ";
      82           0 :     std::cout << std::endl;
      83           0 :   }
      84             : 
      85           0 :   std::string PRINT_PREFIX() const override {
      86           0 :     return std::string("   AFJ(").std::string::append(coupling::filtering::FilterSequence<dim>::_name).std::string::append("): ");
      87             :   }
      88             : 
      89             : private:
      90             :   std::vector<coupling::datastructures::CouplingCell<dim>*> _inputCellVector_secondary;
      91             : 
      92             :   std::vector<coupling::datastructures::CouplingCell<dim>*> _cellVector1_secondary;
      93             :   std::vector<coupling::datastructures::CouplingCell<dim>*> _cellVector2_secondary;
      94             : };
      95             : 
      96             : // inlcude implementation
      97             : #include "coupling/filtering/sequencing/AsymmetricalFilterJunction.cpph"