sortWaves.h

#ifndef PRDC_SORT_WAVES_H
#define PRDC_SORT_WAVES_H
 
/*
* PSCF - Polymer Self-Consistent Field
*
* Copyright 2015 - 2025, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/
 
#include <pscf/math/Sort.h>
#include <pscf/math/IntVec.h>
#include <util/containers/GArray.h>
#include <vector>
 
// Forward declarations
namespace Pscf {
   namespace Prdc {
      template <int D> class UnitCell;
   }
}
 
namespace Pscf {
namespace Prdc {
 
   using namespace Prdc;

   template <int D>
   void sortWaves(UnitCell<D> const & cell,
                  IntVec<D> const & meshDimensions,
                  std::vector< Sort::Item<double> >& items,
                  GArray< Sort::Bunch >& bunches,
                  double epsilon,
                  bool isRealField = true);
 
   // Explicit instantiation declarations
   extern template 
   void sortWaves<1>(UnitCell<1> const &, IntVec<1> const &, 
                     std::vector< Sort::Item<double> >& , 
                     GArray< Sort::Bunch >& , double, bool);
 
   extern template 
   void sortWaves<2>(UnitCell<2> const &, IntVec<2> const &, 
                     std::vector< Sort::Item<double> >& , 
                     GArray< Sort::Bunch >& , double, bool);
 
   extern template 
   void sortWaves<3>(UnitCell<3> const &, IntVec<3> const &, 
                     std::vector< Sort::Item<double> >& , 
                     GArray< Sort::Bunch >& , double, bool);
 
}
}
#endif