/**
	@file SolarSourceModeling.h
	@brief The SolarSourceModeling Class header file.

	@datecreated April 2, 2008

	@version 1.0

	@author
		- Kelly Teague
		- Greg Paxton

	@copyright
		(©) Copyright 2008 by GATS Inc.
		11864 Canon Blvd., Suite 101, Newport News, VA 23606

		All Rights Reserved. No part of this software or publication may be
		reproduced, stored in a retrieval system, or transmitted, in any form
		or by any means, electronic, mechanical, photocopying, recording, or
		otherwise without the prior written permission of GATS Inc.

	@see		SolarSourceModeling
	@bug	None known
	@todo	None known
	@history

*/

#ifndef SOLARSOURCEMODELING_H_
#define SOLARSOURCEMODELING_H_

#include <iostream>
#include <vector>


typedef std::vector<std::vector<double> > Vector2D;

class SolarSourceModeling
{
public:

	SolarSourceModeling(
								const bool mode,
								const double balanceStart,
								const double balanceEnd,
								const std::vector<double> elevationLockdown,
								const std::vector<double> azimuthLockdown,
								const std::vector<double> signal,
								const std::vector<double> time,
								const double calibrationStart,
								const double calibrationEnd,
								const unsigned int offset
							);
				/*!< @brief Constructor for the SolarSourceModeling object. */
	inline void getCoefficients(double *a, double *b, double *c) {*a = _coefA; *b = _coefB; *c = _coefC; } ;
				/*!< @brief Method that records the a, b, and c model values in the passed in variables. */
	double getIntensity(const double elevation, const double azimuth) const;
				/*!< @brief Method for retrieving model intensities from elevation and azimuth input data. */
	virtual ~SolarSourceModeling();
				/*!< @brief Destructor for the SolarSourceModeling object */

	inline bool good(void) { return isTheDataOk; };

private:

	double centerElSolarCoordinates(double hiX, double loX);
	double centerAzSolarCoordinates(double hiY, double loY);
	double sum(std::vector<double> data);
	double sumOfSquares(std::vector<double> data);
	double sumCombined1d(std::vector<double> xData, std::vector<double> yData);
	double getMaxValue(std::vector<double> data);

	void init(void);
	void validateDataValues(void);
	void setCalibrationData();
	void getCruciformSubsetData(void);
	void setCalibrationTimes(void);
	void getCalibrationIndices(std::vector<double> data);
	void setCoefficients(void);
	void getOffsetData(void);

	std::vector<double> getCruciformStart(void);
	double* dot(double X[][3], double *Y);

	int getAzStart(void);
	void getCruciformSubsections(void);
	std::vector<double> normalizeIntensity(void);
	int getIndexOfNearestIntensity(double x0, double y0, std::vector<double> xData, std::vector<double> yData);

	bool isTheDataOk;
	bool sunsetMode;
	unsigned int _offsetData;
	double _balStart, _balEnd, _calStart, _calEnd;
	std::vector<double> _det;
	std::vector<double> _t;

	double _coefA;
	double _coefB;
	double _coefC;

	std::vector<double> _elLockdown;
	std::vector<double> _azLockdown;

	//these are used for normalization of subsection data
	std::vector<double> _azData1;
	std::vector<double> _azData2;
	std::vector<double> _azData3;

	std::vector<double> _elData1;
	std::vector<double> _elData2;
	std::vector<double> _elData3;

	std::vector<double> _intData1;
	std::vector<double> _intData2;
	std::vector<double> _intData3;

	double _relCalStartTime;
	double _relCalEndTime;
	int calibrationStartIndex;
	int calibrationEndIndex;
	
 };
 

#endif /*SOLARSOURCEMODELING_H_*/