#include <iostream>
#include <vector>
#include <blitz/array.h>
//  ---------------------------------------------------------------------------
using namespace std;
using namespace blitz;
//  ---------------------------------------------------------------------------
typedef Array<double,1>  Dvector;
//  ---------------------------------------------------------------------------
#define MAGIC  3
//  ---------------------------------------------------------------------------
class MyClass {
public:
    MyClass() {
        vpDV.resize(MAGIC);
    }
    ~MyClass() {
        for (size_t i=0;i<vpDV.size();i++)
            vpDV[i]=0;
        vpDV.clear();
    }
    void mydisplay() {
        for (size_t i=0;i<vpDV.size();i++) {
            for (int k=0;k<vpDV[i]->extent(firstDim);k++) {
                cout<<(*vpDV[i])(k)<<' ';
            }
            cout<<endl;
        }
    }
public:
    vector<Dvector*> vpDV;
};
//  ---------------------------------------------------------------------------
int main() {
    MyClass obj;
    Dvector DV0(3), DV1(2), DV2(4);
//
    obj.vpDV[0] = &DV0;
    obj.vpDV[1] = &DV1;
    obj.vpDV[2] = &DV2;
//
    obj.mydisplay();
//
    obj.(*vpDV[0])(1) = 0.1;
    // obj.vpDV[0]->(1) = 0.1;
//  The operator-> isn't supported as i expected normally!!!
//
    obj.mydisplay();
//
    return 0;
}
//  ---------------------------------------------------------------------------
#undef MAGIC
//  ---------------------------------------------------------------------------