seams_input.hpp

Go to the documentation of this file.

//-----------------------------------------------------------------------------------
// d-SEAMS - Deferred Structural Elucidation Analysis for Molecular Simulations
//
// Copyright (c) 2018--present d-SEAMS core team
//
// This program is free software: you can redistribute it and/or modify
// it under the terms of the MIT License as published by
// the Open Source Initiative.
//
// A copy of the MIT License is included in the LICENSE file of this repository.
// You should have received a copy of the MIT License along with this program.
// If not, see <https://opensource.org/licenses/MIT>.
//-----------------------------------------------------------------------------------
 
#ifndef __SEAMS_INPUT_H_
#define __SEAMS_INPUT_H_
 
#include <iostream>
#include <memory>
#include <mol_sys.hpp>
#include <ring.hpp>
#include <string>
 
#include "boost/filesystem/operations.hpp"
#include "boost/filesystem/path.hpp"
 
namespace sinp {
 
std::vector<std::string> getInpFileList(std::string inputFolder);
 
molSys::PointCloud<molSys::Point<double>, double>
readLammpsTrj(std::string filename, int targetFrame,
              molSys::PointCloud<molSys::Point<double>, double> *yCloud,
              bool isSlice = false,
              std::array<double, 3> coordLow = std::array<double, 3>{0, 0, 0},
              std::array<double, 3> coordHigh = std::array<double, 3>{0, 0, 0});
 
molSys::PointCloud<molSys::Point<double>, double> readLammpsTrjO(
    std::string filename, int targetFrame,
    molSys::PointCloud<molSys::Point<double>, double> *yCloud, int typeO,
    bool isSlice = false,
    std::array<double, 3> coordLow = std::array<double, 3>{0, 0, 0},
    std::array<double, 3> coordHigh = std::array<double, 3>{0, 0, 0});
 
molSys::PointCloud<molSys::Point<double>, double> readLammpsTrjreduced(
    std::string filename, int targetFrame,
    molSys::PointCloud<molSys::Point<double>, double> *yCloud, int typeI,
    bool isSlice = false,
    std::array<double, 3> coordLow = std::array<double, 3>{0, 0, 0},
    std::array<double, 3> coordHigh = std::array<double, 3>{0, 0, 0});
 
molSys::PointCloud<molSys::Point<double>, double> readXYZ(std::string filename);
 
std::vector<std::vector<int>> readBonds(std::string filename);
 
inline bool atomInSlice(double x, double y, double z,
                        std::array<double, 3> coordLow,
                        std::array<double, 3> coordHigh) {
  int flag = 0; 
 
  if (((x >= coordLow[0]) && (x <= coordHigh[0])) ||
      coordLow[0] == coordHigh[0]) {
    flag++;
  }
  if (((y >= coordLow[1]) && (y <= coordHigh[1])) ||
      coordLow[1] == coordHigh[1]) {
    flag++;
  }
  if (((z >= coordLow[2]) && (z <= coordHigh[2])) ||
      coordLow[2] == coordHigh[2]) {
    flag++;
  }
 
  if (flag == 3) {
    return true;
  } else {
    return false;
  }
}
 
} // namespace sinp
 
#endif