From mboxrd@z Thu Jan  1 00:00:00 1970
Return-Path: <gsl-discuss-return-2533-listarch-gsl-discuss=sourceware.cygnus.com@sources.redhat.com>
Received: (qmail 30777 invoked by alias); 16 Jul 2002 13:25:37 -0000
Mailing-List: contact gsl-discuss-help@sources.redhat.com; run by ezmlm
Precedence: bulk
List-Subscribe: <mailto:gsl-discuss-subscribe@sources.redhat.com>
List-Archive: <http://sources.redhat.com/ml/gsl-discuss/>
List-Post: <mailto:gsl-discuss@sources.redhat.com>
List-Help: <mailto:gsl-discuss-help@sources.redhat.com>, <http://sources.redhat.com/ml/#faqs>
Sender: gsl-discuss-owner@sources.redhat.com
Received: (qmail 30766 invoked from network); 16 Jul 2002 13:25:36 -0000
Received: from unknown (HELO numenor.qualcomm.com) (129.46.51.58)
  by sources.redhat.com with SMTP; 16 Jul 2002 13:25:36 -0000
Received: from sabrina.qualcomm.com (sabrina.qualcomm.com [129.46.61.150])
	by numenor.qualcomm.com (8.12.3/8.12.3/1.0) with ESMTP id g6GDPHgD022174
	(version=TLSv1/SSLv3 cipher=EDH-RSA-DES-CBC3-SHA bits=168 verify=NO);
	Tue, 16 Jul 2002 06:25:17 -0700 (PDT)
Received: from JOHNK1.qualcomm.com (johnk1.qualcomm.com [129.46.32.157])
	by sabrina.qualcomm.com (8.12.3/8.12.3/1.0) with ESMTP id g6GDPDK3011964;
	Tue, 16 Jul 2002 06:25:14 -0700 (PDT)
Message-Id: <5.1.0.14.0.20020716091405.01e18b40@mail1.qualcomm.com>
X-Sender: johnk@mail1.qualcomm.com
Date: Tue, 16 Jul 2002 06:25:00 -0000
To: sliwa@euv-frankfurt-o.de, gsl-discuss@sources.redhat.com
From: John Ketchum <johnk@qualcomm.com>
Subject: Re: Arrays in higher dimensions
In-Reply-To: <3D33E6DD.8060004@euv-frankfurt-o.de>
Mime-Version: 1.0
Content-Type: text/plain; charset="us-ascii"; format=flowed
X-SW-Source: 2002-q3/txt/msg00068.txt.bz2

For 4-d array, with dimensions n1, n2, n3, and n4, something like the 
following:

gsl_matrix ***m4;
m4 = (gsl_matrix***)malloc(n4*sizeof(gsl_matrix**));
for (i=0; i<n4; i++) {
         m4[i] = (gsl_matrix**)malloc(n3*sizeof(gsl_matrix*));
         for (j=0; j<n3; j++) m4[i][j]=gsl_matrix_alloc(n1,n2);
}


This gives you a matrix of pointers to gsl matrices,  which is a 4-d construct.

Make sure you free all of the memory when you are done:

for (i=0; i<n4; i++) {
         for (j=0; j<n3; j++) gsl_matrix_free(m4[i][j]);
         free(m4[i]);
}
free(m4);



At 11:26 AM 7/16/2002 +0200, Przemyslaw Sliwa wrote:
>Hi,
>my problem is - is it possible to create an ayray of gsl_matrices and 
>gsl_vectors within the gsl library?
>I would like to create 3 and 4 dimansional arrays - that means I need one 
>and two dimensional array of
>gsl_matrices and gsl_vectors.
>
>If it is possible - how can I access the elements of such an array.
>Thanks to everyone.
>
>
>Przemyslaw Sliwa
>
>
>Witout stupid mementos.
>

John Ketchum
Qualcomm Inc.
9 Damonmill Square Suite 2A
Concord, MA 01742
Office: 60915
Direct: 781-276-0915
Fax: 781-276-0901
Mobile: 617-312-7023