817 lines
14 KiB
C++
817 lines
14 KiB
C++
/*
|
|
This file is part of SECONDO.
|
|
|
|
Copyright (C) 2013, University in Hagen, Department of Computer Science,
|
|
Database Systems for New Applications.
|
|
|
|
SECONDO is free software; you can redistribute it and/or modify
|
|
it under the terms of the GNU General Public License as published by
|
|
the Free Software Foundation; either version 2 of the License, or
|
|
(at your option) any later version.
|
|
|
|
SECONDO is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
GNU General Public License for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with SECONDO; if not, write to the Free Software
|
|
Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
|
|
|
|
*/
|
|
|
|
#ifndef TILEALGEBRA_MTCELLITERATOR_H
|
|
#define TILEALGEBRA_MTCELLITERATOR_H
|
|
|
|
/*
|
|
TileAlgebra includes
|
|
|
|
*/
|
|
|
|
#include "../grid/mtgrid.h"
|
|
#include "../Index/Index.h"
|
|
|
|
/*
|
|
declaration of namespace TileAlgebra
|
|
|
|
*/
|
|
|
|
namespace TileAlgebra
|
|
{
|
|
|
|
/*
|
|
Template class mtCellIterator represents an iteration class
|
|
for mt datatype cells.
|
|
|
|
author: Dirk Zacher
|
|
|
|
*/
|
|
|
|
template <typename Type>
|
|
class mtCellIterator
|
|
{
|
|
public:
|
|
|
|
/*
|
|
Constructor mtCellIterator initializes some members of the class
|
|
with corresponding parameter values and all the other members
|
|
of the class with default values.
|
|
|
|
author: Dirk Zacher
|
|
parameters: rmt - reference to a mt datatype object
|
|
rX1 - reference to X1 coordinate
|
|
rY1 - reference to Y1 coordinate
|
|
rX2 - reference to X2 coordinate
|
|
rY2 - reference to Y2 coordinate
|
|
rT1 - reference to T1 coordinate
|
|
rT2 - reference to T2 coordinate
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
mtCellIterator(const Type& rmt,
|
|
const double& rX1,
|
|
const double& rY1,
|
|
const double& rX2,
|
|
const double& rY2,
|
|
const double& rT1,
|
|
const double& rT2);
|
|
|
|
/*
|
|
Destructor ~mtCellIterator deinitializes a mtCellIterator object.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
virtual ~mtCellIterator();
|
|
|
|
/*
|
|
Method HasNext checks if mtCellIterator object has a next cell to iterate.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: true, if mtCellIterator object has a next cell to iterate,
|
|
otherwise false
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
bool HasNext();
|
|
|
|
/*
|
|
Method Next returns a pair of last delta value and current delta value
|
|
and calculates next delta values.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: a pair of last delta value and current delta value
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
std::pair<double, double> Next();
|
|
|
|
private:
|
|
|
|
/*
|
|
Method Init initializes mtCellIterator.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
void Init();
|
|
|
|
/*
|
|
Method ComputeNextX calculates next x value.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
void ComputeNextX();
|
|
|
|
/*
|
|
Method ComputeNextY calculates next y value.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
void ComputeNextY();
|
|
|
|
/*
|
|
Method ComputeNextT calculates next time value.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
void ComputeNextT();
|
|
|
|
/*
|
|
Method GetNextValue returns next value of given index.
|
|
|
|
author: Dirk Zacher
|
|
parameters: index - index of next value to return.
|
|
return value: next value of given index
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
double GetNextValue(int index);
|
|
|
|
/*
|
|
Method ComputeNextDelta calculates next delta values.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
void ComputeNextDelta();
|
|
|
|
/*
|
|
Method ShiftX decrements or increments x dimension of current cell index
|
|
dependent on current delta x value.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
void ShiftX();
|
|
|
|
/*
|
|
Method ShiftY decrements or increments y dimension of current cell index
|
|
dependent on current delta y value.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
void ShiftY();
|
|
|
|
/*
|
|
Method ShiftT decrements or increments time dimension of current cell index
|
|
dependent on current delta time value.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
void ShiftT();
|
|
|
|
/*
|
|
Member m_Grid contains the mtgrid object.
|
|
|
|
*/
|
|
|
|
mtgrid m_Grid;
|
|
|
|
/*
|
|
Member m_dX1 contains the X1 coordinate.
|
|
|
|
*/
|
|
|
|
double m_dX1;
|
|
|
|
/*
|
|
Member m_dY1 contains the Y1 coordinate.
|
|
|
|
*/
|
|
|
|
double m_dY1;
|
|
|
|
/*
|
|
Member m_dX2 contains the X2 coordinate.
|
|
|
|
*/
|
|
|
|
double m_dX2;
|
|
|
|
/*
|
|
Member m_dY2 contains the Y2 coordinate.
|
|
|
|
*/
|
|
|
|
double m_dY2;
|
|
|
|
/*
|
|
Member m_dT1 contains the T1 coordinate.
|
|
|
|
*/
|
|
|
|
double m_dT1;
|
|
|
|
/*
|
|
Member m_dT2 contains the T2 coordinate.
|
|
|
|
*/
|
|
|
|
double m_dT2;
|
|
|
|
/*
|
|
Member m_bHasNextX is true if mtCellIterator has a next cell in x dimension,
|
|
otherwise m_bHasNextX is false.
|
|
|
|
*/
|
|
|
|
bool m_bHasNextX;
|
|
|
|
/*
|
|
Member m_bHasNextY is true if mtCellIterator has a next cell in y dimension,
|
|
otherwise m_bHasNextY is false.
|
|
|
|
*/
|
|
|
|
bool m_bHasNextY;
|
|
|
|
/*
|
|
Member m_bHasNextT is true if mtCellIterator has a next cell
|
|
in time dimension, otherwise m_bHasNextT is false.
|
|
|
|
*/
|
|
|
|
bool m_bHasNextT;
|
|
|
|
/*
|
|
Member m_dNextX contains the next value in x dimension.
|
|
|
|
*/
|
|
|
|
double m_dNextX;
|
|
|
|
/*
|
|
Member m_dNextY contains the next value in y dimension.
|
|
|
|
*/
|
|
|
|
double m_dNextY;
|
|
|
|
/*
|
|
Member m_dNextT contains the next value in time dimension.
|
|
|
|
*/
|
|
|
|
double m_dNextT;
|
|
|
|
/*
|
|
Member m_dLastDelta contains the last delta value.
|
|
|
|
*/
|
|
|
|
double m_dLastDelta;
|
|
|
|
/*
|
|
Member m_dCurrentDelta contains the current delta value.
|
|
|
|
*/
|
|
|
|
double m_dCurrentDelta;
|
|
|
|
/*
|
|
Member m_LastCellIndex contains the last 3-dimensional cell index.
|
|
|
|
*/
|
|
|
|
Index<3> m_LastCellIndex;
|
|
|
|
/*
|
|
Member m_CurrentCellIndex contains the current 3-dimensional cell index.
|
|
|
|
*/
|
|
|
|
Index<3> m_CurrentCellIndex;
|
|
|
|
/*
|
|
Member m_dDeltaX contains the delta x value.
|
|
|
|
*/
|
|
|
|
double m_dDeltaX;
|
|
|
|
/*
|
|
Member m_dDeltaY contains the delta y value.
|
|
|
|
*/
|
|
|
|
double m_dDeltaY;
|
|
|
|
/*
|
|
Member m_dDeltaT contains the delta time value.
|
|
|
|
*/
|
|
|
|
double m_dDeltaT;
|
|
};
|
|
|
|
/*
|
|
Constructor mtCellIterator inintializes some members of the class
|
|
with corresponding parameter values and all the other members
|
|
of the class with default values.
|
|
|
|
author: Dirk Zacher
|
|
parameters: rmt - reference to a mt object
|
|
rX1 - reference to X1 coordinate
|
|
rY1 - reference to Y1 coordinate
|
|
rX2 - reference to X2 coordinate
|
|
rY2 - reference to Y2 coordinate
|
|
rT1 - reference to T1 coordinate
|
|
rT2 - reference to T2 coordinate
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
template <typename Type>
|
|
mtCellIterator<Type>::mtCellIterator(const Type& rmt,
|
|
const double& rX1,
|
|
const double& rY1,
|
|
const double& rX2,
|
|
const double& rY2,
|
|
const double& rT1,
|
|
const double& rT2)
|
|
:m_dX1(rX1),
|
|
m_dY1(rY1),
|
|
m_dX2(rX2),
|
|
m_dY2(rY2),
|
|
m_dT1(rT1),
|
|
m_dT2(rT2),
|
|
m_bHasNextX(false),
|
|
m_bHasNextY(false),
|
|
m_bHasNextT(false),
|
|
m_dNextX(0.0),
|
|
m_dNextY(0.0),
|
|
m_dNextT(0.0),
|
|
m_dLastDelta(0.0),
|
|
m_dCurrentDelta(0.0)
|
|
{
|
|
rmt.getgrid(m_Grid);
|
|
|
|
if(rmt.IsValidLocation(m_dX2, m_dY2, m_dT2))
|
|
{
|
|
m_LastCellIndex = rmt.GetLocationIndex(m_dX2, m_dY2, m_dT2);
|
|
}
|
|
|
|
if(rmt.IsValidLocation(m_dX1, m_dY1, m_dT1))
|
|
{
|
|
m_CurrentCellIndex = rmt.GetLocationIndex(m_dX1, m_dY1, m_dT1);
|
|
}
|
|
|
|
Init();
|
|
}
|
|
|
|
/*
|
|
Destructor deinitializes a mtCellIterator object.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
template <typename Type>
|
|
mtCellIterator<Type>::~mtCellIterator()
|
|
{
|
|
|
|
}
|
|
|
|
/*
|
|
Method HasNext checks if mtCellIterator object has a next cell to iterate.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: true, if mtCellIterator object has a next cell to iterate,
|
|
otherwise false
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
template <typename Type>
|
|
bool mtCellIterator<Type>::HasNext()
|
|
{
|
|
return m_dLastDelta < 1.0;
|
|
}
|
|
|
|
/*
|
|
Method Next returns a pair of last delta value and current delta value
|
|
and calculates next delta values.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: a pair of last delta value and current delta value
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
template <typename Type>
|
|
std::pair<double, double> mtCellIterator<Type>::Next()
|
|
{
|
|
std::pair<double, double> retVal(m_dLastDelta, m_dCurrentDelta);
|
|
|
|
ComputeNextDelta();
|
|
|
|
return retVal;
|
|
}
|
|
|
|
/*
|
|
Method Init initializes mtCellIterator.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
template <typename Type>
|
|
void mtCellIterator<Type>::Init()
|
|
{
|
|
if(m_CurrentCellIndex == m_LastCellIndex)
|
|
{
|
|
m_dCurrentDelta = 1.0;
|
|
}
|
|
|
|
else
|
|
{
|
|
m_dDeltaX = m_dX2 - m_dX1;
|
|
m_dDeltaY = m_dY2 - m_dY1;
|
|
m_dDeltaT = m_dT2 - m_dT1;
|
|
|
|
ComputeNextX();
|
|
ComputeNextY();
|
|
ComputeNextT();
|
|
}
|
|
}
|
|
|
|
/*
|
|
Method ComputeNextX calculates next x value.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
template <typename Type>
|
|
void mtCellIterator<Type>::ComputeNextX()
|
|
{
|
|
m_bHasNextX = m_CurrentCellIndex[0] != m_LastCellIndex[0];
|
|
|
|
if(m_bHasNextX)
|
|
{
|
|
double verticalX = 0.0;
|
|
|
|
if(m_dDeltaX < 0.0)
|
|
{
|
|
// left border of cell
|
|
verticalX = m_Grid.GetX() +
|
|
m_CurrentCellIndex[0] * m_Grid.GetLength();
|
|
}
|
|
|
|
else
|
|
{
|
|
// right border of cell
|
|
verticalX = m_Grid.GetX() +
|
|
(m_CurrentCellIndex[0] + 1) * m_Grid.GetLength();
|
|
}
|
|
|
|
m_dNextX = (verticalX - m_dX1)/ m_dDeltaX;
|
|
}
|
|
}
|
|
|
|
/*
|
|
Method ComputeNextY calculates next y value.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
template <typename Type>
|
|
void mtCellIterator<Type>::ComputeNextY()
|
|
{
|
|
m_bHasNextY = m_CurrentCellIndex[1] != m_LastCellIndex[1];
|
|
|
|
if(m_bHasNextY)
|
|
{
|
|
double horizontalY = 0.0;
|
|
|
|
if(m_dDeltaY < 0.0)
|
|
{
|
|
// bottom border of cell
|
|
horizontalY = m_Grid.GetY() +
|
|
m_CurrentCellIndex[1] * m_Grid.GetLength();
|
|
}
|
|
|
|
else
|
|
{
|
|
// right border of cell
|
|
horizontalY = m_Grid.GetY() +
|
|
(m_CurrentCellIndex[1] + 1) * m_Grid.GetLength();
|
|
}
|
|
|
|
m_dNextY = (horizontalY - m_dY1) / m_dDeltaY;
|
|
}
|
|
}
|
|
|
|
/*
|
|
Method ComputeNextT calculates next time value.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
template <typename Type>
|
|
void mtCellIterator<Type>::ComputeNextT()
|
|
{
|
|
m_bHasNextT = m_CurrentCellIndex[2] != m_LastCellIndex[2];
|
|
|
|
if(m_bHasNextT)
|
|
{
|
|
double horizontalT = 0.0;
|
|
|
|
if(m_dDeltaT < 0.0)
|
|
{
|
|
// bottom border of cell
|
|
horizontalT = m_CurrentCellIndex[2] *
|
|
m_Grid.GetDuration().ToDouble();
|
|
}
|
|
|
|
else
|
|
{
|
|
// right border of cell
|
|
horizontalT = (m_CurrentCellIndex[2] + 1) *
|
|
m_Grid.GetDuration().ToDouble();
|
|
}
|
|
|
|
m_dNextT = (horizontalT - m_dT1) / m_dDeltaT;
|
|
}
|
|
}
|
|
|
|
/*
|
|
Method GetNextValue returns next value of given index.
|
|
|
|
author: Dirk Zacher
|
|
parameters: index - index of next value to return.
|
|
return value: next value of given index
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
template <typename Type>
|
|
double mtCellIterator<Type>::GetNextValue(int index)
|
|
{
|
|
double dNextValue = -1.0;
|
|
|
|
switch(index)
|
|
{
|
|
case 0: dNextValue = m_dNextX;
|
|
break;
|
|
case 1: dNextValue = m_dNextY;
|
|
break;
|
|
case 2: dNextValue = m_dNextT;
|
|
break;
|
|
default: assert(false);
|
|
break;
|
|
}
|
|
|
|
return dNextValue;
|
|
}
|
|
|
|
/*
|
|
Method ComputeNextDelta calculates next delta values.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
template <typename Type>
|
|
void mtCellIterator<Type>::ComputeNextDelta()
|
|
{
|
|
m_dLastDelta = m_dCurrentDelta;
|
|
|
|
if(m_CurrentCellIndex == m_LastCellIndex)
|
|
{
|
|
m_dCurrentDelta = 1.0;
|
|
}
|
|
|
|
else
|
|
{
|
|
assert(m_bHasNextX || m_bHasNextY || m_bHasNextT);
|
|
|
|
bool cand[3];
|
|
cand[0] = m_bHasNextX;
|
|
cand[1] = m_bHasNextY;
|
|
cand[2] = m_bHasNextT;
|
|
|
|
double minimum = 0.0;
|
|
bool first = true;
|
|
|
|
for(int i = 0; i < 3; i++)
|
|
{
|
|
if(cand[i])
|
|
{
|
|
if(first)
|
|
{
|
|
first = false;
|
|
minimum = GetNextValue(i);
|
|
}
|
|
|
|
else
|
|
{
|
|
double dv = GetNextValue(i);
|
|
|
|
if(dv < minimum)
|
|
{
|
|
// new minimum found
|
|
minimum = dv;
|
|
|
|
for(int j = 0; j < i; j++)
|
|
{
|
|
cand[j] = false;
|
|
}
|
|
}
|
|
}
|
|
}
|
|
}
|
|
|
|
m_dCurrentDelta = minimum;
|
|
|
|
if(cand[0] == true)
|
|
{
|
|
ShiftX();
|
|
ComputeNextX();
|
|
}
|
|
|
|
if(cand[1] == true)
|
|
{
|
|
ShiftY();
|
|
ComputeNextY();
|
|
}
|
|
|
|
if(cand[2] == true)
|
|
{
|
|
ShiftT();
|
|
ComputeNextT();
|
|
}
|
|
}
|
|
}
|
|
|
|
/*
|
|
Method ShiftX decrements or increments x dimension of current cell index
|
|
dependent on current delta x value.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
template <typename Type>
|
|
void mtCellIterator<Type>::ShiftX()
|
|
{
|
|
if(m_dDeltaX > 0.0)
|
|
{
|
|
m_CurrentCellIndex.Increment(0);
|
|
}
|
|
|
|
else if(m_dDeltaX < 0.0)
|
|
{
|
|
m_CurrentCellIndex.Decrement(0);
|
|
}
|
|
}
|
|
|
|
/*
|
|
Method ShiftY decrements or increments y dimension of current cell index
|
|
dependent on current delta y value.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
template <typename Type>
|
|
void mtCellIterator<Type>::ShiftY()
|
|
{
|
|
if(m_dDeltaY > 0.0)
|
|
{
|
|
m_CurrentCellIndex.Increment(1);
|
|
}
|
|
|
|
else if(m_dDeltaY < 0.0)
|
|
{
|
|
m_CurrentCellIndex.Decrement(1);
|
|
}
|
|
}
|
|
|
|
/*
|
|
Method ShiftT decrements or increments time dimension of current cell index
|
|
dependent on current delta time value.
|
|
|
|
author: Dirk Zacher
|
|
parameters: -
|
|
return value: -
|
|
exceptions: -
|
|
|
|
*/
|
|
|
|
template <typename Type>
|
|
void mtCellIterator<Type>::ShiftT()
|
|
{
|
|
if(m_dDeltaT > 0.0)
|
|
{
|
|
m_CurrentCellIndex.Increment(2);
|
|
}
|
|
|
|
else if(m_dDeltaT < 0.0)
|
|
{
|
|
m_CurrentCellIndex.Decrement(2);
|
|
}
|
|
}
|
|
|
|
}
|
|
|
|
#endif // TILEALGEBRA_MTCELLITERATOR_H
|