fieldHeader.tpp

#ifndef PRDC_FIELD_HEADER_TPP
#define PRDC_FIELD_HEADER_TPP
 
/*
* PSCF - Polymer Self-Consistent Field Theory
*
* Copyright 2016 - 2022, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/
 
#include "fieldHeader.h"
 
#include <util/format/Dbl.h> 
#include <util/format/Int.h> 
#include <iomanip>
 
namespace Pscf {
namespace Prdc {
 
   using namespace Util;
 
   /*
   * Read common part of field header (fortran PSCF format).
   */
   template <int D>
   void readFieldHeader(std::istream& in, 
                        int& ver1, int& ver2, 
                        UnitCell<D>& cell, 
                        std::string& groupName,
                        int& nMonomer)
   {
      std::string label;
 
      // Read format v1 v2
      in >> label;
      UTIL_CHECK(label == "format");
      in >> ver1;
      in >> ver2;
 
      // Read dimension of space 
      in >> label;
      UTIL_CHECK(label == "dim");
      int dim;
      in >> dim;
      UTIL_CHECK(dim == D);
 
      // Read unit cell lattice type, number of parameters, parameters
      // Function declared in UnitCell.h and defined in UnitCell.tpp
      readUnitCellHeader(in, cell);
 
      // Optionally read groupName
      in >> label;
      groupName = "";
      if (label == "group_name") {
         in >> groupName;
         in >> label;
      }
 
      // Read nMonomer value
      UTIL_CHECK(label == "N_monomer");
      in >> nMonomer;
      UTIL_CHECK(nMonomer > 0);
   }
 
   /*
   * Write common part of field header (fortran PSCF format).
   */
   template <int D>
   void writeFieldHeader(std::ostream &out, 
                         int ver1, int ver2,
                         UnitCell<D> const & cell,
                         std::string const & groupName,
                         int nMonomer)
   {
      // Write file format id and dimension of space D
      out << "format " << Int(ver1,3) << " " << Int(ver2,3) <<  std::endl;
      out << "dim" <<  std::endl 
          << "          " << D << std::endl;
 
      // Write unit cell lattice type, number of parameters, parameters
      // Function declared in UnitCell.h and defined in UnitCell.tpp
      writeUnitCellHeader(out, cell); 
 
      // Write group_name if not empty
      if (groupName != "") {
         out << "group_name" << std::endl 
             << "          " << groupName <<  std::endl;
      }
 
      // Write nMonomer, number of monomer types
      out << "N_monomer"  << std::endl 
          << "          " << nMonomer << std::endl;
   }
 
}
}
#endif