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=-0.3 required=5.0 tests=BAYES_00, REPLYTO_WITHOUT_TO_CC autolearn=no autolearn_force=no version=3.4.4 X-Google-Thread: 103376,39ace4c4032cf84b X-Google-NewGroupId: yes X-Google-Attributes: gida07f3367d7,domainid0,public,usenet X-Google-Language: ENGLISH,ASCII-7-bit Path: g2news2.google.com!news4.google.com!feeder.news-service.com!feeds.phibee-telecom.net!de-l.enfer-du-nord.net!feeder1.enfer-du-nord.net!feeder.news-service.com!94.75.214.39.MISMATCH!aioe.org!not-for-mail From: "Nasser M. Abbasi" Newsgroups: comp.lang.ada Subject: Ada Annex G (was Re: question on GtkAda) Date: Sun, 19 Dec 2010 16:50:59 -0800 Organization: Aioe.org NNTP Server Message-ID: References: <4e69ab19-1d1e-4ab9-935b-a5da160dadc4@e16g2000pri.googlegroups.com> Reply-To: nma@12000.org NNTP-Posting-Host: tUYQ4Ty9mMw9Pdc8TJRFQA.user.speranza.aioe.org Mime-Version: 1.0 Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 7bit X-Complaints-To: abuse@aioe.org X-Notice: Filtered by postfilter v. 0.8.2 User-Agent: Mozilla/5.0 (Windows; U; Windows NT 6.1; en-US; rv:1.9.2.13) Gecko/20101207 Thunderbird/3.1.7 Xref: g2news2.google.com comp.lang.ada:17023 Date: 2010-12-19T16:50:59-08:00 List-Id: On 12/18/2010 9:57 AM, Simon Wright wrote: > "Nasser M. Abbasi" writes: > >> Now, given all the above, for scientific applications, I still >> need to be able to access lapack functions from Ada. Ada Annex G >> just does not cut it. Very limited. As well as what other small >> Ada math libraries I could find. >> >> >> So far, I could not find full binding to these. There are 100's of >> functions there, that has been written over span of more than 30 >> years or more: >> >> http://www.netlib.org/lapack/release_notes.html >> >> that one would need to call to do any type of scientific >> work. > > As a matter of interest, what main functionality are you missing? > (bearing in mind that many functions come in 4 variants (s, d, c, z) and > that there are many that are optimized for a particular kind of matrix, > eg tridiagonal). Ok, here is what I found. But I preface this by saying that I am no expert in Ada, nor in linear algebra for that matter, but what I did is look at the lapack functions, and looked at the AnnexG, and determined if that lapack function can be implemented using only functions in the Ada 2005 annexG. I only looked at the real (single) functions listed here http://www.netlib.org/lapack/single/ and went only to the first 58 functions (up to where it says "Available COMPUTATIONAL routines" ) I made a list of the result here (including Matlab and Mathematica functions). http://12000.org/my_notes/lapack_analysis/lapack/index.htm Please correct me if you see any errors, and I will fix the note above. Conclusion on Ada: 1) no support complex types Ax=b, only real, but this does not seem too important. Not sure now. 2) No condition number function 3) NO solver for Constrained Linear Least Squares 4) No Generalized Linear Regression Model 5) No Schur Decomposition 6) No Singular Value Decomposition (SVD) 7) No generalized Singular Value Decomposition 8) No support for finding Eigenvalues/Eigenvectors for the generalized case. It is very well possible that some of the missing functionality above is present is some other Ada libraries over the net, but I was only looking at the annex G. I also, did not try to code any Ada examples to verify any of the above, I was just looking at the manual. It seems to me it would have been more efficient to provide an "offical" bindings to the most stable known math libraries out there, than to try to provide this in an Ada standard library such as the annex? May be a thich binding? But this is not an easy problem. Whcih libraries? which versions? who will make such an "official" binding, etc... --Nasser