From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 588 invoked by alias); 28 Aug 2009 13:57:47 -0000 Received: (qmail 21198 invoked by uid 22791); 27 Aug 2009 14:37:00 -0000 X-SWARE-Spam-Status: No, hits=-2.0 required=5.0 tests=BAYES_00,SARE_MSGID_LONG40,SPF_PASS X-Spam-Check-By: sourceware.org MIME-Version: 1.0 In-Reply-To: References: <48E25CA9.6080306@iki.fi> <498727E5.6080407@iki.fi> <49AA9DB5.6030908@iki.fi> <49FB01D1.30000@iki.fi> <4A7ADFDC.9080408@iki.fi> <4A967114.8080600@iki.fi> From: Jordi Burguet Castell Date: Fri, 28 Aug 2009 13:57:00 -0000 Message-ID: <1d6905490908270736s626985b4ycef9a9e53e4a4bb5@mail.gmail.com> Subject: Re: GSL 2.0 roadmap To: "Robert G. Brown" Cc: GSL Discuss Mailing List Content-Type: text/plain; charset=ISO-8859-1 Content-Transfer-Encoding: quoted-printable Mailing-List: contact gsl-discuss-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: gsl-discuss-owner@sourceware.org X-SW-Source: 2009-q3/txt/msg00024.txt.bz2 Hello Robert, On Thu, Aug 27, 2009 at 7:51 AM, Robert G. Brown wrote: > On the same note, a rather long time ago I developed a set of GSL-ish > "tensor" structures to enable the GSL to handle things beyond matrices > (which involved of course redefining slightly 1st and 2nd rank tensors, > or at least paralleling the matrix and the vector). =A0Tensors are, > obviously useful to all sorts of people doing computations in physics; > not only spatial tensors in e.g. relativistic theory computations but > also angular moment tensors that may not be efficiently rectangular. > > [...] > > I haven't worked with the first two for a while now, but I've still got > the code (of course) and would be happy to de-mothball it if there is > any interest in reconsidering the vector/matrix/tensor interface so it > is more extensible. =A0I think I implemented all the way out to 10th rank > tensors -- possibly overkill even for relativity, but I myself use at > least sixth rank in some of my code, and I'm aware of people who need > eighth rank. You may want to check the extensions "tensor" and "marray" that I also developed a while ago: http://www.gnu.org/software/gsl/#extensions I don't know how much overlap they have with your ideas, the main thing for them is to manipulate tensors of arbitrary rank and do the common operations on them. This may also interest you, quoting from the web: "If you want to add a feature to GSL we recommend that you make it an extension first. We will list it here so that people can try it out. Extensions can be incorporated after they have been tested in real use (see "How to help" for more information)." So maybe you can consider making your code an extension if it is not already? Cheers, Jordi