From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 26596 invoked by alias); 6 Jan 2010 19:57:28 -0000 Received: (qmail 26587 invoked by uid 22791); 6 Jan 2010 19:57:27 -0000 X-SWARE-Spam-Status: No, hits=-10.2 required=5.0 tests=AWL,BAYES_00,RCVD_IN_DNSWL_HI,SPF_PASS X-Spam-Check-By: sourceware.org Received: from proofpoint1.lanl.gov (HELO proofpoint1.lanl.gov) (204.121.3.25) by sourceware.org (qpsmtpd/0.43rc1) with ESMTP; Wed, 06 Jan 2010 19:57:24 +0000 Received: from mailrelay2.lanl.gov (mailrelay2.lanl.gov [128.165.4.103]) by proofpoint1.lanl.gov (8.14.3/8.14.3) with ESMTP id o06JvMRV015008 for ; Wed, 6 Jan 2010 12:57:22 -0700 Received: from alvie-mail.lanl.gov (alvie-mail.lanl.gov [128.165.4.110]) by mailrelay2.lanl.gov (Postfix) with ESMTP id 35F3F15DB84B for ; Wed, 6 Jan 2010 12:57:22 -0700 (MST) Received: from localhost (localhost.localdomain [127.0.0.1]) by alvie-mail.lanl.gov (Postfix) with ESMTP id 347C27D0053 for ; Wed, 6 Jan 2010 12:57:22 -0700 (MST) X-NIE-2-Virus-Scanner: amavisd-new at alvie-mail.lanl.gov Received: from [130.55.124.157] (manticore.lanl.gov [130.55.124.157]) by alvie-mail.lanl.gov (Postfix) with ESMTP id 1BBEE7D004F for ; Wed, 6 Jan 2010 12:57:22 -0700 (MST) Subject: Re: gsl container designs From: Gerard Jungman To: gsl-discuss@sourceware.org In-Reply-To: <4B447C43.2050408@iki.fi> References: <1259110486.3028.69.camel@manticore.lanl.gov> <4B447C43.2050408@iki.fi> Content-Type: text/plain Date: Wed, 06 Jan 2010 19:57:00 -0000 Message-Id: <1262807849.27244.8.camel@manticore.lanl.gov> Mime-Version: 1.0 Content-Transfer-Encoding: 7bit X-Proofpoint-Virus-Version: vendor=fsecure engine=1.12.8161:2.4.5,1.2.40,4.0.166 definitions=2010-01-06_12:2010-01-05,2010-01-06,2010-01-06 signatures=0 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: 2010-q1/txt/msg00004.txt.bz2 On Wed, 2010-01-06 at 14:04 +0200, Tuomo Keskitalo wrote: > Your design 1 looks very interesting. I don't think your designs would > support banded and packed matrices, or would it? (see the last pages of > "Legacy BLAS: C Interface to the Legacy BLAS" at > http://www.netlib.org/blas/blast-forum/) Thanks. It's been a long time since I looked at that document. You are right; I have only considered the dense case, purely from laziness. However, I think the banded and packed cases can be handled similarly. One needs only to define the indexing schemes to layer on top of the (essentially universal) data structs. Finally, add a syntactic layer, to inject appropriate names, analogous to the syntactic layer that declares vectors, matrices, etc. At least, I think it is not conceptually difficult. Maybe I am missing something. -- G. Jungman