diff options
author | Constantin Jucovschi <cj@ubuntu.localdomain> | 2009-04-24 07:20:22 -0400 |
---|---|---|
committer | Constantin Jucovschi <cj@ubuntu.localdomain> | 2009-04-24 07:20:22 -0400 |
commit | 8f27e65bddd7d4b8515ce620fb485fdd78fcdf89 (patch) | |
tree | bd328a4dd4f92d32202241b5e3a7f36177792c5f /raslib/mitera.cc | |
download | rasdaman-upstream-8f27e65bddd7d4b8515ce620fb485fdd78fcdf89.tar.gz rasdaman-upstream-8f27e65bddd7d4b8515ce620fb485fdd78fcdf89.tar.xz rasdaman-upstream-8f27e65bddd7d4b8515ce620fb485fdd78fcdf89.zip |
Initial commitv8.0
Diffstat (limited to 'raslib/mitera.cc')
-rw-r--r-- | raslib/mitera.cc | 141 |
1 files changed, 141 insertions, 0 deletions
diff --git a/raslib/mitera.cc b/raslib/mitera.cc new file mode 100644 index 0000000..b92cb72 --- /dev/null +++ b/raslib/mitera.cc @@ -0,0 +1,141 @@ +/* +* This file is part of rasdaman community. +* +* Rasdaman community 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 3 of the License, or +* (at your option) any later version. +* +* Rasdaman community 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 rasdaman community. If not, see <http://www.gnu.org/licenses/>. +* +* Copyright 2003, 2004, 2005, 2006, 2007, 2008, 2009 Peter Baumann / +rasdaman GmbH. +* +* For more information please see <http://www.rasdaman.org> +* or contact Peter Baumann via <baumann@rasdaman.com>. +/ +/** + * SOURCE: mitera.cc + * + * MODULE: raslib + * CLASS: r_MiterArea + * +*/ + +#include "raslib/mitera.hh" +#include "raslib/minterval.hh" +#include "raslib/rmdebug.hh" + +r_MiterArea::r_MiterArea( const r_Minterval* newIterDom, + const r_Minterval* newImgDom ) throw(r_Error) + : iterDom(newIterDom), imgDom(newImgDom), done(false) +{ + RMDBGENTER( 1, RMDebug::module_raslib, "r_MiterArea", "r_MiterArea()"); + + if (imgDom->dimension() != iterDom->dimension()) + { + //in this case we have an undefined situation + RMInit::logOut << "r_MiterArea::rMiterArea(" << iterDom << ", " << imgDom << ") dimension mismatch" << endl; + throw r_Error(INTERVALSWITHDIFFERENTDIMENSION); + } + + if(!imgDom->is_origin_fixed() || + !imgDom->is_high_fixed()) + { + //in this case we have an undefined situation + RMInit::logOut << "r_MiterArea::rMiterArea(" << iterDom << ", " << imgDom << ") imgDom is opened" << endl; + throw r_Error(INTERVALOPEN); + } + + if(!iterDom->is_origin_fixed() || + !iterDom->is_high_fixed()) + { + //in this case we have an undefined situation + RMInit::logOut << "r_MiterArea::rMiterArea(" << iterDom << ", " << imgDom << ") iterDom is opened" << endl; + throw r_Error(INTERVALOPEN); + } + + + // dimensionality of both iterDom and imgDom + r_Dimension dim = imgDom->dimension(); + // stores the increments + incArrIter = new incArrElem[dim]; + + for(r_Dimension i=0; i<dim; i++ ) { + // used for counting in iteration, initialize with 0 + incArrIter[i].curr = 0; + // how often is the iterDom moved inside the imgDom + incArrIter[i].repeat = (imgDom->get_extent()[i] / iterDom->get_extent()[i]) + + ( imgDom->get_extent()[i] % iterDom->get_extent()[i] != 0 ); + + RMDBGMIDDLE( 4, RMDebug::module_raslib, "r_MiterArea", "repeat dim " << i << ": " << incArrIter[i].repeat ); + } + reset(); + RMDBGEXIT( 1, RMDebug::module_raslib, "r_MiterArea", "r_MiterArea()" ); +} + + +r_MiterArea::~r_MiterArea() +{ + delete [] incArrIter; +} + +void +r_MiterArea::reset() +{ + done = false; + for( int i=0; i<iterDom->dimension(); i++ ) { + incArrIter[i].curr = 0; + } +} + +r_Minterval +r_MiterArea::nextArea() +{ + r_Dimension i = 0; + + if(done) + return retVal; + + r_Minterval currDom(iterDom->dimension()); + + // calculate new result domain here + if(!done) { + for(i=0; i < iterDom->dimension(); i++) { + currDom << r_Sinterval( (*imgDom)[i].low() + incArrIter[i].curr*iterDom->get_extent()[i], + (*imgDom)[i].low() + (incArrIter[i].curr+1)*iterDom->get_extent()[i] + - 1 ); + } + } + retVal = currDom.intersection_with((*imgDom)); + + // increment dimensions + for(i=0; i < iterDom->dimension(); i++) { + incArrIter[i].curr++; + if(incArrIter[i].curr < incArrIter[i].repeat) { + // no overflow in this dimension + break; + } else { + // overflow in this dimension + incArrIter[i].curr = 0; + } + } + if( i == iterDom->dimension() ) { + // overflow in last dimension + done = true; + } + + return retVal; +} + +bool +r_MiterArea::isDone() +{ + return done; +} |