]> git.treefish.org Git - phys/latlib.git/blob - hypercache.cpp
implemented partial virtual equilibration.
[phys/latlib.git] / hypercache.cpp
1 #include "hypercache.h"
2
3 configcache *hypercache::O = NULL;
4 configcache *hypercache::C = NULL;
5 vector<hypercache::para> hypercache::delayedParaAdd;
6 vector<hypercache::para> hypercache::delayedParaSet;
7 string hypercache::activeCFile = "";
8 vector<unsigned long> hypercache::parentConfigs;
9 writeout *hypercache::out = NULL;
10
11 void hypercache::initCache(configcache **cache,
12                            const string& cacheid, const int& nequi, const int& nskip, const string& datadir, 
13                            char **configmem, const int& configMemSize, const int& cachemode, writeout *out_a) {
14   
15   bool allInitBefore = true;
16   
17   if (O == NULL || C == NULL)
18     allInitBefore = false;
19
20   if ( out_a != NULL )
21     out = out_a;
22
23   *cache = new configcache(cacheid, nequi, nskip, datadir, configmem, configMemSize, cachemode, out_a);
24
25   if (O != NULL && C != NULL && !allInitBefore) {
26     for (vector<para>::iterator parit=delayedParaAdd.begin(); parit != delayedParaAdd.end(); ++parit) {
27       O->addPara(parit->parid, parit->val);
28       C->addPara(parit->parid, parit->val);
29     }
30     for (vector<para>::iterator parit=delayedParaSet.begin(); parit != delayedParaSet.end(); ++parit) {
31       O->setPara(parit->parid, parit->val);
32       C->setPara(parit->parid, parit->val);
33     }
34   }
35 }
36
37 void hypercache::addPara(const string& parid, const double& val) {
38   if (O == NULL || C == NULL) {
39     para newpara; 
40     newpara.parid = parid;
41     newpara.val = val;
42     delayedParaAdd.push_back(newpara);
43   }
44   else {
45     O->addPara(parid, val);
46     C->addPara(parid, val);
47   }
48 }
49
50 void hypercache::setPara(const string& parid, const double& val) {
51   if (O == NULL || C == NULL) {
52     para newpara; 
53     newpara.parid = parid;
54     newpara.val = val;
55     delayedParaSet.push_back(newpara);
56   }
57   else {
58     O->setPara(parid, val);
59     C->setPara(parid, val);
60     activeCFile = "";
61     parentConfigs.clear();
62   }
63 }
64
65 void hypercache::finalize() {
66   delete C;
67   delete O;
68 }
69
70 string hypercache::fileOfPath(const string& dressedfile) {
71   return dressedfile.substr(dressedfile.find_last_of("\\/")+1);
72 }
73
74 int hypercache::readC() { 
75   int readret;
76
77   readret = C->readConfig(&parentConfigs);
78
79   if ( readret == -1 )
80     activeCFile = fileOfPath(C->getInFileName());
81   else
82     activeCFile = "";
83
84   return readret; 
85 }
86
87 void hypercache::writeC() {
88   C->writeConfig();
89   activeCFile = fileOfPath( C->getOutFileName().substr( 0, C->getOutFileName().length()-4 ) );
90 }
91
92 void hypercache::writeO() {
93   if ( activeCFile != "" ) {
94     unsigned long afilehash = configcache::hash(activeCFile);
95     O->writeHeader("concurrent_cfile", (char*)(&afilehash), sizeof(unsigned long));
96   }
97
98   O->writeConfig(); 
99 }
100
101 int hypercache::readO() {
102   int readret;
103
104   readret = O->readConfig();
105
106   if ( readret == -1 ) {
107     unsigned long *parentconfig = (unsigned long*)O->getHeader("concurrent_cfile");
108     if ( parentconfig != NULL )
109       addParentConfig(parentconfig);
110   }
111   
112   return readret;
113 }
114
115 void hypercache::addParentConfig(const unsigned long *parentconfig) {
116   for (vector<unsigned long>::iterator parit = parentConfigs.begin(); parit != parentConfigs.end(); ++parit)
117     if ( *parit == *parentconfig )
118       return;
119
120   parentConfigs.push_back(*parentconfig);
121 }