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
	float inMin = S.front().nu;
	float inMax = S.back().nu;

	int nuMin = ceil(inMin);
	int nuMax = 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;

	float 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 < 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(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(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());
	for(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;
}