Hi, Tried this, and it doesn't work. Any way around it? Thanks. vector <vector <vector <vector<double> > > > v0(N, vector<double> (N, vector<double> (N, vector<double> (N, 0.0))));