1 |
mmeineke |
377 |
#ifndef __READWRITE_H__ |
2 |
|
|
#define __READWRITE_H__ |
3 |
|
|
|
4 |
|
|
#include <iostream> |
5 |
|
|
#include <fstream> |
6 |
|
|
|
7 |
gezelter |
829 |
#include <string.h> |
8 |
|
|
#include <stdio.h> |
9 |
|
|
#include <stdlib.h> |
10 |
mmeineke |
804 |
#include <unistd.h> |
11 |
|
|
#include <sys/types.h> |
12 |
|
|
#include <sys/stat.h> |
13 |
mmeineke |
377 |
|
14 |
mmeineke |
646 |
|
15 |
mmeineke |
377 |
#include "Atom.hpp" |
16 |
|
|
#include "SimInfo.hpp" |
17 |
|
|
#include "Thermo.hpp" |
18 |
|
|
|
19 |
|
|
using namespace std; |
20 |
|
|
|
21 |
mmeineke |
804 |
class FilePos{ |
22 |
tim |
837 |
|
23 |
mmeineke |
804 |
public: |
24 |
tim |
837 |
FilePos(){ |
25 |
mmeineke |
804 |
myPos = NULL; |
26 |
|
|
next = NULL; |
27 |
|
|
} |
28 |
tim |
837 |
FilePos( fpos_t* thePos ) { |
29 |
|
|
myPos = thePos; |
30 |
mmeineke |
804 |
next = NULL; |
31 |
|
|
} |
32 |
tim |
837 |
~FilePos(){ |
33 |
mmeineke |
804 |
if( next != NULL ) delete next; |
34 |
|
|
if( myPos != NULL ) delete myPos; |
35 |
|
|
} |
36 |
tim |
837 |
|
37 |
mmeineke |
804 |
void add( fpos_t *thePos ){ |
38 |
tim |
837 |
if( next != NULL ) |
39 |
mmeineke |
804 |
next->add( thePos ); |
40 |
|
|
else |
41 |
|
|
next = new FilePos( thePos ); |
42 |
|
|
} |
43 |
tim |
837 |
|
44 |
mmeineke |
804 |
FilePos &operator=(fpos_t *thePos){ myPos = thePos; return *this; } |
45 |
tim |
837 |
|
46 |
mmeineke |
804 |
void setPos( fpos_t *thePos ){ myPos = thePos; } |
47 |
|
|
fpos_t *getPos( void ){ return myPos; } |
48 |
tim |
837 |
|
49 |
mmeineke |
804 |
FilePos* getNext( void ) { return next; } |
50 |
|
|
private: |
51 |
tim |
837 |
|
52 |
mmeineke |
804 |
fpos_t *myPos; |
53 |
|
|
FilePos* next; |
54 |
tim |
837 |
|
55 |
mmeineke |
804 |
}; |
56 |
mmeineke |
646 |
|
57 |
mmeineke |
377 |
class DumpWriter{ |
58 |
|
|
|
59 |
|
|
public: |
60 |
|
|
DumpWriter( SimInfo* the_entry_plug ); |
61 |
|
|
~DumpWriter(); |
62 |
|
|
|
63 |
|
|
void writeDump( double currentTime ); |
64 |
tim |
929 |
void writeFinal( double currentTime); |
65 |
tim |
934 |
void writeFrame( vector<ofstream*>& outFile, double finalTime ); |
66 |
mmeineke |
377 |
|
67 |
tim |
929 |
#ifdef IS_MPI |
68 |
|
|
void update(); |
69 |
|
|
#endif |
70 |
|
|
|
71 |
mmeineke |
377 |
private: |
72 |
tim |
929 |
|
73 |
|
|
#ifdef IS_MPI |
74 |
|
|
void sortByGlobalIndex(); |
75 |
|
|
#endif |
76 |
|
|
|
77 |
mmeineke |
377 |
SimInfo* entry_plug; |
78 |
tim |
929 |
ofstream dumpFile; |
79 |
tim |
934 |
ofstream finalOut; |
80 |
|
|
|
81 |
|
|
#ifdef IS_MPI |
82 |
tim |
929 |
vector<pair<int, int> > indexArray; |
83 |
tim |
934 |
#endif |
84 |
mmeineke |
377 |
}; |
85 |
|
|
|
86 |
|
|
class StatWriter{ |
87 |
|
|
|
88 |
|
|
public: |
89 |
|
|
StatWriter( SimInfo* the_entry_plug ); |
90 |
|
|
~StatWriter(); |
91 |
|
|
|
92 |
|
|
void writeStat( double currentTime ); |
93 |
tim |
837 |
|
94 |
mmeineke |
377 |
private: |
95 |
tim |
837 |
|
96 |
mmeineke |
377 |
SimInfo* entry_plug; |
97 |
|
|
ofstream outFile; |
98 |
|
|
char outName[500]; |
99 |
|
|
Thermo* tStats; |
100 |
|
|
|
101 |
|
|
}; |
102 |
|
|
|
103 |
|
|
class InitializeFromFile{ |
104 |
|
|
|
105 |
|
|
public: |
106 |
|
|
InitializeFromFile( char *in_name ); |
107 |
|
|
~InitializeFromFile(); |
108 |
|
|
|
109 |
mmeineke |
641 |
void readInit( SimInfo* the_entry_plug ); |
110 |
tim |
837 |
|
111 |
mmeineke |
377 |
private: |
112 |
|
|
char* parseDumpLine(char* line, int atomIndex); |
113 |
tim |
837 |
char* parseCommentLine(char* line, SimInfo* entry_plug); |
114 |
mmeineke |
377 |
FILE *c_in_file; |
115 |
|
|
char c_in_name[500]; |
116 |
mmeineke |
586 |
SimInfo* simnfo; |
117 |
mmeineke |
377 |
}; |
118 |
|
|
|
119 |
gezelter |
637 |
class DumpReader{ |
120 |
mmeineke |
644 |
|
121 |
|
|
public: |
122 |
gezelter |
637 |
DumpReader( char *in_name ); |
123 |
|
|
~DumpReader(); |
124 |
mmeineke |
377 |
|
125 |
gezelter |
637 |
int getNframes(); |
126 |
mmeineke |
646 |
void scanFile( void ); |
127 |
mmeineke |
377 |
|
128 |
mmeineke |
656 |
void getNextFrame() {} |
129 |
|
|
void readFrame(SimInfo* the_simnfo, int whichFrame); |
130 |
|
|
|
131 |
gezelter |
637 |
private: |
132 |
mmeineke |
647 |
|
133 |
chuckv |
657 |
void readSet( int whichFrame ); |
134 |
gezelter |
637 |
char* parseDumpLine(char* line, int atomIndex); |
135 |
mmeineke |
644 |
char* parseCommentLine(char* line, double &time, double boxMat[9] ); |
136 |
mmeineke |
647 |
FILE *inFile; |
137 |
|
|
char inName[500]; |
138 |
|
|
bool isScanned; |
139 |
|
|
int nFrames; |
140 |
tim |
837 |
|
141 |
mmeineke |
804 |
FilePos** frameStart; |
142 |
|
|
FilePos* headFP; |
143 |
mmeineke |
647 |
|
144 |
|
|
SimInfo *simnfo; |
145 |
gezelter |
637 |
}; |
146 |
|
|
|
147 |
|
|
|
148 |
mmeineke |
644 |
|
149 |
mmeineke |
377 |
#endif |