[ander]

Hi

I am going to implement FVM for 2D Euler equations in C++. I wonder what data types I should consider using to store the conserved variables in each node (as well as other variables) One idea is to use c arrays in this fashion:
double U[ni][nj][4].
Another Idea I have is to store them in matrices, in a linalg library like Eigen, but not sure if Eigen supports rank 3 tensors. I have previously written a 2D Euler solver in Matlab using this approach.

Thanks!

[Continuum]

Suggest a class that carries i, j and k. Use k for future 3D if you evolve to that. I would not save in 2D vectors. I would devise a flattening method that works for you structured mesh. You could have a class method that returns the index for the flattened matrix based on ijk input indices.

[requou]

Using the Matrix class from Eigen is the best choice in my opinion.

[ander]

Originally Posted by Continuum

Suggest a class that carries i, j and k. Use k for future 3D if you evolve to that. I would not save in 2D vectors. I would devise a flattening method that works for you structured mesh. You could have a class method that returns the index for the flattened matrix based on ijk input indices.

Thanks. Do you have an example or link to explain how this is done in more detail?