CSimulation.cpp

/*
 * File:   CSimulation.cpp
 * Author: Mark Wittekind and Drew Murray
 *
 * Created on March 6, 2015, 11:33 AM
 */

#include "CSimulation.h"
#include "CObject.h"
#include "SimulationParametersStruct.cpp"

CSimulation::CSimulation(SimulationParameters *SimPar) {
	m_simPar = SimPar;
}

void CSimulation::addObject(CObject* obj){
	m_objects.push_back(obj);
}


CSimulation::~CSimulation() {
}

float CSimulation::calcMag(){
	float magTot = 0.0;
	for (vector<CObject*>::iterator it2 = m_objects.begin(); it2 != m_objects.end(); it2++){
		magTot += (*it2)->calcMag();
	}
	return magTot;
}

float CSimulation::calcEnergy(){
	float enTot = 0.0;
	for (vector<CObject*>::iterator it2 = m_objects.begin(); it2 != m_objects.end(); it2++){
		enTot += (*it2)->calcEnergy();
	}
	return enTot;
}

void CSimulation::run(ostream& pos_stream, ostream& mag_stream, SimulationParameters *SimPar){

	//mystery line to allow virtual functions to work as they should (this variable is never used)
	//CLattice uselessvariable(m_simPar->WIDTH,m_simPar->HEIGHT);

	for (vector<CObject*>::iterator it2 = m_objects.begin(); it2 != m_objects.end(); it2++){
				mag_stream << m_simPar->temperature << "\n";
				(*it2)->print_all(pos_stream, mag_stream, m_simPar);
			}

	for(unsigned int t = 0; t<m_simPar->MAXTIMESTEPS;t++)
	{
		for (vector<CObject*>::iterator it = m_objects.begin(); it != m_objects.end(); it++){
			(*it)->update(m_simPar);
		}
		if((t%SimPar->PRINT_FREQ==m_simPar->PRINT_FREQ-1) || ((t+1)>=(m_simPar->MAXTIMESTEPS)))
		{
			cout << "printing timestep " << t << " of " << SimPar->MAXTIMESTEPS << "\n";
			for (vector<CObject*>::iterator it2 = m_objects.begin(); it2 != m_objects.end(); it2++){
				(*it2)->print_all(pos_stream,mag_stream, m_simPar);
			}
		}
	}
	mag_stream<<"\n";
	pos_stream<<"TEMP CHANGED\n";
		return;
    }