Fix tickmarks plot

Benjamin Renard
1 parent 413a6873
Showing 5 changed files with 103 additions and 138 deletions Show diff stats
src/ParamOutputImpl/Plot/ParameterData.cc
src/ParamOutputImpl/Plot/ParameterData.hh
src/ParamOutputImpl/Plot/PlotOutput.cc
src/ParamOutputImpl/Plot/TickPlot/TickMarkDecorator.cc
src/ParamOutputImpl/Plot/TickPlot/TickMarkDecorator.hh
@@ -14,7 +14,7 @@
 namespace plot {
  
 ParameterData::ParameterData() :
-		_minTime(nan("")), _maxTime(nan("")), _dim1Size(-1), _dim2Size(-1) {
+		_minTime(nan("")), _maxTime(nan("")), _dim1Size(-1), _dim2Size(-1), _paramGapSize(0.) {
 }
  
 ParameterData::~ParameterData() {
@@ -44,7 +44,7 @@ void ParameterData::preAllocate(int nbData)
  * @brief Adds a value by time to parameter
  */
 void ParameterData::addValue(double time, double value,
-		double sampling, double gapSize, bool &gapDetected, AMDA::Common::ParameterIndexComponent& index) {
+		double sampling, bool &gapDetected, AMDA::Common::ParameterIndexComponent& index) {
 	int valueIndex = componentToValuesIndex(index);
 	if (valueIndex >= (int)_values.size())
 	{
@@ -53,7 +53,7 @@ void ParameterData::addValue(double time, double value,
 	// add time if not already in
 	if (_dates.empty() || _dates.back() != time) {
  
-		if (!_dates.empty() && ((time - _dates.back()) > gapSize))
+		if (!_dates.empty() && ((time - _dates.back()) > _paramGapSize))
 		{
 			//gap detected => push a NaN value in data list
 			if (_dates.back() + sampling > time)
@@ -90,7 +90,7 @@ public:
 	/**
 	 * @brief Adds a value by time to parameter
 	 */
-	void addValue(double time_, double value_, double sampling, double gapSize,
+	void addValue(double time_, double value_, double sampling,
 			bool &gapDetected, AMDA::Common::ParameterIndexComponent& index_);
  
 	/**
@@ -209,6 +209,22 @@ public:
 		_dim2Size = dim2Size;
 	}
  
+	/*
+	 * @brief Get param gap size (in seconds)
+	 */
+	double getParamGapSize()
+	{
+		return _paramGapSize;
+	}
+
+	/*
+	 * @brief Set param gap size
+	 */
+	void setParamGapSize(double gapSize)
+	{
+		_paramGapSize = gapSize;
+	}
+
 	/**
 	 * Gets the index of the given time (or just after) in time array,
 	 * returns -1 if not found
@@ -289,6 +305,11 @@ private:
 	 * @brief Dim 2 Size
 	 */
 	int _dim2Size;
+
+	/*
+	 * @brief Param gap size (in seconds)
+	 */
+	double _paramGapSize;
 };
  
  
@@ -623,6 +623,7 @@ void PlotOutput::visit(
 	bool gapDetected = false;
 	//ParamDataScalaireShort
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
 	AMDA::Common::ParameterIndexComponent componentIndex = AMDA::Common::ParameterIndexComponent(-1,-1);
  
 	crtParameterData.preAllocate(_paramDataIndexInfo._nbDataToProcess
@@ -636,7 +637,7 @@ void PlotOutput::visit(
 			pParamData->getTime(index),
 			pParamData->get(index),
 			pParamData->getMinSampling(),
-			gapSize, gapDetected, componentIndex);
+			gapDetected, componentIndex);
 	}
 }
  
@@ -649,6 +650,7 @@ void PlotOutput::visit(
 	bool gapDetected = false;;
 	//ParamDataScalaireFloat
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
 	AMDA::Common::ParameterIndexComponent componentIndex = AMDA::Common::ParameterIndexComponent(-1,-1);
  
 	crtParameterData.preAllocate(_paramDataIndexInfo._nbDataToProcess
@@ -663,7 +665,7 @@ void PlotOutput::visit(
 			pParamData->getTime(index),
 			pParamData->get(index),
 			pParamData->getMinSampling(),
-			gapSize, gapDetected, componentIndex);
+			gapDetected, componentIndex);
 	}
 }
  
@@ -676,6 +678,7 @@ void PlotOutput::visit(
 	bool gapDetected = false;;
 	//ParamDataScalaireDouble
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
 	AMDA::Common::ParameterIndexComponent componentIndex = AMDA::Common::ParameterIndexComponent(-1,-1);
  
 	crtParameterData.preAllocate(_paramDataIndexInfo._nbDataToProcess
@@ -689,7 +692,7 @@ void PlotOutput::visit(
 			pParamData->getTime(index),
 			pParamData->get(index),
 			pParamData->getMinSampling(),
-			gapSize, gapDetected, componentIndex);
+			gapDetected, componentIndex);
 	}
 }
  
@@ -702,6 +705,7 @@ void PlotOutput::visit(
 	bool gapDetected = false;;
 	//ParamDataScalaireLongDouble
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
 	AMDA::Common::ParameterIndexComponent componentIndex = AMDA::Common::ParameterIndexComponent(-1,-1);
  
 	crtParameterData.preAllocate(_paramDataIndexInfo._nbDataToProcess
@@ -715,7 +719,7 @@ void PlotOutput::visit(
 			pParamData->getTime(index),
 			pParamData->get(index),
 			pParamData->getMinSampling(),
-			gapSize, gapDetected, componentIndex);
+			gapDetected, componentIndex);
 	}
 }
  
@@ -728,6 +732,7 @@ void PlotOutput::visit(
 	bool gapDetected = false;;
 	//ParamDataScalaireInt
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
 	AMDA::Common::ParameterIndexComponent componentIndex = AMDA::Common::ParameterIndexComponent(-1,-1);
  
 	crtParameterData.preAllocate(_paramDataIndexInfo._nbDataToProcess
@@ -741,7 +746,7 @@ void PlotOutput::visit(
 			pParamData->getTime(index),
 			pParamData->get(index),
 			pParamData->getMinSampling(),
-			gapSize, gapDetected, componentIndex);
+			gapDetected, componentIndex);
 	}
 }
  
@@ -760,6 +765,7 @@ void PlotOutput::visit(AMDA::Parameters::ParamDataTab1DShort *pParamData) {
 	bool gapDetected = false;;
 	//ParamDataTab1DShort
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
  
 	//init parameter values container
 	if (crtParameterData._indexes.size() == 0)
@@ -777,7 +783,7 @@ void PlotOutput::visit(AMDA::Parameters::ParamDataTab1DShort *pParamData) {
 			crtParameterData.addValue(time,
 				pParamData->get(i)[index.getDim1Index()],
 				pParamData->getMinSampling(),
-				gapSize, gapDetected, index);
+				gapDetected, index);
 		}
 	}
 }
@@ -791,6 +797,7 @@ void PlotOutput::visit(
 	bool gapDetected = false;;
 	//ParamDataTab1DFloat
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
  
 	//init parameter values container
 	if (crtParameterData._indexes.size() == 0)
@@ -808,7 +815,7 @@ void PlotOutput::visit(
 			crtParameterData.addValue(time,
 				pParamData->get(i)[index.getDim1Index()],
 				pParamData->getMinSampling(),
-				gapSize, gapDetected, index);
+				gapDetected, index);
 		}
 	}
 }
@@ -822,6 +829,7 @@ void PlotOutput::visit(
 	bool gapDetected = false;;
 	//ParamDataTab1DDouble
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
  
 	//init parameter values container
 	if (crtParameterData._indexes.size() == 0)
@@ -839,7 +847,7 @@ void PlotOutput::visit(
 			crtParameterData.addValue(time,
 				pParamData->get(i)[index.getDim1Index()],
 				pParamData->getMinSampling(),
-				gapSize, gapDetected, index);
+				gapDetected, index);
 		}
 	}
 }
@@ -853,6 +861,7 @@ void PlotOutput::visit(
 	bool gapDetected = false;;
 	//ParamDataTab1DLongDouble
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
  
 	//init parameter values container
 	if (crtParameterData._indexes.size() == 0)
@@ -870,7 +879,7 @@ void PlotOutput::visit(
 			crtParameterData.addValue(time,
 				pParamData->get(i)[index.getDim1Index()],
 				pParamData->getMinSampling(),
-				gapSize, gapDetected, index);
+				gapDetected, index);
 		}
 	}
 }
@@ -883,6 +892,7 @@ void PlotOutput::visit(AMDA::Parameters::ParamDataTab1DInt *pParamData) {
 	bool gapDetected = false;;
 	//ParamDataTab1DInt
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
  
 	//init parameter values container
 	if (crtParameterData._indexes.size() == 0)
@@ -900,7 +910,7 @@ void PlotOutput::visit(AMDA::Parameters::ParamDataTab1DInt *pParamData) {
 			crtParameterData.addValue(time,
 				pParamData->get(i)[index.getDim1Index()],
 				pParamData->getMinSampling(),
-				gapSize, gapDetected, index);
+				gapDetected, index);
 		}
 	}
 }
@@ -920,6 +930,7 @@ void PlotOutput::visit(AMDA::Parameters::ParamDataTab2DShort *pParamData) {
 	bool gapDetected = false;;
 	//ParamDataTab2DShort
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
  
 	//init parameter values container
 	if (crtParameterData._indexes.size() == 0)
@@ -937,7 +948,7 @@ void PlotOutput::visit(AMDA::Parameters::ParamDataTab2DShort *pParamData) {
 			crtParameterData.addValue(time,
 				pParamData->get(i)[index.getDim1Index()][index.getDim2Index()],
 				pParamData->getMinSampling(),
-				gapSize, gapDetected, index);
+				gapDetected, index);
 		}
 	}
 }
@@ -950,6 +961,7 @@ void PlotOutput::visit(AMDA::Parameters::ParamDataTab2DFloat *pParamData) {
 	bool gapDetected = false;;
 	//ParamDataTab2DFloat
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
  
 	//init parameter values container
 	if (crtParameterData._indexes.size() == 0)
@@ -967,7 +979,7 @@ void PlotOutput::visit(AMDA::Parameters::ParamDataTab2DFloat *pParamData) {
 			crtParameterData.addValue(time,
 				pParamData->get(i)[index.getDim1Index()][index.getDim2Index()],
 				pParamData->getMinSampling(),
-				gapSize, gapDetected, index);
+				gapDetected, index);
 		}
 	}
 }
@@ -980,6 +992,7 @@ void PlotOutput::visit(AMDA::Parameters::ParamDataTab2DDouble *pParamData) {
 	bool gapDetected = false;;
 	//ParamDataTab2DDouble
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
  
 	//init parameter values container
 	if (crtParameterData._indexes.size() == 0)
@@ -997,7 +1010,7 @@ void PlotOutput::visit(AMDA::Parameters::ParamDataTab2DDouble *pParamData) {
 			crtParameterData.addValue(time,
 				pParamData->get(i)[index.getDim1Index()][index.getDim2Index()],
 				pParamData->getMinSampling(),
-				gapSize, gapDetected, index);
+				gapDetected, index);
 		}
 	}
 }
@@ -1010,6 +1023,7 @@ void PlotOutput::visit(AMDA::Parameters::ParamDataTab2DLongDouble *pParamData) {
 	bool gapDetected = false;;
 	//ParamDataTab2DLongDouble
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
  
 	//init parameter values container
 	if (crtParameterData._indexes.size() == 0)
@@ -1027,7 +1041,7 @@ void PlotOutput::visit(AMDA::Parameters::ParamDataTab2DLongDouble *pParamData) {
 			crtParameterData.addValue(time,
 				pParamData->get(i)[index.getDim1Index()][index.getDim2Index()],
 				pParamData->getMinSampling(),
-				gapSize, gapDetected, index);
+				gapDetected, index);
 		}
 	}
 }
@@ -1040,6 +1054,7 @@ void PlotOutput::visit(AMDA::Parameters::ParamDataTab2DInt *pParamData) {
 	bool gapDetected = false;;
 	//ParamDataTab2DInt
 	ParameterData& crtParameterData = _parameterValues[_currentParamId];
+	crtParameterData.setParamGapSize(gapSize);
  
 	//init parameter values container
 	if (crtParameterData._indexes.size() == 0)
@@ -1057,7 +1072,7 @@ void PlotOutput::visit(AMDA::Parameters::ParamDataTab2DInt *pParamData) {
 			crtParameterData.addValue(time,
 				pParamData->get(i)[index.getDim1Index()][index.getDim2Index()],
 				pParamData->getMinSampling(),
-				gapSize, gapDetected, index);
+				gapDetected, index);
 		}
 	}
 }
@@ -55,7 +55,7 @@ TickMarkDecorator&amp; TickMarkDecorator::operator=(const TickMarkDecorator&amp; ref_) {
 	} else {
 		_pTimeInfo->_timeFormat = ref_._pTimeInfo->_timeFormat;
 		_pTimeInfo->_timeTitle = ref_._pTimeInfo->_timeTitle;
-		_pTimeInfo->_min = ref_._pTimeInfo->_min;
+		_pTimeInfo->_isFirstLabel = ref_._pTimeInfo->_isFirstLabel;
 		_pTimeInfo->_maxLen = ref_._pTimeInfo->_maxLen;
 	}
 	_seriesInfoList = ref_._seriesInfoList;
@@ -342,91 +342,6 @@ void TickMarkDecorator::draw(PanelPlotOutput* pplot_, Axis* axis_,
 	restoreColor(pls_, lInitialColor, pplot_->_panel->_page->_mode);
  
 }
-/**
- * @brief Search for each tick time, the closest parameter time in the paramdata_ and store the corresponding value.
- * Selected parameter time is in range ]ticktime-step, ticktime+step[ where step is the space between two tick.
- * If no such time is found for a given tick, associated value is marked as NaN.
- * Note that this algorithm may retrieve same value twice if there is exactly one time between two ticks
- * and it is exactly at the middle.
- * @param seriesdata_ series data to fill
- * @param pAxis_ related axis
- * @param paramdata_ related parameter
- * @param serieIndex_ series index.
- * @return longest label to display for that series (in number of char).
- */
-unsigned int TickMarkDecorator::fillSeriesData(SeriesData& seriesdata_,
-		TimeAxis* pAxis_, ParameterData& paramdata_, AMDA::Common::ParameterIndexComponent serieIndex_) {
-	LOG4CXX_DEBUG(_logger, " filling series data...");
-	// string used to compute formattedvalue max length
-	char formattedValue [40];
-	unsigned int maxlen = 0;
-
-	// time range :
-	double start = pAxis_->getRange().getMin();
-	double stop = pAxis_->getRange().getMax();
-	double step = getMajorTickSpace(pAxis_, start, stop);
-	LOG4CXX_TRACE(_logger,
-			" for all ticks between " << start << " and " << stop << " spaced by " << step);
-
-	int index = -1;
-	int size = paramdata_.getSize();
-
-	for (double tickTime = start; tickTime <= stop; tickTime += step) {
-		// search index closest value to that time :
-		double* ptime = paramdata_.getTimes();
-		int i = index + 1;
-		while (i < size && (*(ptime + i) < tickTime)) {
-			i++;
-		}
-
-		double cur = (i >= size) ? -1 : *(ptime + i);
-		double prev = (i == 0) ? -1 : *(ptime + i - 1);
-
-		if (prev != -1 && prev <= (tickTime - step)) {
-			prev = -1;
-		}
-		if (cur != -1 && cur >= (tickTime + step)) {
-			cur = -1;
-		}
-
-		if (prev == -1) {
-			// only cur can be eligible
-			if (cur == -1) {
-				index = -1; // no suitable value
-			} else {
-				index = i;
-			}
-		} else if (cur == -1) {
-			// only prev can be eligible and it is not -1 due to previous test
-			index = i - 1;
-		} else {
-			// both are eligible, choose closest one :
-			if ((tickTime - prev) > (cur - tickTime)) {
-				index = i;
-			} else {
-				index = i - 1;
-			}
-		}
-
-		// time found, retrieve related value or put a nan one:
-		if (index == -1) {
-			LOG4CXX_TRACE(_logger,
-					"updating series data at time:"<<tickTime<<" with value:(NaN)");
-			seriesdata_[tickTime] = nan("NaN");
-			// restore last index for next loop :
-			index = i;
-		} else {
-			double value = *(paramdata_.getValues(serieIndex_) + index);
-			LOG4CXX_TRACE(_logger,
-					"updating series data at time:"<<tickTime<<" with value :("<<value<<")");
-			// compute formatted value length :
-			const unsigned int len = snprintf(formattedValue, 40, getValuesFormat().c_str(), value);
-			maxlen = std::max((const unsigned int) maxlen, len);
-			seriesdata_[tickTime] = value;
-		}
-	}
-	return std::min((const unsigned int) maxlen, (const unsigned int) 40);
-}
  
 void TickMarkDecorator::installLabelGenerator(PanelPlotOutput* /*pplot_*/, TimeAxis* timeAxis_) {
  
@@ -436,7 +351,7 @@ void TickMarkDecorator::installLabelGenerator(PanelPlotOutput* /*pplot_*/, TimeA
 	// install time format
 	_pTimeInfo->_timeFormat = getPlFormat();
 	_pTimeInfo->_timeTitle = timeAxis_->_legend._text;
-	_pTimeInfo->_min = timeAxis_->getRange().getMin();
+	_pTimeInfo->_isFirstLabel = true;
  
 	unsigned int maxNameLen = _pTimeInfo->_timeTitle.size();
 	// init maximum value len with formatted time :
@@ -473,15 +388,18 @@ void TickMarkDecorator::installLabelGenerator(PanelPlotOutput* /*pplot_*/, TimeA
 			if (pSerie->_name.size() > maxNameLen) {
 				maxNameLen = pSerie->_name.size();
 			}
-			pSerie->_min = timeAxis_->getRange().getMin();
+			pSerie->_isFirstLabel = true;
 			pSerie->_valuesFormat = getValuesFormat();
-			// get serie at index :
-			SeriesData data;
-			_maxValueLen = std::max((const unsigned int) _maxValueLen,
-					fillSeriesData(data, timeAxis_,
-							(*_pParameterValues)[p.getYSeriePropertiesAt(q).getParamId()],
-							q));
-			pSerie->_data = data;
+
+			//Compute max len
+			pSerie->_parameterData = &((*_pParameterValues)[p.getYSeriePropertiesAt(q).getParamId()]);
+			pSerie->_indexComponents = AMDA::Common::ParameterIndexComponent(q.getDim1Index(),q.getDim2Index());
+
+			for (int i = 0; i < pSerie->_parameterData->getSize(); ++i) {
+				const unsigned int len = snprintf(formattedStr, 40, getValuesFormat().c_str(), *pSerie->_parameterData->getValues(pSerie->_indexComponents, i));
+				_maxValueLen = std::max(_maxValueLen, len);
+			}
+
 			_seriesInfoList.push_back(pSerie);
 		}
 	}
@@ -517,24 +435,35 @@ void generateYLabel(PLINT axis, PLFLT value, char *label, PLINT length,
 		SeriesInfo* pInfo = static_cast<SeriesInfo*> (data);
 		// if value equals min, display _name string
 		// otherwise, format value and display it
-		if (value == pInfo->_min) {
+		if (pInfo->_isFirstLabel) {
 			// do not display min value label
 			// but display title :
 			snprintf(label, length, "%*s", pInfo->_maxLen,
 					pInfo->_name.c_str());
-
+			pInfo->_isFirstLabel = false;
 		} else {
 			// display formatted value
 			std::ostringstream labels;
-			// loops over each series and display value at this time
-			auto q = pInfo->_data.find((double) value);
-			if (q != pInfo->_data.end()) {
-				// we found it !
-				labels << boost::format(pInfo->_valuesFormat) % q->second;
-			} else {
-				// not found just display Nan indicator
+
+			//get nearest value
+			double crtTime = (double) value;
+
+			int crtIndex = pInfo->_parameterData->indexOf(crtTime);
+
+			if (crtIndex == -1 && (crtTime - (pInfo->_parameterData->getTimes()[pInfo->_parameterData->getSize()-1]) < pInfo->_parameterData->getParamGapSize()))
+				//Fix last value
+				crtIndex = pInfo->_parameterData->getSize()-1;
+			else if ((crtTime - (pInfo->_parameterData->getTimes()[crtIndex])) > pInfo->_parameterData->getParamGapSize())
+				//Be sure that we are not in a data gap
+				crtIndex = -1;
+
+			if (crtIndex == -1) {
 				labels << "NaN" << std::endl;
 			}
+			else {
+				labels << boost::format(pInfo->_valuesFormat) % (*pInfo->_parameterData->getValues(pInfo->_indexComponents, crtIndex));
+			}
+
 			snprintf(label, length, "%s", labels.str().c_str());
 		}
 	}
@@ -557,11 +486,12 @@ void generateOrbitTimeLabel(PLINT axis, PLFLT value, char *label, PLINT length,
 		// if value equals min, display an empty string
 		// otherwise, format time and display formatted time
  
-		if (value == pInfo->_min) {
+		if (pInfo->_isFirstLabel) {
 			// do not display min value label
 			// but display title :
 			snprintf(label, length, "%*s", pInfo->_maxLen,
 					pInfo->_timeTitle.c_str());
+			pInfo->_isFirstLabel = false;
  
 		} else {
 			// display formatted time
@@ -23,8 +23,6 @@ using std::map;
  
 namespace plot {
  
-typedef std::map<double,double> SeriesData;
-
 struct SeriesInfo{
 	/**
 	 * @brief name of that series (i.e. name of parameter and index if not a scalar parameter).
@@ -42,15 +40,21 @@ struct SeriesInfo{
 	 * @brief maximum len for value. May be useful to justify formatted output.
 	 */
 	unsigned int _maxValueLen;
+
 	/**
-	 * @brief min time index. (used display _name instead of a value).
+	 * @brief to know if it's the first label to draw
 	 */
-	double _min;
+	bool _isFirstLabel;
+
 	/**
-	 * @brief values for that series. Holds for each "time" at a tick ,
-	 *  the related series value.
+	 * @brief link to the parameter data
 	 */
-	SeriesData _data;
+	ParameterData* _parameterData;
+
+	/**
+	 * @brief serie components
+	 */
+	AMDA::Common::ParameterIndexComponent _indexComponents;
 };
  
 typedef std::vector<boost::shared_ptr<SeriesInfo>> SeriesList;
@@ -71,9 +75,9 @@ struct TimeInfo{
 	 */
 	unsigned int _maxLen;
 	/**
-	 * @brief minimum time index. (used display _timeTitle instead of a value).
+	 * @brief to know if it's the first label to draw
 	 */
-	double _min;
+	bool _isFirstLabel;
 };
 typedef boost::shared_ptr<TimeInfo> TimeInfoPtr;
  
@@ -123,11 +127,6 @@ protected:
 private:
  
 	/**
-	 * @brief Calculates tick values (closest value from the tick time)
-	 */
-	unsigned int fillSeriesData(SeriesData& seriesdata_, TimeAxis* pAxis_, ParameterData& paramdata_, AMDA::Common::ParameterIndexComponent indexSerie_);
-
-	/**
 	 * @brief Calculates plot area for date display
 	 */
 	double computeStartDateWidth(PanelPlotOutput* pplot_, TimeAxis* pAxis_);