pscf_pg.cu

/*
* 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 <prdc/crystal/getDimension.h>
#include <rpg/system/System.h>
#include <iostream>
 
namespace Pscf {
namespace Rpg {

   template <int D>
   void run(int argc, char **argv) {
 
      // Construct System object for specific dimension D.
      System<D> system;
 
      // Process command line options
      system.setOptions(argc, argv);
 
      // Read parameters from default parameter file
      system.readParam();
 
      // Read command script to run system
      system.readCommands();
 
   }
 
}
}

int main(int argc, char **argv)
{
 
   // Extract the dimension of space from argument of -d option
   int D = Pscf::Prdc::getDimension(argc, argv);
   std::cout << "dimension    " << D << std::endl;
 
   if (1 == D) {
      Pscf::Rpg::run<1>(argc, argv);
   } else
   if (2 == D) {
      Pscf::Rpg::run<2>(argc, argv);
   } else
   if (3 == D) {
      Pscf::Rpg::run<3>(argc, argv);
   } else {
      std::cout << " Invalid dimension = " << D << std::endl;
   }
 
}