ViewVC Help

View File | Revision Log | Show Annotations | View Changeset | Root Listing

root/OpenMD/trunk/src/io/RestReader.hpp

Comparing trunk/src/io/RestReader.hpp (file contents):
Revision 507 by gezelter, Fri Apr 15 22:04:00 2005 UTC vs.
Revision 1879 by gezelter, Sun Jun 16 15:15:42 2013 UTC

#	Line 1 | Line 1
1	<	/*
2	<	* Copyright (c) 2005 The University of Notre Dame. All Rights Reserved.
1	>	/*
2	>	* Copyright (c) 2005 The University of Notre Dame. All Rights Reserved.
3		*
4		* The University of Notre Dame grants you ("Licensee") a
5		* non-exclusive, royalty free, license to use, modify and
6		* redistribute this software in source and binary code form, provided
7		* that the following conditions are met:
8		*
9	<	* 1. Acknowledgement of the program authors must be made in any
10	<	*    publication of scientific results based in part on use of the
11	<	*    program.  An acceptable form of acknowledgement is citation of
12	<	*    the article in which the program was described (Matthew
13	<	*    A. Meineke, Charles F. Vardeman II, Teng Lin, Christopher
14	<	*    J. Fennell and J. Daniel Gezelter, "OOPSE: An Object-Oriented
15	<	*    Parallel Simulation Engine for Molecular Dynamics,"
16	<	*    J. Comput. Chem. 26, pp. 252-271 (2005))
17	<	*
18	<	* 2. Redistributions of source code must retain the above copyright
9	>	* 1. Redistributions of source code must retain the above copyright
10		*    notice, this list of conditions and the following disclaimer.
11		*
12	<	* 3. Redistributions in binary form must reproduce the above copyright
12	>	* 2. Redistributions in binary form must reproduce the above copyright
13		*    notice, this list of conditions and the following disclaimer in the
14		*    documentation and/or other materials provided with the
15		*    distribution.
#	Line 37 | Line 28
28		* arising out of the use of or inability to use software, even if the
29		* University of Notre Dame has been advised of the possibility of
30		* such damages.
31	<	*/
32	<	#ifndef IO_RESTREADER_HPP
33	<	#define IO_RESTREADER_HPP
34	<
35	<	#define _LARGEFILE_SOURCE64
36	<	#ifndef _FILE_OFFSET_BITS
37	<	#define _FILE_OFFSET_BITS 64
38	<	#endif
39	<	#include <iostream>
40	<	#include <fstream>
31	>	*
32	>	* SUPPORT OPEN SCIENCE!  If you use OpenMD or its source code in your
33	>	* research, please cite the appropriate papers when you publish your
34	>	* work.  Good starting points are:
35	>	*
36	>	* [1]  Meineke, et al., J. Comp. Chem. 26, 252-271 (2005).
37	>	* [2]  Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006).
38	>	* [3]  Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008).
39	>	* [4]  Kuang & Gezelter,  J. Chem. Phys. 133, 164101 (2010).
40	>	* [5]  Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011).
41	>	*/
42	>
43	>	/**
44	>	* @file RestReader.hpp
45	>	*/
46	>
47	>	#ifndef IO_RESTREADER_HPP
48	>	#define IO_RESTREADER_HPP
49	>
50	>	#include <cstdio>
51		#include <string>
51	–	#include <stdio.h>
52	–	#include <stdlib.h>
53	–	#include <unistd.h>
54	–	#include <sys/types.h>
55	–	#include <sys/stat.h>
56	–
52		#include "brains/SimInfo.hpp"
53		#include "primitives/StuntDouble.hpp"
54
55	<	namespace oopse {
56	<
55	>	namespace OpenMD {
56	>
57	>	/**
58	>	* @class RestReader RestReader.hpp "io/RestReader.hpp"
59	>	*/
60	>	class RestReader {
61	>	public:
62	>	RestReader(SimInfo* info, const std::string &filename,
63	>	std::vector<int> stuntDoubleIndex)
64	>	: info_(info), filename_(filename), stuntDoubleIndex_(stuntDoubleIndex) {
65	>	#ifdef IS_MPI
66	>
67	>	if (worldRank == 0) {
68	>	#endif
69	>
70	>	inFile_ = new std::ifstream(filename_.c_str(),
71	>	ifstream::in | ifstream::binary);
72	>
73	>	if (inFile_->fail()) {
74	>	sprintf(painCave.errMsg,
75	>	"RestReader: Cannot open file: %s\n",
76	>	filename_.c_str());
77	>	painCave.isFatal = 1;
78	>	simError();
79	>	}
80	>
81	>	#ifdef IS_MPI
82	>
83	>	}
84	>
85	>	strcpy(checkPointMsg, "Reference file opened for reading successfully.");
86	>	errorCheckPoint();
87	>
88	>	#endif
89	>
90	>	return;
91	>	}
92
93	<	/**
64	<	* @class RestReader RestReader.hpp "io/RestReader.hpp"
65	<	* @for extra file reading when using thermodynamic integration
66	<	*/
67	<	class RestReader {
68	<
69	<	public:
70	<	RestReader( SimInfo* theInfo );
71	<	~RestReader();
72	<
73	<	void readIdealCrystal();
74	<	void readZangle();
75	<	void zeroZangle();
76	<
77	<	private:
78	<	char* parseIdealLine(char* readLine, StuntDouble* sd);
79	<	char *idealName;
80	<	FILE *inAngFile;
81	<	FILE *inIdealFile;
82	<	std::string inAngFileName;
83	<	std::string inIdealFileName;
84	<	bool isScanned;
85	<
86	<	double angleTransfer;
87	<	std::vector<fpos_t*> framePos;
88	<	SimInfo *info_;
89	<	};
93	>	void readReferenceStructure(void);
94
95	<	} //end namespace oopse
95	>	~RestReader(){
96	>	#ifdef IS_MPI
97	>	if (worldRank == 0){
98	>	#endif
99	>
100	>	delete inFile_;
101	>
102	>	#ifdef IS_MPI
103	>	}
104	>
105	>	strcpy(checkPointMsg, "Reference file closed successfully.");
106	>	errorCheckPoint();
107		#endif
108	+	}
109	+
110	+	protected:
111	+	void parseDumpLine(const std::string& line);
112	+	void readStuntDoubles(std::istream& inpuStream);
113	+	void readFrameProperties(std::istream& inputStream);
114	+	void scanFile(void);
115	+	void readSet(void);
116	+	private:
117	+
118	+	SimInfo* info_;
119	+
120	+	std::vector<Vector3d> all_pos_;
121	+	std::vector<int> stuntDoubleIndex_;
122	+
123	+	std::istream* inFile_;
124	+	std::string filename_;
125	+
126	+	std::streampos framePos_;
127	+
128	+	const static int bufferSize = 4096;
129	+	char buffer[bufferSize];
130	+	};
131	+
132	+	}      //end namespace OpenMD
133	+
134	+	#endif //IO_RESTREADER_HPP

Comparing trunk/src/io/RestReader.hpp (property svn:keywords):
Revision 507 by gezelter, Fri Apr 15 22:04:00 2005 UTC vs.
Revision 1879 by gezelter, Sun Jun 16 15:15:42 2013 UTC

#	Line 0 | Line 1
1	+	Author Id Revision Date

Diff Legend

–	Removed lines
+	Added lines
<	Changed lines
>	Changed lines