| 1 | /*************************************************************************** | 
|---|
| 2 | * blitz/vector.h      Declaration of the Vector<P_numtype> class | 
|---|
| 3 | * | 
|---|
| 4 | * $Id: vector.h,v 1.1.1.1 1999-04-09 17:59:00 ansari Exp $ | 
|---|
| 5 | * | 
|---|
| 6 | * Copyright (C) 1997,1998 Todd Veldhuizen <tveldhui@seurat.uwaterloo.ca> | 
|---|
| 7 | * | 
|---|
| 8 | * This program is free software; you can redistribute it and/or | 
|---|
| 9 | * modify it under the terms of the GNU General Public License | 
|---|
| 10 | * as published by the Free Software Foundation; either version 2 | 
|---|
| 11 | * of the License, or (at your option) any later version. | 
|---|
| 12 | * | 
|---|
| 13 | * This program is distributed in the hope that it will be useful, | 
|---|
| 14 | * but WITHOUT ANY WARRANTY; without even the implied warranty of | 
|---|
| 15 | * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the | 
|---|
| 16 | * GNU General Public License for more details. | 
|---|
| 17 | * | 
|---|
| 18 | * Suggestions:          blitz-suggest@cybervision.com | 
|---|
| 19 | * Bugs:                 blitz-bugs@cybervision.com | 
|---|
| 20 | * | 
|---|
| 21 | * For more information, please see the Blitz++ Home Page: | 
|---|
| 22 | *    http://seurat.uwaterloo.ca/blitz/ | 
|---|
| 23 | * | 
|---|
| 24 | *************************************************************************** | 
|---|
| 25 | * $Log: not supported by cvs2svn $ | 
|---|
| 26 | * Revision 1.8  1998/03/14 00:04:47  tveldhui | 
|---|
| 27 | * 0.2-alpha-05 | 
|---|
| 28 | * | 
|---|
| 29 | * Revision 1.7  1997/07/16 14:51:20  tveldhui | 
|---|
| 30 | * Update: Alpha release 0.2 (Arrays) | 
|---|
| 31 | * | 
|---|
| 32 | * Revision 1.6  1997/01/24 14:42:00  tveldhui | 
|---|
| 33 | * Periodic RCS update | 
|---|
| 34 | * | 
|---|
| 35 | * Revision 1.5  1997/01/23 03:28:28  tveldhui | 
|---|
| 36 | * Periodic RCS update | 
|---|
| 37 | * | 
|---|
| 38 | * Revision 1.4  1997/01/13 22:19:58  tveldhui | 
|---|
| 39 | * Periodic RCS update | 
|---|
| 40 | * | 
|---|
| 41 | * Revision 1.3  1996/11/11 17:29:13  tveldhui | 
|---|
| 42 | * Periodic RCS update | 
|---|
| 43 | * | 
|---|
| 44 | * Revision 1.2  1996/10/31 21:06:54  tveldhui | 
|---|
| 45 | * Did away with multiple template parameters. | 
|---|
| 46 | */ | 
|---|
| 47 |  | 
|---|
| 48 | /* | 
|---|
| 49 | * KNOWN BUGS | 
|---|
| 50 | * | 
|---|
| 51 | * 1. operator[](Vector<int>) won't match; compiler complains of no | 
|---|
| 52 | *       suitable copy constructor for VectorPick<T> | 
|---|
| 53 | * 2. Vector<T>(_bz_VecExpr<E>) constructor generates warning | 
|---|
| 54 | * 3. operator+=,-=,..etc.(Random<D>) won't match; compiler complains of | 
|---|
| 55 | *       no suitable copy constructor for _bz_VecExpr(...). | 
|---|
| 56 | */ | 
|---|
| 57 |  | 
|---|
| 58 | #ifndef BZ_VECTOR_H | 
|---|
| 59 | #define BZ_VECTOR_H | 
|---|
| 60 |  | 
|---|
| 61 | #ifndef BZ_BLITZ_H | 
|---|
| 62 | #include <blitz/blitz.h> | 
|---|
| 63 | #endif | 
|---|
| 64 |  | 
|---|
| 65 | #ifndef BZ_MEMBLOCK_H | 
|---|
| 66 | #include <blitz/memblock.h> | 
|---|
| 67 | #endif | 
|---|
| 68 |  | 
|---|
| 69 | #ifndef BZ_RANGE_H | 
|---|
| 70 | #include <blitz/range.h> | 
|---|
| 71 | #endif | 
|---|
| 72 |  | 
|---|
| 73 | #ifndef BZ_LISTINIT_H | 
|---|
| 74 | #include <blitz/listinit.h> | 
|---|
| 75 | #endif | 
|---|
| 76 |  | 
|---|
| 77 | BZ_NAMESPACE(blitz) | 
|---|
| 78 |  | 
|---|
| 79 | // Forward declarations | 
|---|
| 80 | template<class P_numtype> class VectorIter; | 
|---|
| 81 | template<class P_numtype> class VectorIterConst; | 
|---|
| 82 | template<class P_expr>    class _bz_VecExpr; | 
|---|
| 83 | template<class P_numtype> class VectorPick; | 
|---|
| 84 | template<class P_numtype> class Random; | 
|---|
| 85 |  | 
|---|
| 86 | // Declaration of class Vector<P_numtype> | 
|---|
| 87 |  | 
|---|
| 88 | template<class P_numtype> | 
|---|
| 89 | class Vector : protected MemoryBlockReference<P_numtype> { | 
|---|
| 90 |  | 
|---|
| 91 | public: | 
|---|
| 92 | ////////////////////////////////////////////// | 
|---|
| 93 | // Public Types | 
|---|
| 94 | ////////////////////////////////////////////// | 
|---|
| 95 |  | 
|---|
| 96 | typedef P_numtype                  T_numtype; | 
|---|
| 97 | typedef Vector<T_numtype>          T_vector; | 
|---|
| 98 | typedef VectorIter<T_numtype>      T_iterator; | 
|---|
| 99 | typedef VectorIterConst<T_numtype> T_constIterator; | 
|---|
| 100 | typedef VectorPick<T_numtype>      T_pick; | 
|---|
| 101 | typedef Vector<int>                T_indexVector; | 
|---|
| 102 |  | 
|---|
| 103 | ////////////////////////////////////////////// | 
|---|
| 104 | // Constructors                             // | 
|---|
| 105 | ////////////////////////////////////////////// | 
|---|
| 106 | // Most of the constructors are inlined so that | 
|---|
| 107 | // the setting of the stride_ data member will | 
|---|
| 108 | // be visible to the optimizer. | 
|---|
| 109 |  | 
|---|
| 110 | Vector() | 
|---|
| 111 | { | 
|---|
| 112 | length_ = 0; | 
|---|
| 113 | stride_ = 0; | 
|---|
| 114 | } | 
|---|
| 115 |  | 
|---|
| 116 | // This constructor is provided inline because it involves | 
|---|
| 117 | // no memory allocation. | 
|---|
| 118 | Vector(Vector<T_numtype>& vec) | 
|---|
| 119 | : MemoryBlockReference<T_numtype>(vec) | 
|---|
| 120 | { | 
|---|
| 121 | length_ = vec.length_; | 
|---|
| 122 | stride_ = vec.stride_; | 
|---|
| 123 | } | 
|---|
| 124 |  | 
|---|
| 125 | // Slightly unsafe cast-away-const version | 
|---|
| 126 | Vector(const Vector<T_numtype>& vec) | 
|---|
| 127 | : MemoryBlockReference<T_numtype> | 
|---|
| 128 | (const_cast<Vector<T_numtype>& >(vec)) | 
|---|
| 129 | { | 
|---|
| 130 | length_ = vec.length_; | 
|---|
| 131 | stride_ = vec.stride_; | 
|---|
| 132 | } | 
|---|
| 133 |  | 
|---|
| 134 | _bz_explicit Vector(int length) | 
|---|
| 135 | : MemoryBlockReference<T_numtype>(length) | 
|---|
| 136 | { | 
|---|
| 137 | length_ = length; | 
|---|
| 138 | stride_ = 1; | 
|---|
| 139 | } | 
|---|
| 140 |  | 
|---|
| 141 | Vector(Vector<T_numtype>& vec, Range r) | 
|---|
| 142 | : MemoryBlockReference<T_numtype>(vec, | 
|---|
| 143 | r.first() * vec.stride()) | 
|---|
| 144 | { | 
|---|
| 145 | BZPRECONDITION((r.first() >= 0) && (r.first() < vec.length())); | 
|---|
| 146 | BZPRECONDITION((r.last(vec.length()-1) >= 0) | 
|---|
| 147 | && (r.last(vec.length()-1) < vec.length())); | 
|---|
| 148 | length_ = (r.last(vec.length()-1) - r.first()) / r.stride() + 1; | 
|---|
| 149 | stride_ = vec.stride_ * r.stride(); | 
|---|
| 150 | } | 
|---|
| 151 |  | 
|---|
| 152 | Vector(int length, T_numtype initValue) | 
|---|
| 153 | : MemoryBlockReference<T_numtype>(length) | 
|---|
| 154 | { | 
|---|
| 155 | length_ = length; | 
|---|
| 156 | stride_ = 1; | 
|---|
| 157 | (*this) = initValue; | 
|---|
| 158 | } | 
|---|
| 159 |  | 
|---|
| 160 | Vector(int length, T_numtype firstValue, T_numtype delta) | 
|---|
| 161 | : MemoryBlockReference<T_numtype>(length) | 
|---|
| 162 | { | 
|---|
| 163 | length_ = length; | 
|---|
| 164 | stride_ = 1; | 
|---|
| 165 | for (int i=0; i < length; ++i) | 
|---|
| 166 | data_[i] = firstValue + i * delta; | 
|---|
| 167 | } | 
|---|
| 168 |  | 
|---|
| 169 | template<class P_distribution> | 
|---|
| 170 | Vector(int length, Random<P_distribution>& random) | 
|---|
| 171 | : MemoryBlockReference<T_numtype>(length) | 
|---|
| 172 | { | 
|---|
| 173 | length_ = length; | 
|---|
| 174 | stride_ = 1; | 
|---|
| 175 | (*this) = random; | 
|---|
| 176 | } | 
|---|
| 177 |  | 
|---|
| 178 | template<class P_expr> | 
|---|
| 179 | Vector(_bz_VecExpr<P_expr> expr) | 
|---|
| 180 | : MemoryBlockReference<T_numtype>(expr._bz_suggestLength()) | 
|---|
| 181 | { | 
|---|
| 182 | length_ = expr._bz_suggestLength(); | 
|---|
| 183 | stride_ = 1; | 
|---|
| 184 | (*this) = expr; | 
|---|
| 185 | } | 
|---|
| 186 |  | 
|---|
| 187 | // Create a vector view of an already allocated block of memory. | 
|---|
| 188 | // Note that the memory will not be freed when this vector is | 
|---|
| 189 | // destroyed. | 
|---|
| 190 | Vector(int length, T_numtype* _bz_restrict data, int stride = 1) | 
|---|
| 191 | : MemoryBlockReference<T_numtype>(length, data) | 
|---|
| 192 | { | 
|---|
| 193 | length_ = length; | 
|---|
| 194 | stride_ = stride; | 
|---|
| 195 | } | 
|---|
| 196 |  | 
|---|
| 197 | // Create a vector containing a range of numbers | 
|---|
| 198 | Vector(Range r) | 
|---|
| 199 | : MemoryBlockReference<T_numtype>(r._bz_suggestLength()) | 
|---|
| 200 | { | 
|---|
| 201 | length_ = r._bz_suggestLength(); | 
|---|
| 202 | stride_ = 1; | 
|---|
| 203 | (*this) = _bz_VecExpr<Range>(r); | 
|---|
| 204 | } | 
|---|
| 205 |  | 
|---|
| 206 | ////////////////////////////////////////////// | 
|---|
| 207 | // Member functions | 
|---|
| 208 | ////////////////////////////////////////////// | 
|---|
| 209 |  | 
|---|
| 210 | // assertUnitStride() is provided as an optimizing trick.  When | 
|---|
| 211 | // vectors are constructed outside the function scope, the optimizer | 
|---|
| 212 | // is unaware that they have unit stride.  This function sets the | 
|---|
| 213 | // stride to 1 in the local scope so the optimizer can do copy | 
|---|
| 214 | // propagation & dead code elimination.  Obviously, you don't | 
|---|
| 215 | // want to use this routine unless you are certain that the | 
|---|
| 216 | // vectors have unit stride. | 
|---|
| 217 | void            assertUnitStride() | 
|---|
| 218 | { | 
|---|
| 219 | BZPRECONDITION(stride_ == 1); | 
|---|
| 220 | stride_ = 1; | 
|---|
| 221 | } | 
|---|
| 222 |  | 
|---|
| 223 | T_iterator      begin() | 
|---|
| 224 | { return T_iterator(*this); } | 
|---|
| 225 |  | 
|---|
| 226 | T_constIterator begin()  const | 
|---|
| 227 | { return T_constIterator(*this); } | 
|---|
| 228 |  | 
|---|
| 229 | T_vector        copy()   const; | 
|---|
| 230 |  | 
|---|
| 231 | // T_iterator      end(); | 
|---|
| 232 | // T_constIterator end()    const; | 
|---|
| 233 |  | 
|---|
| 234 | T_numtype * _bz_restrict data() | 
|---|
| 235 | { return data_; } | 
|---|
| 236 |  | 
|---|
| 237 | const T_numtype * _bz_restrict data() const | 
|---|
| 238 | { return data_; } | 
|---|
| 239 |  | 
|---|
| 240 | _bz_bool        isUnitStride() const | 
|---|
| 241 | { return stride_ == 1; } | 
|---|
| 242 |  | 
|---|
| 243 | int        length() const | 
|---|
| 244 | { return length_; } | 
|---|
| 245 |  | 
|---|
| 246 | void            makeUnique(); | 
|---|
| 247 |  | 
|---|
| 248 | // int        storageSize() const; | 
|---|
| 249 |  | 
|---|
| 250 | // void            storeToBuffer(void* buffer, int bufferLength) const; | 
|---|
| 251 |  | 
|---|
| 252 | void            reference(T_vector&); | 
|---|
| 253 |  | 
|---|
| 254 | void            resize(int length); | 
|---|
| 255 |  | 
|---|
| 256 | void            resizeAndPreserve(int newLength); | 
|---|
| 257 |  | 
|---|
| 258 | // int             restoreFromBuffer(void* buffer, int bufferLength); | 
|---|
| 259 |  | 
|---|
| 260 | T_vector        reverse() | 
|---|
| 261 | { return T_vector(*this,Range(length()-1,0,-1)); } | 
|---|
| 262 |  | 
|---|
| 263 | int             stride() const | 
|---|
| 264 | { return stride_; } | 
|---|
| 265 |  | 
|---|
| 266 | operator _bz_VecExpr<VectorIterConst<T_numtype> > () const | 
|---|
| 267 | { return _bz_VecExpr<VectorIterConst<T_numtype> >(begin()); } | 
|---|
| 268 |  | 
|---|
| 269 | ///////////////////////////////////////////// | 
|---|
| 270 | // Library-internal member functions | 
|---|
| 271 | // These are undocumented and may change or | 
|---|
| 272 | // disappear in future releases. | 
|---|
| 273 | ///////////////////////////////////////////// | 
|---|
| 274 |  | 
|---|
| 275 | int        _bz_suggestLength() const | 
|---|
| 276 | { return length_; } | 
|---|
| 277 |  | 
|---|
| 278 | _bz_bool        _bz_hasFastAccess() const | 
|---|
| 279 | { return stride_ == 1; } | 
|---|
| 280 |  | 
|---|
| 281 | T_numtype&      _bz_fastAccess(int i) | 
|---|
| 282 | { return data_[i]; } | 
|---|
| 283 |  | 
|---|
| 284 | T_numtype       _bz_fastAccess(int i) const | 
|---|
| 285 | { return data_[i]; } | 
|---|
| 286 |  | 
|---|
| 287 | template<class P_expr, class P_updater> | 
|---|
| 288 | void            _bz_assign(P_expr, P_updater); | 
|---|
| 289 |  | 
|---|
| 290 | _bz_VecExpr<T_constIterator> _bz_asVecExpr() const | 
|---|
| 291 | { return _bz_VecExpr<T_constIterator>(begin()); } | 
|---|
| 292 |  | 
|---|
| 293 | ////////////////////////////////////////////// | 
|---|
| 294 | // Subscripting operators | 
|---|
| 295 | ////////////////////////////////////////////// | 
|---|
| 296 |  | 
|---|
| 297 | // operator()(int) may be used only when the vector has unit | 
|---|
| 298 | // stride.  Otherwise, use operator[]. | 
|---|
| 299 | T_numtype        operator()(int i) const | 
|---|
| 300 | { | 
|---|
| 301 | BZPRECONDITION(i < length_); | 
|---|
| 302 | BZPRECONDITION(stride_ == 1); | 
|---|
| 303 | return data_[i]; | 
|---|
| 304 | } | 
|---|
| 305 |  | 
|---|
| 306 | // operator()(int) may be used only when the vector has unit | 
|---|
| 307 | // stride.  Otherwise, use operator[]. | 
|---|
| 308 | T_numtype& _bz_restrict operator()(int i) | 
|---|
| 309 | { | 
|---|
| 310 | BZPRECONDITION(i < length_); | 
|---|
| 311 | BZPRECONDITION(stride_ == 1); | 
|---|
| 312 | return data_[i]; | 
|---|
| 313 | } | 
|---|
| 314 |  | 
|---|
| 315 | T_numtype        operator[](int i) const | 
|---|
| 316 | { | 
|---|
| 317 | BZPRECONDITION(i < length_); | 
|---|
| 318 | return data_[i * stride_]; | 
|---|
| 319 | } | 
|---|
| 320 |  | 
|---|
| 321 | T_numtype& _bz_restrict operator[](int i) | 
|---|
| 322 | { | 
|---|
| 323 | BZPRECONDITION(i < length_); | 
|---|
| 324 | return data_[i * stride_]; | 
|---|
| 325 | } | 
|---|
| 326 |  | 
|---|
| 327 | T_vector      operator()(Range r) | 
|---|
| 328 | { | 
|---|
| 329 | return T_vector(*this, r); | 
|---|
| 330 | } | 
|---|
| 331 |  | 
|---|
| 332 | T_vector      operator[](Range r) | 
|---|
| 333 | { | 
|---|
| 334 | return T_vector(*this, r); | 
|---|
| 335 | } | 
|---|
| 336 |  | 
|---|
| 337 | T_pick        operator()(T_indexVector i) | 
|---|
| 338 | { | 
|---|
| 339 | return T_pick(*this, i); | 
|---|
| 340 | } | 
|---|
| 341 |  | 
|---|
| 342 | T_pick        operator[](T_indexVector i) | 
|---|
| 343 | { | 
|---|
| 344 | return T_pick(*this, i); | 
|---|
| 345 | } | 
|---|
| 346 |  | 
|---|
| 347 | // T_vector      operator()(difference-equation-expression) | 
|---|
| 348 |  | 
|---|
| 349 | ////////////////////////////////////////////// | 
|---|
| 350 | // Assignment operators | 
|---|
| 351 | ////////////////////////////////////////////// | 
|---|
| 352 |  | 
|---|
| 353 | // Scalar operand | 
|---|
| 354 | ListInitializationSwitch<T_vector,T_iterator> operator=(T_numtype x) | 
|---|
| 355 | { | 
|---|
| 356 | return ListInitializationSwitch<T_vector,T_iterator>(*this, x); | 
|---|
| 357 | } | 
|---|
| 358 |  | 
|---|
| 359 | T_iterator getInitializationIterator() | 
|---|
| 360 | { return begin(); } | 
|---|
| 361 |  | 
|---|
| 362 | T_vector& initialize(T_numtype); | 
|---|
| 363 | T_vector& operator+=(T_numtype); | 
|---|
| 364 | T_vector& operator-=(T_numtype); | 
|---|
| 365 | T_vector& operator*=(T_numtype); | 
|---|
| 366 | T_vector& operator/=(T_numtype); | 
|---|
| 367 | T_vector& operator%=(T_numtype); | 
|---|
| 368 | T_vector& operator^=(T_numtype); | 
|---|
| 369 | T_vector& operator&=(T_numtype); | 
|---|
| 370 | T_vector& operator|=(T_numtype); | 
|---|
| 371 | T_vector& operator>>=(int); | 
|---|
| 372 | T_vector& operator<<=(int); | 
|---|
| 373 |  | 
|---|
| 374 | // Vector operand | 
|---|
| 375 |  | 
|---|
| 376 | template<class P_numtype2> T_vector& operator=(const Vector<P_numtype2> &); | 
|---|
| 377 |  | 
|---|
| 378 | // Specialization uses memcpy instead of element-by-element cast and | 
|---|
| 379 | // copy | 
|---|
| 380 | // NEEDS_WORK -- KCC won't accept this syntax; standard?? | 
|---|
| 381 | // template<> T_vector& operator=(const T_vector&); | 
|---|
| 382 |  | 
|---|
| 383 | template<class P_numtype2> T_vector& operator+=(const Vector<P_numtype2> &); | 
|---|
| 384 | template<class P_numtype2> T_vector& operator-=(const Vector<P_numtype2> &); | 
|---|
| 385 | template<class P_numtype2> T_vector& operator*=(const Vector<P_numtype2> &); | 
|---|
| 386 | template<class P_numtype2> T_vector& operator/=(const Vector<P_numtype2> &); | 
|---|
| 387 | template<class P_numtype2> T_vector& operator%=(const Vector<P_numtype2> &); | 
|---|
| 388 | template<class P_numtype2> T_vector& operator^=(const Vector<P_numtype2> &); | 
|---|
| 389 | template<class P_numtype2> T_vector& operator&=(const Vector<P_numtype2> &); | 
|---|
| 390 | template<class P_numtype2> T_vector& operator|=(const Vector<P_numtype2> &); | 
|---|
| 391 | template<class P_numtype2> T_vector& operator>>=(const Vector<P_numtype2> &); | 
|---|
| 392 | template<class P_numtype2> T_vector& operator<<=(const Vector<P_numtype2> &); | 
|---|
| 393 |  | 
|---|
| 394 | // Vector expression operand | 
|---|
| 395 | template<class P_expr> T_vector& operator=(_bz_VecExpr<P_expr>); | 
|---|
| 396 | template<class P_expr> T_vector& operator+=(_bz_VecExpr<P_expr>); | 
|---|
| 397 | template<class P_expr> T_vector& operator-=(_bz_VecExpr<P_expr>); | 
|---|
| 398 | template<class P_expr> T_vector& operator*=(_bz_VecExpr<P_expr>); | 
|---|
| 399 | template<class P_expr> T_vector& operator/=(_bz_VecExpr<P_expr>); | 
|---|
| 400 | template<class P_expr> T_vector& operator%=(_bz_VecExpr<P_expr>); | 
|---|
| 401 | template<class P_expr> T_vector& operator^=(_bz_VecExpr<P_expr>); | 
|---|
| 402 | template<class P_expr> T_vector& operator&=(_bz_VecExpr<P_expr>); | 
|---|
| 403 | template<class P_expr> T_vector& operator|=(_bz_VecExpr<P_expr>); | 
|---|
| 404 | template<class P_expr> T_vector& operator>>=(_bz_VecExpr<P_expr>); | 
|---|
| 405 | template<class P_expr> T_vector& operator<<=(_bz_VecExpr<P_expr>); | 
|---|
| 406 |  | 
|---|
| 407 | // VectorPick operand | 
|---|
| 408 | template<class P_numtype2> | 
|---|
| 409 | T_vector& operator=(const VectorPick<P_numtype2> &); | 
|---|
| 410 | template<class P_numtype2> | 
|---|
| 411 | T_vector& operator+=(const VectorPick<P_numtype2> &); | 
|---|
| 412 | template<class P_numtype2> | 
|---|
| 413 | T_vector& operator-=(const VectorPick<P_numtype2> &); | 
|---|
| 414 | template<class P_numtype2> | 
|---|
| 415 | T_vector& operator*=(const VectorPick<P_numtype2> &); | 
|---|
| 416 | template<class P_numtype2> | 
|---|
| 417 | T_vector& operator/=(const VectorPick<P_numtype2> &); | 
|---|
| 418 | template<class P_numtype2> | 
|---|
| 419 | T_vector& operator%=(const VectorPick<P_numtype2> &); | 
|---|
| 420 | template<class P_numtype2> | 
|---|
| 421 | T_vector& operator^=(const VectorPick<P_numtype2> &); | 
|---|
| 422 | template<class P_numtype2> | 
|---|
| 423 | T_vector& operator&=(const VectorPick<P_numtype2> &); | 
|---|
| 424 | template<class P_numtype2> | 
|---|
| 425 | T_vector& operator|=(const VectorPick<P_numtype2> &); | 
|---|
| 426 | template<class P_numtype2> | 
|---|
| 427 | T_vector& operator>>=(const VectorPick<P_numtype2> &); | 
|---|
| 428 | template<class P_numtype2> | 
|---|
| 429 | T_vector& operator<<=(const VectorPick<P_numtype2> &); | 
|---|
| 430 |  | 
|---|
| 431 | // Range operand | 
|---|
| 432 | T_vector& operator=(Range); | 
|---|
| 433 | T_vector& operator+=(Range); | 
|---|
| 434 | T_vector& operator-=(Range); | 
|---|
| 435 | T_vector& operator*=(Range); | 
|---|
| 436 | T_vector& operator/=(Range); | 
|---|
| 437 | T_vector& operator%=(Range); | 
|---|
| 438 | T_vector& operator^=(Range); | 
|---|
| 439 | T_vector& operator&=(Range); | 
|---|
| 440 | T_vector& operator|=(Range); | 
|---|
| 441 | T_vector& operator>>=(Range); | 
|---|
| 442 | T_vector& operator<<=(Range); | 
|---|
| 443 |  | 
|---|
| 444 | // Random operand | 
|---|
| 445 | template<class P_distribution> | 
|---|
| 446 | T_vector& operator=(Random<P_distribution>& random); | 
|---|
| 447 | template<class P_distribution> | 
|---|
| 448 | T_vector& operator+=(Random<P_distribution>& random); | 
|---|
| 449 | template<class P_distribution> | 
|---|
| 450 | T_vector& operator-=(Random<P_distribution>& random); | 
|---|
| 451 | template<class P_distribution> | 
|---|
| 452 | T_vector& operator*=(Random<P_distribution>& random); | 
|---|
| 453 | template<class P_distribution> | 
|---|
| 454 | T_vector& operator/=(Random<P_distribution>& random); | 
|---|
| 455 | template<class P_distribution> | 
|---|
| 456 | T_vector& operator%=(Random<P_distribution>& random); | 
|---|
| 457 | template<class P_distribution> | 
|---|
| 458 | T_vector& operator^=(Random<P_distribution>& random); | 
|---|
| 459 | template<class P_distribution> | 
|---|
| 460 | T_vector& operator&=(Random<P_distribution>& random); | 
|---|
| 461 | template<class P_distribution> | 
|---|
| 462 | T_vector& operator|=(Random<P_distribution>& random); | 
|---|
| 463 |  | 
|---|
| 464 | ////////////////////////////////////////////// | 
|---|
| 465 | // Unary operators | 
|---|
| 466 | ////////////////////////////////////////////// | 
|---|
| 467 |  | 
|---|
| 468 | //    T_vector& operator++(); | 
|---|
| 469 | //    void operator++(int); | 
|---|
| 470 | //    T_vector& operator--(); | 
|---|
| 471 | //    void operator--(int); | 
|---|
| 472 |  | 
|---|
| 473 | private: | 
|---|
| 474 | int      length_; | 
|---|
| 475 | int      stride_; | 
|---|
| 476 | }; | 
|---|
| 477 |  | 
|---|
| 478 | // Global I/O functions | 
|---|
| 479 |  | 
|---|
| 480 | template<class P_numtype> | 
|---|
| 481 | ostream& operator<<(ostream& os, const Vector<P_numtype>& x); | 
|---|
| 482 |  | 
|---|
| 483 | template<class P_expr> | 
|---|
| 484 | ostream& operator<<(ostream& os, _bz_VecExpr<P_expr> expr); | 
|---|
| 485 |  | 
|---|
| 486 | BZ_NAMESPACE_END | 
|---|
| 487 |  | 
|---|
| 488 | #include <blitz/veciter.h>          // Iterators | 
|---|
| 489 | #include <blitz/vecpick.h>          // VectorPick | 
|---|
| 490 | #include <blitz/vecexpr.h>          // Expression templates | 
|---|
| 491 | #include <blitz/vecglobs.h>         // Global functions | 
|---|
| 492 | #include <blitz/vector.cc>          // Member functions | 
|---|
| 493 | #include <blitz/vecio.cc>           // IO functions | 
|---|
| 494 |  | 
|---|
| 495 | #endif // BZ_VECTOR_H | 
|---|