C++ 中 std::vector 的内存问题

我在 c++ 中使用 std：： vector 时遇到了内存问题。这是我的代码：

#include <iostream>
#include <vector>
int main () {
std::vector< std::vector<float> > mesh_points_A;
int N=10;
for(int i=0;i<N;i++){ 
for(int j=0;j<N;j++){ 
std::vector<float> xyz;
xyz.push_back(i);
xyz.push_back(j);
xyz.push_back(0.3);
mesh_points_A.push_back(xyz);
}
}
return 0;
}

当我将 N 增加到 10000 或更高时，我的内存不足......但我认为我做错了什么，因为例如，如果我将 python 与 numpy 数组一起使用，这很容易......

提前非常感谢。

编辑：这是原始代码。上面的代码只是一个简化，以更好地说明问题。问题是，是否可以以某种方式创建许多 Surface 对象(在代码中当前有两个(，而不会耗尽内存，同时保持 N=10000。

// classes example compile with c++ -o Surface Surface.cpp -std=c++11
#include <iostream>
#include <vector>
#include <array>
class Surface {
private:
std::vector< std::array<float,3> > mesh_points_A; 
public:
float R;
float z; // z position if the suface
int n_A; //number of area mesh points mesh_points_A.size()
Surface(int nxA, float R , float z); 
};
Surface::Surface(int nxA, float  R,float z)
: z(z)
, R(R)
{
float dxA= 2*R/(nxA*1.0-1.0);
//determine n_A, 
n_A=0;
for(int i=0;i<nxA;i++){ 
float x = -R+i*dxA;
for(int j=0;j<nxA;j++){ 
float y = -R+j*dxA;
if(x*x+y*y<R*R){
n_A+=1;
}
}
}
std::cout<<"Number of area mesh points: "<<n_A<<std::endl;
mesh_points_A.reserve(n_A);
for(int i=0;i<nxA;i++){ 
float x = -R+i*dxA;
for(int j=0;j<nxA;j++){ 
float y = -R+j*dxA;
if(x*x+y*y<R*R){
std::array<float,3> xyz{ {x,y,z} };
mesh_points_A.push_back(xyz);       
}
}
}

}

int main () {
int N= 20000;
Surface s1(N,0.1,0.0);
Surface s2(N,0.1,0.1);
return 0;
}