home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #31 / NN_1992_31.iso / spool / comp / lang / c / 18982 < prev    next >
Encoding:
Internet Message Format  |  1992-12-29  |  2.0 KB
  1. Path: sparky!uunet!olivea!spool.mu.edu!sgiblab!news.kpc.com!kpc!hollasch
  2. From: hollasch@kpc.com (Steve Hollasch)
  3. Newsgroups: comp.lang.c
  4. Subject: Re: Accessing a dynamical 3-D field
  5. Summary: Multi-dimensional array element access
  6. Message-ID: <1992Dec29.190722.20859@kpc.com>
  7. Date: 29 Dec 92 19:07:22 GMT
  8. References: <1992Dec29.142134.25641@Informatik.TU-Muenchen.DE>
  9. Sender: usenet@kpc.com
  10. Organization: Kubota Pacific Computer, Inc.
  11. Lines: 32
  12.  
  13. backscha@Informatik.TU-Muenchen.DE (Martin Backschat) writes:
  14. | 
  15. | In my program I use 3-D fields of doubles with variable bounds.
  16. | This field is implemented as an array of doubles which is allocated
  17. | dynamically with malloc.
  18. | To access for example the [i][j][k] component in the field f of dimension
  19. | U,V,W I have to write   f[i*V*W + j*W + k].
  20. | 
  21. | Is it somehow possible to declare and allocate a 3-D field u (of doubles) with
  22. | variable bounds, allocated by malloc, in such a way that it can be accessed
  23. | with u[i,j,k] (with type casting)?
  24.  
  25.     Quite possibly, but to what advantage?  When the compiler performs a
  26. mutli-dimensional array access, it does this exact computation to yield the
  27. address of the element in question.  There are _some_ CPU's that have
  28. instructions to handle (non-general) 2D arrays, but I don't think there are
  29. many (if any) CPU's that have a single instruction to do a 3D array access.
  30.  
  31.     So, if you WERE able to typecast this to a multi-dimensional array (I
  32. tried to find a way but gave up), the compiler would still generate this
  33. code.  What I would do is to just write a macro, e.g.:
  34.  
  35.         #define FIELD(array,i,j,k) \
  36.             array [((((i) * V) + (j)) * W) + (k)]
  37.  
  38.     So then you could just do "FIELD(f,i,j,k)" to refer to a given element.
  39. Note the way that I computed the index takes one less multiply than the
  40. formula you listed.
  41.  
  42. ______________________________________________________________________________
  43. Steve Hollasch                                   Kubota Pacific Computer, Inc.
  44. hollasch@kpc.com                                 Santa Clara, California
  45.