Work Log by Kevin Beason

void nkToRefl(const char *file){
  ifstream inFile;
  inFile.open(file);
  if (!inFile){
    cerr << "Error opening " << file << "\n";
    return;
  }
  string str;
  getline(inFile, str);
  getline(inFile, str);
  fprintf(stderr, "%8s %30s reflectance = ", file, str.substr(18).c_str());
 
  double eV, n, k, r;
  double h = 6.626068e-34;       // m^2 kg/s
  double c = 299792458;          // m/2
  double ec = 1.60217646e-19;    // C
 
  map reflsMap;
  while (!inFile.eof()){
    inFile >> eV >> n >> k >> r;
    double nrg = eV * ec;          // V = J/C => J = VC
    double wavelen = h*c/nrg;      // nrg = h*c/l  => l=h*c/nrg
    double nm = wavelen*1e9;
    reflsMap[nm] = r;
  }
  inFile.close();
 
  vector wavelens, refls;
  for (map::iterator i=reflsMap.begin(); i!=reflsMap.end(); i++){
    wavelens.push_back(i->first);
    refls.push_back(i->second);
  }
 
  Spectrum lightSpec(6500);                // start out with light
  Spectrum reflSpec(&wavelens[0], &refls[0], wavelens.size());
  lightSpec.scaleSpec(1/4.73223e+06);      // normalize by red channel
  reflSpec.scaleSpec(lightSpec);           // mult light by refl
  SbColor rgb = reflSpec.getRGB();
  fprintf(stderr, "(%5f, %5f, %5f)\n", rgb[0], rgb[1], rgb[2]);
}