My question is about implementation alone.Consider a square domain with regular square, cell centred finitevolumes.  This is for the multiscale finite volume method (Jenny and Lunati)I need to  solve the Poisson equation on each cell of the dual mesh, i.e.,the mesh constructed by joining centroids of each cell. (the dual mesh will look like a translation of the original mesh; 4 cells of the dual mesh make up one cell of primary mesh)I will need to use this to construct the transmissibilities for each cell of the primary mesh.My question is:What kind of date structure should I use that will help me go from primary to dual mesh easilyI am confused what I shld store, the coordinates of the cell centres, or the locations of the faces ? How do I store the data that will tell me which 4 dual cells make up a given primary cell ?It seems to me I will need to do a kind of assembly for my transmissibility matrix as I visit each dual cell and solve a local Poisson problem.
data structures for efficient/easy implementation of finite volume method for 2D Poisson equation
finite volume;poisson;mesh;data structures