From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 19546 invoked by alias); 18 Dec 2003 17:55:50 -0000 Mailing-List: contact gsl-discuss-help@sources.redhat.com; run by ezmlm Precedence: bulk List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sources.redhat.com Received: (qmail 19529 invoked from network); 18 Dec 2003 17:55:49 -0000 Received: from unknown (HELO scumm) (147.156.163.185) by sources.redhat.com with SMTP; 18 Dec 2003 17:55:49 -0000 Received: from (client is using the wrong hostname!) scumm (really [127.0.0.1]) by scumm via in.smtpd with esmtp id (Debian Smail3.2.0.115) Thu, 18 Dec 2003 18:49:46 +0100 (CET) From: Jordi Burguet Castell To: gsl-discuss@sources.redhat.com Subject: Re: Tensor in GSL. Date: Thu, 18 Dec 2003 17:55:00 -0000 User-Agent: KMail/1.5.4 References: In-Reply-To: MIME-Version: 1.0 Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: 7bit Content-Disposition: inline Message-Id: <200312181849.43823.jordi.burguet-castell@cern.ch> X-SW-Source: 2003-q4/txt/msg00136.txt.bz2 > I'm wondering if tensors (i.e. arrays with more than 2 dimensions) will > be implemented in GSL or if someone is already working on this. I need > high dimensional arrays for my own code, and I was considering to > encode such a data structure. Perhaps it could be included in GSL when > (and if) I get it working. Hi, I did a simple library in C++ to do simple tensor algebra. I also started an attempt to gslize it. You can see my post at the list archives, 2003-05-25. I think it is the one Brian refers to. The C++ one worked fine, but for lack of knowledge about GSL & lack of time & maybe laziness, the gsl version was far from finished. In my old post I described it a bit, but let me recall that it can do arbitrary tensors (any rank, any dimension, as far as there is memory) and simple algebra (tensor product, index contraction,...). Also, a drawback I can think of: there are no strides (maybe easy to add, though). I was counting with this xmas to get it finished or at least usable for GSL, but if anyone gets something before based or not on my code, I will be very happy. Jordi PS: Sorry Andrea, you'll get this twice :)