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/point.cc | |
download | rasdaman-upstream-8f27e65bddd7d4b8515ce620fb485fdd78fcdf89.tar.gz rasdaman-upstream-8f27e65bddd7d4b8515ce620fb485fdd78fcdf89.tar.xz rasdaman-upstream-8f27e65bddd7d4b8515ce620fb485fdd78fcdf89.zip |
Initial commitv8.0
Diffstat (limited to 'raslib/point.cc')
-rw-r--r-- | raslib/point.cc | 372 |
1 files changed, 372 insertions, 0 deletions
diff --git a/raslib/point.cc b/raslib/point.cc new file mode 100644 index 0000000..8679dc5 --- /dev/null +++ b/raslib/point.cc @@ -0,0 +1,372 @@ +/* +* 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: point.cc + * + * MODULE: raslib + * CLASS: r_Point + * + * COMMENTS: + * +*/ + +static const char rcsid[] = "@(#)raslib, r_Point: $Id: point.cc,v 1.22 2002/08/28 11:58:13 coman Exp $"; + +#include "point.hh" + +#include <string.h> +#ifdef __VISUALC__ +#include <strstrea.h> +#else +#include <strstream> +#endif + +#include "raslib/rminit.hh" +#include "raslib/error.hh" + +r_Point::r_Point( char* stringRep ) throw (r_Error) + : dimensionality(1), streamInitCnt(0), points(NULL) +{ + char charToken = 0; + r_Range valueToken = 0; + + // for parsing the string + std::istrstream str( stringRep, strlen(stringRep) + 1 ); + + // calculate dimensionality + char* p = stringRep; + while(p = strchr(++p, ',')) + dimensionality++; + + // allocate space for intervals + points = new r_Range[ dimensionality ]; + + str >> charToken; + if(charToken != '[') + { + // error + dimensionality = 0; + delete[] points; + points = NULL; + throw r_Error(NOPOINT); + } + + for( r_Dimension i=0; i<dimensionality; i++ ) + { + str >> valueToken; + points[i] = valueToken; + + if( i < dimensionality-1 ) + { + str >> charToken; + if(charToken != ',') + { + // error + dimensionality = 0; + delete[] points; + points = NULL; + throw r_Error(NOPOINT); + } + } + } +} + + +r_Point::r_Point( r_Dimension dim ) + : dimensionality( dim ), + streamInitCnt(0) +{ + points = new r_Range[ dimensionality ]; + + for( r_Dimension i=0; i< dimensionality; i++ ) + points[i] = 0; +} + + +r_Point& r_Point::operator<<( r_Range newElement ) throw( r_Einit_overflow ) +{ + if( streamInitCnt >= dimensionality ) + { + RMInit::logOut << "r_Point::operator<<(" << newElement << ") already fully initialised" << endl; + throw r_Einit_overflow(); + } + + points[streamInitCnt++] = newElement; + return *this; +} + + +r_Point::r_Point( r_Range p1, r_Range p2 ) + : dimensionality(2), + streamInitCnt(2) +{ + points = new r_Range[dimensionality]; + points[0] = p1; + points[1] = p2; +} + + +r_Point::r_Point( r_Range p1, r_Range p2, r_Range p3 ) + : dimensionality(3), + streamInitCnt(3) +{ + points = new r_Range[dimensionality]; + points[0] = p1; + points[1] = p2; + points[2] = p3; +} + + +r_Point::r_Point( r_Range p1, r_Range p2, r_Range p3, r_Range p4 ) + : dimensionality(4), + streamInitCnt(4) +{ + points = new r_Range[dimensionality]; + points[0] = p1; + points[1] = p2; + points[2] = p3; + points[3] = p4; +} + + +r_Point::r_Point( r_Range p1, r_Range p2, r_Range p3, r_Range p4, r_Range p5 ) + : dimensionality(5), + streamInitCnt(5) +{ + points = new r_Range[dimensionality]; + points[0] = p1; + points[1] = p2; + points[2] = p3; + points[3] = p4; + points[4] = p5; +} + + +r_Point::r_Point() + : dimensionality(0), + streamInitCnt(0), + points(NULL) +{ +} + + +r_Point::r_Point( const r_Point& pt ) + : dimensionality(pt.dimensionality), + streamInitCnt(pt.streamInitCnt), + points(new r_Range[pt.dimensionality]) +{ + for( r_Dimension i=0; i<dimensionality; i++ ) + points[i] = pt[i]; +} + + +r_Point::~r_Point() +{ + if( points ) + { + delete[] points; + points = NULL; + } +} + + +r_Range +r_Point::operator[]( r_Dimension i ) const throw( r_Eindex_violation ) +{ + if( i < 0 || i >= dimensionality ) + { + RMInit::logOut << "r_Point::operator[](" << i << ") const dimension out of bounds (" << dimensionality << ")" << endl; + throw r_Eindex_violation( 0, dimensionality-1, i ); + } + + return points[i]; +} + + +r_Range& +r_Point::operator[]( r_Dimension i ) throw( r_Eindex_violation ) +{ + if( i < 0 || i >= dimensionality ) + { + RMInit::logOut << "r_Point::operator[](" << i << ") dimension out of bounds (" << dimensionality << ")" << endl; + throw r_Eindex_violation( 0, dimensionality-1, i ); + } + + return points[i]; +} + + +const r_Point& +r_Point::operator=( const r_Point& pt ) +{ + if( this != &pt ) + { + if( points && dimensionality != pt.dimension() ) + { + delete[] points; + points = NULL; + } + + dimensionality = pt.dimension(); + streamInitCnt = dimensionality; + + if( !points ) + points = new r_Range[ dimensionality ]; + + for( r_Dimension i=0; i<dimensionality; i++ ) + points[i] = pt[i]; + + } + + return *this; +} + + + +bool +r_Point::operator==( const r_Point& pt ) const +{ + bool returnValue = false; + + if( dimensionality == pt.dimensionality ) + { + returnValue = true; + + for( r_Dimension i=0; i<dimensionality && returnValue ; i++ ) + { + if (points[i] != pt[i]) + { + returnValue = false; + break; + } + } + } + + return returnValue; +} + + + +bool +r_Point::operator!=( const r_Point& pt ) const +{ + return !operator==( pt ); +} + + + +r_Point +r_Point::operator+( const r_Point& pt ) const throw( r_Edim_mismatch ) +{ + if( dimensionality != pt.dimension() ) + { + RMInit::logOut << "r_Point::operator+(" << pt << ") dimensions (" << dimensionality << ") do not match" << endl; + throw r_Edim_mismatch( dimensionality, pt.dimension() ); + } + + r_Point result( dimensionality ); + + for( r_Dimension i=0; i<dimensionality; i++ ) + result[i] = points[i] + pt[i]; + + return result; +} + +r_Point +r_Point::operator-( const r_Point& pt ) const throw( r_Edim_mismatch ) +{ + if( dimensionality != pt.dimension() ) + { + RMInit::logOut << "r_Point::operator-(" << pt << ") dimensions (" << dimensionality << ") do not match" << endl; + throw r_Edim_mismatch( dimensionality, pt.dimension() ); + } + + r_Point result( dimensionality ); + + for( r_Dimension i=0; i<dimensionality; i++ ) + result[i] = points[i] - pt[i]; + + return result; +} + +r_Point +r_Point::operator*( const r_Point& pt ) const throw( r_Edim_mismatch ) +{ + if( dimensionality != pt.dimension() ) + { + RMInit::logOut << "r_Point::operator*(" << pt << ") dimensions (" << dimensionality << ") do not match" << endl; + throw r_Edim_mismatch( dimensionality, pt.dimension() ); + } + + r_Point result( dimensionality ); + + for( r_Dimension i=0; i<dimensionality; i++ ) + result[i] = points[i] * pt[i]; + + return result; +} + +void +r_Point::print_status( std::ostream& s ) const +{ + s << "["; + + if( dimensionality > 0 ) + { + for( r_Dimension i=0; i<dimensionality-1; i++ ) + s << points[i] << ","; + + s << points[dimensionality-1]; + } + + s << "]"; +} + + +char* +r_Point::get_string_representation() const +{ + r_Bytes bufferSize = dimensionality*25+3; // should be enough + + // allocate buffer and initialize string stream + char* buffer = new char[bufferSize]; + std::ostrstream domainStream( buffer, bufferSize ); + + // write into string stream + domainStream << (*this) << std::ends; + + // allocate memory taking the final string + char* returnString = strdup(buffer); + + // delete buffer + delete[] buffer; + + return returnString; +} + +std::ostream& operator<<( std::ostream& s, const r_Point& d ) +{ + d.print_status( s ); + return s; +} + |