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,FREEMAIL_FROM autolearn=unavailable autolearn_force=no version=3.4.4 X-Received: by 10.159.48.153 with SMTP id j25mr10731486uab.43.1494336038263; Tue, 09 May 2017 06:20:38 -0700 (PDT) X-Received: by 10.157.35.104 with SMTP id k37mr1378567otd.14.1494336038219; Tue, 09 May 2017 06:20:38 -0700 (PDT) Path: eternal-september.org!reader01.eternal-september.org!reader02.eternal-september.org!news.eternal-september.org!news.eternal-september.org!feeder.eternal-september.org!2.eu.feeder.erje.net!feeder.erje.net!2.us.feeder.erje.net!weretis.net!feeder6.news.weretis.net!feeder.usenetexpress.com!feeder1.iad1.usenetexpress.com!border1.nntp.dca1.giganews.com!border2.nntp.dca1.giganews.com!nntp.giganews.com!l39no760100qtb.0!news-out.google.com!m134ni535itb.0!nntp.google.com!c26no1784316itd.0!postnews.google.com!glegroupsg2000goo.googlegroups.com!not-for-mail Newsgroups: comp.lang.ada Date: Tue, 9 May 2017 06:20:37 -0700 (PDT) In-Reply-To: Complaints-To: groups-abuse@google.com Injection-Info: glegroupsg2000goo.googlegroups.com; posting-host=185.30.133.97; posting-account=hya6vwoAAADTA0O27Aq3u6Su3lQKpSMz NNTP-Posting-Host: 185.30.133.97 References: User-Agent: G2/1.0 MIME-Version: 1.0 Message-ID: <7543686f-136c-44a6-bf50-64bef4cdc66d@googlegroups.com> Subject: Re: ANN:Ada 2005 Math Extensions 20170427 From: Vincent Injection-Date: Tue, 09 May 2017 13:20:38 +0000 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: quoted-printable Xref: news.eternal-september.org comp.lang.ada:46731 Date: 2017-05-09T06:20:37-07:00 List-Id: Le mardi 9 mai 2017 14:49:03 UTC+2, Simon Wright a =C3=A9crit=C2=A0: > > - AdaLAPACK : https://sourceforge.net/projects/adalapack/ >=20 > This is a binding, and its interface is explicitly in terms of > Fortran-convention objects; for example, >=20 > type Fortran_Real_Matrix is > array (Fortran_Integer range <>, Fortran_Integer range <>) of Real; > pragma Convention (Fortran, Fortran_Real_Matrix); >=20 > ... Thank you very much Simon for your kind answer. I will use the third one, t= he Fortran binding, because I prefer column-major order, for different reas= ons : 1. It is easier to interface with Fortran solvers (obviously), and many exc= ellent solvers are still in Fortran. 2. It is easier to interface with Matlab compatible solvers written in C th= at use sparse matrices in CSC format, since CSC is column-major. 3. It is easier to deal with multidimensionnal arrays that contain matrices= (a very common case in scientific computation). 4. It is consistent with little endian ordering. 5. It is the convention in math to see matrices composed of column vectors.= Using matrices in row-major order forces you to transpose formulas, and it= is a pain. 6. A simple Pragma Convention (Fortran) is enough to transpose to row-major= an array manually entered in the source code. Kind regards, Vincent