X-Git-Url: http://git.treefish.org/~alex/phys/latlib.git/blobdiff_plain/d39c90a61334fa4c953842483bb79f3a86d319d0..20acc11c8054200c02794f58e3cc63851b5627ff:/configcache.cpp diff --git a/configcache.cpp b/configcache.cpp index b5496c7..0cabc9d 100644 --- a/configcache.cpp +++ b/configcache.cpp @@ -5,15 +5,24 @@ #include #include +#define HEADER_READOK 0 +#define HEADER_READERR 1 +#define HEADER_READLAST 2 + configcache::configcache(const string& cacheid, const int& nequi, const int& nskip, const string& datadir, char **configmem, const int& configMemSize, const int& cachemode, - writeout *out_a){ - out = out_a; + ostream *_log){ + log = _log; NEQUI = nequi; NSKIP = nskip; DATADIR = datadir; CACHEID = cacheid; + if ( cacheid.find("_") != -1 ) { + if(log) *log << "CCACHE: Invalid cacheid \"" << cacheid << "\" given. Cacheids must not contain underscores!" << endl << flush; + exit(1); + } + configMem = (char*)malloc(configMemSize); tmpConfig = (char*)malloc(configMemSize); @@ -28,11 +37,11 @@ configcache::configcache(const string& cacheid, const int& nequi, const int& nsk refetchDataFiles = false; } -string configcache::getFileId(const bool& shortid) +string configcache::getFileId(int actnequi, const bool& shortid) { stringstream fileid; - if(!shortid) fileid << CACHEID << "_" << NEQUI << "_" << NSKIP; + if(!shortid) fileid << CACHEID << "_" << actnequi << "_" << NSKIP; for(int ipara=0; iparad_name; - if(isValidInFile(filename)) + if(isValidInFile(filename, &filedesc)) { - inFiles.push_back(filename); + inFiles.push_back(filedesc); } } } } -bool configcache::isValidInFile(const string& infile) +bool configcache::isValidInFile(const string& infile, infiledesc *filedesc) { char *inchar, *inParts; string truncIn, truncOut; + filedesc->filename = infile; + if( infile.size() < 4 ) return false; - if( infile.substr(infile.size()-3) != "dat" ) return false; + if( infile.substr(infile.size()-4) == ".dat" ) + filedesc->extended = false; + else if( infile.substr(infile.size()-4) == "edat" ) + filedesc->extended = true; + else + return false; inchar = new char [infile.size()+1]; strcpy (inchar, infile.c_str()); - + inParts = strtok( inchar, "_" ); for(int iPart=0; inParts!=NULL; iPart++) { @@ -75,9 +92,17 @@ bool configcache::isValidInFile(const string& infile) switch(iPart) { - case 1: if(inParts != CACHEID) return false; break; - case 2: if(atoi(inParts) != NEQUI) return false; break; - case 3: if(atoi(inParts) != NSKIP) return false; break; + case 1: if(inParts != CACHEID) + return false; + break; + case 2: + filedesc->nequi = atoi(inParts); + break; + case 3: + if(atoi(inParts) != NSKIP) + return false; + filedesc->nskip = atoi(inParts); + break; } inParts = strtok( NULL, "_"); } @@ -85,14 +110,77 @@ bool configcache::isValidInFile(const string& infile) delete[] inchar; - if( truncIn.find( getFileId(true) + "_" ) == string::npos ) return false; + if( truncIn.find( getFileId(NEQUI, true) + "_" ) == string::npos ) return false; return true; } -bool configcache::readConfig() +int configcache::readHeader() { - if(DATADIR == "" || MODE == CACHE_MODE_DISABLED) return false; + long unsigned int headersize; + + if( readDataToMem((char *)&headersize, sizeof(long unsigned int)) == sizeof(long unsigned int) && inFile.is_open() ) { + if ( headersize == 0 ) + return HEADER_READLAST; + + pair newHeader; + + if( readDataToMem((char *)&newHeader.first, sizeof(unsigned long)) == sizeof(unsigned long) && inFile.is_open() ) { + newHeader.second = malloc(headersize); + + if( readDataToMem((char *)newHeader.second, headersize) == headersize && inFile.is_open() ) { + headerStore.push_back(newHeader); + return HEADER_READOK; + } + else { + if(log) *log << "CCACHE: Could not read heade-data! Closing dat-file: " << openFileDesc.filename << endl << flush; + inFile.close(); + return HEADER_READERR; + } + } + else { + if(log) *log << "CCACHE: Could not read headerid-hash! Closing dat-file: " << openFileDesc.filename << endl << flush; + inFile.close(); + return HEADER_READERR; + } + } + else { + if(log) *log << "CCACHE: Could not read header size. Closing dat-file: " << openFileDesc.filename << endl << flush; + inFile.close(); + return HEADER_READERR; + } +} + +bool configcache::readAllHeaders() +{ + int readHeaderStatus; + + deleteHeaderStore(); + + do { + readHeaderStatus = readHeader(); + } + while ( readHeaderStatus == HEADER_READOK ); + + if ( readHeaderStatus == HEADER_READLAST ) return true; + else if ( readHeaderStatus == HEADER_READERR ) return false; +} + +void * configcache::getHeader(const string& headerid) { + for (vector< pair >::iterator headerStoreIt = headerStore.begin(); headerStoreIt != headerStore.end(); ++headerStoreIt) + if ( headerStoreIt->first == hash(headerid) ) + return headerStoreIt->second; + + return NULL; +} + +/* returns number of equilibration-steps left with new read configuration + or zero if no new configuration was read */ +void configcache::readConfig(bool& readnewconfig, int& nequileft, vector *excludeFileHashes) +{ + readnewconfig = false; + + if(DATADIR == "" || MODE == CACHE_MODE_DISABLED) return; if(refetchDataFiles){ refetchDataFiles = false; @@ -101,56 +189,116 @@ bool configcache::readConfig() while(true) { - if( (!inFile.is_open()) && inFiles.size() == 0 ) return false; + vector::iterator inFileIt = getNextInfile(excludeFileHashes); + int iDidVirtualSkips; - while( (!inFile.is_open()) && inFiles.size() > 0 ) - { - if(out) *out->log << "Opening dat-file: " << inFiles.back() << endl; + if( (!inFile.is_open()) && inFileIt == inFiles.end() ) { + if (readnewconfig) + nequileft = nequileft_internal; + return; + } - inFile.open( (DATADIR + "/" + inFiles.back()).c_str(), std::ios::binary ); - inFiles.pop_back(); + while( (!inFile.is_open()) && inFiles.size() > 0 ) { + openFileDesc = *inFileIt; - if( !inFile.is_open() ) continue; + if (openFileDesc.nequi < NEQUI) + doVirtualEquilibration = true; + else + doVirtualEquilibration = false; - inBuffer = new boost::iostreams::filtering_istreambuf; - inBuffer->push( boost::iostreams::bzip2_decompressor() ); - inBuffer->push(inFile); - } + firstUsedConfig = true; + + if(log) *log << "CCACHE: Opening dat-file: " << inFileIt->filename << endl << flush; + inFile.open( (DATADIR + "/" + inFileIt->filename).c_str(), std::ios::binary ); + + inFiles.erase(inFileIt); + + if( !inFile.is_open() ) continue; + + inBuffer = new boost::iostreams::filtering_istreambuf; + inBuffer->push( boost::iostreams::bzip2_decompressor() ); + inBuffer->push(inFile); + } if( inFile.is_open() ) { - int readturn = -1; - try { readturn = boost::iostreams::read(*inBuffer, tmpConfig, configSize); } - catch (std::exception const& ex) { - if(out) *out->log << "Caught exception: " << ex.what() << endl; - inFile.close(); + if (doVirtualEquilibration) { + if(log) *log << "CCACHE: Trying virtual equilibration." << endl << flush; + doVirtualEquilibration = false; + for (iDidVirtualSkips=0; iDidVirtualSkips < (NEQUI-openFileDesc.nequi)/openFileDesc.nskip; iDidVirtualSkips++) { + if( readFullBlock(tmpConfig, configSize) != configSize || ! inFile.is_open() ) + break; + else if ( (NEQUI-openFileDesc.nequi) - (iDidVirtualSkips+1)*openFileDesc.nskip < nequileft_internal ) { + memcpy(configMem, tmpConfig, configSize); + nequileft_internal = NEQUI - openFileDesc.nequi - (iDidVirtualSkips+1)*openFileDesc.nskip; + readnewconfig = true; + firstUsedConfig = false; + } + } } - - if( readturn == configSize && inFile.is_open() ) + + if( readFullBlock(tmpConfig, configSize) == configSize && inFile.is_open() ) { memcpy(configMem, tmpConfig, configSize); - return true; + readnewconfig = true; + if (firstUsedConfig) { + firstUsedConfig = false; + if (openFileDesc.nequi < NEQUI) + nequileft_internal = NEQUI - openFileDesc.nequi - iDidVirtualSkips*openFileDesc.nskip; + else + nequileft_internal = NEQUI - openFileDesc.nequi; + } + nequileft_internal -= openFileDesc.nskip; + nequileft = nequileft_internal; + return; } - else inFile.close(); + else { + if(log) *log << "CCACHE: Could not read configuration. Closing dat-file: " << openFileDesc.filename << endl << flush; + inFile.close(); + } } } } -void configcache::writeConfig() +void configcache::openOutFile(int actnequi) { + time_t secstamp = time(NULL); + + outFileName.str(""); + outFileName << DATADIR << "/" << secstamp << "_" << getFileId(actnequi) << "_.edat.tmp"; + + outFile.open( outFileName.str().c_str(), std::ios::binary ); + + outBuffer = new boost::iostreams::filtering_ostreambuf; + outBuffer->push(boost::iostreams::bzip2_compressor()); + outBuffer->push(outFile); +} + +void configcache::writeHeader(const string& headerid, const char *header, long unsigned int size, int actnequi) { + unsigned long headeridhash; + if( DATADIR == "" || MODE < 2 ) return; - if(!outFile.is_open()){ - time_t secstamp = time(NULL); + if(!outFile.is_open()) + openOutFile(actnequi); - outFileName.str(""); - outFileName << DATADIR << "/" << secstamp << "_" << getFileId() << "_.dat.tmp"; - outFile.open( outFileName.str().c_str(), std::ios::binary ); + headeridhash = hash(headerid); - outBuffer = new boost::iostreams::filtering_ostreambuf; - outBuffer->push(boost::iostreams::bzip2_compressor()); - outBuffer->push(outFile); - } + boost::iostreams::write(*outBuffer, (char*)&size, sizeof(long unsigned int)); + boost::iostreams::write(*outBuffer, (char*)&headeridhash, sizeof(unsigned long)); + boost::iostreams::write(*outBuffer, header, size); +} + +void configcache::writeConfig(int actnequi) +{ + long unsigned int zeroheader=0; + + if ( DATADIR == "" || MODE < 2 ) return; + + if ( ! outFile.is_open() ) + openOutFile(actnequi); + + boost::iostreams::write(*outBuffer, (char*)&zeroheader, sizeof(long unsigned int)); boost::iostreams::write(*outBuffer, configMem, configSize); } @@ -175,6 +323,7 @@ void configcache::setPara(const string& parid, const double& value){ inFiles.clear(); refetchDataFiles = true; + nequileft_internal = NEQUI; } configcache::~configcache() @@ -198,3 +347,74 @@ void configcache::finishOutFile() rename( outFileName.str().c_str(), outFileName.str().substr(0, outFileName.str().size()-4).c_str() ); } } + +int configcache::readFullBlock(char *tmpData, long unsigned int dataSize) +{ + /* try to read header */ + if ( openFileDesc.extended ) + if ( ! readAllHeaders() ) + return -1; + + /* read data */ + return readDataToMem(tmpData, dataSize); +} + +int configcache::readDataToMem(char *tmpData, long unsigned int dataSize) +{ + int readturn = -1; + + if ( dataSize == 0 ) return 0; + + try { readturn = boost::iostreams::read(*inBuffer, tmpData, dataSize); } + catch(boost::iostreams::bzip2_error& error) { + if(log) *log << "CCACHE: Caught bzip2 exception with error code: " << error.error() << endl << flush; + inFile.close(); + } + catch (std::exception const& ex) { + if(log) *log << "CCACHE: Caught exception: " << ex.what() << endl << flush; + inFile.close(); + } + catch( ... ) { + if(log) *log << "CCACHE: Caught unknown exception while reading." << endl << flush; + inFile.close(); + } + + return readturn; +} + +unsigned long configcache::hash(const string& str) +{ + unsigned long hash = 5381; + + for(string::const_iterator it=str.begin();it!=str.end();it++) + hash = ((hash << 5) + hash) + *it; /* hash * 33 + character */ + + return hash; +} + +void configcache::deleteHeaderStore() +{ + while ( headerStore.size() > 0 ) { + free(headerStore.back().second); + headerStore.pop_back(); + } +} + +vector::iterator configcache::getNextInfile(vector *excludeFileHashes) { + for (vector::iterator init = inFiles.begin(); init != inFiles.end(); ++init) { + if (excludeFileHashes != NULL) { + bool excludethisfile = false; + + for (vector::iterator exit = excludeFileHashes->begin(); exit != excludeFileHashes->end(); ++exit) + if ( *exit == hash(init->filename) ) { + excludethisfile = true; + break; + } + + if (excludethisfile) + continue; + } + return init; + } + return inFiles.end(); +}