Line data    Source code

       1             : // This file is part of the Mamico proJunctionn conditions of distribution
       2             : // and use, please see the copyright notice in Mamico's main folder, or at
       3             : // www5.in.tum.de/mamic::o
       4             : 
       5             : // Member functions of coupling::filtering::AsymmetricalFilterJunction
       6             : 
       7             : template <unsigned int dim>
       8           0 : int coupling::filtering::AsymmetricalFilterJunction<dim>::loadFiltersFromXML(
       9             :     tinyxml2::XMLElement* sequenceNode) { // in our case "junctionNode" would be more precise
      10           0 :   if (!sequenceNode)
      11             :     return EXIT_FAILURE;
      12             : 
      13             :   // null if junction is empty
      14           0 :   tinyxml2::XMLElement* currJunctorNode = sequenceNode->FirstChildElement();
      15             : 
      16             :   // These four are not to be confused with member variables of similar names.
      17             :   // Watch out for the "_" in front of member variables.
      18           0 :   std::vector<coupling::datastructures::CouplingCell<dim>*> inputCellVector_primary = coupling::filtering::FilterSequence<dim>::_inputCellVector;
      19           0 :   std::vector<coupling::datastructures::CouplingCell<dim>*> outputCellVector_primary = coupling::filtering::FilterSequence<dim>::_cellVector2;
      20             : 
      21           0 :   std::vector<coupling::datastructures::CouplingCell<dim>*> inputCellVector_secondary = _inputCellVector_secondary;
      22             :   // Note: To prevent counter intuitive design, secondary output has been
      23             :   // disabled, i.e. all junctors in an AFJ take the same secondary input.
      24             :   // There
      25             :   // is a chance that this will be re-enabled in the future, if it is needed for
      26             :   // a certain junctor to work.
      27             :   // std::vector<coupling::datastructures::CouplingCell<dim>* >
      28             :   // outputCellVector_secondary = _cellVector2_secondary;
      29             : 
      30             : #ifdef DEBUG_FILTER_PIPELINE
      31             :   std::cout << PRINT_PREFIX() << "Initializing junctor objects of asymmetrical junction." << std::endl;
      32             : #endif
      33             : 
      34             :   // Keep track of this for each indivial input partition.
      35             :   bool firstIteration = {true};
      36             : 
      37           0 :   while (currJunctorNode) {
      38             : 
      39             :     // All Junctors must implement one of the {X,Y}-Junctor interfaces, each of
      40             :     // which implement FilterInterface.
      41           0 :     coupling::filtering::FilterInterface<dim>* newjunctor = nullptr;
      42             : 
      43             :     // ###############################################################
      44             :     //  This is where you have to define how to handle each junctor individually.
      45             : 
      46             :     // Write To File Junction
      47             :     // TODO: documentation about parametrization in .xml
      48           0 :     if (std::strcmp(currJunctorNode->Value(), "write-to-file-junction") == 0) {
      49             : 
      50             :       // Check "overwrite" property for both underlying WriteToFile instances.
      51           0 :       std::array<bool, 2> overwrite = {false};
      52           0 :       if (currJunctorNode->Attribute("wtf1-overwrite"))
      53           0 :         if (std::strcmp(currJunctorNode->Attribute("wtf1-overwrite"), "true") == 0)
      54           0 :           overwrite[0] = true;
      55           0 :       if (currJunctorNode->Attribute("wtf2-overwrite"))
      56           0 :         if (std::strcmp(currJunctorNode->Attribute("wtf2-overwrite"), "true") == 0)
      57           0 :           overwrite[1] = true;
      58             : 
      59             :       // Check "one-cell-only" property for both underlying WriteToFile
      60             :       // instances.
      61             :       // TODO: replace std::couts with proper exceptions
      62           0 :       std::array<int, 2> oco = {-1, -1};
      63           0 :       tarch::configuration::ParseConfiguration::readIntOptional(oco[0], currJunctorNode, "wtf1-one-cell-only");
      64           0 :       if (oco[0] < -1 || oco[0] >= (int)inputCellVector_primary.size()) {
      65           0 :         std::cout << "ERROR: XML config file references undefined cell index " << oco[0] << " for primary write-to-file in write-to-file-junctor." << std::endl;
      66           0 :         exit(EXIT_FAILURE);
      67             :       }
      68           0 :       tarch::configuration::ParseConfiguration::readIntOptional(oco[1], currJunctorNode, "wtf2-one-cell-only");
      69           0 :       if (oco[1] < -1 || oco[1] >= (int)inputCellVector_secondary.size()) {
      70           0 :         std::cout << "ERROR: XML config file references undefined cell index " << oco[1] << " for secondary write-to-file in write-to-file-junctor."
      71           0 :                   << std::endl;
      72           0 :         exit(EXIT_FAILURE);
      73             :       }
      74             : 
      75           0 :       newjunctor = new coupling::filtering::WriteToFileJunctor<dim>(
      76           0 :           inputCellVector_primary, outputCellVector_primary, inputCellVector_secondary, this->_filteredValues,
      77             :           {currJunctorNode->Attribute("wtf1-location"), currJunctorNode->Attribute("wtf2-location")}, overwrite, oco);
      78             :     }
      79             : 
      80             :     // ###############################################################
      81             : 
      82             :     // Unknown junctor identifier
      83             :     else {
      84             :       std::cout << "ERROR: Asymmetrical Filter-Junction does not recognize "
      85             :                    "junctor called "
      86           0 :                 << currJunctorNode->Value() << std::endl;
      87           0 :       return EXIT_FAILURE;
      88             :     }
      89             : 
      90           0 :     if (firstIteration) {
      91           0 :       inputCellVector_primary = coupling::filtering::FilterSequence<dim>::_cellVector1;
      92           0 :       inputCellVector_secondary = _cellVector1_secondary;
      93             :       firstIteration = false;
      94             :     }
      95             : 
      96           0 :     inputCellVector_primary.swap(outputCellVector_primary);
      97             : 
      98             :     // See comment above regarding secondary output of junctors.
      99             :     // inputCellVector_secondary.swap(outputCellVector_secondary);
     100             : 
     101           0 :     coupling::filtering::FilterSequence<dim>::_filters.push_back(newjunctor);
     102           0 :     currJunctorNode = currJunctorNode->NextSiblingElement();
     103             :   }
     104             : 
     105             : #ifdef DEBUG_FILTER_JUNCTION_ASYM
     106           0 :   std::cout << PRINT_PREFIX();
     107           0 :   this->printOutputCellVector();
     108             : #endif
     109             : 
     110             :   return 0;
     111           0 : }