mersenneforum.org  

Go Back   mersenneforum.org > Extra Stuff > Programming > ALGOL 68

Reply
 
Thread Tools
Old 2015-07-06, 20:20   #12
Brian Gladman
 
Brian Gladman's Avatar
 
May 2008
Worcester, United Kingdom

523 Posts
Default

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.
Brian Gladman is offline   Reply With Quote
Old 2015-07-07, 01:03   #13
ewmayer
2ω=0
 
ewmayer's Avatar
 
Sep 2002
Rep├║blica de California

2B8616 Posts
Default

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
Ax = cBx to one of standard form A'x = cx 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.
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!
ewmayer is offline   Reply With Quote
Reply

Thread Tools


All times are UTC. The time now is 06:58.

Thu May 28 06:58:16 UTC 2020 up 64 days, 4:31, 1 user, load averages: 1.67, 1.80, 1.72

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2020, Jelsoft Enterprises Ltd.

This forum has received and complied with 0 (zero) government requests for information.

Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation.
A copy of the license is included in the FAQ.