FunctorWriter.cc

#include "FunctorWriter.hh"
#include <fstream> 
#include <map>
#include "PdfBase.hh" 

void writeToFile (PdfBase* pdf, const char* fname) {
  PdfBase::parCont params;
  pdf->getParameters(params); 

  std::ofstream writer;
  writer.open(fname);

  for (PdfBase::parIter p = params.begin(); p != params.end(); ++p) {
    writer << (*p)->name << " " 
	   << (*p)->value << " " 
	   << (*p)->error << " "
	   << (*p)->numbins << " "
	   << (*p)->lowerlimit << " "
	   << (*p)->upperlimit 
	   << std::endl; 
  }

  writer.close(); 
}


void readFromFile (PdfBase* pdf, const char* fname) {
  PdfBase::parCont params;
  pdf->getParameters(params); 

  std::map<string, Variable*> tempMap;
  for (PdfBase::parIter p = params.begin(); p != params.end(); ++p) {
    tempMap[(*p)->name] = (*p); 
  }  

  std::ifstream reader;
  reader.open(fname);
  std::string buffer;
  char discard[1000]; 
  int numSet = 0; 
  while (true) {
    reader >> buffer;
    if (reader.eof()) break; 
    Variable* var = tempMap[buffer];
    if (var) {
      reader >> var->value
	     >> var->error
	     >> var->numbins 
	     >> var->lowerlimit
	     >> var->upperlimit; 
      if (++numSet == tempMap.size()) break; 
    }
    else {
      reader.getline(discard, 1000); 
    }
  }

  reader.close(); 
}

void readListOfNumbers (thrust::host_vector<fptype>& target, const char* fname) {
  std::ifstream reader;
  reader.open(fname);
  fptype buffer = 0; 
  while (true) {
    reader >> buffer;
    if (reader.eof()) break;  
    target.push_back(buffer);
  }
  reader.close(); 
}

void writeListOfNumbers (thrust::host_vector<fptype>& target, const char* fname) {
  std::ofstream writer;
  writer.open(fname);

  for (thrust::host_vector<fptype>::iterator t = target.begin(); t != target.end(); ++t) {
    writer << (*t) << " "; 
  }

  writer.close(); 
}