![]() |
Thanks for the encouragement Paul.
I have taken a quick look at compiling Algol 68 Genie with Visual Studio 2015 and, although I will have to overcome some *nix file permission issues, the initial results are encouraging. I am especially pleased to find that the GSL and pthreads libraries are used since I already build these on Visual Studio. I am not a member of the LinkedIn group but I'll look into joining this. |
IIRC most of the longstanding components - 'longstanding' implies HLL implementation, as opposed to architecture-targeted BLAS - of the Linpack and Eispack libraries were written in Algol ... yep, here is a snip from the top-of-file comment in my local copy of the 1983-dated Eispack comqr.f file, of which I made much use for my PhD work, along with a fair bit of custom code, e.g. a routine for high-accuracy reduction of the generalized linear eigenproblem
[b]Ax[/b] = c[b]Bx[/b] to one of standard form [b]A'x[/b] = c[b]x[/b] with A' of full rank. (The rank deficiency is a near-ubiquitous result of the imposition of boundary conditions in numerical discretization of a continuous differential-operator eigenvalue problem.) The custom-reduction + Eispack-COMQR algorithm is far more accurate than the generalized-eigensystem routines in LAPACK package, or at least was as of the late 1990s, last time I used the code and compared to a free-trial copy of Lapack). But I digress: [code]C THIS SUBROUTINE IS A TRANSLATION OF A UNITARY ANALOGUE OF THE C ALGOL PROCEDURE COMLR, NUM. MATH. 12, 369-376(1968) BY MARTIN C AND WILKINSON. C HANDBOOK FOR AUTO. COMP., VOL.II-LINEAR ALGEBRA, 396-403(1971). C THE UNITARY ANALOGUE SUBSTITUTES THE QR ALGORITHM OF FRANCIS C (COMP. JOUR. 4, 332-345(1962)) FOR THE LR ALGORITHM. C C THIS SUBROUTINE FINDS THE EIGENVALUES OF A COMPLEX C UPPER HESSENBERG MATRIX BY THE QR METHOD.[/code] This is one of those things I wish I had time to play with, but alas my free coding time is fully occupied with code for real-world high-performance applications. But have fun, y'all! |
| All times are UTC. The time now is 11:54. |
Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, Jelsoft Enterprises Ltd.