[369] | 1 | #ifndef __XIOS_ARRAY_HPP__ |
---|
| 2 | #define __XIOS_ARRAY_HPP__ |
---|
| 3 | |
---|
| 4 | #define BZ_COLUMN_MAJOR_ARRAY |
---|
| 5 | #include <blitz/array.h> |
---|
| 6 | #include "buffer_in.hpp" |
---|
| 7 | #include "buffer_out.hpp" |
---|
| 8 | #include "message.hpp" |
---|
[464] | 9 | #include <cmath> |
---|
[369] | 10 | |
---|
[677] | 11 | using namespace blitz; |
---|
[464] | 12 | BZ_DECLARE_FUNCTION(round) |
---|
[677] | 13 | |
---|
[369] | 14 | namespace xios |
---|
| 15 | { |
---|
[677] | 16 | template <typename T_numtype,int N_rank> |
---|
| 17 | class CArray : public Array<T_numtype,N_rank>, public virtual CBaseType |
---|
| 18 | { |
---|
| 19 | private: |
---|
| 20 | bool initialized; |
---|
[369] | 21 | |
---|
[677] | 22 | public: |
---|
| 23 | typedef typename Array<T_numtype,N_rank>::T_default_storage T_default_storage; |
---|
| 24 | using Array<T_numtype,N_rank>::operator =; |
---|
[369] | 25 | |
---|
[677] | 26 | template<typename T_expr> explicit CArray(_bz_ArrayExpr<T_expr> expr) |
---|
| 27 | : Array<T_numtype,N_rank>(expr) |
---|
| 28 | , initialized(true) |
---|
| 29 | {} |
---|
[369] | 30 | |
---|
[677] | 31 | CArray(GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 32 | : Array<T_numtype,N_rank>(storage) |
---|
| 33 | , initialized(false) |
---|
| 34 | {} |
---|
[369] | 35 | |
---|
[677] | 36 | explicit CArray(int length0, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 37 | : Array<T_numtype,N_rank>(length0, storage) |
---|
| 38 | , initialized(true) |
---|
| 39 | {} |
---|
[369] | 40 | |
---|
[677] | 41 | CArray(int length0, int length1, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 42 | : Array<T_numtype,N_rank>(length0, length1, storage) |
---|
| 43 | , initialized(true) |
---|
| 44 | {} |
---|
[369] | 45 | |
---|
[677] | 46 | CArray(int length0, int length1, int length2, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 47 | : Array<T_numtype,N_rank>(length0, length1, length2, storage) |
---|
| 48 | , initialized(true) |
---|
| 49 | {} |
---|
[369] | 50 | |
---|
[677] | 51 | CArray(int length0, int length1, int length2, int length3, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 52 | : Array<T_numtype,N_rank>(length0, length1, length2, length3, storage) |
---|
| 53 | , initialized(true) |
---|
| 54 | {} |
---|
[369] | 55 | |
---|
[677] | 56 | CArray(int length0, int length1, int length2, int length3, int length4, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 57 | : Array<T_numtype,N_rank>(length0,length1, length2, length3, length4, storage) |
---|
| 58 | , initialized(true) |
---|
| 59 | {} |
---|
[369] | 60 | |
---|
[677] | 61 | CArray(int length0, int length1, int length2, int length3, int length4, int length5, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 62 | : Array<T_numtype,N_rank>(length0, length1, length2, length3, length4, length5, storage) |
---|
| 63 | , initialized(true) |
---|
| 64 | {} |
---|
[369] | 65 | |
---|
[677] | 66 | CArray(int length0, int length1, int length2, int length3, int length4, int length5, int length6, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 67 | : Array<T_numtype,N_rank>(length0, length1, length2, length3, length4, length5, length6, storage) |
---|
| 68 | , initialized(true) |
---|
| 69 | {} |
---|
[369] | 70 | |
---|
[677] | 71 | CArray(int length0, int length1, int length2, int length3, int length4, int length5, int length6, int length7, |
---|
| 72 | GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 73 | : Array<T_numtype,N_rank>(length0, length1, length2, length3, length4, length5, length6, length7, storage) |
---|
| 74 | , initialized(true) |
---|
| 75 | {} |
---|
[369] | 76 | |
---|
[677] | 77 | CArray(int length0, int length1, int length2, int length3, int length4, int length5, int length6, |
---|
| 78 | int length7, int length8, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 79 | : Array<T_numtype,N_rank>(length0, length1, length2, length3, length4, length5, length6, length7, length8, storage) |
---|
| 80 | , initialized(true) |
---|
| 81 | {} |
---|
[369] | 82 | |
---|
[677] | 83 | CArray(int length0, int length1, int length2, int length3, int length4, |
---|
| 84 | int length5, int length6, int length7, int length8, int length9, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 85 | : Array<T_numtype,N_rank>(length0, length1, length2, length3, length4, length5, length6, length7, length8, length9, storage) |
---|
| 86 | , initialized(true) |
---|
| 87 | {} |
---|
[369] | 88 | |
---|
[677] | 89 | CArray(int length0, int length1, int length2, int length3, int length4, int length5, int length6, |
---|
| 90 | int length7, int length8, int length9, int length10, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 91 | : Array<T_numtype,N_rank>(length0, length1, length2, length3, length4, length5, length6, length7, length8, length9, length10, storage) |
---|
| 92 | , initialized(true) |
---|
| 93 | {} |
---|
[369] | 94 | |
---|
[677] | 95 | CArray(T_numtype* restrict dataFirst, TinyVector<int, N_rank> shape, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 96 | : Array<T_numtype,N_rank>(dataFirst, shape, storage) |
---|
| 97 | , initialized(true) |
---|
| 98 | {} |
---|
[369] | 99 | |
---|
[677] | 100 | CArray(T_numtype* restrict dataFirst, TinyVector<int, N_rank> shape, TinyVector<diffType, N_rank> stride, |
---|
| 101 | GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 102 | : Array<T_numtype,N_rank>(dataFirst, shape, stride, storage) |
---|
| 103 | , initialized(true) |
---|
| 104 | {} |
---|
[369] | 105 | |
---|
[677] | 106 | CArray(T_numtype* restrict dataFirst, TinyVector<int, N_rank> shape, preexistingMemoryPolicy deletionPolicy, |
---|
| 107 | GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 108 | : Array<T_numtype,N_rank>(dataFirst, shape, deletionPolicy, storage) |
---|
| 109 | , initialized(true) |
---|
| 110 | {} |
---|
[369] | 111 | |
---|
[677] | 112 | CArray(T_numtype* restrict dataFirst, TinyVector<int, N_rank> shape, TinyVector<diffType, N_rank> stride, |
---|
| 113 | preexistingMemoryPolicy deletionPolicy, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 114 | : Array<T_numtype,N_rank>(dataFirst, shape, stride, deletionPolicy, storage) |
---|
| 115 | , initialized(true) |
---|
| 116 | {} |
---|
[369] | 117 | |
---|
[677] | 118 | CArray(const TinyVector<int, N_rank>& extent, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 119 | : Array<T_numtype,N_rank>(extent, storage) |
---|
| 120 | , initialized(true) |
---|
| 121 | {} |
---|
[369] | 122 | |
---|
[677] | 123 | CArray(const TinyVector<int, N_rank>& lbounds, const TinyVector<int, N_rank>& extent, |
---|
| 124 | const GeneralArrayStorage<N_rank>& storage) |
---|
| 125 | : Array<T_numtype,N_rank>(lbounds, extent, storage) |
---|
| 126 | , initialized(true) |
---|
| 127 | {} |
---|
[369] | 128 | |
---|
[677] | 129 | CArray(Range r0, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 130 | : Array<T_numtype,N_rank>(r0, storage) |
---|
| 131 | , initialized(true) |
---|
| 132 | {} |
---|
[369] | 133 | |
---|
[677] | 134 | CArray(Range r0, Range r1, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 135 | : Array<T_numtype,N_rank>(r0, r1, storage) |
---|
| 136 | , initialized(true) |
---|
| 137 | {} |
---|
[369] | 138 | |
---|
[677] | 139 | CArray(Range r0, Range r1, Range r2, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 140 | : Array<T_numtype,N_rank>(r0, r1, r2, storage) |
---|
| 141 | , initialized(true) |
---|
| 142 | {} |
---|
[369] | 143 | |
---|
[677] | 144 | CArray(Range r0, Range r1, Range r2, Range r3, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 145 | : Array<T_numtype,N_rank>(r0, r1, r2, r3, storage) |
---|
| 146 | , initialized(true) |
---|
| 147 | {} |
---|
[369] | 148 | |
---|
[677] | 149 | CArray(Range r0, Range r1, Range r2, Range r3, Range r4, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 150 | : Array<T_numtype,N_rank>(r0, r1, r2, r3, r4, storage) |
---|
| 151 | , initialized(true) |
---|
| 152 | {} |
---|
[369] | 153 | |
---|
[677] | 154 | CArray(Range r0, Range r1, Range r2, Range r3, Range r4, Range r5, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 155 | : Array<T_numtype,N_rank>(r0, r1, r2, r3, r4, r5, storage) |
---|
| 156 | , initialized(true) |
---|
| 157 | {} |
---|
[369] | 158 | |
---|
[677] | 159 | CArray(Range r0, Range r1, Range r2, Range r3, Range r4, Range r5, Range r6, |
---|
| 160 | GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 161 | : Array<T_numtype,N_rank>(r0, r1, r2, r3, r4, r5, r6, storage) |
---|
| 162 | , initialized(true) |
---|
| 163 | {} |
---|
[369] | 164 | |
---|
[677] | 165 | CArray(Range r0, Range r1, Range r2, Range r3, Range r4, Range r5, Range r6, Range r7, |
---|
| 166 | GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 167 | : Array<T_numtype,N_rank>(r0, r1, r2, r3, r4, r5, r6, r7, storage) |
---|
| 168 | , initialized(true) |
---|
| 169 | {} |
---|
[369] | 170 | |
---|
[677] | 171 | CArray(Range r0, Range r1, Range r2, Range r3, Range r4, Range r5, |
---|
| 172 | Range r6, Range r7, Range r8, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 173 | : Array<T_numtype,N_rank>(r0, r1, r2, r3, r4, r5, r6, r7, r8, storage) |
---|
| 174 | , initialized(true) |
---|
| 175 | {} |
---|
[369] | 176 | |
---|
[677] | 177 | CArray(Range r0, Range r1, Range r2, Range r3, Range r4, Range r5, |
---|
| 178 | Range r6, Range r7, Range r8, Range r9, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 179 | : Array<T_numtype,N_rank>(r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, storage) |
---|
| 180 | , initialized(true) |
---|
| 181 | {} |
---|
[369] | 182 | |
---|
[677] | 183 | CArray(Range r0, Range r1, Range r2, Range r3, Range r4, Range r5, Range r6, Range r7, |
---|
| 184 | Range r8, Range r9, Range r10, GeneralArrayStorage<N_rank> storage = T_default_storage()) |
---|
| 185 | : Array<T_numtype,N_rank>(r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10, storage) |
---|
| 186 | , initialized(true) |
---|
| 187 | {} |
---|
[369] | 188 | |
---|
[677] | 189 | CArray(const CArray<T_numtype, N_rank>& array) |
---|
| 190 | : Array<T_numtype,N_rank>(array) |
---|
| 191 | , initialized(array.initialized) |
---|
| 192 | {} |
---|
[369] | 193 | |
---|
[677] | 194 | CArray(const Array<T_numtype, N_rank>& array) |
---|
| 195 | : Array<T_numtype,N_rank>(array) |
---|
| 196 | , initialized(true) |
---|
| 197 | {} |
---|
[369] | 198 | |
---|
[677] | 199 | CArray(const TinyVector<int,N_rank-1>& shape, int lastExtent, const GeneralArrayStorage<N_rank>& storage) |
---|
| 200 | : Array<T_numtype,N_rank>(shape, lastExtent, storage) |
---|
| 201 | , initialized(true) |
---|
| 202 | {} |
---|
[369] | 203 | |
---|
[677] | 204 | CArray(Array<T_numtype, N_rank>& array, Range r0) |
---|
| 205 | : Array<T_numtype,N_rank>(array, r0) |
---|
| 206 | , initialized(true) |
---|
| 207 | {} |
---|
[369] | 208 | |
---|
[677] | 209 | CArray(Array<T_numtype, N_rank>& array, Range r0, Range r1) |
---|
| 210 | : Array<T_numtype,N_rank>(array, r0, r1) |
---|
| 211 | , initialized(true) |
---|
| 212 | {} |
---|
[369] | 213 | |
---|
[677] | 214 | CArray(Array<T_numtype, N_rank>& array, Range r0, Range r1, Range r2) |
---|
| 215 | : Array<T_numtype,N_rank>( array, r0, r1, r2) |
---|
| 216 | , initialized(true) |
---|
| 217 | {} |
---|
[369] | 218 | |
---|
[677] | 219 | CArray(Array<T_numtype, N_rank>& array, Range r0, Range r1, Range r2, Range r3) |
---|
| 220 | : Array<T_numtype,N_rank>(array, r0, r1, r2, r3) |
---|
| 221 | , initialized(true) |
---|
| 222 | {} |
---|
[369] | 223 | |
---|
[677] | 224 | CArray(Array<T_numtype, N_rank>& array, Range r0, Range r1, Range r2, |
---|
| 225 | Range r3, Range r4) |
---|
| 226 | : Array<T_numtype,N_rank>(array, r0, r1, r2, r3, r4) |
---|
| 227 | , initialized(true) |
---|
| 228 | {} |
---|
[369] | 229 | |
---|
[677] | 230 | CArray(Array<T_numtype, N_rank>& array, Range r0, Range r1, Range r2, |
---|
| 231 | Range r3, Range r4, Range r5) |
---|
| 232 | : Array<T_numtype,N_rank>( array, r0, r1, r2, r3, r4, r5) |
---|
| 233 | , initialized(true) |
---|
| 234 | {} |
---|
[369] | 235 | |
---|
[677] | 236 | CArray(Array<T_numtype, N_rank>& array, Range r0, Range r1, Range r2, Range r3, |
---|
| 237 | Range r4, Range r5, Range r6) |
---|
| 238 | : Array<T_numtype,N_rank>( array, r0, r1, r2, r3, r4, r5, r6) |
---|
| 239 | , initialized(true) |
---|
| 240 | {} |
---|
[369] | 241 | |
---|
[677] | 242 | CArray(Array<T_numtype, N_rank>& array, Range r0, Range r1, Range r2, Range r3, Range r4, |
---|
| 243 | Range r5, Range r6, Range r7) |
---|
| 244 | : Array<T_numtype,N_rank>(array, r0, r1, r2, r3, r4, r5, r6, r7) |
---|
| 245 | , initialized(true) |
---|
| 246 | {} |
---|
[369] | 247 | |
---|
[677] | 248 | CArray(Array<T_numtype, N_rank>& array, Range r0, Range r1, Range r2, Range r3, Range r4, Range r5, |
---|
| 249 | Range r6, Range r7, Range r8) |
---|
| 250 | : Array<T_numtype,N_rank>(array, r0, r1, r2, r3, r4, r5, r6, r7, r8) |
---|
| 251 | , initialized(true) |
---|
| 252 | {} |
---|
[369] | 253 | |
---|
[677] | 254 | CArray(Array<T_numtype, N_rank>& array, Range r0, Range r1, Range r2, Range r3, Range r4, Range r5, |
---|
| 255 | Range r6, Range r7, Range r8, Range r9) |
---|
| 256 | : Array<T_numtype,N_rank>(array, r0, r1, r2, r3, r4, r5, r6, r7, r8, r9) |
---|
| 257 | , initialized(true) |
---|
| 258 | {} |
---|
[369] | 259 | |
---|
[677] | 260 | CArray(Array<T_numtype, N_rank>& array, Range r0, Range r1, Range r2, Range r3, Range r4, Range r5, Range r6, |
---|
| 261 | Range r7, Range r8, Range r9, Range r10) |
---|
| 262 | : Array<T_numtype,N_rank>(array, r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10) |
---|
| 263 | , initialized(true) |
---|
| 264 | {} |
---|
| 265 | |
---|
| 266 | CArray(Array<T_numtype, N_rank>& array, const RectDomain<N_rank>& subdomain) |
---|
| 267 | : Array<T_numtype,N_rank>(array, subdomain) |
---|
| 268 | , initialized(true) |
---|
| 269 | {} |
---|
| 270 | |
---|
| 271 | CArray(Array<T_numtype, N_rank>& array, const StridedDomain<N_rank>& subdomain) |
---|
| 272 | : Array<T_numtype,N_rank>(array, subdomain) |
---|
| 273 | , initialized(true) |
---|
| 274 | {} |
---|
| 275 | |
---|
| 276 | template<int N_rank2, typename R0, typename R1, typename R2, typename R3, typename R4, typename R5, |
---|
| 277 | typename R6, typename R7, typename R8, typename R9, typename R10> |
---|
| 278 | CArray(Array<T_numtype,N_rank2>& array, R0 r0, R1 r1, R2 r2, R3 r3, R4 r4, R5 r5, R6 r6, R7 r7, R8 r8, R9 r9, R10 r10) |
---|
| 279 | : Array<T_numtype,N_rank>(array, r0,r1, r2, r3, r4, r5, r6, r7, r8, r9, r10) |
---|
| 280 | , initialized(true) |
---|
| 281 | {} |
---|
| 282 | |
---|
| 283 | virtual ~CArray() {} |
---|
| 284 | |
---|
| 285 | CArray<T_numtype, N_rank> copy() const |
---|
| 286 | { |
---|
| 287 | CArray<T_numtype, N_rank> copy = Array<T_numtype, N_rank>::copy(); |
---|
| 288 | copy.initialized = initialized; |
---|
| 289 | return copy; |
---|
| 290 | } |
---|
| 291 | |
---|
| 292 | void reference(const CArray<T_numtype, N_rank>& array) |
---|
| 293 | { |
---|
| 294 | Array<T_numtype,N_rank>::reference(array); |
---|
| 295 | initialized = array.initialized; |
---|
| 296 | } |
---|
| 297 | |
---|
| 298 | void reference(const Array<T_numtype,N_rank>& array) |
---|
| 299 | { |
---|
| 300 | Array<T_numtype,N_rank>::reference(array); |
---|
| 301 | initialized = true; |
---|
| 302 | } |
---|
| 303 | |
---|
[1158] | 304 | bool operator== (const CArray<T_numtype,N_rank>& array) |
---|
| 305 | { |
---|
| 306 | size_t nbThis = this->numElements(); |
---|
| 307 | size_t nbArr = array.numElements(); |
---|
| 308 | if (nbThis != nbArr) return false; |
---|
| 309 | typename Array<T_numtype,N_rank>::const_iterator itx=array.begin(), itxe=array.end(), ity=this->begin() ; |
---|
| 310 | for(;itx!=itxe;++itx,++ity) if (*itx!=*ity) return false ; |
---|
| 311 | return true; |
---|
| 312 | } |
---|
| 313 | |
---|
| 314 | bool operator== (const Array<T_numtype,N_rank>& array) |
---|
| 315 | { |
---|
| 316 | return ((*this) == (dynamic_cast<const CArray<T_numtype,N_rank>& >(array))); |
---|
| 317 | } |
---|
| 318 | |
---|
[677] | 319 | void resize(int extent) |
---|
| 320 | { |
---|
| 321 | Array<T_numtype,N_rank>::resize(extent); |
---|
| 322 | initialized = true; |
---|
| 323 | } |
---|
| 324 | void resize(int extent1, int extent2) |
---|
| 325 | { |
---|
| 326 | Array<T_numtype,N_rank>::resize(extent1, extent2); |
---|
| 327 | initialized = true; |
---|
| 328 | } |
---|
| 329 | void resize(int extent1, int extent2, int extent3) |
---|
| 330 | { |
---|
| 331 | Array<T_numtype,N_rank>::resize(extent1, extent2, extent3); |
---|
| 332 | initialized = true; |
---|
| 333 | } |
---|
| 334 | void resize(int extent1, int extent2, int extent3, int extent4) |
---|
| 335 | { |
---|
| 336 | Array<T_numtype,N_rank>::resize(extent1, extent2, extent3, extent4); |
---|
| 337 | initialized = true; |
---|
| 338 | } |
---|
| 339 | void resize(int extent1, int extent2, int extent3, int extent4, int extent5) |
---|
| 340 | { |
---|
| 341 | Array<T_numtype,N_rank>::resize(extent1, extent2, extent3, extent4, extent5); |
---|
| 342 | initialized = true; |
---|
| 343 | } |
---|
| 344 | void resize(int extent1, int extent2, int extent3, int extent4, int extent5, int extent6) |
---|
| 345 | { |
---|
| 346 | Array<T_numtype,N_rank>::resize(extent1, extent2, extent3, extent4, extent5, extent6); |
---|
| 347 | initialized = true; |
---|
| 348 | } |
---|
| 349 | void resize(int extent1, int extent2, int extent3, int extent4, int extent5, |
---|
| 350 | int extent6, int extent7) |
---|
| 351 | { |
---|
| 352 | Array<T_numtype,N_rank>::resize(extent1, extent2, extent3, extent4, extent5, |
---|
| 353 | extent6, extent7); |
---|
| 354 | initialized = true; |
---|
| 355 | } |
---|
| 356 | void resize(int extent1, int extent2, int extent3, int extent4, int extent5, |
---|
| 357 | int extent6, int extent7, int extent8) |
---|
| 358 | { |
---|
| 359 | Array<T_numtype,N_rank>::resize(extent1, extent2, extent3, extent4, extent5, |
---|
| 360 | extent6, extent7, extent8); |
---|
| 361 | initialized = true; |
---|
| 362 | } |
---|
| 363 | void resize(int extent1, int extent2, int extent3, int extent4, int extent5, |
---|
| 364 | int extent6, int extent7, int extent8, int extent9) |
---|
| 365 | { |
---|
| 366 | Array<T_numtype,N_rank>::resize(extent1, extent2, extent3, extent4, extent5, |
---|
| 367 | extent6, extent7, extent8, extent9); |
---|
| 368 | initialized = true; |
---|
| 369 | } |
---|
| 370 | void resize(int extent1, int extent2, int extent3, int extent4, int extent5, |
---|
| 371 | int extent6, int extent7, int extent8, int extent9, int extent10) |
---|
| 372 | { |
---|
| 373 | Array<T_numtype,N_rank>::resize(extent1, extent2, extent3, extent4, extent5, |
---|
| 374 | extent6, extent7, extent8, extent9, extent10); |
---|
| 375 | initialized = true; |
---|
| 376 | } |
---|
| 377 | void resize(int extent1, int extent2, int extent3, int extent4, int extent5, int extent6, |
---|
| 378 | int extent7, int extent8, int extent9, int extent10, int extent11) |
---|
| 379 | { |
---|
| 380 | Array<T_numtype,N_rank>::resize(extent1, extent2, extent3, extent4, extent5, extent6, |
---|
| 381 | extent7, extent8, extent9, extent10, extent11); |
---|
| 382 | initialized = true; |
---|
| 383 | } |
---|
| 384 | |
---|
| 385 | void resize(Range r1) |
---|
| 386 | { |
---|
| 387 | Array<T_numtype,N_rank>::resize(r1); |
---|
| 388 | initialized = true; |
---|
| 389 | } |
---|
| 390 | void resize(Range r1, Range r2) |
---|
| 391 | { |
---|
| 392 | Array<T_numtype,N_rank>::resize(r1, r2); |
---|
| 393 | initialized = true; |
---|
| 394 | } |
---|
| 395 | void resize(Range r1, Range r2, Range r3) |
---|
| 396 | { |
---|
| 397 | Array<T_numtype,N_rank>::resize(r1, r2, r3); |
---|
| 398 | initialized = true; |
---|
| 399 | } |
---|
| 400 | void resize(Range r1, Range r2, Range r3, Range r4) |
---|
| 401 | { |
---|
| 402 | Array<T_numtype,N_rank>::resize(r1, r2, r3, r4); |
---|
| 403 | initialized = true; |
---|
| 404 | } |
---|
| 405 | void resize(Range r1, Range r2, Range r3, Range r4, Range r5) |
---|
| 406 | { |
---|
| 407 | Array<T_numtype,N_rank>::resize(r1, r2, r3, r4, r5); |
---|
| 408 | initialized = true; |
---|
| 409 | } |
---|
| 410 | void resize(Range r1, Range r2, Range r3, Range r4, Range r5, Range r6) |
---|
| 411 | { |
---|
| 412 | Array<T_numtype,N_rank>::resize(r1, r2, r3, r4, r5, r6); |
---|
| 413 | initialized = true; |
---|
| 414 | } |
---|
| 415 | void resize(Range r1, Range r2, Range r3, Range r4, Range r5, |
---|
| 416 | Range r6, Range r7) |
---|
| 417 | { |
---|
| 418 | Array<T_numtype,N_rank>::resize(r1, r2, r3, r4, r5, |
---|
| 419 | r6, r7); |
---|
| 420 | initialized = true; |
---|
| 421 | } |
---|
| 422 | void resize(Range r1, Range r2, Range r3, Range r4, Range r5, |
---|
| 423 | Range r6, Range r7, Range r8) |
---|
| 424 | { |
---|
| 425 | Array<T_numtype,N_rank>::resize(r1, r2, r3, r4, r5, |
---|
| 426 | r6, r7, r8); |
---|
| 427 | initialized = true; |
---|
| 428 | } |
---|
| 429 | void resize(Range r1, Range r2, Range r3, Range r4, Range r5, |
---|
| 430 | Range r6, Range r7, Range r8, Range r9) |
---|
| 431 | { |
---|
| 432 | Array<T_numtype,N_rank>::resize(r1, r2, r3, r4, r5, |
---|
| 433 | r6, r7, r8, r9); |
---|
| 434 | initialized = true; |
---|
| 435 | } |
---|
| 436 | void resize(Range r1, Range r2, Range r3, Range r4, Range r5, |
---|
| 437 | Range r6, Range r7, Range r8, Range r9, Range r10) |
---|
| 438 | { |
---|
| 439 | Array<T_numtype,N_rank>::resize(r1, r2, r3, r4, r5, |
---|
| 440 | r6, r7, r8, r9, r10); |
---|
| 441 | initialized = true; |
---|
| 442 | } |
---|
| 443 | void resize(Range r1, Range r2, Range r3, Range r4, Range r5, Range r6, |
---|
| 444 | Range r7, Range r8, Range r9, Range r10, Range r11) |
---|
| 445 | { |
---|
| 446 | Array<T_numtype,N_rank>::resize(r1, r2, r3, r4, r5, r6, |
---|
| 447 | r7, r8, r9, r10, r11); |
---|
| 448 | initialized = true; |
---|
| 449 | } |
---|
| 450 | |
---|
| 451 | void resize(const TinyVector<int,N_rank>& extent) |
---|
| 452 | { |
---|
| 453 | Array<T_numtype,N_rank>::resize(extent); |
---|
| 454 | initialized = true; |
---|
| 455 | } |
---|
| 456 | |
---|
| 457 | void resizeAndPreserve(const TinyVector<int,N_rank>& extent) |
---|
| 458 | { |
---|
| 459 | Array<T_numtype,N_rank>::resizeAndPreserve(extent); |
---|
| 460 | initialized = true; |
---|
| 461 | } |
---|
| 462 | |
---|
| 463 | void resizeAndPreserve(int extent) |
---|
| 464 | { |
---|
| 465 | Array<T_numtype,N_rank>::resizeAndPreserve(extent); |
---|
| 466 | initialized = true; |
---|
| 467 | } |
---|
| 468 | void resizeAndPreserve(int extent1, int extent2) |
---|
| 469 | { |
---|
| 470 | Array<T_numtype,N_rank>::resizeAndPreserve(extent1, extent2); |
---|
| 471 | initialized = true; |
---|
| 472 | } |
---|
| 473 | void resizeAndPreserve(int extent1, int extent2, int extent3) |
---|
| 474 | { |
---|
| 475 | Array<T_numtype,N_rank>::resizeAndPreserve(extent1, extent2, extent3); |
---|
| 476 | initialized = true; |
---|
| 477 | } |
---|
| 478 | void resizeAndPreserve(int extent1, int extent2, int extent3, int extent4) |
---|
| 479 | { |
---|
| 480 | Array<T_numtype,N_rank>::resizeAndPreserve(extent1, extent2, extent3, extent4); |
---|
| 481 | initialized = true; |
---|
| 482 | } |
---|
| 483 | void resizeAndPreserve(int extent1, int extent2, int extent3, int extent4, int extent5) |
---|
| 484 | { |
---|
| 485 | Array<T_numtype,N_rank>::resizeAndPreserve(extent1, extent2, extent3, extent4, extent5); |
---|
| 486 | initialized = true; |
---|
| 487 | } |
---|
| 488 | void resizeAndPreserve(int extent1, int extent2, int extent3, int extent4, |
---|
| 489 | int extent5, int extent6) |
---|
| 490 | { |
---|
| 491 | Array<T_numtype,N_rank>::resizeAndPreserve(extent1, extent2, extent3, extent4, |
---|
| 492 | extent5, extent6); |
---|
| 493 | initialized = true; |
---|
| 494 | } |
---|
| 495 | void resizeAndPreserve(int extent1, int extent2, int extent3, int extent4, |
---|
| 496 | int extent5, int extent6, int extent7) |
---|
| 497 | { |
---|
| 498 | Array<T_numtype,N_rank>::resizeAndPreserve(extent1, extent2, extent3, extent4, |
---|
| 499 | extent5, extent6, extent7); |
---|
| 500 | initialized = true; |
---|
| 501 | } |
---|
| 502 | void resizeAndPreserve(int extent1, int extent2, int extent3, int extent4, |
---|
| 503 | int extent5, int extent6, int extent7, int extent8) |
---|
| 504 | { |
---|
| 505 | Array<T_numtype,N_rank>::resizeAndPreserve(extent1, extent2, extent3, extent4, |
---|
| 506 | extent5, extent6, extent7, extent8); |
---|
| 507 | initialized = true; |
---|
| 508 | } |
---|
| 509 | void resizeAndPreserve(int extent1, int extent2, int extent3, int extent4, int extent5, |
---|
| 510 | int extent6, int extent7, int extent8, int extent9) |
---|
| 511 | { |
---|
| 512 | Array<T_numtype,N_rank>::resizeAndPreserve(extent1, extent2, extent3, extent4, extent5, |
---|
| 513 | extent6, extent7, extent8, extent9); |
---|
| 514 | initialized = true; |
---|
| 515 | } |
---|
| 516 | void resizeAndPreserve(int extent1, int extent2, int extent3, int extent4, int extent5, |
---|
| 517 | int extent6, int extent7, int extent8, int extent9, int extent10) |
---|
| 518 | { |
---|
| 519 | Array<T_numtype,N_rank>::resizeAndPreserve(extent1, extent2, extent3, extent4, extent5, |
---|
| 520 | extent6, extent7, extent8, extent9, extent10); |
---|
| 521 | initialized = true; |
---|
| 522 | } |
---|
| 523 | void resizeAndPreserve(int extent1, int extent2, int extent3, int extent4, int extent5, int extent6, |
---|
| 524 | int extent7, int extent8, int extent9, int extent10, int extent11) |
---|
| 525 | { |
---|
| 526 | Array<T_numtype,N_rank>::resizeAndPreserve(extent1, extent2, extent3, extent4, extent5, extent6, |
---|
| 527 | extent7, extent8, extent9, extent10, extent11); |
---|
| 528 | initialized = true; |
---|
| 529 | } |
---|
| 530 | |
---|
| 531 | virtual void fromString(const string& str) { istringstream iss(str); iss >> *this; initialized = true; } |
---|
| 532 | virtual string toString(void) const { ostringstream oss; oss << *this; return oss.str(); } |
---|
| 533 | virtual void reset(void) { this->free(); initialized = false; } |
---|
| 534 | virtual bool isEmpty(void) const { return !initialized; } |
---|
[731] | 535 | virtual size_t size(void) const { return size(this->numElements()); } |
---|
| 536 | static size_t size(sizeType numElements) { return (N_rank + 1) * sizeof(int) + sizeof(size_t) + numElements * sizeof(T_numtype); } |
---|
[677] | 537 | |
---|
| 538 | virtual CBaseType* clone(void) const { return new CArray(*this); } |
---|
| 539 | |
---|
| 540 | virtual bool toBuffer(CBufferOut& buffer) const |
---|
| 541 | { |
---|
| 542 | bool ret; |
---|
| 543 | ret = buffer.put(this->dimensions()); |
---|
| 544 | ret &= buffer.put(this->shape().data(), this->dimensions()); |
---|
| 545 | ret &= buffer.put(this->numElements()); |
---|
| 546 | ret &= buffer.put(this->dataFirst(), this->numElements()); |
---|
| 547 | return ret; |
---|
| 548 | } |
---|
| 549 | |
---|
| 550 | virtual bool fromBuffer(CBufferIn& buffer) |
---|
| 551 | { |
---|
| 552 | bool ret; |
---|
| 553 | int numDim; |
---|
| 554 | TinyVector<int,N_rank> vect; |
---|
| 555 | size_t ne; |
---|
| 556 | |
---|
| 557 | ret = buffer.get(numDim); |
---|
| 558 | ret &= buffer.get(vect.data(), N_rank); |
---|
| 559 | this->resize(vect); |
---|
| 560 | ret &= buffer.get(ne); |
---|
| 561 | ret &= buffer.get(this->dataFirst(), ne); |
---|
| 562 | |
---|
| 563 | initialized = true; |
---|
| 564 | |
---|
| 565 | return ret; |
---|
| 566 | } |
---|
| 567 | }; |
---|
| 568 | |
---|
[1158] | 569 | |
---|
| 570 | #define macro(NRANK)\ |
---|
| 571 | \ |
---|
| 572 | template <>\ |
---|
| 573 | inline size_t CArray<StdString,NRANK>::size(void) const\ |
---|
| 574 | {\ |
---|
| 575 | size_t size=(1 + NRANK) * sizeof(int) ;\ |
---|
| 576 | Array<StdString,NRANK>::const_iterator it, itb=this->begin(), ite=this->end() ;\ |
---|
| 577 | \ |
---|
| 578 | for(it=itb;it!=ite;++it)\ |
---|
| 579 | {\ |
---|
| 580 | size+= sizeof(size_t) ;\ |
---|
| 581 | size+= (*it).size();\ |
---|
| 582 | } \ |
---|
| 583 | return size ;\ |
---|
| 584 | }\ |
---|
| 585 | \ |
---|
| 586 | /* for array string this function is an evaluation of maximum size of the array, considering stringArrayLen is the maximum size of the strings*/ \ |
---|
| 587 | template <>\ |
---|
| 588 | inline size_t CArray<StdString,NRANK>::size(sizeType numElements)\ |
---|
| 589 | {\ |
---|
| 590 | return (NRANK + 1) * sizeof(int) + numElements * stringArrayLen;\ |
---|
| 591 | }\ |
---|
| 592 | \ |
---|
| 593 | template <>\ |
---|
| 594 | inline bool CArray<StdString,NRANK>::toBuffer(CBufferOut& buffer) const\ |
---|
| 595 | {\ |
---|
| 596 | bool ret;\ |
---|
| 597 | ret = buffer.put(this->dimensions());\ |
---|
| 598 | ret &= buffer.put(this->shape().data(), this->dimensions());\ |
---|
| 599 | \ |
---|
| 600 | Array<StdString,NRANK>::const_iterator it, itb=this->begin(), ite=this->end() ;\ |
---|
| 601 | \ |
---|
| 602 | for(it=itb;it!=ite;++it)\ |
---|
| 603 | {\ |
---|
| 604 | ret &= buffer.put((*it).size()) ;\ |
---|
| 605 | ret &= buffer.put((*it).data(),(*it).size());\ |
---|
| 606 | } \ |
---|
| 607 | return ret;\ |
---|
| 608 | }\ |
---|
| 609 | \ |
---|
| 610 | template <>\ |
---|
| 611 | inline bool CArray<StdString,NRANK>::fromBuffer(CBufferIn& buffer)\ |
---|
| 612 | {\ |
---|
| 613 | bool ret;\ |
---|
| 614 | int numDim;\ |
---|
| 615 | TinyVector<int,NRANK> vect;\ |
---|
| 616 | size_t ne;\ |
---|
| 617 | \ |
---|
| 618 | ret = buffer.get(numDim);\ |
---|
| 619 | ret &= buffer.get(vect.data(), NRANK);\ |
---|
| 620 | this->resize(vect);\ |
---|
| 621 | \ |
---|
| 622 | Array<StdString,NRANK>::iterator it, itb=this->begin(), ite=this->end() ;\ |
---|
| 623 | for(it=itb;it!=ite;++it)\ |
---|
| 624 | {\ |
---|
| 625 | ret &= buffer.get(ne) ;\ |
---|
| 626 | char* str = new char[ne] ;\ |
---|
| 627 | ret &= buffer.get(str, ne);\ |
---|
| 628 | *it = string(str,ne) ;\ |
---|
| 629 | delete [] str ;\ |
---|
| 630 | }\ |
---|
| 631 | initialized = true;\ |
---|
| 632 | return ret;\ |
---|
| 633 | } |
---|
| 634 | macro(1) |
---|
| 635 | macro(2) |
---|
| 636 | macro(3) |
---|
| 637 | macro(4) |
---|
| 638 | macro(5) |
---|
| 639 | macro(6) |
---|
| 640 | macro(7) |
---|
| 641 | |
---|
| 642 | #undef macro |
---|
| 643 | |
---|
[369] | 644 | template <typename T_numtype,int N_rank> inline CBufferOut& operator<<(CBufferOut& buffer, const CArray<T_numtype,N_rank>& array) |
---|
| 645 | { |
---|
[677] | 646 | if (!array.toBuffer(buffer)) ERROR("template <typename T_numtype,int N_rank> inline CBufferOut& operator<<(CBufferOut& buffer, const CArray& array)", |
---|
[680] | 647 | << "Not enough free space in buffer to queue the array."); |
---|
[677] | 648 | return buffer; |
---|
[369] | 649 | } |
---|
| 650 | |
---|
| 651 | template <typename T_numtype,int N_rank> inline CBufferIn& operator>>(CBufferIn& buffer, CArray<T_numtype, N_rank>& array) |
---|
| 652 | { |
---|
| 653 | if (!array.fromBuffer(buffer)) ERROR("template <typename T_numtype,int N_rank> inline CBufferIn& operator>>(CBufferIn& buffer, CArray& array)", |
---|
[680] | 654 | << "Not enough data in buffer to unqueue the array."); |
---|
[677] | 655 | return buffer; |
---|
[369] | 656 | } |
---|
[677] | 657 | |
---|
[369] | 658 | template <typename T_numtype,int N_rank> inline CMessage& operator<<(CMessage& msg, const CArray<T_numtype, N_rank>& array) |
---|
| 659 | { |
---|
[677] | 660 | msg.push(array); |
---|
| 661 | return msg; |
---|
[369] | 662 | } |
---|
| 663 | |
---|
| 664 | template <typename T_numtype,int N_rank> inline CMessage& operator<<(CMessage& msg, CArray<T_numtype, N_rank>& array) |
---|
| 665 | { |
---|
[677] | 666 | msg.push(array); |
---|
| 667 | return msg; |
---|
[369] | 668 | } |
---|
| 669 | } |
---|
| 670 | |
---|
| 671 | #endif |
---|