| 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 | 
 | 
|---|