TweetFollow Us on Twitter

LS FORTRAN 3.3
Volume9
Number11
Column TagJörg's Folder

Language Systems FORTRAN 3.3

Performance test & FORTRAN to C conversion

By Jörg Langowski, MacTech Magazine Regular Contributing Author

Note: Source code files accompanying article are located on MacTech CD-ROM or source code disks.

Since Absoft has just brought out a new Fortran compiler (see July column, MacFORTRAN 3.2), the competition wasn’t sleeping: as I told you in my last column, Language Systems has brought out version 3.3 of their Fortran. It’s now installed (I got it the day I sent off the last column) and here are the Benchmark figures. Note: All benchmark tests were run at optimization level 3 for the LS Fortran compiler and at level -O (i.e. basic optimizations but no loop unrolling and subroutine folding) for the Absoft compiler. First, the results of the Linpack program for LS Fortran 3.3:

Linpack performance, single precision

LSF 3.3 MacIIx 0.125 MFlops

Q700, 68030 code 1.33 MFlops

Q700, 68040 code 1.33 MFlops

These performances are exactly the same as for version 3.0.1. It is interesting that it makes no difference at all whether 68030 or 68040 code is generated; it seems that the Absoft compiler is still making better use of the 68040 specifics. For example, Absoft MacFortran v3.3 runs the same code at 1.30 MFlops on a Quadra 700 when 68030 code is generated (thus at the same speed as LS Fortran) and at 1.61 MFlops for native 68040 code, i.e. 25% faster.

Language Systems has added a faster transcendental function library for the 68040, however. Since the Linpack does not use transcendental functions, this feature is not tested. Another known program, the Savage benchmark, tests accuracy and speed of some built-in math functions (see listing). The results, comparing the latest versions of the LS and Absoft compilers:

Savage Benchmark (execution times and results)

LSF 3.3 MacIIx 3.60 s 2500.000

Q700, 68030 code 3.02 s 2500.000

Q700, 68040 code 1.80 s 2500.000

Absoft MacIIx 3.52 s 2499.886

extended precision 3.65 s 2500.000

Q700, 68030 code 3.23 s 2499.886

Q700, 68040 code 1.37 s 2477.244

Q700, double precision

68040 code 1.43 s 2500.000

One thing that is immediately evident is the shabby performance of the 68040 on 68030 code; the Quadra 700 is no faster here than a Mac IIx! But this comes as no surprise: the 68040 FPU has no transcendental functions and EXP, LOG, etc. must be emulated in software. The 68040 math libraries of both LS and Absoft Fortran are much faster than the Apple emulation. For Absoft on the 68040, you have to take care that you compile with the ‘double precision transcendentals’ option, on the 68040, or with the extended precision option on the 68030, otherwise precision really breaks down. On the 68030, two compilers have identical performance speed- and precision-wise on the Savage benchmark. On the 68040, Absoft again has the speed edge over LS Fortran by about 30%.

So it is fair to say that Absoft’s compiler still produces the fastest-running Fortran code for the 68040; however, when you port programs from other systems, you should be very careful in checking the numerical accuracy against known results. Or play it safe in any case by using default double or extended precision, which doesn’t slow down execution too much.

To save the honor of the Macintosh, one should say that on so-called ‘bigger’ systems you have similar problems: for instance, on a VAX or on a Silicon Graphics workstation you have to declare A double precision to get the correct answer of 2500.000.

It seems at the moment that the Fortran compilers for the Macintosh are at a stage where they have converged to their respective limits of performance, and both are very similar speed- and precision-wise except for Absoft’s higher speed on a 68040. On the other hand, LS Fortran is still the only compiler that creates code to run on systems without an FPU. I guess what one should expect next from Language Systems is that they work on their Macintosh interface, for instance it would be nice to have multiple output windows connected to Fortran units.

Fortran to C

For those of you who have inherited some thousand lines of Fortran code and want to port it to the Mac, there may still be another way to go than using one of those excellent Fortran compilers. Even though Fortran is much better than its reputation with some people (you just have to look at the code provided by Language Systems of Absoft for full Macintosh applications), you might prefer to put in the work of converting your whole program into C, or C++. It is true that you may can gain a lot in readability and maintainability this way; but that mountain of work rewriting the code

Well, if you want to give up your Fortran nostalgia and like the fast compile/test cycle that Think C offers (there is no Think Fortran), here is a program for you: f2c, a public domain Fortran-to-C translator that has been around on Unix systems for quite a while. Some time ago, Basil Duval from Switzerland took the time to convert this program to Think C and put the result into the public domain Info-Mac archives at Stanford. His README file explains what he did:

“IMPORTANT--- This file and the others are available via anonymous ftp on elpp1.epfl.ch.

f2c Conversion on the MAC

Basil P. Duval Jan 1992, EPFL/CRPP Bassenges 1015 Lausanne, Switzerland

(This software is submitted free of change and may not be sold)

Email: Duval@elpp1.epfl.ch

Subject

This note describes a port of f2c to the THINK C Mac environment. The program f2c and the libraries libI77 and libF77 were converted and fortran programs have been successfully translated and linked to form working applications. In this first port, the default THINK C standard "Glass Window" screen is used to perform user input and output.

Version converted

The version taken from the server research.att.edu the 9 jan 1992 and has the 30 december 1991 time stamp.

Environment

The main problem with the port from UNIX to the Mac has been the 32k limit on global variables present up to version 5 of THINK C. Version 4 of THINK C allows string variables to be placed in a separate section. by simply requesting this option in the project options, the global variables were reduced from 56k to 31.6k which is acceptable in the Mac environment. Thus the program could be linked and tested. This enabled other problems of the port to surface. The main problem is due to the confusion with short and long integers. Up to THINK C 4, short was always defined as 2 bytes. In particular this causes several bugs when using the "strn..." functions which expect a (long) integer and as the code arrived were only given a 2 byte (short). The last problem was to provide a simple user interface. THINK C provides a console library which emulates the UNIX style command line via the "ccommand" routine which has been used to provide a rudimentary, but serviceable, interface. A list of the files that were changed is provided and the details of each change are provided at the end of this note.

(details omitted-JL)

Enclosed is a complete translation of the source using THINK C 5.02 together with a test file temp.f and a test project. The application has been compiled and linked for a Mac with 68030 and 68882 FPU unit but with the projects this could easily be changed. I spend the time since you mailed me getting the latest version of the f2c and putting in the bits needed to make it fly on the Mac, and so this version is labelled May 1992.

Good luck. The C source code is not supposed to be easy to edit, but it does get you quite a long way from the Fortran.

Basil Duval 2 June 1992

PS:Language Systems Fortran together with MPW 3.2 is a great combination if you decide to stay in fortran. THINK C is however a great environment for the Mac....”

The ready-to-use f2c application will be on the source code disk; for space reasons we cannot distribute the full system. If you have Internet access, you can however get the files via ftp either from the address given in the README file (see above), or from the Info-Mac archives at sumex-aim.stanford.edu. The conversion of the Savage benchmark is enclosed as an example in Listing 2. Of course you could have done this better by hand - however when you’re dealing with 4000 lines of code, that’s another story. I’ll try my hands on one of the big programs I have here until the next column and tell you my experiences - but first I have to get familiar with Think C, which I have never seriously used, believe it or not. Until then.

Listing 1: Savage Benchmark
C SAVAGE BENCHMARK
 REAL*4 A
C
CSECNDS is a real*4 function that returns 
Cseconds and fractional seconds.  The value is
Cmodified by subtracting the supplied argument.
CIt acts as in the VMS FORTRAN Manual.
C
 real*4 SECNDS
 real*4 time
 write(*,*) 'Starting '
 time = SECNDS(0.0)
 DO 200 J=1,10
 A=1.0D0
 DO 100 I=1,2499
 A=TAN(ATAN(EXP(LOG(SQRT(A*A))))) + 1.0
100CONTINUE
200CONTINUE
 time = SECNDS(time)
 WRITE (*,10) A
10 FORMAT(' A=',F17.12)
 write(*,*) time,' seconds'
 PAUSE
 END
 
 REAL FUNCTION SECNDS(time)
 real time
 n = LONG(362)
 SECNDS = FLOAT(n)/60.0 - time
 END
Listing 2: Savage Benchmark, converted to C by the f2c translator

/*  -- translated by f2c (version of 23 May 1992  14:18:33).
   You must link the resulting object file with the libraries:
 -lF77 -lI77 -lm -lc   (in that order)
*/

#include "f2c.h"

/* Table of constant values */

static integer c__9 = 9;
static integer c__1 = 1;
static real c_b4 = (float)0.;
static integer c__4 = 4;

/* !MP inlines.f */
/* SAVAGE BENCHMARK */
/* Main program */ MAIN__()
{
    /* Format strings */
    static char fmt_10[] = "(\002 A=\002,f17.12)";

    /* Builtin functions */
    integer s_wsle(), do_lio(), e_wsle();
    double sqrt(), log(), exp(), atan(), tan();
    integer s_wsfe(), do_fio(), e_wsfe();

    /* Local variables */
    static real time, a;
    static integer i, j;
    extern doublereal secnds_();

    /* Fortran I/O blocks */
    static cilist io___1 = { 0, 6, 0, 0, 0 };
    static cilist io___6 = { 0, 6, 0, fmt_10, 0 };
    static cilist io___7 = { 0, 6, 0, 0, 0 };


/* SECNDS is a real*4 function that returns */
/* seconds and fractional seconds.  The value is */
/* modified by subtracting the supplied argument. */
/* It acts as in the VMS FORTRAN Manual. */

    s_wsle(&io___1);
    do_lio(&c__9, &c__1, "Starting\311", 9L);
    e_wsle();
    time = secnds_(&c_b4);
    for (j = 1; j <= 10; ++j) {
      a = (float)1.;
      for (i = 1; i <= 2499; ++i) {
        a = tan(atan(exp(log(sqrt(a * a))))) + (float)1.;
/* L100: */
      }
/* L200: */
    }
    time = secnds_(&time);
    s_wsfe(&io___6);
    do_fio(&c__1, (char *)&a, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsle(&io___7);
    do_lio(&c__4, &c__1, (char *)&time, (ftnlen)sizeof(real));
    do_lio(&c__9, &c__1, " seconds", 8L);
    e_wsle();
} /* MAIN__ */

 

Community Search:
MacTech Search:

Software Updates via MacUpdate

FontExplorer X Pro 7.1.3 - Font manageme...
FontExplorer X Pro is optimized for professional use; it's the solution that gives you the power you need to manage all your fonts. Now you can more easily manage, activate and organize your... Read more
DiskCatalogMaker 8.2.5 - Catalog your di...
DiskCatalogMaker is a simple disk management tool which catalogs disks. Simple, light-weight, and fast Finder-like intuitive look and feel Super-fast search algorithm Can compress catalog data for... Read more
Skim 1.5.12 - PDF reader and note-taker...
Skim is a PDF reader and note-taker for OS X. It is designed to help you read and annotate scientific papers in PDF, but is also great for viewing any PDF file. Skim includes many features and has a... Read more
rekordbox 6.1.0.0030 - Professional DJ m...
rekordbox is the best way of preparing and managing your tracks, be it at home, in the studio, or even on the plane! It allows you to import music from other music-management software using the... Read more
iExplorer 4.4.0 - View and transfer file...
iExplorer is an iPhone browser for Mac lets you view the files on your iOS device. By using a drag and drop interface, you can quickly copy files and folders between your Mac and your iPhone or... Read more
OmniGraffle 7.17.5 - Create diagrams, fl...
OmniGraffle helps you draw beautiful diagrams, family trees, flow charts, org charts, layouts, and (mathematically speaking) any other directed or non-directed graphs. We've had people use Graffle to... Read more
Apple Configurator 2.13.1 - Configure an...
Apple Configurator makes it easy to deploy iPad, iPhone, iPod touch, and Apple TV devices in your school or business. Use Apple Configurator to quickly configure large numbers of devices connected to... Read more
OmniGraffle Pro 7.17.5 - Create diagrams...
OmniGraffle Pro helps you draw beautiful diagrams, family trees, flow charts, org charts, layouts, and (mathematically speaking) any other directed or non-directed graphs. We've had people use... Read more
FoldersSynchronizer 5.1.6 - Synchronize...
FoldersSynchronizer is a popular and useful utility that synchronizes and backs-up files, folders, disks and boot disks. On each session you can apply special options like Timers, Multiple Folders,... Read more
Adobe After Effects 17.1.4 - Create prof...
After Effects is available as part of Adobe Creative Cloud for $52.99/month (or $20.99/month for a single app license). The new, more connected After Effects can make the impossible possible. Get... Read more

Latest Forum Discussions

See All

Infinity Mechs is an upcoming idle game...
Indie developer SkullStar studio has announced an upcoming idle mech game called Infinity Mechs. It draws inspiration from the mobile game Iron Saga and has been officially licensed by Game Duchy. It's set to launch for both iOS and Android on... | Read more »
PUBG Mobile Lite's latest update se...
PUBG Mobile Lite, the streamlined version of the popular battle royale that's designed to work on less powerful devices, sees the return of a popular game variant today, Survive Till Dawn mode. It arrives as part of the 0.19.0 content update. [... | Read more »
Matchy Catch, Jyamma Games’ new hyper-ca...
Matchy Catch is a new hyper-casual puzzler from Jyamma Games, the Italian studio behind the Pong-inspired puzzle-adventure Hi-Ball Rush. It’s only the developer’s second game for iOS and Android devices, but it promises to be every bit as fun and... | Read more »
Among Us! Imposter Guide - How to be a s...
Among Us! continues to be getting a lot of play in these parts, and since our first guide we've learned a thing or two about the game. This is especially true regarding the imposter role, as its a relatively rare opportunity that we've now put... | Read more »
Paladin's Story is an upcoming fant...
Paladin's Story is an upcoming fantasy RPG with an off-kilter sense of humour that's heading for iOS and Android. It will officially launch for both on September 16th though the game is already available on Google Play in Early Access. [Read more... | Read more »
Among Us! Guide - Tips for the uninitiat...
| Read more »
A Pretty Odd Bunny is a stealth-platform...
A Pretty Odd Bunny is a stealth-platformer from two-man team AJ Ordaz and René Rivera. It follows the story of a red-eyed rabbit who is allergic to carrots and instead has a penchant for devouring pigs. It's available now for Android devices. [... | Read more »
Apple Arcade: Ranked - Top 25 [Updated 9...
In case you missed it, I am on a quest to rank every Apple Arcade game there is. [Read more] | Read more »
The 5 Best Mobile Games Like Tony Hawk...
Tony Hawk’s Pro Skater 1+2 dropped last week, meaning you can get remastered versions of two of the most iconic “sports” games ever made and experience some of the highest watermarks set in games of generations’ past. [Read more] | Read more »
Bouncing Box is a challenging platformer...
In platformers, we often spend a large amount of time mindless destroying boxes by jumping onto or into them, whatever makes them splinter apart. Bouncing Box from developer The K Brothers aims to give those poor destructible cubes a chance to be... | Read more »

Price Scanner via MacPrices.net

Price drop! Get a 44mm Apple Watch Series 5 G...
Amazon has dropped their price on the 44mm Apple Watch Series 5 GPS + Cellular by $100 to $429 shipped. That’s $100 off Apple’s original MSRP for this model. For the latest prices and sales, see our... Read more
Verizon offers $200 discount on new Apple Wat...
Verizon will take up to $200 off the purchase of a new GPS + Cellular Apple Watch Series 6 or Apple Watch SE with select trade-in and the purchase of a new iPhone with service. The fine print: “Get... Read more
Verizon offers $250 discount on new 8th gener...
Verizon will take up to $250 off the price of an 8th generation 2020 Apple Cellular iPad with select trade-ins and a new iPhone purchase. Plus get Apple News+ free for 6 months. The fine print: “Save... Read more
Apple’s Implementation Of COVID-19 Exposure...
NEWS: 09.18.20 – The latest effort by Apple to embed exposure notifications for COVID-19 contact tracing right into its mobile operating system has some iPhone users weary of being exposed to... Read more
Here’s how to get a 16″ MacBook Pro for $300...
B&H Photo has new 16″ MacBook Pros on sale today for $250-$300 off Apple’s MSRP, starting at $2099. Expedited shipping is free to many addresses in the US: – 2019 16″ 2.6GHz 6-Core MacBook Pro... Read more
Apple has Certified Refurbished 16″ MacBook P...
Apple has Certified Refurbished 2019 16″ MacBook Pros available for up to $420 off the cost of new models, starting at $2039. Each model features a new outer case, shipping is free, and an Apple 1-... Read more
Price drops! Apple reseller B&H drops App...
B&H Photo has dropped prices on Apple Watch Series 5 models by $50-$70 off Apple’s original MSRP. Shipping is free. These are the same Apple Watch models sold by Apple in their retail and online... Read more
Apple extends their 2020 Back-to-School promo...
As part of their Back to School promotion for 2020, Apple will include one free pair Apple AirPods (with charging case) with the purchase of a MacBook Air, MacBook Pro, iMac, iMac Pro (Mac Pro and... Read more
Apple 7.9″ iPad minis are on sale today for $...
Amazon has new 7.9″ 64GB WiFi iPad minis on sale today for $50 off Apple’s MSRP, each including free shipping. Prices start at $349. These are the same iPad minis sold by Apple in their retail and... Read more
Lowest price anywhere! New 16″ 6-core MacBook...
Expercom has the Silver 16″ 6-core MacBook Pro back on sale for a limited time for $2079 shipped. Their price is $320 off Apple’s MSRP for this model, and it’s the cheapest price currently available... Read more

Jobs Board

Security Officer ($23.00/Hourly) - *Apple*...
**Security Officer \($23\.00/Hourly\) \- Apple Store** **Description** About NMS Built on a culture of safety and integrity, NMSdelivers award\-winning, integrated Read more
Platform - Workplace Eng - *Apple* Enterpri...
MORE ABOUT THIS JOB We are looking for an Apple Platform Engineer who will bring a unique engineering skill set, support, clarity, organization and above all else, Read more
*Apple* Certified Repair Technician - Utah S...
…selected candidate will work in the USU Campus Store Tech Department as an Apple Certified Repair Technician and floor associate. This position is for both summer Read more
Senior Data Engineer - *Apple* - Theorem, L...
Job Summary Apple is seeking an experienced, detail-minded data engineeringconsultant to join our worldwide business development and strategy team. If you are Read more
Desktop Support Technician - A+, *Apple* -...
Desktop Support Technician - A+, Apple **Ref No.:** 20-01604 **Location:** Miami, Florida **e** at http://www.excell.com/ **X** at http://www.excell.com/ **cell** at Read more
All contents are Copyright 1984-2011 by Xplain Corporation. All rights reserved. Theme designed by Icreon.