Pauli Virtanen (el 2007-11-23 a les 01:18:18 +0000) va dir:: >[...] > Let me first clarify the reasons why I chose to describe the data as an > array of dimensions > > (2, 2, N, N, P, Q) > > In the problem I solved, the quantity of interest is a 2 x 2 tensor which > is a function defined on a N x N grid of spatial coordinates. This needs > to be solved for different control parameters (say, p and q), and for > each combination of them I get the solution at once for all grid points > --- and in Fortran order. >[...]
Dropping by the very interesting discussion, if 2x2 tensors is the basic
element you want to store (varying with a series of paramenters), and
since it doesn't make much sense for the chunkshape to split a tensor in
the middle, you could have replaced::
f.createCArray(f.root, 'huge_array', tables.Float64Atom(), (2,2,N,N,50,50))
by::
f.createCArray(..., tables.Float64Atom(shape=(2,2)), (N,N,50,50))
This would have had the side-effect of placing the largest dimensions
first and using a C order. It would be a pity, since the disussion may
not have been spurred! ;)
::
Ivan Vilata i Balaguer >qo< http://www.carabos.com/
Cárabos Coop. V. V V Enjoy Data
""
signature.asc
Description: Digital signature
------------------------------------------------------------------------- This SF.net email is sponsored by: Microsoft Defy all challenges. Microsoft(R) Visual Studio 2005. http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/
_______________________________________________ Pytables-users mailing list Pytables-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/pytables-users
