Re: [Matlisp-users] help with gross array hack needed

[Raymond T. <ray...@er...>]

>>>>> "rif" == rif  <ri...@MI...> writes:

    rif> Is there any way, in CMUCL, using some implementation-specific
    rif> construct, to create an "aliased array"?  Suppose I create a
    rif> double-float array of size 10:

    rif> (make-array 10 :element-type 'double-float)

    rif> I know I can take it's address, with #'system:vector-sap.

    rif> What I want to do is manipulate that address and somehow "invert"
    rif> vector-sap to create a new double-float array whose first element is
    rif> (for example) the third element of the original array.

    rif> I need this because I need to pass matlisp a subset of an array.  More
    rif> specifically, I want to allocate an n by n+1 array, so that I can
    rif> store half of each of two n by n symmetric arrays in n(n+1) space
    rif> rather than 2n^2 space.  This is a common technique for using
    rif> BLAS/LAPACK, but to get at it, I need to be able to give a Fortran
    rif> matlisp routine a way to get at the square n by n matrix that "starts
    rif> with the second column" of the n by n+1 matrix, and the way the
    rif> interface is defined, it wants actual arrays, not pointers.  If I can
    rif> somehow create the appropriate array, I don't have to go mucking about
    rif> in the internals of matlisp's def-fortran-routine and so forth.

Marc has already mentioned displaced arrays.  This should work, but
depends on what you are really trying to do.  What matlisp function
are you trying to call?

Ray