| 1 | /*************************************************************************** | 
|---|
| 2 | * blitz/array/map.h      Declaration of the ArrayIndexMapping class | 
|---|
| 3 | * | 
|---|
| 4 | * $Id: map.h,v 1.1.1.1 1999-04-09 17:59:03 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.2  1998/03/14 00:04:47  tveldhui | 
|---|
| 27 | * 0.2-alpha-05 | 
|---|
| 28 | * | 
|---|
| 29 | * Revision 1.1  1997/07/16 14:51:20  tveldhui | 
|---|
| 30 | * Update: Alpha release 0.2 (Arrays) | 
|---|
| 31 | * | 
|---|
| 32 | */ | 
|---|
| 33 |  | 
|---|
| 34 | /* | 
|---|
| 35 | * ArrayIndexMapping is used to implement tensor array notation.  For | 
|---|
| 36 | * example: | 
|---|
| 37 | * | 
|---|
| 38 | * Array<float, 2> A, B; | 
|---|
| 39 | * firstIndex i; | 
|---|
| 40 | * secondIndex j; | 
|---|
| 41 | * thirdIndex k; | 
|---|
| 42 | * Array<float, 3> C = A(i,j) * B(j,k); | 
|---|
| 43 | * | 
|---|
| 44 | * For expression templates purposes, something like B(j,k) is represented | 
|---|
| 45 | * by an instance of class ArrayIndexMapping.  This class maps an array onto | 
|---|
| 46 | * the destination array coordinate system, e.g. B(j,k) -> C(i,j,k) | 
|---|
| 47 | */ | 
|---|
| 48 |  | 
|---|
| 49 | #ifndef BZ_ARRAYMAP_H | 
|---|
| 50 | #define BZ_ARRAYMAP_H | 
|---|
| 51 |  | 
|---|
| 52 | #ifndef BZ_ARRAY_H | 
|---|
| 53 | #error <blitz/array/map.h> must be included via <blitz/array.h> | 
|---|
| 54 | #endif | 
|---|
| 55 |  | 
|---|
| 56 | BZ_NAMESPACE(blitz) | 
|---|
| 57 |  | 
|---|
| 58 | /* | 
|---|
| 59 | * _bz_doArrayIndexMapping is a helper class.  It is specialized for | 
|---|
| 60 | * ranks 1, 2, 3, ..., 11. | 
|---|
| 61 | */ | 
|---|
| 62 |  | 
|---|
| 63 | template<int N_rank> | 
|---|
| 64 | struct _bz_doArrayIndexMapping { | 
|---|
| 65 | template<class T_numtype, int N_destRank> | 
|---|
| 66 | static T_numtype map(const Array<T_numtype, N_rank>&, | 
|---|
| 67 | const TinyVector<int,N_destRank>&, int, int, int, int, int, int, | 
|---|
| 68 | int, int, int, int, int) | 
|---|
| 69 | { | 
|---|
| 70 | // If you try to use an array index mapping on an array with | 
|---|
| 71 | // rank greater than 11, then you'll get a precondition failure | 
|---|
| 72 | // here. | 
|---|
| 73 | BZPRECONDITION(0); | 
|---|
| 74 | } | 
|---|
| 75 | }; | 
|---|
| 76 |  | 
|---|
| 77 | template<> | 
|---|
| 78 | struct _bz_doArrayIndexMapping<1> { | 
|---|
| 79 | template<class T_numtype, int N_destRank> | 
|---|
| 80 | static T_numtype map(const Array<T_numtype, 1>& array, | 
|---|
| 81 | const TinyVector<int,N_destRank>& index, int i0, int, int, int, int, | 
|---|
| 82 | int, int, int, int, int, int) | 
|---|
| 83 | { | 
|---|
| 84 | return array(index[i0]); | 
|---|
| 85 | } | 
|---|
| 86 | }; | 
|---|
| 87 |  | 
|---|
| 88 |  | 
|---|
| 89 | template<> | 
|---|
| 90 | struct _bz_doArrayIndexMapping<2> { | 
|---|
| 91 | template<class T_numtype, int N_destRank> | 
|---|
| 92 | static T_numtype map(const Array<T_numtype, 2>& array, | 
|---|
| 93 | const TinyVector<int,N_destRank>& index, int i0, int i1, int, | 
|---|
| 94 | int, int, int, int, int, int, int, int) | 
|---|
| 95 | { | 
|---|
| 96 | return array(index[i0], index[i1]); | 
|---|
| 97 | } | 
|---|
| 98 | }; | 
|---|
| 99 |  | 
|---|
| 100 | template<> | 
|---|
| 101 | struct _bz_doArrayIndexMapping<3> { | 
|---|
| 102 | template<class T_numtype, int N_destRank> | 
|---|
| 103 | static T_numtype map(const Array<T_numtype, 3>& array, | 
|---|
| 104 | const TinyVector<int,N_destRank>& index, int i0, int i1, int i2, | 
|---|
| 105 | int, int, int, int, int, int, int, int) | 
|---|
| 106 | { | 
|---|
| 107 | return array(index[i0], index[i1], index[i2]); | 
|---|
| 108 | } | 
|---|
| 109 | }; | 
|---|
| 110 |  | 
|---|
| 111 | template<> | 
|---|
| 112 | struct _bz_doArrayIndexMapping<4> { | 
|---|
| 113 | template<class T_numtype, int N_destRank> | 
|---|
| 114 | static T_numtype map(const Array<T_numtype, 4>& array, | 
|---|
| 115 | const TinyVector<int,N_destRank>& index, int i0, int i1, int i2, | 
|---|
| 116 | int i3, int, int, int, int, int, int, int) | 
|---|
| 117 | { | 
|---|
| 118 | return array(index[i0], index[i1], index[i2], index[i3]); | 
|---|
| 119 | } | 
|---|
| 120 | }; | 
|---|
| 121 |  | 
|---|
| 122 | template<> | 
|---|
| 123 | struct _bz_doArrayIndexMapping<5> { | 
|---|
| 124 | template<class T_numtype, int N_destRank> | 
|---|
| 125 | static T_numtype map(const Array<T_numtype, 5>& array, | 
|---|
| 126 | const TinyVector<int,N_destRank>& index, int i0, int i1, int i2, | 
|---|
| 127 | int i3, int i4, int, int, int, int, int, int) | 
|---|
| 128 | { | 
|---|
| 129 | return array(index[i0], index[i1], index[i2], index[i3], index[i4]); | 
|---|
| 130 | } | 
|---|
| 131 | }; | 
|---|
| 132 |  | 
|---|
| 133 | template<> | 
|---|
| 134 | struct _bz_doArrayIndexMapping<6> { | 
|---|
| 135 | template<class T_numtype, int N_destRank> | 
|---|
| 136 | static T_numtype map(const Array<T_numtype, 6>& array, | 
|---|
| 137 | const TinyVector<int,N_destRank>& index, int i0, int i1, int i2, | 
|---|
| 138 | int i3, int i4, int i5, int, int, int, int, int) | 
|---|
| 139 | { | 
|---|
| 140 | return array(index[i0], index[i1], index[i2], index[i3], index[i4], | 
|---|
| 141 | index[i5]); | 
|---|
| 142 | } | 
|---|
| 143 | }; | 
|---|
| 144 |  | 
|---|
| 145 | template<> | 
|---|
| 146 | struct _bz_doArrayIndexMapping<7> { | 
|---|
| 147 | template<class T_numtype, int N_destRank> | 
|---|
| 148 | static T_numtype map(const Array<T_numtype, 7>& array, | 
|---|
| 149 | const TinyVector<int,N_destRank>& index, int i0, int i1, int i2, | 
|---|
| 150 | int i3, int i4, int i5, int i6, int, int, int, int) | 
|---|
| 151 | { | 
|---|
| 152 | return array(index[i0], index[i1], index[i2], index[i3], index[i4], | 
|---|
| 153 | index[i5], index[i6]); | 
|---|
| 154 | } | 
|---|
| 155 | }; | 
|---|
| 156 |  | 
|---|
| 157 | template<> | 
|---|
| 158 | struct _bz_doArrayIndexMapping<8> { | 
|---|
| 159 | template<class T_numtype, int N_destRank> | 
|---|
| 160 | static T_numtype map(const Array<T_numtype, 8>& array, | 
|---|
| 161 | const TinyVector<int,N_destRank>& index, int i0, int i1, int i2, | 
|---|
| 162 | int i3, int i4, int i5, int i6, int i7, int, int, int) | 
|---|
| 163 | { | 
|---|
| 164 | return array(index[i0], index[i1], index[i2], index[i3], index[i4], | 
|---|
| 165 | index[i5], index[i6], index[i7]); | 
|---|
| 166 | } | 
|---|
| 167 | }; | 
|---|
| 168 |  | 
|---|
| 169 | template<> | 
|---|
| 170 | struct _bz_doArrayIndexMapping<9> { | 
|---|
| 171 | template<class T_numtype, int N_destRank> | 
|---|
| 172 | static T_numtype map(const Array<T_numtype, 9>& array, | 
|---|
| 173 | const TinyVector<int,N_destRank>& index, int i0, int i1, int i2, | 
|---|
| 174 | int i3, int i4, int i5, int i6, int i7, int i8, int, int) | 
|---|
| 175 | { | 
|---|
| 176 | return array(index[i0], index[i1], index[i2], index[i3], index[i4], | 
|---|
| 177 | index[i5], index[i6], index[i7], index[i8]); | 
|---|
| 178 | } | 
|---|
| 179 | }; | 
|---|
| 180 |  | 
|---|
| 181 | template<> | 
|---|
| 182 | struct _bz_doArrayIndexMapping<10> { | 
|---|
| 183 | template<class T_numtype, int N_destRank> | 
|---|
| 184 | static T_numtype map(const Array<T_numtype, 10>& array, | 
|---|
| 185 | const TinyVector<int,N_destRank>& index, int i0, int i1, int i2, | 
|---|
| 186 | int i3, int i4, int i5, int i6, int i7, int i8, int i9, int) | 
|---|
| 187 | { | 
|---|
| 188 | return array(index[i0], index[i1], index[i2], index[i3], index[i4], | 
|---|
| 189 | index[i5], index[i6], index[i7], index[i8], index[i9]); | 
|---|
| 190 | } | 
|---|
| 191 | }; | 
|---|
| 192 |  | 
|---|
| 193 | template<> | 
|---|
| 194 | struct _bz_doArrayIndexMapping<11> { | 
|---|
| 195 | template<class T_numtype, int N_destRank> | 
|---|
| 196 | static T_numtype map(const Array<T_numtype, 11>& array, | 
|---|
| 197 | const TinyVector<int,N_destRank>& index, int i0, int i1, int i2, | 
|---|
| 198 | int i3, int i4, int i5, int i6, int i7, int i8, int i9, int i10) | 
|---|
| 199 | { | 
|---|
| 200 | return array(index[i0], index[i1], index[i2], index[i3], index[i4], | 
|---|
| 201 | index[i5], index[i6], index[i7], index[i8], index[i9], | 
|---|
| 202 | index[i10]); | 
|---|
| 203 | } | 
|---|
| 204 | }; | 
|---|
| 205 |  | 
|---|
| 206 | template<class P_numtype, int N_rank, int N_map0, int N_map1=0, int N_map2=0, | 
|---|
| 207 | int N_map3=0, int N_map4=0, int N_map5=0, int N_map6=0, int N_map7=0, | 
|---|
| 208 | int N_map8=0, int N_map9=0, int N_map10=0> | 
|---|
| 209 | class ArrayIndexMapping { | 
|---|
| 210 | public: | 
|---|
| 211 | typedef P_numtype T_numtype; | 
|---|
| 212 | typedef const Array<T_numtype,N_rank>& T_ctorArg1; | 
|---|
| 213 | typedef int                            T_ctorArg2;    // dummy | 
|---|
| 214 |  | 
|---|
| 215 | /* | 
|---|
| 216 | * This enum block finds the maximum of the N_map0, N_map1, ..., N_map10 | 
|---|
| 217 | * parameters and stores it in maxRank10.  The rank of the expression is | 
|---|
| 218 | * then maxRank10 + 1, since the IndexPlaceholders start at 0 rather than | 
|---|
| 219 | * 1. | 
|---|
| 220 | */ | 
|---|
| 221 | enum { | 
|---|
| 222 | maxRank1 = (N_map0 > N_map1) ? N_map0 : N_map1, | 
|---|
| 223 | maxRank2 = (N_map2 > maxRank1) ? N_map2 : maxRank1, | 
|---|
| 224 | maxRank3 = (N_map3 > maxRank2) ? N_map3 : maxRank2, | 
|---|
| 225 | maxRank4 = (N_map4 > maxRank3) ? N_map4 : maxRank3, | 
|---|
| 226 | maxRank5 = (N_map5 > maxRank4) ? N_map5 : maxRank4, | 
|---|
| 227 | maxRank6 = (N_map6 > maxRank5) ? N_map6 : maxRank5, | 
|---|
| 228 | maxRank7 = (N_map7 > maxRank6) ? N_map7 : maxRank6, | 
|---|
| 229 | maxRank8 = (N_map8 > maxRank7) ? N_map8 : maxRank7, | 
|---|
| 230 | maxRank9 = (N_map9 > maxRank8) ? N_map9 : maxRank8, | 
|---|
| 231 | maxRank10 = (N_map10 > maxRank9) ? N_map10 : maxRank9 | 
|---|
| 232 | }; | 
|---|
| 233 |  | 
|---|
| 234 | enum { numArrayOperands = 1, numIndexPlaceholders = 1, | 
|---|
| 235 | rank = maxRank10 + 1 }; | 
|---|
| 236 |  | 
|---|
| 237 | ArrayIndexMapping(const Array<T_numtype, N_rank>& array) | 
|---|
| 238 | : array_(array) | 
|---|
| 239 | { | 
|---|
| 240 | } | 
|---|
| 241 |  | 
|---|
| 242 | ArrayIndexMapping(const ArrayIndexMapping<T_numtype,N_rank,N_map0, | 
|---|
| 243 | N_map1,N_map2,N_map3,N_map4,N_map5,N_map6,N_map7,N_map8,N_map9, | 
|---|
| 244 | N_map10>& z) | 
|---|
| 245 | : array_(z.array_) | 
|---|
| 246 | { | 
|---|
| 247 | } | 
|---|
| 248 |  | 
|---|
| 249 | #ifdef BZ_ARRAY_EXPR_PASS_INDEX_BY_VALUE | 
|---|
| 250 | template<int N_destRank> | 
|---|
| 251 | T_numtype operator()(TinyVector<int, N_destRank> i) | 
|---|
| 252 | { | 
|---|
| 253 | return _bz_doArrayIndexMapping<N_rank>::map(array_, i, | 
|---|
| 254 | N_map0, N_map1, N_map2, N_map3, N_map4, N_map5, N_map6, | 
|---|
| 255 | N_map7, N_map8, N_map9, N_map10); | 
|---|
| 256 | } | 
|---|
| 257 | #else | 
|---|
| 258 | template<int N_destRank> | 
|---|
| 259 | T_numtype operator()(const TinyVector<int, N_destRank>& i) | 
|---|
| 260 | { | 
|---|
| 261 | return _bz_doArrayIndexMapping<N_rank>::map(array_, i, | 
|---|
| 262 | N_map0, N_map1, N_map2, N_map3, N_map4, N_map5, N_map6, | 
|---|
| 263 | N_map7, N_map8, N_map9, N_map10); | 
|---|
| 264 | } | 
|---|
| 265 | #endif | 
|---|
| 266 |  | 
|---|
| 267 | int lbound(int rank) | 
|---|
| 268 | { | 
|---|
| 269 | if (N_map0 == rank) | 
|---|
| 270 | return array_.lbound(0); | 
|---|
| 271 | else if ((N_map1 == rank) && (N_rank > 1)) | 
|---|
| 272 | return array_.lbound(1); | 
|---|
| 273 | else if ((N_map2 == rank) && (N_rank > 2)) | 
|---|
| 274 | return array_.lbound(2); | 
|---|
| 275 | else if ((N_map3 == rank) && (N_rank > 3)) | 
|---|
| 276 | return array_.lbound(3); | 
|---|
| 277 | else if ((N_map4 == rank) && (N_rank > 4)) | 
|---|
| 278 | return array_.lbound(4); | 
|---|
| 279 | else if ((N_map5 == rank) && (N_rank > 5)) | 
|---|
| 280 | return array_.lbound(5); | 
|---|
| 281 | else if ((N_map6 == rank) && (N_rank > 6)) | 
|---|
| 282 | return array_.lbound(6); | 
|---|
| 283 | else if ((N_map7 == rank) && (N_rank > 7)) | 
|---|
| 284 | return array_.lbound(7); | 
|---|
| 285 | else if ((N_map8 == rank) && (N_rank > 8)) | 
|---|
| 286 | return array_.lbound(8); | 
|---|
| 287 | else if ((N_map9 == rank) && (N_rank > 9)) | 
|---|
| 288 | return array_.lbound(9); | 
|---|
| 289 | else if ((N_map10 == rank) && (N_rank > 10)) | 
|---|
| 290 | return array_.lbound(10); | 
|---|
| 291 | else | 
|---|
| 292 | return INT_MIN;   // tiny(int()); | 
|---|
| 293 | } | 
|---|
| 294 |  | 
|---|
| 295 | int ubound(int rank) | 
|---|
| 296 | { | 
|---|
| 297 | if (N_map0 == rank) | 
|---|
| 298 | return array_.ubound(0); | 
|---|
| 299 | else if ((N_map1 == rank) && (N_rank > 1)) | 
|---|
| 300 | return array_.ubound(1); | 
|---|
| 301 | else if ((N_map2 == rank) && (N_rank > 2)) | 
|---|
| 302 | return array_.ubound(2); | 
|---|
| 303 | else if ((N_map3 == rank) && (N_rank > 3)) | 
|---|
| 304 | return array_.ubound(3); | 
|---|
| 305 | else if ((N_map4 == rank) && (N_rank > 4)) | 
|---|
| 306 | return array_.ubound(4); | 
|---|
| 307 | else if ((N_map5 == rank) && (N_rank > 5)) | 
|---|
| 308 | return array_.ubound(5); | 
|---|
| 309 | else if ((N_map6 == rank) && (N_rank > 6)) | 
|---|
| 310 | return array_.ubound(6); | 
|---|
| 311 | else if ((N_map7 == rank) && (N_rank > 7)) | 
|---|
| 312 | return array_.ubound(7); | 
|---|
| 313 | else if ((N_map8 == rank) && (N_rank > 8)) | 
|---|
| 314 | return array_.ubound(8); | 
|---|
| 315 | else if ((N_map9 == rank) && (N_rank > 9)) | 
|---|
| 316 | return array_.ubound(9); | 
|---|
| 317 | else if ((N_map10 == rank) && (N_rank > 10)) | 
|---|
| 318 | return array_.ubound(10); | 
|---|
| 319 | else | 
|---|
| 320 | return INT_MAX;   // huge(int()); | 
|---|
| 321 | } | 
|---|
| 322 |  | 
|---|
| 323 | // If you have a precondition failure on this routine, it means | 
|---|
| 324 | // you are trying to use stack iteration mode on an expression | 
|---|
| 325 | // which contains an index placeholder.  You must use index | 
|---|
| 326 | // iteration mode instead. | 
|---|
| 327 | int operator*() | 
|---|
| 328 | { | 
|---|
| 329 | BZPRECONDITION(0); | 
|---|
| 330 | return 0; | 
|---|
| 331 | } | 
|---|
| 332 |  | 
|---|
| 333 | // See operator*() note | 
|---|
| 334 | void push(int) | 
|---|
| 335 | { | 
|---|
| 336 | BZPRECONDITION(0); | 
|---|
| 337 | } | 
|---|
| 338 |  | 
|---|
| 339 | // See operator*() note | 
|---|
| 340 | void pop(int) | 
|---|
| 341 | { | 
|---|
| 342 | BZPRECONDITION(0); | 
|---|
| 343 | } | 
|---|
| 344 |  | 
|---|
| 345 | // See operator*() note | 
|---|
| 346 | void advance() | 
|---|
| 347 | { | 
|---|
| 348 | BZPRECONDITION(0); | 
|---|
| 349 | } | 
|---|
| 350 |  | 
|---|
| 351 | // See operator*() note | 
|---|
| 352 | void advance(int) | 
|---|
| 353 | { | 
|---|
| 354 | BZPRECONDITION(0); | 
|---|
| 355 | } | 
|---|
| 356 |  | 
|---|
| 357 | // See operator*() note | 
|---|
| 358 | void loadStride(int) | 
|---|
| 359 | { | 
|---|
| 360 | BZPRECONDITION(0); | 
|---|
| 361 | } | 
|---|
| 362 |  | 
|---|
| 363 | _bz_bool isUnitStride(int rank) const | 
|---|
| 364 | { | 
|---|
| 365 | BZPRECONDITION(0); | 
|---|
| 366 | return false; | 
|---|
| 367 | } | 
|---|
| 368 |  | 
|---|
| 369 | void advanceUnitStride() | 
|---|
| 370 | { | 
|---|
| 371 | BZPRECONDITION(0); | 
|---|
| 372 | } | 
|---|
| 373 |  | 
|---|
| 374 | _bz_bool canCollapse(int,int) const | 
|---|
| 375 | {   BZPRECONDITION(0);  return _bz_false; } | 
|---|
| 376 |  | 
|---|
| 377 | T_numtype operator[](int) | 
|---|
| 378 | { | 
|---|
| 379 | BZPRECONDITION(0); | 
|---|
| 380 | return T_numtype(); | 
|---|
| 381 | } | 
|---|
| 382 |  | 
|---|
| 383 | T_numtype fastRead(int) | 
|---|
| 384 | { | 
|---|
| 385 | BZPRECONDITION(0); | 
|---|
| 386 | return T_numtype(); | 
|---|
| 387 | } | 
|---|
| 388 |  | 
|---|
| 389 | int suggestStride(int) const | 
|---|
| 390 | { | 
|---|
| 391 | BZPRECONDITION(0); | 
|---|
| 392 | return 0; | 
|---|
| 393 | } | 
|---|
| 394 |  | 
|---|
| 395 | _bz_bool isStride(int,int) const | 
|---|
| 396 | { | 
|---|
| 397 | BZPRECONDITION(0); | 
|---|
| 398 | return _bz_true; | 
|---|
| 399 | } | 
|---|
| 400 |  | 
|---|
| 401 | template<int N_rank2> | 
|---|
| 402 | void moveTo(const TinyVector<int,N_rank2>& i) | 
|---|
| 403 | { | 
|---|
| 404 | BZPRECONDITION(0); | 
|---|
| 405 | return ; | 
|---|
| 406 | } | 
|---|
| 407 |  | 
|---|
| 408 | void prettyPrint(string& str, prettyPrintFormat& format) const | 
|---|
| 409 | { | 
|---|
| 410 | // NEEDS_WORK-- do real formatting for reductions | 
|---|
| 411 | str += "map[NEEDS_WORK]"; | 
|---|
| 412 | } | 
|---|
| 413 |  | 
|---|
| 414 | template<class T_shape> | 
|---|
| 415 | _bz_bool shapeCheck(const T_shape& shape) const | 
|---|
| 416 | { | 
|---|
| 417 | // NEEDS_WORK-- do a real shape check (tricky) | 
|---|
| 418 | return _bz_true; | 
|---|
| 419 | } | 
|---|
| 420 |  | 
|---|
| 421 | private: | 
|---|
| 422 | ArrayIndexMapping() { } | 
|---|
| 423 |  | 
|---|
| 424 | const Array<T_numtype, N_rank>& array_; | 
|---|
| 425 | }; | 
|---|
| 426 |  | 
|---|
| 427 | BZ_NAMESPACE_END | 
|---|
| 428 |  | 
|---|
| 429 | #endif // BZ_ARRAYMAP_H | 
|---|
| 430 |  | 
|---|