Jump to content

Fortran

From Wikipedia, the free encyclopedia

Fortran
ParadigmMulti-paradigm:structured,imperative(procedural,object-oriented),generic,array
Designed byJohn Backus
DeveloperJohn BackusandIBM
First appeared1957;67 years ago(1957)
Stable release
Fortran 2023 (ISO/IEC 1539:2023) / November 17, 2023;10 months ago(2023-11-17)
Typing disciplinestrong,static,manifest
Filename extensions.f90,.f,.for
Websitefortran-lang.org
Majorimplementations
Absoft,Cray,GFortran,G95,IBMXL Fortran,Intel,Hitachi,Lahey/Fujitsu,Numerical Algorithms Group,Open Watcom,PathScale,PGI,Silverfrost,Oracle Solaris Studio,others
Influenced by
Speedcoding
Influenced
ALGOL 58,BASIC,C,Chapel,[1]CMS-2,DOPE,Fortress,MATLAB,PL/I,PACT I,MUMPS,IDL,Ratfor,SAKO[2]

Fortran(/ˈfɔːrtræn/;formerlyFORTRAN) is athird generation,compiled,imperativeprogramming languagethat is especially suited tonumeric computationandscientific computing.

Fortran was originally developed byIBM.[3]It first compiled correctly in 1958.[4]Fortrancomputer programshave been written to support scientific and engineering applications, such asnumerical weather prediction,finite element analysis,computational fluid dynamics,geophysics,computational physics,crystallographyandcomputational chemistry.It is a popular language forhigh-performance computing[5]and is used for programs that benchmark and rank the world'sfastest supercomputers.[6][7]

TheIBM Blue Gene/Psupercomputer installation in 2007 at the Argonne Leadership Angela Yang Computing Facility located in theArgonne National Laboratory,inLemont,Illinois,US

Fortran has evolved through numerous versions and dialects. In 1966, theAmerican National Standards Institute(ANSI) developed a standard for Fortran to limit proliferation of compilers using slightly different syntax.[8]Successive versions have added support for a character data type (Fortran 77),structured programming,array programming,modular programming,generic programming(Fortran 90),parallel computing(Fortran 95),object-oriented programming(Fortran 2003), andconcurrent programming(Fortran 2008).

FORTRAN and COBOL genealogy tree

Since April 2024, Fortran has ranked among the top ten languages in theTIOBE index,a measure of the popularity of programming languages.[9]

Naming

[edit]

The first manual for FORTRAN describes it as aFormula Translating System,and printed the name withsmall caps,Fortran.[10]: p.2 [11]Other sources suggest the name stands forFormula Translator,[12]orFormula Translation.[13]

Early IBM computers did not supportlowercaseletters, and the names of versions of the language through FORTRAN 77 were usually spelled in all-uppercase.[14]FORTRAN 77 was the last version in which the Fortran character set included only uppercase letters.[15]

The official languagestandardsfor Fortran have referred to the language as "Fortran" with initial caps since Fortran 90.[citation needed]

Origins

[edit]
John Backus (1924–2007) proposed the FORTRAN project in December 1953 and received the A.M. Turing Award in 1977.

In late 1953,John W. Backussubmitted a proposal to his superiors atIBMto develop a more practical alternative toassembly languagefor programming theirIBM 704mainframe computer.[11]: 69 Backus' historic FORTRAN team consisted of programmers Richard Goldberg, Sheldon F. Best, Harlan Herrick, Peter Sheridan,Roy Nutt,Robert Nelson, Irving Ziller, Harold Stern,Lois Haibt,andDavid Sayre.[16]Its concepts included easier entry of equations into a computer, an idea developed byJ. Halcombe Laningand demonstrated in theLaning and Zierler systemof 1952.[17]

The Fortran Automatic Coding System for theIBM 704(October 15, 1956), the first programmer's reference manual for Fortran[10]

A draft specification forThe IBM Mathematical Formula Translating Systemwas completed by November 1954.[11]: 71 The first manual for FORTRAN appeared in October 1956,[10][11]: 72 with the first FORTRANcompilerdelivered in April 1957.[11]: 75 Fortran produced efficient enough code forassembly languageprogrammers to accept ahigh-level programming languagereplacement.[18]

John Backus said during a 1979 interview withThink,the IBM employee magazine, "Much of my work has come from being lazy. I didn't like writing programs, and so, when I was working on theIBM 701,writing programs for computing missile trajectories, I started work on a programming system to make it easier to write programs. "[19]

The language was widely adopted by scientists for writing numerically intensive programs, which encouraged compiler writers to produce compilers that could generate faster and more efficient code. The inclusion of acomplex number data typein the language made Fortran especially suited to technical applications such as electrical engineering.[20]

By 1960, versions of FORTRAN were available for theIBM 709,650,1620,and7090computers. Significantly, the increasing popularity of FORTRAN spurred competing computer manufacturers to provide FORTRAN compilers for their machines, so that by 1963 over 40 FORTRAN compilers existed.

FORTRAN was provided for theIBM 1401computer by an innovative 63-phase compiler that ran entirely in itscore memoryof only 8000 (six-bit) characters. The compiler could be run from tape, or from a 2200-card deck; it used no further tape or disk storage. It kept the program in memory and loadedoverlaysthat gradually transformed it, in place, into executable form, as described by Haines.[21] This article was reprinted, edited, in both editions ofAnatomy of a Compiler[22]and in the IBM manual "Fortran Specifications and Operating Procedures, IBM 1401".[23]The executable form was not entirelymachine language;rather, floating-point arithmetic, sub-scripting, input/output, and function references were interpreted, precedingUCSD PascalP-codeby two decades.GOTRAN,a simplified, interpreted version of FORTRAN I (with only 12 statements not 32) for "load and go" operation was available (at least for the earlyIBM 1620computer).[24]Modern Fortran, and almost all later versions, are fully compiled, as done for other high-performance languages.

The development of Fortran paralleled theearly evolution of compiler technology,and many advances in the theory and design ofcompilerswere specifically motivated by the need to generate efficient code for Fortran programs.

FORTRAN

[edit]

The initial release of FORTRAN for the IBM 704[10]contained 32statements,including:

  • DIMENSIONandEQUIVALENCEstatements
  • Assignment statements
  • Three-wayarithmeticIFstatement, which passed control to one of three locations in the program depending on whether the result of the arithmetic expression was negative, zero, or positive
  • Control statements for checking exceptions (IF ACCUMULATOR OVERFLOW,IF QUOTIENT OVERFLOW,andIF DIVIDE CHECK); and control statements for manipulatingsense switches and sense lights(IF (SENSE SWITCH),IF (SENSE LIGHT),andSENSE LIGHT)
  • GO TO,computedGO TO,ASSIGN,and assignedGO TO
  • DOloops
  • Formatted I/O:FORMAT,READ,READ INPUT TAPE,WRITE OUTPUT TAPE,PRINT,andPUNCH
  • Unformatted I/O:READ TAPE,READ DRUM,WRITE TAPE,andWRITE DRUM
  • Other I/O:END FILE,REWIND,andBACKSPACE
  • PAUSE,STOP,andCONTINUE
  • FREQUENCYstatement (for providingoptimizationhints to the compiler).

The arithmeticIFstatement was reminiscent of (but not readily implementable by) a three-way comparison instruction (CAS—Compare Accumulator with Storage) available on the 704. The statement provided the only way to compare numbers—by testing their difference, with an attendant risk of overflow. This deficiency was later overcome by "logical" facilities introduced in FORTRAN IV.

TheFREQUENCYstatement was used originally (and optionally) to give branch probabilities for the three branch cases of the arithmeticIFstatement. It could also be used to suggest how many iterations aDOloop might run. The first FORTRAN compiler used this weighting to performat compile timeaMonte Carlo simulationof the generated code, the results of which were used to optimize the placement of basic blocks in memory—a very sophisticated optimization for its time. The Monte Carlo technique is documented in Backus et al.'s paper on this original implementation,The FORTRAN Automatic Coding System:

The fundamental unit of program is thebasic block;a basic block is a stretch of program which has one entry point and one exit point. The purpose of section 4 is to prepare for section 5 a table of predecessors (PRED table) which enumerates the basic blocks and lists for every basic block each of the basic blocks which can be its immediate predecessor in flow, together with the absolute frequency of each such basic block link. This table is obtained by running the program once in Monte-Carlo fashion, in which the outcome of conditional transfers arising out of IF-type statements and computed GO TO's is determined by a random number generator suitably weighted according to whatever FREQUENCY statements have been provided.[16]

The first FORTRAN compiler reported diagnostic information by halting the program when an error was found and outputting an error code on its console. That code could be looked up by the programmer in an error messages table in the operator's manual, providing them with a brief description of the problem.[10]: p.19–20 [25]Later, an error-handling subroutine to handle user errors such as division by zero, developed by NASA,[26]was incorporated, informing users of which line of code contained the error.

Fixed layout and punched cards

[edit]
FORTRAN code on apunched card,showing the specialized uses of columns 1–5, 6 and 73–80
A reproduction of a FORTRAN coding form, printed on paper and intended to be used by programmers to prepare programs for punching onto cards bykeypunchoperators. Now obsolete.

Before the development of disk files, text editors and terminals, programs were most often entered on akeypunchkeyboard onto 80-columnpunched cards,one line to a card. The resulting deck of cards would be fed into a card reader to be compiled. Punched card codes included no lower-case letters or many special characters, and special versions of the IBM 026keypunchwere offered that would correctly print the re-purposed special characters used in FORTRAN.

Reflecting punched card input practice, Fortran programs were originally written in a fixed-column format, with the first 72 columns read into twelve 36-bit words.

A letter "C" in column 1 caused the entire card to be treated as a comment and ignored by the compiler. Otherwise, the columns of the card were divided into four fields:

  • 1 to 5 were the label field: a sequence of digits here was taken as a label for use in DO or control statements such as GO TO and IF, or to identify a FORMAT statement referred to in a WRITE or READ statement. Leading zeros are ignored and 0 is not a valid label number.
  • 6 was a continuation field: a character other than a blank or a zero here caused the card to be taken as a continuation of the statement on the prior card. The continuation cards were usually numbered 1, 2,etc.and the starting card might therefore have zero in its continuation column—which is not a continuation of its preceding card.
  • 7 to 72 served as the statement field.
  • 73 to 80 were ignored (the IBM 704'scard readeronly used 72 columns).[27]

Columns 73 to 80 could therefore be used for identification information, such as punching a sequence number or text, which could be used to re-order cards if a stack of cards was dropped; though in practice this was reserved for stable, production programs. AnIBM 519could be used to copy a program deck and add sequence numbers. Some early compilers, e.g., the IBM 650's, had additional restrictions due to limitations on their card readers.[28]Keypunchescould be programmed to tab to column 7 and skip out after column 72. Later compilers relaxed most fixed-format restrictions, and the requirement was eliminated in the Fortran 90 standard.

Within the statement field,whitespace characters(blanks) were ignored outside a text literal. This allowed omitting spaces between tokens for brevity or including spaces within identifiers for clarity. For example,AVG OF Xwas a valid identifier, equivalent toAVGOFX,and101010DO101I=1,101was a valid statement, equivalent to10101DO101I=1,101because the zero in column 6 is treated as if it were a space (!), while101010DO101I=1.101was instead10101DO101I=1.101,the assignment of 1.101 to a variable calledDO101I.Note the slight visual difference between a comma and a period.

Hollerith strings,originally allowed only in FORMAT and DATA statements, were prefixed by a character count and the letter H (e.g.,26HTHIS IS ALPHANUMERIC DATA.), allowing blanks to be retained within the character string. Miscounts were a problem.

Evolution

[edit]
Timeline of Fortran language
Year Informal name ANSI Standard ISO/IEC Standard
1957 FORTRAN
1958 FORTRAN II
1958 FORTRAN III
1961 FORTRAN IV
1966 FORTRAN 66 X3.9-1966
1978 FORTRAN 77 X3.9-1978
1991 Fortran 90 X3.198-1992 1539:1991
1997 Fortran 95 1539-1:1997
2004 Fortran 2003 1539-1:2004
2010 Fortran 2008 1539-1:2010
2018 Fortran 2018 1539-1:2018
2023 Fortran 2023 1539-1:2023

FORTRAN II

[edit]

IBM'sFORTRAN IIappeared in 1958. The main enhancement was to supportprocedural programmingby allowing user-written subroutines and functions which returned values with parameters passed byreference.The COMMON statement provided a way for subroutines to access common (orglobal) variables. Six new statements were introduced:[29]

  • SUBROUTINE,FUNCTION,andEND
  • CALLandRETURN
  • COMMON

Over the next few years, FORTRAN II added support for theDOUBLE PRECISIONandCOMPLEXdata types.

Early FORTRAN compilers supported norecursionin subroutines. Early computer architectures supported no concept of a stack, and when they did directly support subroutine calls, the return location was often stored in one fixed location adjacent to the subroutine code (e.g. theIBM 1130) or a specific machine register (IBM 360et seq), which only allows recursion if a stack is maintained by software and the return address is stored on the stack before the call is made and restored after the call returns. Although not specified in FORTRAN 77, many F77 compilers supported recursion as an option, and theBurroughs mainframes,designed with recursion built-in, did so by default. It became a standard in Fortran 90 via the new keyword RECURSIVE.[30]

Simple FORTRAN II program

[edit]

This program, forHeron's formula,reads data on a tape reel containing three 5-digit integers A, B, and C as input. There are no "type" declarations available: variables whose name starts with I, J, K, L, M, or N are "fixed-point" (i.e. integers), otherwise floating-point. Since integers are to be processed in this example, the names of the variables start with the letter "I". The name of a variable must start with a letter and can continue with both letters and digits, up to a limit of six characters in FORTRAN II. If A, B, and C cannot represent the sides of a triangle in plane geometry, then the program's execution will end with an error code of "STOP 1". Otherwise, an output line will be printed showing the input values for A, B, and C, followed by the computed AREA of the triangle as a floating-point number occupying ten spaces along the line of output and showing 2 digits after the decimal point, the.2 in F10.2 of the FORMAT statement with label 601.

C AREA OF A TRIANGLE WITH A STANDARD SQUARE ROOT FUNCTION
C INPUT - TAPE READER UNIT 5, INTEGER INPUT
C OUTPUT - LINE PRINTER UNIT 6, REAL OUTPUT
C INPUT ERROR DISPLAY ERROR OUTPUT CODE 1 IN JOB CONTROL LISTING
READINPUTTAPE5,501,IA,IB,IC
501FORMAT(3I5)
C IA, IB, AND IC MAY NOT BE NEGATIVE OR ZERO
C FURTHERMORE, THE SUM OF TWO SIDES OF A TRIANGLE
C MUST BE GREATER THAN THE THIRD SIDE, SO WE CHECK FOR THAT, TOO
IF(IA)777,777,701
701IF(IB)777,777,702
702IF(IC)777,777,703
703IF(IA+IB-IC)777,777,704
704IF(IA+IC-IB)777,777,705
705IF(IB+IC-IA)777,777,799
777STOP1
C USING HERON'S FORMULA WE CALCULATE THE
C AREA OF THE TRIANGLE
799S=FLOATF(IA+IB+IC)/2.0
AREA=SQRTF(S*(S-FLOATF(IA))*(S-FLOATF(IB))*
+(S-FLOATF(IC)))
WRITEOUTPUTTAPE6,601,IA,IB,IC,AREA
601FORMAT(4HA=,I5,5HB=,I5,5HC=,I5,8HAREA=,F10.2,
+13HSQUAREUNITS)
STOP
END

FORTRAN III

[edit]

IBM also developed aFORTRAN IIIin 1958 that allowed forinline assemblycode among other features; however, this version was never released as a product. Like the 704 FORTRAN and FORTRAN II, FORTRAN III included machine-dependent features that made code written in it unportable from machine to machine.[11]: 76 Early versions of FORTRAN provided by other vendors suffered from the same disadvantage.

FORTRAN IV

[edit]

IBM began development of FORTRAN IV starting in 1961, as a result of customer demands.FORTRAN IVremoved the machine-dependent features of FORTRAN II (such asREAD INPUT TAPE), while adding new features such as aLOGICALdata type,logicalBoolean expressionsand thelogical IF statementas an alternative to thearithmetic IF statement.FORTRAN IV was eventually released in 1962, first for theIBM 7030( "Stretch" ) computer, followed by versions for theIBM 7090,IBM 7094,and later for theIBM 1401in 1966.[31]

By 1965, FORTRAN IV was supposed to be compliant with thestandardbeing developed by theAmerican Standards AssociationX3.4.3 FORTRAN Working Group.[32]

Between 1966 and 1968, IBM offered several FORTRAN IV compilers for itsSystem/360,each named by letters that indicated the minimum amount of memory the compiler needed to run. [33] The letters (F, G, H) matched the codes used with System/360 model numbers to indicate memory size, each letter increment being a factor of two larger:[34]: p. 5 

  • 1966: FORTRAN IV F for DOS/360 (64K bytes)
  • 1966: FORTRAN IV G for OS/360 (128K bytes)
  • 1968: FORTRAN IV H for OS/360 (256K bytes)

Digital Equipment Corporationmaintained DECSYSTEM-10 Fortran IV (F40) forPDP-10from 1967 to 1975.[35]Compilers were also available for theUNIVAC 1100 seriesand theControl Data6000 seriesand7000 seriessystems.[36]

At about this time FORTRAN IV had started to become an important educational tool and implementations such as the University of Waterloo's WATFOR andWATFIVwere created to simplify the complex compile and link processes of earlier compilers.

In the FORTRAN IV programming environment of the era, except for that used on Control Data Corporation (CDC) systems, only one instruction was placed per line. The CDC version allowed for multiple instructions per line if separated by a $ (dollar) character. The FORTRANsheetwas divided into four fields, as described above.

Two compilers of the time, IBM "G" and UNIVAC, allowed comments to be written on the same line as instructions, separated by a special character: "master space": V (perforations 7 and 8) for UNIVAC and perforations 12/11/0/7/8/9 (hexadecimal FF) for IBM. These comments were not to be inserted in the middle of continuation cards.[35][36]

FORTRAN 66

[edit]

Perhaps the most significant development in the early history of FORTRAN was the decision by theAmerican Standards Association(nowAmerican National Standards Institute(ANSI)) to form a committee sponsored by theBusiness Equipment Manufacturers Association(BEMA) to develop anAmerican Standard Fortran.The resulting two standards, approved in March 1966, defined two languages,FORTRAN(based on FORTRAN IV, which had served as a de facto standard), andBasic FORTRAN(based on FORTRAN II, but stripped of its machine-dependent features). The FORTRAN defined by the first standard, officially denoted X3.9-1966, became known asFORTRAN 66(although many continued to term it FORTRAN IV, the language on which the standard was largely based). FORTRAN 66 effectively became the first industry-standard version of FORTRAN. FORTRAN 66 included:

  • Main program,SUBROUTINE,FUNCTION,andBLOCK DATAprogram units
  • INTEGER,REAL,DOUBLE PRECISION,COMPLEX,andLOGICALdata types
  • COMMON,DIMENSION,andEQUIVALENCEstatements
  • DATAstatement for specifying initial values
  • IntrinsicandEXTERNAL(e.g., library) functions
  • Assignment statement
  • GO TO,computedGO TO,assignedGO TO,andASSIGNstatements
  • LogicalIFand arithmetic (three-way)IFstatements
  • DOloop statement
  • READ,WRITE,BACKSPACE,REWIND,andENDFILEstatements for sequential I/O
  • FORMATstatement and assigned format
  • CALL,RETURN,PAUSE,andSTOPstatements
  • Hollerith constantsinDATAandFORMATstatements, and as arguments to procedures
  • Identifiers of up to six characters in length
  • Comment lines
  • ENDline

The above Fortran II version of the Heron program needs several modifications to compile as a Fortran 66 program. Modifications include using the more machine independent versions of theREADandWRITEstatements, and removal of the unneededFLOATFtype conversion functions. Though not required, the arithmeticIFstatements can be re-written to use logicalIFstatements and expressions in a more structured fashion.

C AREA OF A TRIANGLE WITH A STANDARD SQUARE ROOT FUNCTION
C INPUT - TAPE READER UNIT 5, INTEGER INPUT
C OUTPUT - LINE PRINTER UNIT 6, REAL OUTPUT
C INPUT ERROR DISPLAY ERROR OUTPUT CODE 1 IN JOB CONTROL LISTING
READ(5,501)IA,IB,IC
501FORMAT(3I5)
C
C IA, IB, AND IC MAY NOT BE NEGATIVE OR ZERO
C FURTHERMORE, THE SUM OF TWO SIDES OF A TRIANGLE
C MUST BE GREATER THAN THE THIRD SIDE, SO WE CHECK FOR THAT, TOO
IF(IA.GT.0.AND.IB.GT.0.AND.IC.GT.0)GOTO10
WRITE(6,602)
602FORMAT(42HIA,IB,ANDICMUSTBEGREATERTHANZERO.)
STOP1
10CONTINUE
C
IF(IA+IB-IC.GT.0
+.AND.IA+IC-IB.GT.0
+.AND.IB+IC-IA.GT.0)GOTO20
WRITE(6,603)
603FORMAT(50HSUMOFTWOSIDESMUSTBEGREATERTHANTHIRDSIDE.)
STOP1
20CONTINUE
C
C USING HERON'S FORMULA WE CALCULATE THE
C AREA OF THE TRIANGLE
S=(IA+IB+IC)/2.0
AREA=SQRT(S*(S-IA)*(S-IB)*(S-IC))
WRITE(6,601)IA,IB,IC,AREA
601FORMAT(4HA=,I5,5HB=,I5,5HC=,I5,8HAREA=,F10.2,
+13HSQUAREUNITS)
STOP
END

FORTRAN 77

[edit]
FORTRAN-77 program with compiler output, written on aCDC175atRWTH Aachen University,Germany, in 1987
4.3 BSDfor theDigital Equipment Corporation(DEC)VAX,displaying themanualfor FORTRAN 77 (f77) compiler

After the release of the FORTRAN 66 standard, compiler vendors introduced several extensions toStandard Fortran,prompting ANSI committee X3J3 in 1969 to begin work on revising the 1966 standard, under sponsorship ofCBEMA,the Computer Business Equipment Manufacturers Association (formerly BEMA). Final drafts of this revised standard circulated in 1977, leading to formal approval of the new FORTRAN standard in April 1978. The new standard, calledFORTRAN 77and officially denoted X3.9-1978, added a number of significant features to address many of the shortcomings of FORTRAN 66:

  • BlockIFandEND IFstatements, with optionalELSE IFandELSEclauses, to provide improved language support forstructured programming
  • DOloop extensions, including parameter expressions, negative increments, and zero trip counts
  • OPEN,CLOSE,andINQUIREstatements for improved I/O capability
  • Direct-access file I/O
  • IMPLICITstatement, to override implicit conventions that undeclared variables areINTEGERif their name begins withI,J,K,L,M,orN(andREALotherwise)
  • CHARACTERdata type, replacing Hollerith strings with vastly expanded facilities for character input and output and processing of character-based data
  • PARAMETERstatement for specifying constants
  • SAVEstatement for persistent local variables
  • Generic names for intrinsic functions (e.g.SQRTalso accepts arguments of other types, such asCOMPLEXorREAL*16).
  • A set of intrinsics (LGE, LGT, LLE, LLT) forlexicalcomparison of strings, based upon theASCIIcollating sequence.(These ASCII functions were demanded by theU.S. Department of Defense,in their conditional approval vote.[citation needed])
  • A maximum of seven dimensions in arrays, rather than three. Allowed subscript expressions were also generalized.

In this revision of the standard, a number of features were removed or altered in a manner that might invalidate formerly standard-conforming programs. (Removal was the only allowable alternative to X3J3 at that time, since the concept of "deprecation"was not yet available for ANSI standards.) While most of the 24 items in the conflict list (see Appendix A2 of X3.9-1978) addressed loopholes or pathological cases permitted by the prior standard but rarely used, a small number of specific capabilities were deliberately removed, such as:

  • Hollerith constantsandHollerithdata, such asGREET=12HHELLOTHERE!
  • Reading into an H edit (Hollerith field) descriptor in a FORMAT specification
  • Overinde xing of array bounds by subscripts
    DIMENSIONA(10,5)
    Y=A(11,1)
    
  • Transfer of control out of and back into the range of a DO loop (also known as "Extended Range" )

A Fortran 77 version of the Heron program requires no modifications to the Fortran 66 version. However this example demonstrates additional cleanup of the I/O statements, including using list-directed I/O, and replacing the Hollerith edit descriptors in theFORMATstatements with quoted strings. It also uses structuredIFandEND IFstatements, rather thanGOTO/CONTINUE.

PROGRAMHERON
C AREA OF A TRIANGLE WITH A STANDARD SQUARE ROOT FUNCTION
C INPUT - DEFAULT STANDARD INPUT UNIT, INTEGER INPUT
C OUTPUT - DEFAULT STANDARD OUTPUT UNIT, REAL OUTPUT
C INPUT ERROR DISPLAY ERROR OUTPUT CODE 1 IN JOB CONTROL LISTING
READ(*,*)IA,IB,IC
C
C IA, IB, AND IC MAY NOT BE NEGATIVE OR ZERO
C FURTHERMORE, THE SUM OF TWO SIDES OF A TRIANGLE
C MUST BE GREATER THAN THE THIRD SIDE, SO WE CHECK FOR THAT, TOO
IF(IA.LE.0.OR.IB.LE.0.OR.IC.LE.0)THEN
WRITE(*,*)'IA, IB, and IC must be greater than zero.'
STOP1
END IF
C
IF(IA+IB-IC.LE.0
+.OR.IA+IC-IB.LE.0
+.OR.IB+IC-IA.LE.0)THEN
WRITE(*,*)'Sum of two sides must be greater than third side.'
STOP1
END IF
C
C USING HERON'S FORMULA WE CALCULATE THE
C AREA OF THE TRIANGLE
S=(IA+IB+IC)/2.0
AREA=SQRT(S*(S-IA)*(S-IB)*(S-IC))
WRITE(*,601)IA,IB,IC,AREA
601FORMAT('A= ',I5,' B= ',I5,' C= ',I5,' AREA= ',F10.2,
+' square units')
STOP
END

Transition to ANSI Standard Fortran

[edit]

The development of a revised standard to succeed FORTRAN 77 would be repeatedly delayed as the standardization process struggled to keep up with rapid changes in computing and programming practice. In the meantime, as the "Standard FORTRAN" for nearly fifteen years, FORTRAN 77 would become the historically most important dialect.

An important practical extension to FORTRAN 77 was the release of MIL-STD-1753 in 1978.[37]This specification, developed by theU.S. Department of Defense,standardized a number of features implemented by most FORTRAN 77 compilers but not included in the ANSI FORTRAN 77 standard. These features would eventually be incorporated into the Fortran 90 standard.

TheIEEE1003.9POSIXStandard, released in 1991, provided a simple means for FORTRAN 77 programmers to issue POSIX system calls.[38]Over 100 calls were defined in the document – allowing access to POSIX-compatible process control, signal handling, file system control, device control, procedure pointing, and stream I/O in a portable manner.

Fortran 90

[edit]

The much-delayed successor to FORTRAN 77, informally known asFortran 90(and prior to that,Fortran 8X), was finally released as ISO/IEC standard 1539:1991 in 1991 and an ANSI Standard in 1992. In addition to changing the official spelling from FORTRAN to Fortran, this major revision added many new features to reflect the significant changes in programming practice that had evolved since the 1978 standard:

  • Free-form source inputremoved the need to skip the first six character positions before entering statements.
  • Lowercase Fortran keywords
  • Identifiers up to 31 characters in length (In the previous standard, it was only six characters).
  • Inline comments
  • Ability to operate on arrays (or array sections) as a whole, thus greatly simplifying math and engineering computations.
    • whole, partial and masked array assignment statements and array expressions, such asX(1:N)=R(1:N)*COS(A(1:N))
    • WHEREstatement for selective array assignment
    • array-valued constants and expressions,
    • user-defined array-valued functions and array constructors.
  • RECURSIVEprocedures
  • Modules,to group relatedproceduresand data together, and make them available to other program units, including the capability to limit the accessibility to only specific parts of the module.
  • A vastly improved argument-passing mechanism, allowinginterfacesto be checked at compile time
  • User-written interfaces for generic procedures
  • Operator overloading
  • Derived (structured) data types
  • New data type declaration syntax, to specify the data type and other attributes of variables
  • Dynamic memory allocationby means of theALLOCATABLEattribute and theALLOCATEandDEALLOCATEstatements
  • POINTERattribute, pointer assignment, andNULLIFYstatement to facilitate the creation and manipulation of dynamicdata structures
  • Structured looping constructs, with anEND DOstatement for loop termination, andEXITandCYCLEstatements for terminating normalDOloop iterations in an orderly way
  • SELECT CASE,CASE,...CASE DEFAULT,END SELECTconstruct formulti-way selection
  • Portable specification of numerical precision under the user's control
  • New and enhanced intrinsic procedures.

Obsolescence and deletions

[edit]

Unlike the prior revision, Fortran 90 removed no features.[39]Any standard-conforming FORTRAN 77 program was also standard-conforming under Fortran 90, and either standard should have been usable to define its behavior.

A small set of features were identified as "obsolescent" and were expected to be removed in a future standard. All of the functionalities of these early-version features can be performed by newer Fortran features. Some are kept to simplify porting of old programs but many were deleted in Fortran 95.

Obsolescence and deletions
Obsolescent feature Current Status
Arithmetic IF-statement Obsolescent in F90, deleted in F2018
Non-integer DO parameters or control variables Obsolescent in F90, deleted in F95
Shared DO-loop termination or termination with a statement other than END DO or CONTINUE Obsolescent in F90, deleted in F2018
Branching to END IF from outside a block Obsolescent in F90, deleted in F95
PAUSE statement Obsolescent in F90, deleted in F95
ASSIGN statement and assigned GO TO statement Obsolescent in F90, deleted in F95
Assigned statement numbers and FORMAT specifiers Obsolescent in F90, deleted in F95
H edit descriptor Obsolescent in F90, deleted in F95
Vertical format control Deleted in F2003
Alternate return Obsolescent in F90
Computed GO TO statement Obsolescent in F90
Statement functions Obsolescent in F90
DATA statements among executable statements Obsolescent in F90
Assumed length character functions Obsolescent in F90
Fixed form source code Obsolescent in F90
CHARACTER* form of CHARACTER declaration Obsolescent in F90
ENTRY statements Obsolescent in F2008
Label form of DO statement Obsolescent in F2018
COMMON and EQUIVALENCE statements, and the BLOCK DATA program unit Obsolescent in F2018
Specific names for intrinsic function Obsolescent in F2018
FORALL construct and statement Obsolescent in F2018

"Hello, World!" example

[edit]
programhelloworld
print*,"Hello, World!"
end programhelloworld

Fortran 95

[edit]

Fortran 95,published officially as ISO/IEC 1539-1:1997, was a minor revision, mostly to resolve some outstanding issues from the Fortran 90 standard. Nevertheless, Fortran 95 also added a number of extensions, notably from theHigh Performance Fortranspecification:

  • FORALLand nestedWHEREconstructs to aid vectorization
  • User-definedPUREandELEMENTALprocedures
  • Default initialization of derived type components, including pointer initialization
  • Expanded the ability to use initialization expressions for data objects
  • Initialization of pointers toNULL()
  • Clearly defined thatALLOCATABLEarrays are automatically deallocated when they go out of scope.

A number of intrinsic functions were extended (for example adimargument was added to themaxlocintrinsic).

Several features noted in Fortran 90 to be "obsolescent" were removed from Fortran 95:

  • DOstatements usingREALandDOUBLE PRECISIONindex variables
  • Branching to anEND IFstatement from outside its block
  • PAUSEstatement
  • ASSIGNand assignedGO TOstatement, and assigned format specifiers
  • HHollerith edit descriptor.

An important supplement to Fortran 95 was theISO technical reportTR-15581: Enhanced Data Type Facilities,informally known as theAllocatable TR.This specification defined enhanced use ofALLOCATABLEarrays, prior to the availability of fully Fortran 2003-compliant Fortran compilers. Such uses includeALLOCATABLEarrays as derived type components, in procedure dummy argument lists, and as function return values. (ALLOCATABLEarrays are preferable toPOINTER-based arrays becauseALLOCATABLEarrays are guaranteed by Fortran 95 to be deallocated automatically when they go out of scope, eliminating the possibility ofmemory leakage.In addition, elements of allocatable arrays are contiguous, andaliasingis not an issue for optimization of array references, allowing compilers to generate faster code than in the case of pointers.[40])

Another important supplement to Fortran 95 was theISOtechnical reportTR-15580: Floating-point exception handling,informally known as theIEEE TR.This specification defined support forIEEE floating-point arithmeticandfloating-pointexception handling.

Conditional compilation and varying length strings

[edit]

In addition to the mandatory "Base language" (defined in ISO/IEC 1539-1: 1997), the Fortran 95 language also included two optional modules:

  • Varying length character strings (ISO/IEC 1539-2: 2000)
  • Conditional compilation (ISO/IEC 1539-3: 1998)

which, together, compose the multi-part International Standard (ISO/IEC 1539).

According to the standards developers, "the optional parts describe self-contained features which have been requested by a substantial body of users and/or implementors, but which are not deemed to be of sufficient generality for them to be required in all standard-conforming Fortran compilers." Nevertheless, if a standard-conforming Fortran does provide such options, then they "must be provided in accordance with the description of those facilities in the appropriate Part of the Standard".

Modern Fortran

[edit]

The language defined by the twenty-first century standards, in particular because of its incorporation ofobject-oriented programmingsupport and subsequentlyCoarray Fortran,is often referred to as 'Modern Fortran', and the term is increasingly used in the literature.[41]

Fortran 2003

[edit]

Fortran 2003,officially published as ISO/IEC 1539-1:2004, was a major revision introducing many new features.[42]A comprehensive summary of the new features of Fortran 2003 is available at the Fortran Working Group (ISO/IEC JTC1/SC22/WG5) official Web site.[43]

From that article, the major enhancements for this revision include:

  • Derived type enhancements: parameterized derived types, improved control of accessibility, improved structure constructors, and finalizers
  • Object-oriented programmingsupport: type extension andinheritance,polymorphism,dynamic type allocation, and type-bound procedures, providing complete support forabstract data types
  • Data manipulation enhancements: allocatable components (incorporating TR 15581), deferred type parameters,VOLATILEattribute,explicit type specification in array constructors and allocate statements, pointer enhancements, extended initialization expressions, and enhanced intrinsic procedures
  • Input/output enhancements:asynchronoustransfer, stream access, user specified transfer operations for derived types, user specified control of rounding during format conversions, named constants for preconnected units, theFLUSHstatement, regularization of keywords, and access to error messages
  • Procedure pointers
  • Support forIEEE floating-point arithmeticandfloating-pointexception handling(incorporating TR 15580)
  • Interoperability with the C programming language
  • Support for international usage: access toISO 106464-byte characters and choice of decimal or comma in numeric formatted input/output
  • Enhanced integration with the host operating system: access tocommand-linearguments,environment variables,and processor error messages

An important supplement to Fortran 2003 was theISO technical reportTR-19767: Enhanced module facilities in Fortran.This report providedsub-modules,which make Fortran modules more similar toModula-2modules. They are similar toAdaprivate child sub-units. This allows the specification and implementation of a module to be expressed in separate program units, which improves packaging of large libraries, allows preservation of trade secrets while publishing definitive interfaces, and prevents compilation cascades.

Fortran 2008

[edit]

ISO/IEC 1539-1:2010, informally known as Fortran 2008, was approved in September 2010.[44][45]As with Fortran 95, this is a minor upgrade, incorporating clarifications and corrections to Fortran 2003, as well as introducing some new capabilities. The new capabilities include:

  • Sub-modules – additional structuring facilities for modules; supersedes ISO/IEC TR 19767:2005
  • Coarray Fortran– a parallel execution model
  • TheDO CONCURRENTconstruct – for loop iterations with no interdependencies
  • The CONTIGUOUS attribute – to specify storage layout restrictions
  • TheBLOCK construct– can contain declarations of objects with construct scope
  • Recursive allocatable components– as an alternative to recursive pointers in derived types

The Final Draft international Standard (FDIS) is available as document N1830.[46]

A supplement to Fortran 2008 is theInternational Organization for Standardization(ISO) Technical Specification (TS) 29113 onFurther Interoperability of Fortran with C,[47][48]which has been submitted to ISO in May 2012 for approval. The specification adds support for accessing the array descriptor from C and allows ignoring the type and rank of arguments.

Fortran 2018

[edit]

The Fortran 2018 revision of the language was earlier referred to as Fortran 2015.[49]It was a significant revision and was released on November 28, 2018.[50]

Fortran 2018 incorporates two previously published Technical Specifications:

  • ISO/IEC TS 29113:2012 Further Interoperability with C[51]
  • ISO/IEC TS 18508:2015 Additional Parallel Features in Fortran[52]

Additional changes and new features include support for ISO/IEC/IEEE 60559:2011 (the version of theIEEE floating-point standardbefore the latest minor revision IEEE 754–2019), hexadecimal input/output, IMPLICIT NONE enhancements and other changes.[53][54][55][56]

Fortran 2018 deleted the arithmetic IF statement. It also deleted non-block DO constructs - loops which do not end with an END DO or CONTINUE statement. These had been an obsolescent part of the language since Fortran 90.

New obsolescences are: COMMON and EQUIVALENCE statements and the BLOCK DATA program unit, labelled DO loops, specific names for intrinsic functions, and the FORALL statement and construct.

Fortran 2023

[edit]

Fortran 2023 (ISO/IEC 1539-1:2023) was published in November 2023, and can be purchased from the ISO.[57] Fortran 2023 is a minor extension of Fortran 2018 that focuses on correcting errors and omissions in Fortran 2018. It also adds some small features, including anenumerated typecapability.

Language features

[edit]

A full description of the Fortran language features brought by Fortran 95 is covered in the related article,Fortran 95 language features.The language versions defined by later standards are often referred to collectively as 'Modern Fortran' and are described in the literature.

Science and engineering

[edit]
General relativistic magnetohydrodynamic Fortran simulation of black hole accretion using the BHAC code with cartesian adaptive mesh
Flow around a cylinder computed in Fortran with OpenCL (University of Bristol, UK)
Velocity and sea surface temperature in the oceans, computed with the NEMO Fortran code (Nucleus for European Modeling of the Ocean,https:// nemo-ocean.eu) in theBarcelona Supercomputing Center(2020)

Although a 1968 journal article by the authors ofBASICalready described FORTRAN as "old-fashioned",[58]programs have been written in Fortran for many decades and there is a vast body of Fortran software in daily use throughout the scientific and engineering communities.[59]Jay Pasachoffwrote in 1984 that "physics and astronomy students simply have to learn FORTRAN. So much exists in FORTRAN that it seems unlikely that scientists will change toPascal,Modula-2,or whatever. "[60]In 1993,Cecil E. Leithcalled FORTRAN the "mother tongue of scientific computing", adding that its replacement by any other possible language "may remain a forlorn hope".[61]

It is the primary language for some of the most intensivesuper-computingtasks, such as inastronomy,climate modeling,computational chemistry,computational economics,computational fluid dynamics,computational physics,data analysis,[62]hydrological modeling,numerical linear algebra and numerical libraries (LAPACK,IMSLandNAG),optimization,satellite simulation,structural engineering,andweather prediction.[63]Many of the floating-point benchmarks to gauge the performance of new computer processors, such as the floating-point components of theSPECbenchmarks (e.g.,CFP2006,CFP2017) are written in Fortran. Math algorithms are well documented inNumerical Recipes.

Apart from this, more modern codes in computational science generally use large program libraries, such asMETISfor graph partitioning,PETScorTrilinosfor linear algebra capabilities,deal.IIorFEniCSfor mesh and finite element support, and other generic libraries. Since the early 2000s, many of the widely used support libraries have also been implemented inCand more recently, inC++.On the other hand, high-level languages such as theWolfram Language,MATLAB,Python,andRhave become popular in particular areas of computational science. Consequently, a growing fraction of scientific programs are also written in such higher-level scripting languages. For this reason,facilities for inter-operation with Cwere added to Fortran 2003 and enhanced by the ISO/IEC technical specification 29113, which was incorporated into Fortran 2018 to allow more flexible interoperation with other programming languages.

Portability

[edit]

Portabilitywas a problem in the early days because there was no agreed upon standard—not even IBM's reference manual—and computer companies vied to differentiate their offerings from others by providing incompatible features. Standards have improved portability. The 1966 standard provided a referencesyntaxand semantics, but vendors continued to provide incompatible extensions. Although careful programmers were coming to realize that use of incompatible extensions caused expensive portability problems, and were therefore using programs such asThe PFORT Verifier,[64][65]it was not until after the 1977 standard, when the National Bureau of Standards (nowNIST) publishedFIPS PUB 69,that processors purchased by the U.S. Government were required to diagnose extensions of the standard. Rather than offer two processors, essentially every compiler eventually had at least an option to diagnose extensions.[66][67]

Incompatible extensions were not the only portability problem. For numerical calculations, it is important to take account of the characteristics of the arithmetic. This was addressed by Fox et al. in the context of the 1966 standard by thePORTlibrary.[65]The ideas therein became widely used, and were eventually incorporated into the 1990 standard by way of intrinsic inquiry functions. The widespread (now almost universal) adoption of theIEEE 754standard for binary floating-point arithmetic has essentially removed this problem.

Access to the computing environment (e.g., the program's command line, environment variables, textual explanation of error conditions) remained a problem until it was addressed by the 2003 standard.

Large collections of library software that could be described as being loosely related to engineering and scientific calculations, such as graphics libraries, have been written in C, and therefore access to them presented a portability problem. This has been addressed by incorporation of C interoperability into the 2003 standard.

It is now possible (and relatively easy) to write an entirely portable program in Fortran, even without recourse to apreprocessor.

Obsolete variants

[edit]

Until the Fortran 66 standard was developed, each compiler supported its own variant of Fortran. Some were more divergent from the mainstream than others.

The first Fortran compiler set a high standard of efficiency for compiled code. This goal made it difficult to create a compiler so it was usually done by the computer manufacturers to support hardware sales. This left an important niche: compilers that were fast and provided good diagnostics for the programmer (often a student). Examples include Watfor, Watfiv, PUFFT, and on a smaller scale, FORGO, Wits Fortran, and Kingston Fortran 2.

Fortran 5was marketed byData GeneralCorp from the early 1970s to the early 1980s, for theNova,Eclipse,andMVline of computers. It had an optimizing compiler that was quite good for minicomputers of its time. The language most closely resembles FORTRAN 66.

FORTRAN Vwas distributed byControl Data Corporationin 1968 for theCDC 6600series. The language was based upon FORTRAN IV.[68]

Univac also offered a compiler for the 1100 series known as FORTRAN V. A spinoff of Univac Fortran V was Athena FORTRAN.

Specific variantsproduced by the vendors of high-performance scientific computers (e.g.,Burroughs,Control Data Corporation(CDC),Cray,Honeywell,IBM,Texas Instruments,andUNIVAC) added extensions to Fortran to take advantage of special hardware features such asinstruction cache,CPUpipelines,and vector arrays. For example, one of IBM's FORTRAN compilers (H Extended IUP) had a level of optimization which reordered themachine codeinstructionsto keep multiple internal arithmetic units busy simultaneously. Another example isCFD,a special variant of FORTRAN designed specifically for theILLIAC IVsupercomputer, running atNASA'sAmes Research Center. IBM Research Labs also developed an extended FORTRAN-based language calledVECTRANfor processing vectors and matrices.

Object-Oriented Fortranwas an object-oriented extension of Fortran, in which data items can be grouped into objects, which can be instantiated and executed in parallel. It was available for Sun, Iris,iPSC,and nCUBE, but is no longer supported.

Such machine-specific extensions have either disappeared over time or have had elements incorporated into the main standards. The major remaining extension isOpenMP,which is a cross-platform extension for shared memory programming. One new extension, Coarray Fortran, is intended to support parallel programming.

FOR TRANSITwas the name of a reduced version of the IBM 704 FORTRAN language, which was implemented for the IBM 650, using a translator program developed at Carnegie in the late 1950s.[69]The following comment appears in the IBM Reference Manual (FOR TRANSIT Automatic Coding SystemC28-4038, Copyright 1957, 1959 by IBM):

The FORTRAN system was designed for a more complex machine than the 650, and consequently some of the 32 statements found in the FORTRAN Programmer's Reference Manual are not acceptable to the FOR TRANSIT system. In addition, certain restrictions to the FORTRAN language have been added. However, none of these restrictions make a source program written for FOR TRANSIT incompatible with the FORTRAN system for the 704.

The permissible statements were:

  • Arithmetic assignment statements, e.g.,a = b
  • GO to n
  • GO TO (n1,n2,..., nm), i
  • IF (a) n1,n2,n3
  • PAUSE
  • STOP
  • DO n i = m1, m2
  • CONTINUE
  • END
  • READ n, list
  • PUNCH n, list
  • DIMENSION V, V, V,...
  • EQUIVALENCE (a,b,c), (d,c),...

Up to ten subroutines could be used in one program.

FOR TRANSIT statements were limited to columns 7 through 56, only. Punched cards were used for input and output on the IBM 650. Three passes were required to translate source code to the "IT" language, then to compile the IT statements into SOAP assembly language, and finally to produce the object program, which could then be loaded into the machine to run the program (using punched cards for data input, and outputting results onto punched cards).

Two versions existed for the 650s with a 2000 word memory drum: FOR TRANSIT I (S) and FOR TRANSIT II, the latter for machines equipped with inde xing registers and automatic floating-point decimal (bi-quinary) arithmetic. Appendix A of the manual included wiring diagrams for theIBM 533card reader/punchcontrol panel.

Fortran-based languages

[edit]

Prior to FORTRAN 77, manypreprocessorswere commonly used to provide a friendlier language, with the advantage that the preprocessed code could be compiled on any machine with a standard FORTRAN compiler.[70]These preprocessors would typically supportstructured programming,variable names longer than six characters, additional data types,conditional compilation,and evenmacrocapabilities. Popular preprocessors includedEFL,FLECS,iftran,MORTRAN,SFtran,S-Fortran,Ratfor,andRatfiv.EFL, Ratfor and Ratfiv, for example, implementedC-like languages, outputting preprocessed code in standard FORTRAN 66. ThePFORTpreprocessor was often used to verify that code conformed to a portable subset of the language. Despite advances in the Fortran language, preprocessors continue to be used for conditional compilation and macro substitution.

One of the earliest versions of FORTRAN, introduced in the '60s, was popularly used in colleges and universities. Developed, supported, and distributed by theUniversity of Waterloo,WATFORwas based largely on FORTRAN IV. A student using WATFOR could submit their batch FORTRAN job and, if there were no syntax errors, the program would move straight to execution. This simplification allowed students to concentrate on their program's syntax and semantics, or execution logic flow, rather than dealing with submissionJob Control Language(JCL), the compile/link-edit/execution successive process(es), or other complexities of the mainframe/minicomputer environment. A down side to this simplified environment was that WATFOR was not a good choice for programmers needing the expanded abilities of their host processor(s), e.g., WATFOR typically had very limited access to I/O devices. WATFOR was succeeded byWATFIVand its later versions.

program;s=0i=1,n;s=s+1;stopi;s='s'Stop

(line programming)

LRLTRANwas developed at theLawrence Radiation Laboratoryto provide support for vector arithmetic and dynamic storage, among other extensions to support systems programming. The distribution included theLivermore Time Sharing System(LTSS) operating system.

The Fortran-95 Standard includes an optionalPart 3which defines an optionalconditional compilationcapability. This capability is often referred to as "CoCo".

Many Fortran compilers have integrated subsets of theC preprocessorinto their systems.

SIMSCRIPTis an application specific Fortran preprocessor for modeling and simulating large discrete systems.

TheF programming languagewas designed to be a clean subset of Fortran 95 that attempted to remove the redundant, unstructured, and deprecated features of Fortran, such as theEQUIVALENCEstatement. F retains the array features added in Fortran 90, and removes control statements that were made obsolete by structured programming constructs added to both FORTRAN 77 and Fortran 90. F is described by its creators as "a compiled, structured, array programming language especially well suited to education and scientific computing".[71]Essential Lahey Fortran 90 (ELF90) was a similar subset.

Lahey and Fujitsu teamed up to create Fortran for the Microsoft.NET Framework.[72]Silverfrost FTN95 is also capable of creating.NET code.[73]

Code examples

[edit]

The following program illustrates dynamic memory allocation and array-based operations, two features introduced with Fortran 90. Particularly noteworthy is the absence ofDOloops andIF/THENstatements in manipulating the array; mathematical operations are applied to the array as a whole. Also apparent is the use of descriptive variable names and general code formatting that conform with contemporary programming style. This example computes an average over data entered interactively.

programaverage

!Read in some numbers and take the average
!As written, if there are no data points, an average of zero is returned
!While this may not be desired behavior, it keeps this example simple

implicit none

real,allocatable::points(:)
integer::number_of_points
real::average_points,positive_average,negative_average
average_points=0.
positive_average=0.
negative_average=0.
write(*,*)"Input number of points to average:"
read(*,*)number_of_points

allocate(points(number_of_points))

write(*,*)"Enter the points to average:"
read(*,*)points

!Take the average by summing points and dividing by number_of_points
if(number_of_points>0)average_points=sum(points)/number_of_points

!Now form average over positive and negative points only
if(count(points>0.)>0)positive_average=sum(points,points>0.)/count(points>0.)
if(count(points<0.)>0)negative_average=sum(points,points<0.)/count(points<0.)

!Print result to terminal stdout unit 6
write(*,'(a,g12.4)')'Average = ',average_points
write(*,'(a,g12.4)')'Average of positive points = ',positive_average
write(*,'(a,g12.4)')'Average of negative points = ',negative_average
deallocate(points)!free memory

end programaverage

Humor

[edit]

During the same FORTRAN standards committee meeting at which the name "FORTRAN 77" was chosen, a satirical technical proposal was incorporated into the official distribution bearing the title "Letter OConsidered Harmful".This proposal purported to address the confusion that sometimes arises between the letter" O "and the numeral zero, by eliminating the letter from allowable variable names. However, the method proposed was to eliminate the letter from the character set entirely (thereby retaining 48 as the number of lexical characters, which the colon had increased to 49). This was considered beneficial in that it would promote structured programming, by making it impossible to use the notoriousGO TOstatement as before. (TroublesomeFORMATstatements would also be eliminated.) It was noted that this "might invalidate some existing programs" but that most of these "probably were non-conforming, anyway".[74][unreliable source?][75]

When X3J3 debated whether the minimum trip count for a DO loop should be zero or one in Fortran 77, Loren Meissner suggested a minimum trip count of two—reasoning(tongue-in-cheek)that if it were less than two, then there would be no reason for a loop.

When assumed-length arrays were being added, there was a dispute as to the appropriate character to separate upper and lower bounds. In a comment examining these arguments, Walt Brainerd penned an article entitled "Astronomy vs. Gastroenterology" because some proponents had suggested using the star or asterisk ( "*" ), while others favored the colon ( ":" ).[citation needed]

Variable names beginning with the letters I–N have a default type of integer, while variables starting with any other letters defaulted to real, although programmers could override the defaults with an explicit declaration.[76]This led to the joke: "In FORTRAN, GOD is REAL (unless declared INTEGER)."

See also

[edit]

References

[edit]
  1. ^"Chapel spec (Acknowledgements)"(PDF).Cray Inc. October 1, 2015.Archived(PDF)from the original on February 5, 2016.RetrievedJanuary 14,2016.
  2. ^Report of a Visit to Discuss Common Programming Languages in Czechoslowakia and Poland, 1963, John A. Gosden (Editor), Roger E. Gay, John L. Jones, Jack N. Merner, Christopher J. Shaw
  3. ^John Backus."The history of FORTRAN I, II and III"(PDF).Softwarepreservation.org.Archived(PDF)from the original on July 26, 2007.RetrievedNovember 19,2014.
  4. ^Wilson, Leslie B. (2001).Comparative Programming Languages, Third Edition.Addison-Wesley. p. 16.ISBN0-201-71012-9.The manual for Fortran I was released in 1956, but it was 1958 before successful compilers were running programs correctly.
  5. ^Loh, Eugene (June 18, 2010)."The Ideal HPC Programming Language".ACM Queue.8(6): 30–38.doi:10.1145/1810226.1820518.
  6. ^"HPL – A Portable Implementation of the High-Performance Linpack Benchmark for Distributed-Memory Computers".RetrievedFebruary 21,2015.
  7. ^"Q13. What are the benchmarks?".Overview – CPU 2017.SPEC.RetrievedNovember 13,2019.
  8. ^Wilson, Leslie B. (2001).Comparative Programming Languages, Third Edition.Addison-Wesley. p. 18.ISBN0-201-71012-9.Another problem was that there was no standard for Fortran and so slightly different versions... would likely fail when used with a different compiler.
  9. ^TIOBE Software BV (May 2024)."TIOBE Index".TIOBE.TIOBE.RetrievedMay 6,2024.
  10. ^abcdeBackus, John Warner;Beeber, R. J.; Best, Sheldon F.;Goldberg, Richard;Herrick, Harlan L.; Hughes, R. A.; Mitchell, L. B.; Nelson, Robert A.;Nutt, Roy;Sayre, David;Sheridan, Peter B.; Stern, Harold; Ziller, Irving (October 15, 1956).Sayre, David(ed.).The FORTRAN Automatic Coding System for the IBM 704 EDPM: Programmer's Reference Manual(PDF).New York, USA: Applied Science Division and Programming Research Department,International Business Machines Corporation.pp. 2, 19–20.Archived(PDF)from the original on July 4, 2022.RetrievedJuly 4,2022.(2+51+1 pages)
  11. ^abcdefBackus, John(October–December 1998)."The History of Fortran I, II, and III"(PDF).IEEE Annals of the History of Computing.20(4): 68–78.doi:10.1109/85.728232.Archived(PDF)from the original on March 3, 2016.RetrievedJune 17,2020.[1][2]
  12. ^"FORTRAN: The Pioneering Programming Language".IBM.March 7, 2012.RetrievedJuly 19,2017.
  13. ^"FORTRAN".Encyclopedia Britannica.June 14, 2021.RetrievedJuly 19,2021.
  14. ^Chapman, Stephen J. (2018).Fortran for Scientists and Engineers(Fourth ed.).New York:McGraw-Hill Education.p. 13.ISBN978-0-07-338589-1.
  15. ^The "Fortran character set" defined by the FORTRAN 77 standard was the minimal character set that standard-compliant compilers were required to support; in practice, many FORTRAN 77 compilers supported the fullASCIIcharacter set.
  16. ^abJ. W. Backus;R. J. Beeber; S. Best; R. Goldberg;L. M. Haibt;H. L. Herrick; R. A. Nelson;D. Sayre;P. B. Sheridan; H. Stern; L. Ziller; R. A. Hughes;R. Nutt(February 1957).The FORTRAN Automatic Coding System(PDF).Western Joint Computer Conference. pp. 188–198.doi:10.1145/1455567.1455599.
  17. ^Mindell, David (2008).Digital Apollo.Cambridge, MA: MIT Press. p. 99.ISBN978-0-262-13497-2.JSTORj.ctt5hhn02.OCLC228031739.
  18. ^Padua, David (January–February 2000)."The Fortran I Compiler"(PDF).Computing in Science & Engineering.2(1): 70–75.Bibcode:2000CSE.....2a..70P.doi:10.1109/5992.814661.Archived fromthe original(PDF)on June 17, 2020.The Fortran I compiler was the first major project in code optimization. It tackled problems of crucial importance whose general solution was an important research focus in compiler technology for several decades. Many classical techniques for compiler analysis and optimization can trace their origins and inspiration to the Fortran I compiler.
  19. ^Brian Bergstein (May 20, 2007)."Fortran creator John Backus dies".MSNBC.RetrievedOctober 29,2018.
  20. ^"FORTRAN – The First Programming Language for Numeric Calculations".SciHi Blog.October 15, 2019.
  21. ^Haines, L. H. (1965). "Serial compilation and the 1401 FORTRAN compiler".IBM Systems Journal.4(1): 73–80.doi:10.1147/sj.41.0073.
  22. ^Lee, John A. N. (1967).Anatomy of a Compiler.Van Nostrand Reinhold.
  23. ^Fortran Specifications and Operating Procedures, IBM 1401(PDF).IBM. C24-1455-2.Archived(PDF)from the original on September 23, 2017.
  24. ^"Subject: GOTRAN ON THE IBM 1316 DISK STORAGE PACK (Serial No. K7402)".Archived fromthe originalon February 21, 2006.RetrievedMarch 3,2006.
  25. ^The FORTRAN Automatic Coding System for the IBM 704 EDPM: Preliminary Operator's Manual(PDF).Programming Research Department, International Business Machines Corporation. April 8, 1957. pp. 6–37.Archived(PDF)from the original on February 26, 2014.
  26. ^Armstead, Betty Jo (January 21, 2015)."My Years at NASA"(PDF).Denver Museum of Nature & Science.Archived fromthe original(PDF)on December 24, 2019.RetrievedJune 15,2019.
  27. ^Reference Manual, IBM 7090 Data Processing System(PDF).1961. A22-6528-3.Archived(PDF)from the original on December 1, 2008.
  28. ^Fortran II General Information Manual(PDF).1963.Archived(PDF)from the original on April 26, 2005.RetrievedNovember 19,2014.
  29. ^Reference Manual, FORTRAN II for the IBM 704 Data Processing System(PDF).1958. C28-6000-2.Archived(PDF)from the original on October 30, 2005.
  30. ^"Ibibilio.org".Ibiblio.org.RetrievedSeptember 15,2014.
  31. ^Fortran IV Language Specifications, Program Specifications, and Operating Procedures, IBM 1401, 1440, and 1460(PDF).IBM. April 1966. C24-3322-2.Archived(PDF)from the original on September 23, 2017.
  32. ^McCracken, Daniel D. (1965)."Preface".A Guide to FORTRAN IV Programming.New York: Wiley. p. v.ISBN978-0-471-58281-6.
  33. ^"List of FORTRAN Implementations 1957 – 1967".BCS Fortran Specialist Group.IEEE Annals. 2017.RetrievedOctober 17,2017.
  34. ^IBM System/360 Model 50 Functional Characteristics(PDF).IBM. 1967. A22-6898-1.Archived(PDF)from the original on October 29, 2021 – via bitsavers.
  35. ^ab"DECSYSTEM-10 FORTRAN IV (F40) Programmers Reference Manual"(PDF).Github.Digital Equipment Corporation.RetrievedApril 15,2022.
  36. ^ab"FORTRAN IV".WorldCat.org.RetrievedDecember 10,2023.
  37. ^MIL-STD 1753 - FORTRAN, DoD Supplement To American National Standard X3.9-1978.United States Government Printing Office.November 9, 1978.RetrievedApril 21,2024.
  38. ^IEEE 1003.9-1992 – IEEE Standard for InformationTechnology – POSIX(R) FORTRAN 77 Language Interfaces – Part 1: Binding for System Application Program Interface (API).IEEE.RetrievedNovember 24,2018.
  39. ^Appendix B.1
  40. ^"Fortran 95 Reference".Gnu.Org.RetrievedMay 10,2014.
  41. ^Lionel, Steve (December 30, 2013)."Doctor Fortran in" It's a Modern Fortran World "".Intel (was DEC) Fortran forum discussion group.Intel.RetrievedApril 11,2022.
  42. ^"Fortran 2003– Last Working Draft".Gnu.Org.RetrievedMay 10,2014.
  43. ^"WG5 Completes Processing Fortran 2003 and the TR".nag.co.uk.May 14, 2004. Archived fromthe originalon August 5, 2004.RetrievedApril 3,2023..It may also be downloaded as a PDF file at"The New Features of Fortran 2003"(PDF).Archived(PDF)from the original on January 15, 2018.RetrievedApril 3,2023.
  44. ^"N1836, Summary of Voting/Table of Replies on ISO/IEC FDIS 1539-1, Information technology – Programming languages – Fortran – Part 1: Base language"(PDF).Archived(PDF)from the original on January 15, 2018.
  45. ^"Fortran 2008 – Last Working Draft"(PDF).Gnu.Org.Archived(PDF)from the original on May 12, 2014.RetrievedMay 10,2014.
  46. ^N1830, Information technology – Programming languages – Fortran – Part 1: Base language[3][dead link]
  47. ^ISO/IEC TS 29113:2012 – Information technology – Further interoperability of Fortran with C.
  48. ^"Draft of the Technical Specification (TS) 29113"(PDF).wg5-fortran.org.Archived(PDF)from the original on January 15, 2018.
  49. ^"Doctor Fortran in" Eighteen is the new Fifteen "".Software.intel.RetrievedNovember 20,2017.
  50. ^"Fortran 2018".ISO.RetrievedNovember 30,2018.
  51. ^"Further Interoperability with C"(PDF).ISO.Archived(PDF)from the original on December 1, 2017.RetrievedNovember 20,2017.
  52. ^"Additional Parallel Features in Fortran".ISO.RetrievedNovember 20,2017.
  53. ^"The New Features of Fortran 2015".ISO.RetrievedJune 23,2017.
  54. ^"Doctor Fortran in" One Door Closes "".Software.intel.RetrievedSeptember 21,2015.
  55. ^"Doctor Fortran Goes Dutch: Fortran 2015".Software.intel.RetrievedNovember 19,2014.
  56. ^Fortran 2018 Interpretation Document(PDF).October 9, 2018.
  57. ^"Fortran 2023".wg5-fortran.org.RetrievedNovember 30,2023.
  58. ^Kemeny, John G.; Kurtz, Thomas E. (October 11, 1968)."Dartmouth Time-Sharing".Science.162(3850): 223–228.Bibcode:1968Sci...162..223K.doi:10.1126/science.162.3850.223.PMID5675464.
  59. ^Phillips, Lee (May 8, 2014)."Scientific computing's future: Can any coding language top a 1950s behemoth?".Ars Technica.RetrievedMay 8,2014.
  60. ^Pasachoff, Jay M.(April 1984)."Scientists: FORTRAN vs. Modula-2".BYTE(letter). p. 404.RetrievedFebruary 6,2015.
  61. ^Galperin, Boris (1993). "26".Large Eddy Simulation of Complex Engineering and Geophysical Flows.London: Cambridgey. p. 573.ISBN978-0-521-43009-8.
  62. ^Brandt, Siegmund (1970).Statistical and Computational Methods in Data Analysis.North-Holland Publishing Company. pp. Preface.ISBN0-7204-2033-4.
  63. ^Curcic, Milan (November 24, 2020).Modern Fortran: Building efficient parallel applications.Simon and Schuster.ISBN978-1-61729-528-7.
  64. ^Gaffney, P W; Wooten, J W (May 1, 1980). Methods to ensure the standardization of FORTRAN software (Report). Oak Ridge National Laboratory.OSTI5361454.PFORT... Library...
  65. ^abP. A. Fox (1977). "Port — A portable mathematical subroutine library".A portable mathematical subroutine library.Lecture Notes in Computer Science. Vol. 57. pp. 165–177.doi:10.1007/3-540-08446-0_42.ISBN978-3-540-08446-4.PORT... written in (PFORT).. ANS Fortran
  66. ^Whitten, Douglas E.; Demaine, Paul A. D. (1975). "A machine and configuration independent Fortran: Portable Fortran {PFortran}".IEEE Transactions on Software Engineering.SE-1(1). Institute of Electrical and Electronics Engineers (IEEE): 111–124.doi:10.1109/tse.1975.6312825.ISSN0098-5589.S2CID16485156.
  67. ^"Portability Issues"... discusses.. portability of.. Fortran
  68. ^Healy, MJR (1968)."Towards FORTRAN VI".Advanced scientific Fortran by CDC.CDC. pp. 169–172. Archived fromthe originalon July 5, 2009.RetrievedApril 10,2009.
  69. ^Perlis, A. J.; Smith, J. W.; Van Zoeren, H. R.Internal Translator (IT) A Compiler for the IBM 650(PDF).Computation Center, Carnegie Institute of Technology.Archived(PDF)from the original on October 2, 2012.
  70. ^Meissner, Loren (September 1, 1975)."On extending Fortran control structures to facilitate structured programming".ACM SIGPLAN Notices.acm.org.
  71. ^"F Programming Language Homepage".Fortran. Archived fromthe originalon January 9, 2015.RetrievedNovember 19,2014.
  72. ^"Fortran for.NET Language System".Archived fromthe originalon October 18, 2014.
  73. ^"FTN95: Fortran 95 for Windows".Silverfrost.RetrievedNovember 19,2014.
  74. ^X3J3 post-meeting distribution for meeting held at Brookhaven National Laboratory in November 1976.
  75. ^"The obliteration of O".Computer Weekly.March 3, 1977.
  76. ^Rules for Data Typing (FORTRAN 77 Language Reference).RetrievedSeptember 29,2016.

Further reading

[edit]
Language standards
Related standards
Other reference material
Books
Articles
[edit]