Vector3D.h

Go to the documentation of this file.

#pragma once 

#include <vector>

template<typename T>
struct Vector3D {
    int xsize = 0;
    int ysize = 0;
    int zsize = 0;
    std::vector<T> value; 
    
    Vector3D() { }
    
    Vector3D(int x, int y, int z) { 
        resize(x,y,z);
    }
    
    Vector3D(int x, int y, int z, T b) { 
        resize(x,y,z);
        fill(b);
    }
    
    
    void fill(T x){
        std::fill(value.begin(), value.end(), x);
    }
    
    void resize(const int x, const int y, const int z) {
        xsize = x; ysize=y; zsize=z;
        value.resize(x*y*z);
    }
    
    void reserve(const int x, const int y, const int z) {
        xsize = x; ysize=y; zsize=z;
        value.reserve(x*y*z);
    }
    
    T& at(const int x, const int y, const int z) {
        return value.at((x*ysize + y)*zsize + z);
    }
        
    T& operator()(const int x, const int y, const int z) {
        // NOTE: this indexing makes iteration over z the fastest
        return value.at((x*ysize + y)*zsize + z);
    }
    
    template<typename X>
    void operator[](X x) {
        assert(false && "**** Cannot use [] with Vector2d, use .at()");
    }
};