From mboxrd@z Thu Jan 1 00:00:00 1970 X-Spam-Checker-Version: SpamAssassin 3.4.4 (2020-01-24) on polar.synack.me X-Spam-Level: X-Spam-Status: No, score=-1.9 required=5.0 tests=BAYES_00 autolearn=ham autolearn_force=no version=3.4.4 X-Google-Thread: 103376,26e047a7a0709881 X-Google-NewGroupId: yes X-Google-Attributes: gida07f3367d7,domainid0,public,usenet X-Google-Language: ENGLISH,UTF8 Received: by 10.68.219.170 with SMTP id pp10mr20805303pbc.1.1341826555549; Mon, 09 Jul 2012 02:35:55 -0700 (PDT) Path: l9ni11215pbj.0!nntp.google.com!news1.google.com!volia.net!news2.volia.net!feed-A.news.volia.net!news.musoftware.de!wum.musoftware.de!news.swapon.de!eternal-september.org!feeder.eternal-september.org!mx04.eternal-september.org!.POSTED!not-for-mail From: Simon Wright Newsgroups: comp.lang.ada Subject: Re: How to initialize a matrix of one column only? Date: Mon, 09 Jul 2012 10:35:53 +0100 Organization: A noiseless patient Spider Message-ID: References: <4d3d481d-a179-44f6-a81d-04a0531641ad@y3g2000pbc.googlegroups.com> Mime-Version: 1.0 Injection-Info: mx04.eternal-september.org; posting-host="xdR+M82wpvuVwH9XwXKBMw"; logging-data="17344"; mail-complaints-to="abuse@eternal-september.org"; posting-account="U2FsdGVkX1/8wWqzcCazq3Usf5SJY+vDVIQzc7EHUyc=" User-Agent: Gnus/5.13 (Gnus v5.13) Emacs/24.1 (darwin) Cancel-Lock: sha1:weOT6o/ZzzPYcvoU+1BB4kn/ZMw= sha1:z+E8pw/a9ctWW9GIwEEELIII/38= Content-Type: text/plain; charset=utf-8 Content-Transfer-Encoding: 8bit Date: 2012-07-09T10:35:53+01:00 List-Id: Jerry writes: > On Jul 7, 5:08 pm, "Nasser M. Abbasi" wrote: > >> where the labase package simply introduces a type for >> Fortran_Real_Matrix as follows >> >> ---------------------------- >>    type Fortran_Real_Matrix is array (Fortran_Integer range <>, >>                                       Fortran_Integer range <>) >>           of Real; >>      pragma Convention (Fortran, Fortran_Real_Matrix); > > I wish these bindings would be rewritten to use Annex G.3 Real_Vector > and Real_Matrix (etc. for complex) so that a single set of types can > be used in a program that uses the G.3 stuff anyway. The problem with that is that the BLAS is a Fortran library, so expects Fortran conventions for matrices. The Ada_BLAS bindings are thin, so don't do the transposition required . The old GNAT Ada.Numerics.Generic_{Real,Complex}_Arrays did the transposition explicitly. There was a discussion about this (and performance aspects) on this group: http://groups.google.com/group/comp.lang.ada/browse_frm/thread/5395b8f41548bf11/37bcbb133b8d5a62?lnk=gst&q=sauvage#37bcbb133b8d5a62