optics / imie

  #include "globals.h"
  #include <sstream>
  #include <fstream>
  #include <iostream>
  using namespace std;
  
  vector<SpecPair> LoadSpectrum(string filename)
  {
  	//load a spectrum from a file and resample to 2wn intervals
  
  	//create the spectrum
  	vector<SpecPair> S;
  
  	//open the file
  	ifstream inFile(filename.c_str());
  
  
  	SpecPair temp;
  	while(!inFile.eof()){
  		inFile>>temp.nu;
  		inFile>>temp.A;
  		S.push_back(temp);
  	}
  
  	//compute the minimum and maximum input wavenumbers

  	double inMin = S.front().nu;
  	double inMax = S.back().nu;

  	int nuMin = (int)ceil(inMin);
  	int nuMax = (int)floor(inMax);

  
  	//make sure both are either even or odd
  	if(nuMin % 2 != nuMax % 2)
  		nuMax--;
  
  	//compute the number of values in the resampled spectrum
  	int nVals = (nuMax - nuMin)/2 + 1;
  
  	//allocate space for the spectrum
  	vector<SpecPair> outSpec;
  

  	double nu, highVal, lowVal, a;

  	int j=1;
  	for(int i=0; i<nVals; i++)
  	{
  		nu = nuMin + i * 2;
  		temp.nu = nu;
  
  		//handle the boundary cases
  		if(nu < inMin || nu > inMax)
  			temp.A = 0.0;
  		else
  		{
  			//move to the correct position in the input array

  			while(j < (int)S.size()-1 && S[j].nu <= nu)

  				j++;
  
  			lowVal = S[j-1].nu;
  			highVal = S[j].nu;
  			a = (nu - lowVal)/(highVal - lowVal);
  			temp.A = S[j-1].A * (1.0 - a) + S[j].A * a;
  		}
  		outSpec.push_back(temp);
  	}
  
  
  
  	return outSpec;
  }
  
  vector<SpecPair> SetReferenceSpectrum(char* text)
  {
  	stringstream inString(text);
  
  	//create the spectrum
  	vector<SpecPair> S;
  
  	SpecPair temp;
  	while(!inString.eof()){
  		inString>>temp.nu;
  		inString>>temp.A;
  		S.push_back(temp);
  	}
  
  	return S;
  }
  
  void SaveK(string fileName)
  {
  	ofstream outFile(fileName.c_str());

  	for(unsigned int i=0; i<EtaK.size(); i++)

  	{
  		outFile<<EtaK[i].nu<<"          ";
  		outFile<<EtaK[i].A<<endl;
  	}
  	outFile.close();
  }
  
  void SaveN(string fileName)
  {
  	ofstream outFile(fileName.c_str());

  	for(unsigned int i=0; i<EtaN.size(); i++)

  	{
  		outFile<<EtaN[i].nu<<"          ";
  		outFile<<EtaN[i].A<<endl;
  	}
  	outFile.close();
  }
  
  void SaveSimulation(string fileName)
  {
  	ofstream outFile(fileName.c_str());

  	outFile.precision(10);

  	for(unsigned int i=0; i<SimSpectrum.size(); i++)

  	{
  		outFile<<SimSpectrum[i].nu<<"          ";
  		outFile<<SimSpectrum[i].A<<endl;
  	}
  	outFile.close();
  }
  
  void SaveState()
  {
  	ofstream outFile("main.prj");
  	//Window Parameters
  	outFile<<nuMin<<endl;
  	outFile<<nuMax<<endl;
  	outFile<<aMin<<endl;
  	outFile<<aMax<<endl;
  	outFile<<dNu<<endl;
  
  	//material parameters
  	outFile<<radius<<endl;
  	outFile<<baseIR<<endl;
  	outFile<<cA<<endl;
  
  	//optical parameters
  	outFile<<cNAi<<endl;
  	outFile<<cNAo<<endl;
  	outFile<<oNAi<<endl;
  	outFile<<oNAo<<endl;
  
  	outFile.close();
  
  }
  
  void LoadState()
  {
  	ifstream inFile("main.prj");
  	//Window Parameters
  	inFile>>nuMin;
  	inFile>>nuMax;
  	inFile>>aMin;
  	inFile>>aMax;
  	inFile>>dNu;
  
  	//material parameters
  	inFile>>radius;
  	inFile>>baseIR;
  	inFile>>cA;
  
  	//optical parameters
  	inFile>>cNAi;
  	inFile>>cNAo;
  	inFile>>oNAi;
  	inFile>>oNAo;
  
  	inFile.close();
  
  }
  
  void SetDefaults()
  {
  
  	nuMin = 800;
  	nuMax = 4000;
  	dNu = 2;
  
  	aMin = 0;
  	aMax = 1;
  
  
  	//material parameters
  	radius = 4.0f;
  	baseIR = 1.49f;
  	cA = 1.0;
  	vector<SpecPair> KMaterial;
  	vector<SpecPair> NMaterial;
  
  	//optical parameters
  	cNAi = 0.0;
  	cNAo = 0.6;
  	oNAi = 0.0;
  	oNAo = 0.6;
  }