1 \input texinfo @c -*-texinfo-*- 2 @c %**start of header 3 @setfilename gfortran.info 4 @set copyrights-gfortran 1999-2022 5 6 @include gcc-common.texi 7 8 @settitle The GNU Fortran Compiler 9 10 @c Create a separate index for command line options 11 @defcodeindex op 12 @c Merge the standard indexes into a single one. 13 @syncodeindex fn cp 14 @syncodeindex vr cp 15 @syncodeindex ky cp 16 @syncodeindex pg cp 17 @syncodeindex tp cp 18 19 @c TODO: The following "Part" definitions are included here temporarily 20 @c until they are incorporated into the official Texinfo distribution. 21 @c They borrow heavily from Texinfo's \unnchapentry definitions. 22 23 @tex 24 \gdef\part#1#2{% 25 \pchapsepmacro 26 \gdef\thischapter{} 27 \begingroup 28 \vglue\titlepagetopglue 29 \titlefonts \rm 30 \leftline{Part #1:@* #2} 31 \vskip4pt \hrule height 4pt width \hsize \vskip4pt 32 \endgroup 33 \writetocentry{part}{#2}{#1} 34 } 35 \gdef\blankpart{% 36 \writetocentry{blankpart}{}{} 37 } 38 % Part TOC-entry definition for summary contents. 39 \gdef\dosmallpartentry#1#2#3#4{% 40 \vskip .5\baselineskip plus.2\baselineskip 41 \begingroup 42 \let\rm=\bf \rm 43 \tocentry{Part #2: #1}{\doshortpageno\bgroup#4\egroup} 44 \endgroup 45 } 46 \gdef\dosmallblankpartentry#1#2#3#4{% 47 \vskip .5\baselineskip plus.2\baselineskip 48 } 49 % Part TOC-entry definition for regular contents. This has to be 50 % equated to an existing entry to not cause problems when the PDF 51 % outline is created. 52 \gdef\dopartentry#1#2#3#4{% 53 \unnchapentry{Part #2: #1}{}{#3}{#4} 54 } 55 \gdef\doblankpartentry#1#2#3#4{} 56 @end tex 57 58 @c %**end of header 59 60 @c Use with @@smallbook. 61 62 @c %** start of document 63 64 @c Cause even numbered pages to be printed on the left hand side of 65 @c the page and odd numbered pages to be printed on the right hand 66 @c side of the page. Using this, you can print on both sides of a 67 @c sheet of paper and have the text on the same part of the sheet. 68 69 @c The text on right hand pages is pushed towards the right hand 70 @c margin and the text on left hand pages is pushed toward the left 71 @c hand margin. 72 @c (To provide the reverse effect, set bindingoffset to -0.75in.) 73 74 @c @tex 75 @c \global\bindingoffset=0.75in 76 @c \global\normaloffset =0.75in 77 @c @end tex 78 79 @copying 80 Copyright @copyright{} @value{copyrights-gfortran} Free Software Foundation, Inc. 81 82 Permission is granted to copy, distribute and/or modify this document 83 under the terms of the GNU Free Documentation License, Version 1.3 or 84 any later version published by the Free Software Foundation; with the 85 Invariant Sections being ``Funding Free Software'', the Front-Cover 86 Texts being (a) (see below), and with the Back-Cover Texts being (b) 87 (see below). A copy of the license is included in the section entitled 88 ``GNU Free Documentation License''. 89 90 (a) The FSF's Front-Cover Text is: 91 92 A GNU Manual 93 94 (b) The FSF's Back-Cover Text is: 95 96 You have freedom to copy and modify this GNU Manual, like GNU 97 software. Copies published by the Free Software Foundation raise 98 funds for GNU development. 99 @end copying 100 101 @ifinfo 102 @dircategory Software development 103 @direntry 104 * gfortran: (gfortran). The GNU Fortran Compiler. 105 @end direntry 106 This file documents the use and the internals of 107 the GNU Fortran compiler, (@command{gfortran}). 108 109 Published by the Free Software Foundation 110 51 Franklin Street, Fifth Floor 111 Boston, MA 02110-1301 USA 112 113 @insertcopying 114 @end ifinfo 115 116 117 @setchapternewpage odd 118 @titlepage 119 @title Using GNU Fortran 120 @versionsubtitle 121 @author The @t{gfortran} team 122 @page 123 @vskip 0pt plus 1filll 124 Published by the Free Software Foundation@* 125 51 Franklin Street, Fifth Floor@* 126 Boston, MA 02110-1301, USA@* 127 @c Last printed ??ber, 19??.@* 128 @c Printed copies are available for $? each.@* 129 @c ISBN ??? 130 @sp 1 131 @insertcopying 132 @end titlepage 133 134 @c TODO: The following "Part" definitions are included here temporarily 135 @c until they are incorporated into the official Texinfo distribution. 136 137 @tex 138 \global\let\partentry=\dosmallpartentry 139 \global\let\blankpartentry=\dosmallblankpartentry 140 @end tex 141 @summarycontents 142 143 @tex 144 \global\let\partentry=\dopartentry 145 \global\let\blankpartentry=\doblankpartentry 146 @end tex 147 @contents 148 149 @page 150 151 @c --------------------------------------------------------------------- 152 @c TexInfo table of contents. 153 @c --------------------------------------------------------------------- 154 155 @ifnottex 156 @node Top 157 @top Introduction 158 @cindex Introduction 159 160 This manual documents the use of @command{gfortran}, 161 the GNU Fortran compiler. You can find in this manual how to invoke 162 @command{gfortran}, as well as its features and incompatibilities. 163 164 @ifset DEVELOPMENT 165 @emph{Warning:} This document, and the compiler it describes, are still 166 under development. While efforts are made to keep it up-to-date, it might 167 not accurately reflect the status of the most recent GNU Fortran compiler. 168 @end ifset 169 170 @comment 171 @comment When you add a new menu item, please keep the right hand 172 @comment aligned to the same column. Do not use tabs. This provides 173 @comment better formatting. 174 @comment 175 @menu 176 * Introduction:: 177 178 Part I: Invoking GNU Fortran 179 * Invoking GNU Fortran:: Command options supported by @command{gfortran}. 180 * Runtime:: Influencing runtime behavior with environment variables. 181 182 Part II: Language Reference 183 * Compiler Characteristics:: User-visible implementation details. 184 * Extensions:: Language extensions implemented by GNU Fortran. 185 * Mixed-Language Programming:: Interoperability with C 186 * Coarray Programming:: 187 * Intrinsic Procedures:: Intrinsic procedures supported by GNU Fortran. 188 * Intrinsic Modules:: Intrinsic modules supported by GNU Fortran. 189 190 * Contributing:: How you can help. 191 * Copying:: GNU General Public License says 192 how you can copy and share GNU Fortran. 193 * GNU Free Documentation License:: 194 How you can copy and share this manual. 195 * Funding:: How to help assure continued work for free software. 196 * Option Index:: Index of command line options 197 * Keyword Index:: Index of concepts 198 @end menu 199 @end ifnottex 200 201 @c --------------------------------------------------------------------- 202 @c Introduction 203 @c --------------------------------------------------------------------- 204 205 @node Introduction 206 @chapter Introduction 207 208 @c The following duplicates the text on the TexInfo table of contents. 209 @iftex 210 This manual documents the use of @command{gfortran}, the GNU Fortran 211 compiler. You can find in this manual how to invoke @command{gfortran}, 212 as well as its features and incompatibilities. 213 214 @ifset DEVELOPMENT 215 @emph{Warning:} This document, and the compiler it describes, are still 216 under development. While efforts are made to keep it up-to-date, it 217 might not accurately reflect the status of the most recent GNU Fortran 218 compiler. 219 @end ifset 220 @end iftex 221 222 @menu 223 * About GNU Fortran:: What you should know about the GNU Fortran compiler. 224 * GNU Fortran and GCC:: You can compile Fortran, C, or other programs. 225 * Standards:: Standards supported by GNU Fortran. 226 @end menu 227 228 229 @c --------------------------------------------------------------------- 230 @c About GNU Fortran 231 @c --------------------------------------------------------------------- 232 233 @node About GNU Fortran 234 @section About GNU Fortran 235 236 The GNU Fortran compiler is the successor to @command{g77}, the 237 Fortran 77 front end included in GCC prior to version 4 (released in 238 2005). While it is backward-compatible with most @command{g77} 239 extensions and command-line options, @command{gfortran} is a completely new 240 implemention designed to support more modern dialects of Fortran. 241 GNU Fortran implements the Fortran 77, 90 and 95 standards 242 completely, most of the Fortran 2003 and 2008 standards, and some 243 features from the 2018 standard. It also implements several extensions 244 including OpenMP and OpenACC support for parallel programming. 245 246 The GNU Fortran compiler passes the 247 @uref{http://www.fortran-2000.com/ArnaudRecipes/fcvs21_f95.html, 248 NIST Fortran 77 Test Suite}, and produces acceptable results on the 249 @uref{https://www.netlib.org/lapack/faq.html#1.21, LAPACK Test Suite}. 250 It also provides respectable performance on 251 the @uref{https://polyhedron.com/?page_id=175, 252 Polyhedron Fortran compiler benchmarks} and the 253 @uref{https://www.netlib.org/benchmark/livermore, 254 Livermore Fortran Kernels test}. It has been used to compile a number of 255 large real-world programs, including 256 @uref{http://hirlam.org/, the HARMONIE and HIRLAM weather forecasting code} and 257 @uref{https://github.com/dylan-jayatilaka/tonto, 258 the Tonto quantum chemistry package}; see 259 @url{https://gcc.gnu.org/@/wiki/@/GfortranApps} for an extended list. 260 261 GNU Fortran provides the following functionality: 262 263 @itemize @bullet 264 @item 265 Read a program, stored in a file and containing @dfn{source code} 266 instructions written in Fortran 77. 267 268 @item 269 Translate the program into instructions a computer 270 can carry out more quickly than it takes to translate the 271 original Fortran instructions. 272 The result after compilation of a program is 273 @dfn{machine code}, 274 which is efficiently translated and processed 275 by a machine such as your computer. 276 Humans usually are not as good writing machine code 277 as they are at writing Fortran (or C++, Ada, or Java), 278 because it is easy to make tiny mistakes writing machine code. 279 280 @item 281 Provide information about the reasons why 282 the compiler may be unable to create a binary from the source code, 283 for example if the source code is flawed. 284 The Fortran language standards require that the compiler can point out 285 mistakes in your code. 286 An incorrect usage of the language causes an @dfn{error message}. 287 288 The compiler also attempts to diagnose cases where your 289 program contains a correct usage of the language, 290 but instructs the computer to do something questionable. 291 This kind of diagnostic message is called a @dfn{warning message}. 292 293 @item 294 Provide optional information about the translation passes 295 from the source code to machine code. 296 This can help you to find the cause of 297 certain bugs which may not be obvious in the source code, 298 but may be more easily found at a lower level compiler output. 299 It also helps developers to find bugs in the compiler itself. 300 301 @item 302 Provide information in the generated machine code that can 303 make it easier to find bugs in the program (using a debugging tool, 304 called a @dfn{debugger}, such as the GNU Debugger @command{gdb}). 305 306 @item 307 Locate and gather machine code already generated to 308 perform actions requested by statements in the program. 309 This machine code is organized into @dfn{modules} and is located 310 and @dfn{linked} to the user program. 311 @end itemize 312 313 The GNU Fortran compiler consists of several components: 314 315 @itemize @bullet 316 @item 317 A version of the @command{gcc} command 318 (which also might be installed as the system's @command{cc} command) 319 that also understands and accepts Fortran source code. 320 The @command{gcc} command is the @dfn{driver} program for 321 all the languages in the GNU Compiler Collection (GCC); 322 With @command{gcc}, 323 you can compile the source code of any language for 324 which a front end is available in GCC. 325 326 @item 327 The @command{gfortran} command itself, 328 which also might be installed as the 329 system's @command{f95} command. 330 @command{gfortran} is just another driver program, 331 but specifically for the Fortran compiler only. 332 The primary difference between the @command{gcc} and @command{gfortran} 333 commands is that the latter automatically links the correct libraries 334 to your program. 335 336 @item 337 A collection of run-time libraries. 338 These libraries contain the machine code needed to support 339 capabilities of the Fortran language that are not directly 340 provided by the machine code generated by the 341 @command{gfortran} compilation phase, 342 such as intrinsic functions and subroutines, 343 and routines for interaction with files and the operating system. 344 @c and mechanisms to spawn, 345 @c unleash and pause threads in parallelized code. 346 347 @item 348 The Fortran compiler itself, (@command{f951}). 349 This is the GNU Fortran parser and code generator, 350 linked to and interfaced with the GCC backend library. 351 @command{f951} ``translates'' the source code to 352 assembler code. You would typically not use this 353 program directly; 354 instead, the @command{gcc} or @command{gfortran} driver 355 programs call it for you. 356 @end itemize 357 358 359 @c --------------------------------------------------------------------- 360 @c GNU Fortran and GCC 361 @c --------------------------------------------------------------------- 362 363 @node GNU Fortran and GCC 364 @section GNU Fortran and GCC 365 @cindex GNU Compiler Collection 366 @cindex GCC 367 368 GNU Fortran is a part of GCC, the @dfn{GNU Compiler Collection}. GCC 369 consists of a collection of front ends for various languages, which 370 translate the source code into a language-independent form called 371 @dfn{GENERIC}. This is then processed by a common middle end which 372 provides optimization, and then passed to one of a collection of back 373 ends which generate code for different computer architectures and 374 operating systems. 375 376 Functionally, this is implemented with a driver program (@command{gcc}) 377 which provides the command-line interface for the compiler. It calls 378 the relevant compiler front-end program (e.g., @command{f951} for 379 Fortran) for each file in the source code, and then calls the assembler 380 and linker as appropriate to produce the compiled output. In a copy of 381 GCC that has been compiled with Fortran language support enabled, 382 @command{gcc} recognizes files with @file{.f}, @file{.for}, @file{.ftn}, 383 @file{.f90}, @file{.f95}, @file{.f03} and @file{.f08} extensions as 384 Fortran source code, and compiles it accordingly. A @command{gfortran} 385 driver program is also provided, which is identical to @command{gcc} 386 except that it automatically links the Fortran runtime libraries into the 387 compiled program. 388 389 Source files with @file{.f}, @file{.for}, @file{.fpp}, @file{.ftn}, @file{.F}, 390 @file{.FOR}, @file{.FPP}, and @file{.FTN} extensions are treated as fixed form. 391 Source files with @file{.f90}, @file{.f95}, @file{.f03}, @file{.f08}, 392 @file{.F90}, @file{.F95}, @file{.F03} and @file{.F08} extensions are 393 treated as free form. The capitalized versions of either form are run 394 through preprocessing. Source files with the lower case @file{.fpp} 395 extension are also run through preprocessing. 396 397 This manual specifically documents the Fortran front end, which handles 398 the programming language's syntax and semantics. The aspects of GCC 399 that relate to the optimization passes and the back-end code generation 400 are documented in the GCC manual; see 401 @ref{Top,,Introduction,gcc,Using the GNU Compiler Collection (GCC)}. 402 The two manuals together provide a complete reference for the GNU 403 Fortran compiler. 404 405 @c --------------------------------------------------------------------- 406 @c Standards 407 @c --------------------------------------------------------------------- 408 409 @node Standards 410 @section Standards 411 @cindex Standards 412 413 @menu 414 * Fortran 95 status:: 415 * Fortran 2003 status:: 416 * Fortran 2008 status:: 417 * Fortran 2018 status:: 418 @end menu 419 420 Fortran is developed by the Working Group 5 of Sub-Committee 22 of the 421 Joint Technical Committee 1 of the International Organization for 422 Standardization and the International Electrotechnical Commission (IEC). 423 This group is known as @uref{http://www.nag.co.uk/sc22wg5/, WG5}. 424 Official Fortran standard documents are available for purchase 425 from ISO; a collection of free documents (typically final drafts) are 426 also available on the @uref{https://gcc.gnu.org/wiki/GFortranStandards, wiki}. 427 428 The GNU Fortran compiler implements ISO/IEC 1539:1997 (Fortran 95). 429 As such, it can also compile essentially all standard-compliant 430 Fortran 90 and Fortran 77 programs. It also supports the ISO/IEC 431 TR-15581 enhancements to allocatable arrays. 432 433 GNU Fortran also supports almost all of ISO/IEC 1539-1:2004 434 (Fortran 2003) and ISO/IEC 1539-1:2010 (Fortran 2008). 435 It has partial support for features introduced in ISO/IEC 436 1539:2018 (Fortran 2018), the most recent version of the Fortran 437 language standard, including full support for the Technical Specification 438 @code{Further Interoperability of Fortran with C} (ISO/IEC TS 29113:2012). 439 More details on support for these standards can be 440 found in the following sections of the documentation. 441 442 Additionally, the GNU Fortran compilers supports the OpenMP specification 443 (version 4.5 and partial support of the features of the 5.0 version, 444 @url{https://openmp.org/@/openmp-specifications/}). 445 There also is support for the OpenACC specification (targeting 446 version 2.6, @uref{https://www.openacc.org/}). See 447 @uref{https://gcc.gnu.org/wiki/OpenACC} for more information. 448 449 @node Fortran 95 status 450 @subsection Fortran 95 status 451 @cindex Varying length strings 452 @cindex strings, varying length 453 @cindex conditional compilation 454 455 The Fortran 95 standard specifies in Part 2 (ISO/IEC 1539-2:2000) 456 varying length character strings. While GNU Fortran currently does not 457 support such strings directly, there exist two Fortran implementations 458 for them, which work with GNU Fortran. They can be found at 459 @uref{https://www.fortran.com/@/iso_varying_string.f95} and at 460 @uref{ftp://ftp.nag.co.uk/@/sc22wg5/@/ISO_VARYING_STRING/}. 461 462 Deferred-length character strings of Fortran 2003 supports part of 463 the features of @code{ISO_VARYING_STRING} and should be considered as 464 replacement. (Namely, allocatable or pointers of the type 465 @code{character(len=:)}.) 466 467 Part 3 of the Fortran 95 standard (ISO/IEC 1539-3:1998) defines 468 Conditional Compilation, which is not widely used and not directly 469 supported by the GNU Fortran compiler. You can use the program coco 470 to preprocess such files (@uref{http://www.daniellnagle.com/coco.html}). 471 472 @node Fortran 2003 status 473 @subsection Fortran 2003 status 474 475 GNU Fortran implements the Fortran 2003 (ISO/IEC 1539-1:2004) standard 476 except for finalization support, which is incomplete. 477 See the 478 @uref{https://gcc.gnu.org/wiki/Fortran2003, wiki page} for a full list 479 of new features introduced by Fortran 2003 and their implementation status. 480 481 @node Fortran 2008 status 482 @subsection Fortran 2008 status 483 484 The GNU Fortran compiler supports almost all features of Fortran 2008; 485 the @uref{https://gcc.gnu.org/wiki/Fortran2008Status, wiki} 486 has some information about the current implementation status. 487 In particular, the following are not yet supported: 488 489 @itemize @bullet 490 @item 491 @code{DO CONCURRENT} and @code{FORALL} do not recognize a 492 type-spec in the loop header. 493 494 @item 495 The change to permit any constant expression in subscripts and 496 nested implied-do limits in a @code{DATA} statement has not been implemented. 497 @end itemize 498 499 500 @node Fortran 2018 status 501 @subsection Fortran 2018 status 502 503 Fortran 2018 (ISO/IEC 1539:2018) is the most recent version 504 of the Fortran language standard. GNU Fortran implements some of the 505 new features of this standard: 506 507 @itemize @bullet 508 @item 509 All Fortran 2018 features derived from ISO/IEC TS 29113:2012, 510 ``Further Interoperability of Fortran with C'', are supported by GNU Fortran. 511 This includes assumed-type and assumed-rank objects and 512 the @code{SELECT RANK} construct as well as the parts relating to 513 @code{BIND(C)} functions. 514 See also @ref{Further Interoperability of Fortran with C}. 515 516 @item 517 GNU Fortran supports a subset of features derived from ISO/IEC TS 18508:2015, 518 ``Additional Parallel Features in Fortran'': 519 520 @itemize @bullet 521 @item 522 The new atomic ADD, CAS, FETCH and ADD/OR/XOR, OR and XOR intrinsics. 523 524 @item 525 The @code{CO_MIN} and @code{CO_MAX} and @code{SUM} reduction intrinsics, 526 and the @code{CO_BROADCAST} and @code{CO_REDUCE} intrinsic, except that those 527 do not support polymorphic types or types with allocatable, pointer or 528 polymorphic components. 529 530 @item 531 Events (@code{EVENT POST}, @code{EVENT WAIT}, @code{EVENT_QUERY}). 532 533 @item 534 Failed images (@code{FAIL IMAGE}, @code{IMAGE_STATUS}, 535 @code{FAILED_IMAGES}, @code{STOPPED_IMAGES}). 536 537 @end itemize 538 539 @item 540 An @code{ERROR STOP} statement is permitted in a @code{PURE} 541 procedure. 542 543 @item 544 GNU Fortran supports the @code{IMPLICIT NONE} statement with an 545 @code{implicit-none-spec-list}. 546 547 @item 548 The behavior of the @code{INQUIRE} statement with the @code{RECL=} 549 specifier now conforms to Fortran 2018. 550 551 @end itemize 552 553 554 @c ===================================================================== 555 @c PART I: INVOCATION REFERENCE 556 @c ===================================================================== 557 558 @tex 559 \part{I}{Invoking GNU Fortran} 560 @end tex 561 562 @c --------------------------------------------------------------------- 563 @c Compiler Options 564 @c --------------------------------------------------------------------- 565 566 @include invoke.texi 567 568 569 @c --------------------------------------------------------------------- 570 @c Runtime 571 @c --------------------------------------------------------------------- 572 573 @node Runtime 574 @chapter Runtime: Influencing runtime behavior with environment variables 575 @cindex environment variable 576 577 The behavior of the @command{gfortran} can be influenced by 578 environment variables. 579 580 Malformed environment variables are silently ignored. 581 582 @menu 583 * TMPDIR:: Directory for scratch files 584 * GFORTRAN_STDIN_UNIT:: Unit number for standard input 585 * GFORTRAN_STDOUT_UNIT:: Unit number for standard output 586 * GFORTRAN_STDERR_UNIT:: Unit number for standard error 587 * GFORTRAN_UNBUFFERED_ALL:: Do not buffer I/O for all units 588 * GFORTRAN_UNBUFFERED_PRECONNECTED:: Do not buffer I/O for preconnected units. 589 * GFORTRAN_SHOW_LOCUS:: Show location for runtime errors 590 * GFORTRAN_OPTIONAL_PLUS:: Print leading + where permitted 591 * GFORTRAN_LIST_SEPARATOR:: Separator for list output 592 * GFORTRAN_CONVERT_UNIT:: Set conversion for unformatted I/O 593 * GFORTRAN_ERROR_BACKTRACE:: Show backtrace on run-time errors 594 * GFORTRAN_FORMATTED_BUFFER_SIZE:: Buffer size for formatted files 595 * GFORTRAN_UNFORMATTED_BUFFER_SIZE:: Buffer size for unformatted files 596 @end menu 597 598 @node TMPDIR 599 @section @env{TMPDIR}---Directory for scratch files 600 601 When opening a file with @code{STATUS='SCRATCH'}, GNU Fortran tries to 602 create the file in one of the potential directories by testing each 603 directory in the order below. 604 605 @enumerate 606 @item 607 The environment variable @env{TMPDIR}, if it exists. 608 609 @item 610 On the MinGW target, the directory returned by the @code{GetTempPath} 611 function. Alternatively, on the Cygwin target, the @env{TMP} and 612 @env{TEMP} environment variables, if they exist, in that order. 613 614 @item 615 The @code{P_tmpdir} macro if it is defined, otherwise the directory 616 @file{/tmp}. 617 @end enumerate 618 619 @node GFORTRAN_STDIN_UNIT 620 @section @env{GFORTRAN_STDIN_UNIT}---Unit number for standard input 621 622 This environment variable can be used to select the unit number 623 preconnected to standard input. This must be a positive integer. 624 The default value is 5. 625 626 @node GFORTRAN_STDOUT_UNIT 627 @section @env{GFORTRAN_STDOUT_UNIT}---Unit number for standard output 628 629 This environment variable can be used to select the unit number 630 preconnected to standard output. This must be a positive integer. 631 The default value is 6. 632 633 @node GFORTRAN_STDERR_UNIT 634 @section @env{GFORTRAN_STDERR_UNIT}---Unit number for standard error 635 636 This environment variable can be used to select the unit number 637 preconnected to standard error. This must be a positive integer. 638 The default value is 0. 639 640 @node GFORTRAN_UNBUFFERED_ALL 641 @section @env{GFORTRAN_UNBUFFERED_ALL}---Do not buffer I/O on all units 642 643 This environment variable controls whether all I/O is unbuffered. If 644 the first letter is @samp{y}, @samp{Y} or @samp{1}, all I/O is 645 unbuffered. This will slow down small sequential reads and writes. If 646 the first letter is @samp{n}, @samp{N} or @samp{0}, I/O is buffered. 647 This is the default. 648 649 @node GFORTRAN_UNBUFFERED_PRECONNECTED 650 @section @env{GFORTRAN_UNBUFFERED_PRECONNECTED}---Do not buffer I/O on preconnected units 651 652 The environment variable named @env{GFORTRAN_UNBUFFERED_PRECONNECTED} controls 653 whether I/O on a preconnected unit (i.e.@: STDOUT or STDERR) is unbuffered. If 654 the first letter is @samp{y}, @samp{Y} or @samp{1}, I/O is unbuffered. This 655 will slow down small sequential reads and writes. If the first letter 656 is @samp{n}, @samp{N} or @samp{0}, I/O is buffered. This is the default. 657 658 @node GFORTRAN_SHOW_LOCUS 659 @section @env{GFORTRAN_SHOW_LOCUS}---Show location for runtime errors 660 661 If the first letter is @samp{y}, @samp{Y} or @samp{1}, filename and 662 line numbers for runtime errors are printed. If the first letter is 663 @samp{n}, @samp{N} or @samp{0}, do not print filename and line numbers 664 for runtime errors. The default is to print the location. 665 666 @node GFORTRAN_OPTIONAL_PLUS 667 @section @env{GFORTRAN_OPTIONAL_PLUS}---Print leading + where permitted 668 669 If the first letter is @samp{y}, @samp{Y} or @samp{1}, 670 a plus sign is printed 671 where permitted by the Fortran standard. If the first letter 672 is @samp{n}, @samp{N} or @samp{0}, a plus sign is not printed 673 in most cases. Default is not to print plus signs. 674 675 @node GFORTRAN_LIST_SEPARATOR 676 @section @env{GFORTRAN_LIST_SEPARATOR}---Separator for list output 677 678 This environment variable specifies the separator when writing 679 list-directed output. It may contain any number of spaces and 680 at most one comma. If you specify this on the command line, 681 be sure to quote spaces, as in 682 @smallexample 683 $ GFORTRAN_LIST_SEPARATOR=' , ' ./a.out 684 @end smallexample 685 when @command{a.out} is the compiled Fortran program that you want to run. 686 Default is a single space. 687 688 @node GFORTRAN_CONVERT_UNIT 689 @section @env{GFORTRAN_CONVERT_UNIT}---Set conversion for unformatted I/O 690 691 By setting the @env{GFORTRAN_CONVERT_UNIT} variable, it is possible 692 to change the representation of data for unformatted files. 693 The syntax for the @env{GFORTRAN_CONVERT_UNIT} variable for 694 most systems is: 695 @smallexample 696 GFORTRAN_CONVERT_UNIT: mode | mode ';' exception | exception ; 697 mode: 'native' | 'swap' | 'big_endian' | 'little_endian' ; 698 exception: mode ':' unit_list | unit_list ; 699 unit_list: unit_spec | unit_list unit_spec ; 700 unit_spec: INTEGER | INTEGER '-' INTEGER ; 701 @end smallexample 702 The variable consists of an optional default mode, followed by 703 a list of optional exceptions, which are separated by semicolons 704 from the preceding default and each other. Each exception consists 705 of a format and a comma-separated list of units. Valid values for 706 the modes are the same as for the @code{CONVERT} specifier: 707 708 @itemize @w{} 709 @item @code{NATIVE} Use the native format. This is the default. 710 @item @code{SWAP} Swap between little- and big-endian. 711 @item @code{LITTLE_ENDIAN} Use the little-endian format 712 for unformatted files. 713 @item @code{BIG_ENDIAN} Use the big-endian format for unformatted files. 714 @end itemize 715 For POWER systems which support @option{-mabi=ieeelongdouble}, 716 there are additional options, which can be combined with the 717 others with commas. Those are 718 @itemize @w{} 719 @item @code{R16_IEEE} Use IEEE 128-bit format for @code{REAL(KIND=16)}. 720 @item @code{R16_IBM} Use IBM @code{long double} format for 721 @code{REAL(KIND=16)}. 722 @end itemize 723 A missing mode for an exception is taken to mean @code{BIG_ENDIAN}. 724 Examples of values for @env{GFORTRAN_CONVERT_UNIT} are: 725 @itemize @w{} 726 @item @code{'big_endian'} Do all unformatted I/O in big_endian mode. 727 @item @code{'little_endian;native:10-20,25'} Do all unformatted I/O 728 in little_endian mode, except for units 10 to 20 and 25, which are in 729 native format. 730 @item @code{'10-20'} Units 10 to 20 are big-endian, the rest is native. 731 @item @code{'big_endian,r16_ibm'} Do all unformatted I/O in big-endian 732 mode and use IBM long double for output of @code{REAL(KIND=16)} values. 733 @end itemize 734 735 Setting the environment variables should be done on the command 736 line or via the @command{export} 737 command for @command{sh}-compatible shells and via @command{setenv} 738 for @command{csh}-compatible shells. 739 740 Example for @command{sh}: 741 @smallexample 742 $ gfortran foo.f90 743 $ GFORTRAN_CONVERT_UNIT='big_endian;native:10-20' ./a.out 744 @end smallexample 745 746 Example code for @command{csh}: 747 @smallexample 748 % gfortran foo.f90 749 % setenv GFORTRAN_CONVERT_UNIT 'big_endian;native:10-20' 750 % ./a.out 751 @end smallexample 752 753 Using anything but the native representation for unformatted data 754 carries a significant speed overhead. If speed in this area matters 755 to you, it is best if you use this only for data that needs to be 756 portable. 757 758 @xref{CONVERT specifier}, for an alternative way to specify the 759 data representation for unformatted files. @xref{Runtime Options}, for 760 setting a default data representation for the whole program. The 761 @code{CONVERT} specifier overrides the @option{-fconvert} compile options. 762 763 @emph{Note that the values specified via the GFORTRAN_CONVERT_UNIT 764 environment variable will override the CONVERT specifier in the 765 open statement}. This is to give control over data formats to 766 users who do not have the source code of their program available. 767 768 @node GFORTRAN_ERROR_BACKTRACE 769 @section @env{GFORTRAN_ERROR_BACKTRACE}---Show backtrace on run-time errors 770 771 If the @env{GFORTRAN_ERROR_BACKTRACE} variable is set to @samp{y}, 772 @samp{Y} or @samp{1} (only the first letter is relevant) then a 773 backtrace is printed when a serious run-time error occurs. To disable 774 the backtracing, set the variable to @samp{n}, @samp{N}, @samp{0}. 775 Default is to print a backtrace unless the @option{-fno-backtrace} 776 compile option was used. 777 778 @node GFORTRAN_FORMATTED_BUFFER_SIZE 779 @section @env{GFORTRAN_FORMATTED_BUFFER_SIZE}---Set buffer size for formatted I/O 780 781 The @env{GFORTRAN_FORMATTED_BUFFER_SIZE} environment variable 782 specifies buffer size in bytes to be used for formatted output. 783 The default value is 8192. 784 785 @node GFORTRAN_UNFORMATTED_BUFFER_SIZE 786 @section @env{GFORTRAN_UNFORMATTED_BUFFER_SIZE}---Set buffer size for unformatted I/O 787 788 The @env{GFORTRAN_UNFORMATTED_BUFFER_SIZE} environment variable 789 specifies buffer size in bytes to be used for unformatted output. 790 The default value is 131072. 791 792 @c ===================================================================== 793 @c PART II: LANGUAGE REFERENCE 794 @c ===================================================================== 795 796 @tex 797 \part{II}{Language Reference} 798 @end tex 799 800 801 802 @c --------------------------------------------------------------------- 803 @c Compiler Characteristics 804 @c --------------------------------------------------------------------- 805 806 @node Compiler Characteristics 807 @chapter Compiler Characteristics 808 809 This chapter describes certain characteristics of the GNU Fortran 810 compiler, that are not specified by the Fortran standard, but which 811 might in some way or another become visible to the programmer. 812 813 @menu 814 * KIND Type Parameters:: 815 * Internal representation of LOGICAL variables:: 816 * Evaluation of logical expressions:: 817 * MAX and MIN intrinsics with REAL NaN arguments:: 818 * Thread-safety of the runtime library:: 819 * Data consistency and durability:: 820 * Files opened without an explicit ACTION= specifier:: 821 * File operations on symbolic links:: 822 * File format of unformatted sequential files:: 823 * Asynchronous I/O:: 824 @end menu 825 826 827 @node KIND Type Parameters 828 @section KIND Type Parameters 829 @cindex kind 830 831 The @code{KIND} type parameters supported by GNU Fortran for the primitive 832 data types are: 833 834 @table @code 835 836 @item INTEGER 837 1, 2, 4, 8*, 16*, default: 4** 838 839 @item LOGICAL 840 1, 2, 4, 8*, 16*, default: 4** 841 842 @item REAL 843 4, 8, 10*, 16*, default: 4*** 844 845 @item COMPLEX 846 4, 8, 10*, 16*, default: 4*** 847 848 @item DOUBLE PRECISION 849 4, 8, 10*, 16*, default: 8*** 850 851 @item CHARACTER 852 1, 4, default: 1 853 854 @end table 855 856 @noindent 857 * not available on all systems @* 858 ** unless @option{-fdefault-integer-8} is used @* 859 *** unless @option{-fdefault-real-8} is used (see @ref{Fortran Dialect Options}) 860 861 @noindent 862 The @code{KIND} value matches the storage size in bytes, except for 863 @code{COMPLEX} where the storage size is twice as much (or both real and 864 imaginary part are a real value of the given size). It is recommended to use 865 the @ref{SELECTED_CHAR_KIND}, @ref{SELECTED_INT_KIND} and 866 @ref{SELECTED_REAL_KIND} intrinsics or the @code{INT8}, @code{INT16}, 867 @code{INT32}, @code{INT64}, @code{REAL32}, @code{REAL64}, and @code{REAL128} 868 parameters of the @code{ISO_FORTRAN_ENV} module instead of the concrete values. 869 The available kind parameters can be found in the constant arrays 870 @code{CHARACTER_KINDS}, @code{INTEGER_KINDS}, @code{LOGICAL_KINDS} and 871 @code{REAL_KINDS} in the @ref{ISO_FORTRAN_ENV} module. For C interoperability, 872 the kind parameters of the @ref{ISO_C_BINDING} module should be used. 873 874 875 @node Internal representation of LOGICAL variables 876 @section Internal representation of LOGICAL variables 877 @cindex logical, variable representation 878 879 The Fortran standard does not specify how variables of @code{LOGICAL} 880 type are represented, beyond requiring that @code{LOGICAL} variables 881 of default kind have the same storage size as default @code{INTEGER} 882 and @code{REAL} variables. The GNU Fortran internal representation is 883 as follows. 884 885 A @code{LOGICAL(KIND=N)} variable is represented as an 886 @code{INTEGER(KIND=N)} variable, however, with only two permissible 887 values: @code{1} for @code{.TRUE.} and @code{0} for 888 @code{.FALSE.}. Any other integer value results in undefined behavior. 889 890 See also @ref{Argument passing conventions} and @ref{Interoperability with C}. 891 892 893 @node Evaluation of logical expressions 894 @section Evaluation of logical expressions 895 896 The Fortran standard does not require the compiler to evaluate all parts of an 897 expression, if they do not contribute to the final result. For logical 898 expressions with @code{.AND.} or @code{.OR.} operators, in particular, GNU 899 Fortran will optimize out function calls (even to impure functions) if the 900 result of the expression can be established without them. However, since not 901 all compilers do that, and such an optimization can potentially modify the 902 program flow and subsequent results, GNU Fortran throws warnings for such 903 situations with the @option{-Wfunction-elimination} flag. 904 905 906 @node MAX and MIN intrinsics with REAL NaN arguments 907 @section MAX and MIN intrinsics with REAL NaN arguments 908 @cindex MAX, MIN, NaN 909 910 The Fortran standard does not specify what the result of the 911 @code{MAX} and @code{MIN} intrinsics are if one of the arguments is a 912 @code{NaN}. Accordingly, the GNU Fortran compiler does not specify 913 that either, as this allows for faster and more compact code to be 914 generated. If the programmer wishes to take some specific action in 915 case one of the arguments is a @code{NaN}, it is necessary to 916 explicitly test the arguments before calling @code{MAX} or @code{MIN}, 917 e.g. with the @code{IEEE_IS_NAN} function from the intrinsic module 918 @code{IEEE_ARITHMETIC}. 919 920 921 @node Thread-safety of the runtime library 922 @section Thread-safety of the runtime library 923 @cindex thread-safety, threads 924 925 GNU Fortran can be used in programs with multiple threads, e.g.@: by 926 using OpenMP, by calling OS thread handling functions via the 927 @code{ISO_C_BINDING} facility, or by GNU Fortran compiled library code 928 being called from a multi-threaded program. 929 930 The GNU Fortran runtime library, (@code{libgfortran}), supports being 931 called concurrently from multiple threads with the following 932 exceptions. 933 934 During library initialization, the C @code{getenv} function is used, 935 which need not be thread-safe. Similarly, the @code{getenv} 936 function is used to implement the @code{GET_ENVIRONMENT_VARIABLE} and 937 @code{GETENV} intrinsics. It is the responsibility of the user to 938 ensure that the environment is not being updated concurrently when any 939 of these actions are taking place. 940 941 The @code{EXECUTE_COMMAND_LINE} and @code{SYSTEM} intrinsics are 942 implemented with the @code{system} function, which need not be 943 thread-safe. It is the responsibility of the user to ensure that 944 @code{system} is not called concurrently. 945 946 For platforms not supporting thread-safe POSIX functions, further 947 functionality might not be thread-safe. For details, please consult 948 the documentation for your operating system. 949 950 The GNU Fortran runtime library uses various C library functions that 951 depend on the locale, such as @code{strtod} and @code{snprintf}. In 952 order to work correctly in locale-aware programs that set the locale 953 using @code{setlocale}, the locale is reset to the default ``C'' 954 locale while executing a formatted @code{READ} or @code{WRITE} 955 statement. On targets supporting the POSIX 2008 per-thread locale 956 functions (e.g. @code{newlocale}, @code{uselocale}, 957 @code{freelocale}), these are used and thus the global locale set 958 using @code{setlocale} or the per-thread locales in other threads are 959 not affected. However, on targets lacking this functionality, the 960 global LC_NUMERIC locale is set to ``C'' during the formatted I/O. 961 Thus, on such targets it's not safe to call @code{setlocale} 962 concurrently from another thread while a Fortran formatted I/O 963 operation is in progress. Also, other threads doing something 964 dependent on the LC_NUMERIC locale might not work correctly if a 965 formatted I/O operation is in progress in another thread. 966 967 @node Data consistency and durability 968 @section Data consistency and durability 969 @cindex consistency, durability 970 971 This section contains a brief overview of data and metadata 972 consistency and durability issues when doing I/O. 973 974 With respect to durability, GNU Fortran makes no effort to ensure that 975 data is committed to stable storage. If this is required, the GNU 976 Fortran programmer can use the intrinsic @code{FNUM} to retrieve the 977 low level file descriptor corresponding to an open Fortran unit. Then, 978 using e.g. the @code{ISO_C_BINDING} feature, one can call the 979 underlying system call to flush dirty data to stable storage, such as 980 @code{fsync} on POSIX, @code{_commit} on MingW, or @code{fcntl(fd, 981 F_FULLSYNC, 0)} on Mac OS X. The following example shows how to call 982 fsync: 983 984 @smallexample 985 ! Declare the interface for POSIX fsync function 986 interface 987 function fsync (fd) bind(c,name="fsync") 988 use iso_c_binding, only: c_int 989 integer(c_int), value :: fd 990 integer(c_int) :: fsync 991 end function fsync 992 end interface 993 994 ! Variable declaration 995 integer :: ret 996 997 ! Opening unit 10 998 open (10,file="foo") 999 1000 ! ... 1001 ! Perform I/O on unit 10 1002 ! ... 1003 1004 ! Flush and sync 1005 flush(10) 1006 ret = fsync(fnum(10)) 1007 1008 ! Handle possible error 1009 if (ret /= 0) stop "Error calling FSYNC" 1010 @end smallexample 1011 1012 With respect to consistency, for regular files GNU Fortran uses 1013 buffered I/O in order to improve performance. This buffer is flushed 1014 automatically when full and in some other situations, e.g. when 1015 closing a unit. It can also be explicitly flushed with the 1016 @code{FLUSH} statement. Also, the buffering can be turned off with the 1017 @code{GFORTRAN_UNBUFFERED_ALL} and 1018 @code{GFORTRAN_UNBUFFERED_PRECONNECTED} environment variables. Special 1019 files, such as terminals and pipes, are always unbuffered. Sometimes, 1020 however, further things may need to be done in order to allow other 1021 processes to see data that GNU Fortran has written, as follows. 1022 1023 The Windows platform supports a relaxed metadata consistency model, 1024 where file metadata is written to the directory lazily. This means 1025 that, for instance, the @code{dir} command can show a stale size for a 1026 file. One can force a directory metadata update by closing the unit, 1027 or by calling @code{_commit} on the file descriptor. Note, though, 1028 that @code{_commit} will force all dirty data to stable storage, which 1029 is often a very slow operation. 1030 1031 The Network File System (NFS) implements a relaxed consistency model 1032 called open-to-close consistency. Closing a file forces dirty data and 1033 metadata to be flushed to the server, and opening a file forces the 1034 client to contact the server in order to revalidate cached 1035 data. @code{fsync} will also force a flush of dirty data and metadata 1036 to the server. Similar to @code{open} and @code{close}, acquiring and 1037 releasing @code{fcntl} file locks, if the server supports them, will 1038 also force cache validation and flushing dirty data and metadata. 1039 1040 1041 @node Files opened without an explicit ACTION= specifier 1042 @section Files opened without an explicit ACTION= specifier 1043 @cindex open, action 1044 1045 The Fortran standard says that if an @code{OPEN} statement is executed 1046 without an explicit @code{ACTION=} specifier, the default value is 1047 processor dependent. GNU Fortran behaves as follows: 1048 1049 @enumerate 1050 @item Attempt to open the file with @code{ACTION='READWRITE'} 1051 @item If that fails, try to open with @code{ACTION='READ'} 1052 @item If that fails, try to open with @code{ACTION='WRITE'} 1053 @item If that fails, generate an error 1054 @end enumerate 1055 1056 1057 @node File operations on symbolic links 1058 @section File operations on symbolic links 1059 @cindex file, symbolic link 1060 1061 This section documents the behavior of GNU Fortran for file operations on 1062 symbolic links, on systems that support them. 1063 1064 @itemize 1065 1066 @item Results of INQUIRE statements of the ``inquire by file'' form will 1067 relate to the target of the symbolic link. For example, 1068 @code{INQUIRE(FILE="foo",EXIST=ex)} will set @var{ex} to @var{.true.} if 1069 @var{foo} is a symbolic link pointing to an existing file, and @var{.false.} 1070 if @var{foo} points to an non-existing file (``dangling'' symbolic link). 1071 1072 @item Using the @code{OPEN} statement with a @code{STATUS="NEW"} specifier 1073 on a symbolic link will result in an error condition, whether the symbolic 1074 link points to an existing target or is dangling. 1075 1076 @item If a symbolic link was connected, using the @code{CLOSE} statement 1077 with a @code{STATUS="DELETE"} specifier will cause the symbolic link itself 1078 to be deleted, not its target. 1079 1080 @end itemize 1081 1082 @node File format of unformatted sequential files 1083 @section File format of unformatted sequential files 1084 @cindex file, unformatted sequential 1085 @cindex unformatted sequential 1086 @cindex sequential, unformatted 1087 @cindex record marker 1088 @cindex subrecord 1089 1090 Unformatted sequential files are stored as logical records using 1091 record markers. Each logical record consists of one of more 1092 subrecords. 1093 1094 Each subrecord consists of a leading record marker, the data written 1095 by the user program, and a trailing record marker. The record markers 1096 are four-byte integers by default, and eight-byte integers if the 1097 @option{-fmax-subrecord-length=8} option (which exists for backwards 1098 compability only) is in effect. 1099 1100 The representation of the record markers is that of unformatted files 1101 given with the @option{-fconvert} option, the @ref{CONVERT specifier} 1102 in an open statement or the @ref{GFORTRAN_CONVERT_UNIT} environment 1103 variable. 1104 1105 The maximum number of bytes of user data in a subrecord is 2147483639 1106 (2 GiB - 9) for a four-byte record marker. This limit can be lowered 1107 with the @option{-fmax-subrecord-length} option, although this is 1108 rarely useful. If the length of a logical record exceeds this limit, 1109 the data is distributed among several subrecords. 1110 1111 The absolute of the number stored in the record markers is the number 1112 of bytes of user data in the corresponding subrecord. If the leading 1113 record marker of a subrecord contains a negative number, another 1114 subrecord follows the current one. If the trailing record marker 1115 contains a negative number, then there is a preceding subrecord. 1116 1117 In the most simple case, with only one subrecord per logical record, 1118 both record markers contain the number of bytes of user data in the 1119 record. 1120 1121 The format for unformatted sequential data can be duplicated using 1122 unformatted stream, as shown in the example program for an unformatted 1123 record containing a single subrecord: 1124 1125 @smallexample 1126 program main 1127 use iso_fortran_env, only: int32 1128 implicit none 1129 integer(int32) :: i 1130 real, dimension(10) :: a, b 1131 call random_number(a) 1132 open (10,file='test.dat',form='unformatted',access='stream') 1133 inquire (iolength=i) a 1134 write (10) i, a, i 1135 close (10) 1136 open (10,file='test.dat',form='unformatted') 1137 read (10) b 1138 if (all (a == b)) print *,'success!' 1139 end program main 1140 @end smallexample 1141 1142 @node Asynchronous I/O 1143 @section Asynchronous I/O 1144 @cindex input/output, asynchronous 1145 @cindex asynchronous I/O 1146 1147 Asynchronous I/O is supported if the program is linked against the 1148 POSIX thread library. If that is not the case, all I/O is performed 1149 as synchronous. On systems which do not support pthread condition 1150 variables, such as AIX, I/O is also performed as synchronous. 1151 1152 On some systems, such as Darwin or Solaris, the POSIX thread library 1153 is always linked in, so asynchronous I/O is always performed. On other 1154 sytems, such as Linux, it is necessary to specify @option{-pthread}, 1155 @option{-lpthread} or @option{-fopenmp} during the linking step. 1156 1157 @c --------------------------------------------------------------------- 1158 @c Extensions 1159 @c --------------------------------------------------------------------- 1160 1161 @c Maybe this chapter should be merged with the 'Standards' section, 1162 @c whenever that is written :-) 1163 1164 @node Extensions 1165 @chapter Extensions 1166 @cindex extensions 1167 1168 The two sections below detail the extensions to standard Fortran that are 1169 implemented in GNU Fortran, as well as some of the popular or 1170 historically important extensions that are not (or not yet) implemented. 1171 For the latter case, we explain the alternatives available to GNU Fortran 1172 users, including replacement by standard-conforming code or GNU 1173 extensions. 1174 1175 @menu 1176 * Extensions implemented in GNU Fortran:: 1177 * Extensions not implemented in GNU Fortran:: 1178 @end menu 1179 1180 1181 @node Extensions implemented in GNU Fortran 1182 @section Extensions implemented in GNU Fortran 1183 @cindex extensions, implemented 1184 1185 GNU Fortran implements a number of extensions over standard Fortran. 1186 This chapter contains information on their syntax and meaning. There 1187 are currently two categories of GNU Fortran extensions, those that 1188 provide functionality beyond that provided by any standard, and those 1189 that are supported by GNU Fortran purely for backward compatibility 1190 with legacy compilers. By default, @option{-std=gnu} allows the 1191 compiler to accept both types of extensions, but to warn about the use 1192 of the latter. Specifying either @option{-std=f95}, 1193 @option{-std=f2003}, @option{-std=f2008}, or @option{-std=f2018} 1194 disables both types of extensions, and @option{-std=legacy} allows 1195 both without warning. The special compile flag @option{-fdec} enables 1196 additional compatibility extensions along with those enabled by 1197 @option{-std=legacy}. 1198 1199 @menu 1200 * Old-style kind specifications:: 1201 * Old-style variable initialization:: 1202 * Extensions to namelist:: 1203 * X format descriptor without count field:: 1204 * Commas in FORMAT specifications:: 1205 * Missing period in FORMAT specifications:: 1206 * Default widths for F@comma{} G and I format descriptors:: 1207 * I/O item lists:: 1208 * @code{Q} exponent-letter:: 1209 * BOZ literal constants:: 1210 * Real array indices:: 1211 * Unary operators:: 1212 * Implicitly convert LOGICAL and INTEGER values:: 1213 * Hollerith constants support:: 1214 * Character conversion:: 1215 * Cray pointers:: 1216 * CONVERT specifier:: 1217 * OpenMP:: 1218 * OpenACC:: 1219 * Argument list functions:: 1220 * Read/Write after EOF marker:: 1221 * STRUCTURE and RECORD:: 1222 * UNION and MAP:: 1223 * Type variants for integer intrinsics:: 1224 * AUTOMATIC and STATIC attributes:: 1225 * Extended math intrinsics:: 1226 * Form feed as whitespace:: 1227 * TYPE as an alias for PRINT:: 1228 * %LOC as an rvalue:: 1229 * .XOR. operator:: 1230 * Bitwise logical operators:: 1231 * Extended I/O specifiers:: 1232 * Legacy PARAMETER statements:: 1233 * Default exponents:: 1234 @end menu 1235 1236 @node Old-style kind specifications 1237 @subsection Old-style kind specifications 1238 @cindex kind, old-style 1239 1240 GNU Fortran allows old-style kind specifications in declarations. These 1241 look like: 1242 @smallexample 1243 TYPESPEC*size x,y,z 1244 @end smallexample 1245 @noindent 1246 where @code{TYPESPEC} is a basic type (@code{INTEGER}, @code{REAL}, 1247 etc.), and where @code{size} is a byte count corresponding to the 1248 storage size of a valid kind for that type. (For @code{COMPLEX} 1249 variables, @code{size} is the total size of the real and imaginary 1250 parts.) The statement then declares @code{x}, @code{y} and @code{z} to 1251 be of type @code{TYPESPEC} with the appropriate kind. This is 1252 equivalent to the standard-conforming declaration 1253 @smallexample 1254 TYPESPEC(k) x,y,z 1255 @end smallexample 1256 @noindent 1257 where @code{k} is the kind parameter suitable for the intended precision. As 1258 kind parameters are implementation-dependent, use the @code{KIND}, 1259 @code{SELECTED_INT_KIND} and @code{SELECTED_REAL_KIND} intrinsics to retrieve 1260 the correct value, for instance @code{REAL*8 x} can be replaced by: 1261 @smallexample 1262 INTEGER, PARAMETER :: dbl = KIND(1.0d0) 1263 REAL(KIND=dbl) :: x 1264 @end smallexample 1265 1266 @node Old-style variable initialization 1267 @subsection Old-style variable initialization 1268 1269 GNU Fortran allows old-style initialization of variables of the 1270 form: 1271 @smallexample 1272 INTEGER i/1/,j/2/ 1273 REAL x(2,2) /3*0.,1./ 1274 @end smallexample 1275 The syntax for the initializers is as for the @code{DATA} statement, but 1276 unlike in a @code{DATA} statement, an initializer only applies to the 1277 variable immediately preceding the initialization. In other words, 1278 something like @code{INTEGER I,J/2,3/} is not valid. This style of 1279 initialization is only allowed in declarations without double colons 1280 (@code{::}); the double colons were introduced in Fortran 90, which also 1281 introduced a standard syntax for initializing variables in type 1282 declarations. 1283 1284 Examples of standard-conforming code equivalent to the above example 1285 are: 1286 @smallexample 1287 ! Fortran 90 1288 INTEGER :: i = 1, j = 2 1289 REAL :: x(2,2) = RESHAPE((/0.,0.,0.,1./),SHAPE(x)) 1290 ! Fortran 77 1291 INTEGER i, j 1292 REAL x(2,2) 1293 DATA i/1/, j/2/, x/3*0.,1./ 1294 @end smallexample 1295 1296 Note that variables which are explicitly initialized in declarations 1297 or in @code{DATA} statements automatically acquire the @code{SAVE} 1298 attribute. 1299 1300 @node Extensions to namelist 1301 @subsection Extensions to namelist 1302 @cindex Namelist 1303 1304 GNU Fortran fully supports the Fortran 95 standard for namelist I/O 1305 including array qualifiers, substrings and fully qualified derived types. 1306 The output from a namelist write is compatible with namelist read. The 1307 output has all names in upper case and indentation to column 1 after the 1308 namelist name. Two extensions are permitted: 1309 1310 Old-style use of @samp{$} instead of @samp{&} 1311 @smallexample 1312 $MYNML 1313 X(:)%Y(2) = 1.0 2.0 3.0 1314 CH(1:4) = "abcd" 1315 $END 1316 @end smallexample 1317 1318 It should be noted that the default terminator is @samp{/} rather than 1319 @samp{&END}. 1320 1321 Querying of the namelist when inputting from stdin. After at least 1322 one space, entering @samp{?} sends to stdout the namelist name and the names of 1323 the variables in the namelist: 1324 @smallexample 1325 ? 1326 1327 &mynml 1328 x 1329 x%y 1330 ch 1331 &end 1332 @end smallexample 1333 1334 Entering @samp{=?} outputs the namelist to stdout, as if 1335 @code{WRITE(*,NML = mynml)} had been called: 1336 @smallexample 1337 =? 1338 1339 &MYNML 1340 X(1)%Y= 0.000000 , 1.000000 , 0.000000 , 1341 X(2)%Y= 0.000000 , 2.000000 , 0.000000 , 1342 X(3)%Y= 0.000000 , 3.000000 , 0.000000 , 1343 CH=abcd, / 1344 @end smallexample 1345 1346 To aid this dialog, when input is from stdin, errors send their 1347 messages to stderr and execution continues, even if @code{IOSTAT} is set. 1348 1349 @code{PRINT} namelist is permitted. This causes an error if 1350 @option{-std=f95} is used. 1351 @smallexample 1352 PROGRAM test_print 1353 REAL, dimension (4) :: x = (/1.0, 2.0, 3.0, 4.0/) 1354 NAMELIST /mynml/ x 1355 PRINT mynml 1356 END PROGRAM test_print 1357 @end smallexample 1358 1359 Expanded namelist reads are permitted. This causes an error if 1360 @option{-std=f95} is used. In the following example, the first element 1361 of the array will be given the value 0.00 and the two succeeding 1362 elements will be given the values 1.00 and 2.00. 1363 @smallexample 1364 &MYNML 1365 X(1,1) = 0.00 , 1.00 , 2.00 1366 / 1367 @end smallexample 1368 1369 When writing a namelist, if no @code{DELIM=} is specified, by default a 1370 double quote is used to delimit character strings. If -std=F95, F2003, 1371 or F2008, etc, the delim status is set to 'none'. Defaulting to 1372 quotes ensures that namelists with character strings can be subsequently 1373 read back in accurately. 1374 1375 @node X format descriptor without count field 1376 @subsection @code{X} format descriptor without count field 1377 1378 To support legacy codes, GNU Fortran permits the count field of the 1379 @code{X} edit descriptor in @code{FORMAT} statements to be omitted. 1380 When omitted, the count is implicitly assumed to be one. 1381 1382 @smallexample 1383 PRINT 10, 2, 3 1384 10 FORMAT (I1, X, I1) 1385 @end smallexample 1386 1387 @node Commas in FORMAT specifications 1388 @subsection Commas in @code{FORMAT} specifications 1389 1390 To support legacy codes, GNU Fortran allows the comma separator 1391 to be omitted immediately before and after character string edit 1392 descriptors in @code{FORMAT} statements. A comma with no following format 1393 decriptor is permited if the @option{-fdec-blank-format-item} is given on 1394 the command line. This is considered non-conforming code and is 1395 discouraged. 1396 1397 @smallexample 1398 PRINT 10, 2, 3 1399 10 FORMAT ('FOO='I1' BAR='I2) 1400 print 20, 5, 6 1401 20 FORMAT (I3, I3,) 1402 @end smallexample 1403 1404 1405 @node Missing period in FORMAT specifications 1406 @subsection Missing period in @code{FORMAT} specifications 1407 1408 To support legacy codes, GNU Fortran allows missing periods in format 1409 specifications if and only if @option{-std=legacy} is given on the 1410 command line. This is considered non-conforming code and is 1411 discouraged. 1412 1413 @smallexample 1414 REAL :: value 1415 READ(*,10) value 1416 10 FORMAT ('F4') 1417 @end smallexample 1418 1419 @node Default widths for F@comma{} G and I format descriptors 1420 @subsection Default widths for @code{F}, @code{G} and @code{I} format descriptors 1421 1422 To support legacy codes, GNU Fortran allows width to be omitted from format 1423 specifications if and only if @option{-fdec-format-defaults} is given on the 1424 command line. Default widths will be used. This is considered non-conforming 1425 code and is discouraged. 1426 1427 @smallexample 1428 REAL :: value1 1429 INTEGER :: value2 1430 WRITE(*,10) value1, value1, value2 1431 10 FORMAT ('F, G, I') 1432 @end smallexample 1433 1434 1435 @node I/O item lists 1436 @subsection I/O item lists 1437 @cindex I/O item lists 1438 1439 To support legacy codes, GNU Fortran allows the input item list 1440 of the @code{READ} statement, and the output item lists of the 1441 @code{WRITE} and @code{PRINT} statements, to start with a comma. 1442 1443 @node @code{Q} exponent-letter 1444 @subsection @code{Q} exponent-letter 1445 @cindex @code{Q} exponent-letter 1446 1447 GNU Fortran accepts real literal constants with an exponent-letter 1448 of @code{Q}, for example, @code{1.23Q45}. The constant is interpreted 1449 as a @code{REAL(16)} entity on targets that support this type. If 1450 the target does not support @code{REAL(16)} but has a @code{REAL(10)} 1451 type, then the real-literal-constant will be interpreted as a 1452 @code{REAL(10)} entity. In the absence of @code{REAL(16)} and 1453 @code{REAL(10)}, an error will occur. 1454 1455 @node BOZ literal constants 1456 @subsection BOZ literal constants 1457 @cindex BOZ literal constants 1458 1459 Besides decimal constants, Fortran also supports binary (@code{b}), 1460 octal (@code{o}) and hexadecimal (@code{z}) integer constants. The 1461 syntax is: @samp{prefix quote digits quote}, where the prefix is 1462 either @code{b}, @code{o} or @code{z}, quote is either @code{'} or 1463 @code{"} and the digits are @code{0} or @code{1} for binary, 1464 between @code{0} and @code{7} for octal, and between @code{0} and 1465 @code{F} for hexadecimal. (Example: @code{b'01011101'}.) 1466 1467 Up to Fortran 95, BOZ literal constants were only allowed to initialize 1468 integer variables in DATA statements. Since Fortran 2003 BOZ literal 1469 constants are also allowed as actual arguments to the @code{REAL}, 1470 @code{DBLE}, @code{INT} and @code{CMPLX} intrinsic functions. 1471 The BOZ literal constant is simply a string of bits, which is padded 1472 or truncated as needed, during conversion to a numeric type. The 1473 Fortran standard states that the treatment of the sign bit is processor 1474 dependent. Gfortran interprets the sign bit as a user would expect. 1475 1476 As a deprecated extension, GNU Fortran allows hexadecimal BOZ literal 1477 constants to be specified using the @code{X} prefix. That the BOZ literal 1478 constant can also be specified by adding a suffix to the string, for 1479 example, @code{Z'ABC'} and @code{'ABC'X} are equivalent. Additionally, 1480 as extension, BOZ literals are permitted in some contexts outside of 1481 @code{DATA} and the intrinsic functions listed in the Fortran standard. 1482 Use @option{-fallow-invalid-boz} to enable the extension. 1483 1484 @node Real array indices 1485 @subsection Real array indices 1486 @cindex array, indices of type real 1487 1488 As an extension, GNU Fortran allows the use of @code{REAL} expressions 1489 or variables as array indices. 1490 1491 @node Unary operators 1492 @subsection Unary operators 1493 @cindex operators, unary 1494 1495 As an extension, GNU Fortran allows unary plus and unary minus operators 1496 to appear as the second operand of binary arithmetic operators without 1497 the need for parenthesis. 1498 1499 @smallexample 1500 X = Y * -Z 1501 @end smallexample 1502 1503 @node Implicitly convert LOGICAL and INTEGER values 1504 @subsection Implicitly convert @code{LOGICAL} and @code{INTEGER} values 1505 @cindex conversion, to integer 1506 @cindex conversion, to logical 1507 1508 As an extension for backwards compatibility with other compilers, GNU 1509 Fortran allows the implicit conversion of @code{LOGICAL} values to 1510 @code{INTEGER} values and vice versa. When converting from a 1511 @code{LOGICAL} to an @code{INTEGER}, @code{.FALSE.} is interpreted as 1512 zero, and @code{.TRUE.} is interpreted as one. When converting from 1513 @code{INTEGER} to @code{LOGICAL}, the value zero is interpreted as 1514 @code{.FALSE.} and any nonzero value is interpreted as @code{.TRUE.}. 1515 1516 @smallexample 1517 LOGICAL :: l 1518 l = 1 1519 @end smallexample 1520 @smallexample 1521 INTEGER :: i 1522 i = .TRUE. 1523 @end smallexample 1524 1525 However, there is no implicit conversion of @code{INTEGER} values in 1526 @code{if}-statements, nor of @code{LOGICAL} or @code{INTEGER} values 1527 in I/O operations. 1528 1529 @node Hollerith constants support 1530 @subsection Hollerith constants support 1531 @cindex Hollerith constants 1532 1533 GNU Fortran supports Hollerith constants in assignments, @code{DATA} 1534 statements, function and subroutine arguments. A Hollerith constant is 1535 written as a string of characters preceded by an integer constant 1536 indicating the character count, and the letter @code{H} or 1537 @code{h}, and stored in bytewise fashion in a numeric (@code{INTEGER}, 1538 @code{REAL}, or @code{COMPLEX}), @code{LOGICAL} or @code{CHARACTER} variable. 1539 The constant will be padded with spaces or truncated to fit the size of 1540 the variable in which it is stored. 1541 1542 Examples of valid uses of Hollerith constants: 1543 @smallexample 1544 complex*16 x(2) 1545 data x /16Habcdefghijklmnop, 16Hqrstuvwxyz012345/ 1546 x(1) = 16HABCDEFGHIJKLMNOP 1547 call foo (4h abc) 1548 @end smallexample 1549 1550 Examples of Hollerith constants: 1551 @smallexample 1552 integer*4 a 1553 a = 0H ! Invalid, at least one character is needed. 1554 a = 4HAB12 ! Valid 1555 a = 8H12345678 ! Valid, but the Hollerith constant will be truncated. 1556 a = 3Hxyz ! Valid, but the Hollerith constant will be padded. 1557 @end smallexample 1558 1559 In general, Hollerith constants were used to provide a rudimentary 1560 facility for handling character strings in early Fortran compilers, 1561 prior to the introduction of @code{CHARACTER} variables in Fortran 77; 1562 in those cases, the standard-compliant equivalent is to convert the 1563 program to use proper character strings. On occasion, there may be a 1564 case where the intent is specifically to initialize a numeric variable 1565 with a given byte sequence. In these cases, the same result can be 1566 obtained by using the @code{TRANSFER} statement, as in this example. 1567 @smallexample 1568 integer(kind=4) :: a 1569 a = transfer ("abcd", a) ! equivalent to: a = 4Habcd 1570 @end smallexample 1571 1572 The use of the @option{-fdec} option extends support of Hollerith constants 1573 to comparisons: 1574 @smallexample 1575 integer*4 a 1576 a = 4hABCD 1577 if (a .ne. 4habcd) then 1578 write(*,*) "no match" 1579 end if 1580 @end smallexample 1581 1582 Supported types are numeric (@code{INTEGER}, @code{REAL}, or @code{COMPLEX}), 1583 and @code{CHARACTER}. 1584 1585 @node Character conversion 1586 @subsection Character conversion 1587 @cindex conversion, to character 1588 1589 Allowing character literals to be used in a similar way to Hollerith constants 1590 is a non-standard extension. This feature is enabled using 1591 -fdec-char-conversions and only applies to character literals of @code{kind=1}. 1592 1593 Character literals can be used in @code{DATA} statements and assignments with 1594 numeric (@code{INTEGER}, @code{REAL}, or @code{COMPLEX}) or @code{LOGICAL} 1595 variables. Like Hollerith constants they are copied byte-wise fashion. The 1596 constant will be padded with spaces or truncated to fit the size of the 1597 variable in which it is stored. 1598 1599 Examples: 1600 @smallexample 1601 integer*4 x 1602 data x / 'abcd' / 1603 1604 x = 'A' ! Will be padded. 1605 x = 'ab1234' ! Will be truncated. 1606 @end smallexample 1607 1608 1609 @node Cray pointers 1610 @subsection Cray pointers 1611 @cindex pointer, Cray 1612 1613 Cray pointers are part of a non-standard extension that provides a 1614 C-like pointer in Fortran. This is accomplished through a pair of 1615 variables: an integer "pointer" that holds a memory address, and a 1616 "pointee" that is used to dereference the pointer. 1617 1618 Pointer/pointee pairs are declared in statements of the form: 1619 @smallexample 1620 pointer ( <pointer> , <pointee> ) 1621 @end smallexample 1622 or, 1623 @smallexample 1624 pointer ( <pointer1> , <pointee1> ), ( <pointer2> , <pointee2> ), ... 1625 @end smallexample 1626 The pointer is an integer that is intended to hold a memory address. 1627 The pointee may be an array or scalar. 1628 If an assumed-size array is permitted within the scoping unit, a 1629 pointee can be an assumed-size array. 1630 That is, the last dimension may be left unspecified by using a @code{*} 1631 in place of a value. A pointee cannot be an assumed shape array. 1632 No space is allocated for the pointee. 1633 1634 The pointee may have its type declared before or after the pointer 1635 statement, and its array specification (if any) may be declared 1636 before, during, or after the pointer statement. The pointer may be 1637 declared as an integer prior to the pointer statement. However, some 1638 machines have default integer sizes that are different than the size 1639 of a pointer, and so the following code is not portable: 1640 @smallexample 1641 integer ipt 1642 pointer (ipt, iarr) 1643 @end smallexample 1644 If a pointer is declared with a kind that is too small, the compiler 1645 will issue a warning; the resulting binary will probably not work 1646 correctly, because the memory addresses stored in the pointers may be 1647 truncated. It is safer to omit the first line of the above example; 1648 if explicit declaration of ipt's type is omitted, then the compiler 1649 will ensure that ipt is an integer variable large enough to hold a 1650 pointer. 1651 1652 Pointer arithmetic is valid with Cray pointers, but it is not the same 1653 as C pointer arithmetic. Cray pointers are just ordinary integers, so 1654 the user is responsible for determining how many bytes to add to a 1655 pointer in order to increment it. Consider the following example: 1656 @smallexample 1657 real target(10) 1658 real pointee(10) 1659 pointer (ipt, pointee) 1660 ipt = loc (target) 1661 ipt = ipt + 1 1662 @end smallexample 1663 The last statement does not set @code{ipt} to the address of 1664 @code{target(1)}, as it would in C pointer arithmetic. Adding @code{1} 1665 to @code{ipt} just adds one byte to the address stored in @code{ipt}. 1666 1667 Any expression involving the pointee will be translated to use the 1668 value stored in the pointer as the base address. 1669 1670 To get the address of elements, this extension provides an intrinsic 1671 function @code{LOC()}. The @code{LOC()} function is equivalent to the 1672 @code{&} operator in C, except the address is cast to an integer type: 1673 @smallexample 1674 real ar(10) 1675 pointer(ipt, arpte(10)) 1676 real arpte 1677 ipt = loc(ar) ! Makes arpte is an alias for ar 1678 arpte(1) = 1.0 ! Sets ar(1) to 1.0 1679 @end smallexample 1680 The pointer can also be set by a call to the @code{MALLOC} intrinsic 1681 (see @ref{MALLOC}). 1682 1683 Cray pointees often are used to alias an existing variable. For 1684 example: 1685 @smallexample 1686 integer target(10) 1687 integer iarr(10) 1688 pointer (ipt, iarr) 1689 ipt = loc(target) 1690 @end smallexample 1691 As long as @code{ipt} remains unchanged, @code{iarr} is now an alias for 1692 @code{target}. The optimizer, however, will not detect this aliasing, so 1693 it is unsafe to use @code{iarr} and @code{target} simultaneously. Using 1694 a pointee in any way that violates the Fortran aliasing rules or 1695 assumptions is illegal. It is the user's responsibility to avoid doing 1696 this; the compiler works under the assumption that no such aliasing 1697 occurs. 1698 1699 Cray pointers will work correctly when there is no aliasing (i.e., when 1700 they are used to access a dynamically allocated block of memory), and 1701 also in any routine where a pointee is used, but any variable with which 1702 it shares storage is not used. Code that violates these rules may not 1703 run as the user intends. This is not a bug in the optimizer; any code 1704 that violates the aliasing rules is illegal. (Note that this is not 1705 unique to GNU Fortran; any Fortran compiler that supports Cray pointers 1706 will ``incorrectly'' optimize code with illegal aliasing.) 1707 1708 There are a number of restrictions on the attributes that can be applied 1709 to Cray pointers and pointees. Pointees may not have the 1710 @code{ALLOCATABLE}, @code{INTENT}, @code{OPTIONAL}, @code{DUMMY}, 1711 @code{TARGET}, @code{INTRINSIC}, or @code{POINTER} attributes. Pointers 1712 may not have the @code{DIMENSION}, @code{POINTER}, @code{TARGET}, 1713 @code{ALLOCATABLE}, @code{EXTERNAL}, or @code{INTRINSIC} attributes, nor 1714 may they be function results. Pointees may not occur in more than one 1715 pointer statement. A pointee cannot be a pointer. Pointees cannot occur 1716 in equivalence, common, or data statements. 1717 1718 A Cray pointer may also point to a function or a subroutine. For 1719 example, the following excerpt is valid: 1720 @smallexample 1721 implicit none 1722 external sub 1723 pointer (subptr,subpte) 1724 external subpte 1725 subptr = loc(sub) 1726 call subpte() 1727 [...] 1728 subroutine sub 1729 [...] 1730 end subroutine sub 1731 @end smallexample 1732 1733 A pointer may be modified during the course of a program, and this 1734 will change the location to which the pointee refers. However, when 1735 pointees are passed as arguments, they are treated as ordinary 1736 variables in the invoked function. Subsequent changes to the pointer 1737 will not change the base address of the array that was passed. 1738 1739 @node CONVERT specifier 1740 @subsection @code{CONVERT} specifier 1741 @cindex @code{CONVERT} specifier 1742 1743 GNU Fortran allows the conversion of unformatted data between little- 1744 and big-endian representation to facilitate moving of data 1745 between different systems. The conversion can be indicated with 1746 the @code{CONVERT} specifier on the @code{OPEN} statement. 1747 @xref{GFORTRAN_CONVERT_UNIT}, for an alternative way of specifying 1748 the data format via an environment variable. 1749 1750 Valid values for @code{CONVERT} on most systems are: 1751 @itemize @w{} 1752 @item @code{CONVERT='NATIVE'} Use the native format. This is the default. 1753 @item @code{CONVERT='SWAP'} Swap between little- and big-endian. 1754 @item @code{CONVERT='LITTLE_ENDIAN'} Use the little-endian representation 1755 for unformatted files. 1756 @item @code{CONVERT='BIG_ENDIAN'} Use the big-endian representation for 1757 unformatted files. 1758 @end itemize 1759 On POWER systems which support @option{-mabi=ieeelongdouble}, 1760 there are additional options, which can be combined with the others 1761 with commas. Those are 1762 @itemize @w{} 1763 @item @code{CONVERT='R16_IEEE'} Use IEEE 128-bit format for 1764 @code{REAL(KIND=16)}. 1765 @item @code{CONVERT='R16_IBM'} Use IBM @code{long double} format for 1766 real@code{REAL(KIND=16)}. 1767 @end itemize 1768 1769 Using the option could look like this: 1770 @smallexample 1771 open(file='big.dat',form='unformatted',access='sequential', & 1772 convert='big_endian') 1773 @end smallexample 1774 1775 The value of the conversion can be queried by using 1776 @code{INQUIRE(CONVERT=ch)}. The values returned are 1777 @code{'BIG_ENDIAN'} and @code{'LITTLE_ENDIAN'}. 1778 1779 @code{CONVERT} works between big- and little-endian for 1780 @code{INTEGER} values of all supported kinds and for @code{REAL} 1781 on IEEE systems of kinds 4 and 8. Conversion between different 1782 ``extended double'' types on different architectures such as 1783 m68k and x86_64, which GNU Fortran 1784 supports as @code{REAL(KIND=10)} and @code{REAL(KIND=16)}, will 1785 probably not work. 1786 1787 @emph{Note that the values specified via the GFORTRAN_CONVERT_UNIT 1788 environment variable will override the CONVERT specifier in the 1789 open statement}. This is to give control over data formats to 1790 users who do not have the source code of their program available. 1791 1792 Using anything but the native representation for unformatted data 1793 carries a significant speed overhead. If speed in this area matters 1794 to you, it is best if you use this only for data that needs to be 1795 portable. 1796 1797 @node OpenMP 1798 @subsection OpenMP 1799 @cindex OpenMP 1800 1801 OpenMP (Open Multi-Processing) is an application programming 1802 interface (API) that supports multi-platform shared memory 1803 multiprocessing programming in C/C++ and Fortran on many 1804 architectures, including Unix and Microsoft Windows platforms. 1805 It consists of a set of compiler directives, library routines, 1806 and environment variables that influence run-time behavior. 1807 1808 GNU Fortran strives to be compatible to the 1809 @uref{https://openmp.org/wp/openmp-specifications/, 1810 OpenMP Application Program Interface v4.5}. 1811 1812 To enable the processing of the OpenMP directive @code{!$omp} in 1813 free-form source code; the @code{c$omp}, @code{*$omp} and @code{!$omp} 1814 directives in fixed form; the @code{!$} conditional compilation sentinels 1815 in free form; and the @code{c$}, @code{*$} and @code{!$} sentinels 1816 in fixed form, @command{gfortran} needs to be invoked with the 1817 @option{-fopenmp}. This also arranges for automatic linking of the 1818 GNU Offloading and Multi Processing Runtime Library 1819 @ref{Top,,libgomp,libgomp,GNU Offloading and Multi Processing Runtime 1820 Library}. 1821 1822 The OpenMP Fortran runtime library routines are provided both in a 1823 form of a Fortran 90 module named @code{omp_lib} and in a form of 1824 a Fortran @code{include} file named @file{omp_lib.h}. 1825 1826 An example of a parallelized loop taken from Appendix A.1 of 1827 the OpenMP Application Program Interface v2.5: 1828 @smallexample 1829 SUBROUTINE A1(N, A, B) 1830 INTEGER I, N 1831 REAL B(N), A(N) 1832 !$OMP PARALLEL DO !I is private by default 1833 DO I=2,N 1834 B(I) = (A(I) + A(I-1)) / 2.0 1835 ENDDO 1836 !$OMP END PARALLEL DO 1837 END SUBROUTINE A1 1838 @end smallexample 1839 1840 Please note: 1841 @itemize 1842 @item 1843 @option{-fopenmp} implies @option{-frecursive}, i.e., all local arrays 1844 will be allocated on the stack. When porting existing code to OpenMP, 1845 this may lead to surprising results, especially to segmentation faults 1846 if the stacksize is limited. 1847 1848 @item 1849 On glibc-based systems, OpenMP enabled applications cannot be statically 1850 linked due to limitations of the underlying pthreads-implementation. It 1851 might be possible to get a working solution if 1852 @command{-Wl,--whole-archive -lpthread -Wl,--no-whole-archive} is added 1853 to the command line. However, this is not supported by @command{gcc} and 1854 thus not recommended. 1855 @end itemize 1856 1857 @node OpenACC 1858 @subsection OpenACC 1859 @cindex OpenACC 1860 1861 OpenACC is an application programming interface (API) that supports 1862 offloading of code to accelerator devices. It consists of a set of 1863 compiler directives, library routines, and environment variables that 1864 influence run-time behavior. 1865 1866 GNU Fortran strives to be compatible to the 1867 @uref{https://www.openacc.org/, OpenACC Application Programming 1868 Interface v2.6}. 1869 1870 To enable the processing of the OpenACC directive @code{!$acc} in 1871 free-form source code; the @code{c$acc}, @code{*$acc} and @code{!$acc} 1872 directives in fixed form; the @code{!$} conditional compilation 1873 sentinels in free form; and the @code{c$}, @code{*$} and @code{!$} 1874 sentinels in fixed form, @command{gfortran} needs to be invoked with 1875 the @option{-fopenacc}. This also arranges for automatic linking of 1876 the GNU Offloading and Multi Processing Runtime Library 1877 @ref{Top,,libgomp,libgomp,GNU Offloading and Multi Processing Runtime 1878 Library}. 1879 1880 The OpenACC Fortran runtime library routines are provided both in a 1881 form of a Fortran 90 module named @code{openacc} and in a form of a 1882 Fortran @code{include} file named @file{openacc_lib.h}. 1883 1884 @node Argument list functions 1885 @subsection Argument list functions @code{%VAL}, @code{%REF} and @code{%LOC} 1886 @cindex argument list functions 1887 @cindex @code{%VAL} 1888 @cindex @code{%REF} 1889 @cindex @code{%LOC} 1890 1891 GNU Fortran supports argument list functions @code{%VAL}, @code{%REF} 1892 and @code{%LOC} statements, for backward compatibility with g77. 1893 It is recommended that these should be used only for code that is 1894 accessing facilities outside of GNU Fortran, such as operating system 1895 or windowing facilities. It is best to constrain such uses to isolated 1896 portions of a program--portions that deal specifically and exclusively 1897 with low-level, system-dependent facilities. Such portions might well 1898 provide a portable interface for use by the program as a whole, but are 1899 themselves not portable, and should be thoroughly tested each time they 1900 are rebuilt using a new compiler or version of a compiler. 1901 1902 @code{%VAL} passes a scalar argument by value, @code{%REF} passes it by 1903 reference and @code{%LOC} passes its memory location. Since gfortran 1904 already passes scalar arguments by reference, @code{%REF} is in effect 1905 a do-nothing. @code{%LOC} has the same effect as a Fortran pointer. 1906 1907 An example of passing an argument by value to a C subroutine foo.: 1908 @smallexample 1909 C 1910 C prototype void foo_ (float x); 1911 C 1912 external foo 1913 real*4 x 1914 x = 3.14159 1915 call foo (%VAL (x)) 1916 end 1917 @end smallexample 1918 1919 For details refer to the g77 manual 1920 @uref{https://gcc.gnu.org/@/onlinedocs/@/gcc-3.4.6/@/g77/@/index.html#Top}. 1921 1922 Also, @code{c_by_val.f} and its partner @code{c_by_val.c} of the 1923 GNU Fortran testsuite are worth a look. 1924 1925 @node Read/Write after EOF marker 1926 @subsection Read/Write after EOF marker 1927 @cindex @code{EOF} 1928 @cindex @code{BACKSPACE} 1929 @cindex @code{REWIND} 1930 1931 Some legacy codes rely on allowing @code{READ} or @code{WRITE} after the 1932 EOF file marker in order to find the end of a file. GNU Fortran normally 1933 rejects these codes with a run-time error message and suggests the user 1934 consider @code{BACKSPACE} or @code{REWIND} to properly position 1935 the file before the EOF marker. As an extension, the run-time error may 1936 be disabled using -std=legacy. 1937 1938 1939 @node STRUCTURE and RECORD 1940 @subsection @code{STRUCTURE} and @code{RECORD} 1941 @cindex @code{STRUCTURE} 1942 @cindex @code{RECORD} 1943 1944 Record structures are a pre-Fortran-90 vendor extension to create 1945 user-defined aggregate data types. Support for record structures in GNU 1946 Fortran can be enabled with the @option{-fdec-structure} compile flag. 1947 If you have a choice, you should instead use Fortran 90's ``derived types'', 1948 which have a different syntax. 1949 1950 In many cases, record structures can easily be converted to derived types. 1951 To convert, replace @code{STRUCTURE /}@var{structure-name}@code{/} 1952 by @code{TYPE} @var{type-name}. Additionally, replace 1953 @code{RECORD /}@var{structure-name}@code{/} by 1954 @code{TYPE(}@var{type-name}@code{)}. Finally, in the component access, 1955 replace the period (@code{.}) by the percent sign (@code{%}). 1956 1957 Here is an example of code using the non portable record structure syntax: 1958 1959 @example 1960 ! Declaring a structure named ``item'' and containing three fields: 1961 ! an integer ID, an description string and a floating-point price. 1962 STRUCTURE /item/ 1963 INTEGER id 1964 CHARACTER(LEN=200) description 1965 REAL price 1966 END STRUCTURE 1967 1968 ! Define two variables, an single record of type ``item'' 1969 ! named ``pear'', and an array of items named ``store_catalog'' 1970 RECORD /item/ pear, store_catalog(100) 1971 1972 ! We can directly access the fields of both variables 1973 pear.id = 92316 1974 pear.description = "juicy D'Anjou pear" 1975 pear.price = 0.15 1976 store_catalog(7).id = 7831 1977 store_catalog(7).description = "milk bottle" 1978 store_catalog(7).price = 1.2 1979 1980 ! We can also manipulate the whole structure 1981 store_catalog(12) = pear 1982 print *, store_catalog(12) 1983 @end example 1984 1985 @noindent 1986 This code can easily be rewritten in the Fortran 90 syntax as following: 1987 1988 @example 1989 ! ``STRUCTURE /name/ ... END STRUCTURE'' becomes 1990 ! ``TYPE name ... END TYPE'' 1991 TYPE item 1992 INTEGER id 1993 CHARACTER(LEN=200) description 1994 REAL price 1995 END TYPE 1996 1997 ! ``RECORD /name/ variable'' becomes ``TYPE(name) variable'' 1998 TYPE(item) pear, store_catalog(100) 1999 2000 ! Instead of using a dot (.) to access fields of a record, the 2001 ! standard syntax uses a percent sign (%) 2002 pear%id = 92316 2003 pear%description = "juicy D'Anjou pear" 2004 pear%price = 0.15 2005 store_catalog(7)%id = 7831 2006 store_catalog(7)%description = "milk bottle" 2007 store_catalog(7)%price = 1.2 2008 2009 ! Assignments of a whole variable do not change 2010 store_catalog(12) = pear 2011 print *, store_catalog(12) 2012 @end example 2013 2014 @noindent 2015 GNU Fortran implements STRUCTURES like derived types with the following 2016 rules and exceptions: 2017 2018 @itemize @bullet 2019 @item Structures act like derived types with the @code{SEQUENCE} attribute. 2020 Otherwise they may contain no specifiers. 2021 2022 @item Structures may contain a special field with the name @code{%FILL}. 2023 This will create an anonymous component which cannot be accessed but occupies 2024 space just as if a component of the same type was declared in its place, useful 2025 for alignment purposes. As an example, the following structure will consist 2026 of at least sixteen bytes: 2027 2028 @smallexample 2029 structure /padded/ 2030 character(4) start 2031 character(8) %FILL 2032 character(4) end 2033 end structure 2034 @end smallexample 2035 2036 @item Structures may share names with other symbols. For example, the following 2037 is invalid for derived types, but valid for structures: 2038 2039 @smallexample 2040 structure /header/ 2041 ! ... 2042 end structure 2043 record /header/ header 2044 @end smallexample 2045 2046 @item Structure types may be declared nested within another parent structure. 2047 The syntax is: 2048 @smallexample 2049 structure /type-name/ 2050 ... 2051 structure [/<type-name>/] <field-list> 2052 ... 2053 @end smallexample 2054 2055 The type name may be ommitted, in which case the structure type itself is 2056 anonymous, and other structures of the same type cannot be instantiated. The 2057 following shows some examples: 2058 2059 @example 2060 structure /appointment/ 2061 ! nested structure definition: app_time is an array of two 'time' 2062 structure /time/ app_time (2) 2063 integer(1) hour, minute 2064 end structure 2065 character(10) memo 2066 end structure 2067 2068 ! The 'time' structure is still usable 2069 record /time/ now 2070 now = time(5, 30) 2071 2072 ... 2073 2074 structure /appointment/ 2075 ! anonymous nested structure definition 2076 structure start, end 2077 integer(1) hour, minute 2078 end structure 2079 character(10) memo 2080 end structure 2081 @end example 2082 2083 @item Structures may contain @code{UNION} blocks. For more detail see the 2084 section on @ref{UNION and MAP}. 2085 2086 @item Structures support old-style initialization of components, like 2087 those described in @ref{Old-style variable initialization}. For array 2088 initializers, an initializer may contain a repeat specification of the form 2089 @code{<literal-integer> * <constant-initializer>}. The value of the integer 2090 indicates the number of times to repeat the constant initializer when expanding 2091 the initializer list. 2092 @end itemize 2093 2094 @node UNION and MAP 2095 @subsection @code{UNION} and @code{MAP} 2096 @cindex @code{UNION} 2097 @cindex @code{MAP} 2098 2099 Unions are an old vendor extension which were commonly used with the 2100 non-standard @ref{STRUCTURE and RECORD} extensions. Use of @code{UNION} and 2101 @code{MAP} is automatically enabled with @option{-fdec-structure}. 2102 2103 A @code{UNION} declaration occurs within a structure; within the definition of 2104 each union is a number of @code{MAP} blocks. Each @code{MAP} shares storage 2105 with its sibling maps (in the same union), and the size of the union is the 2106 size of the largest map within it, just as with unions in C. The major 2107 difference is that component references do not indicate which union or map the 2108 component is in (the compiler gets to figure that out). 2109 2110 Here is a small example: 2111 @smallexample 2112 structure /myunion/ 2113 union 2114 map 2115 character(2) w0, w1, w2 2116 end map 2117 map 2118 character(6) long 2119 end map 2120 end union 2121 end structure 2122 2123 record /myunion/ rec 2124 ! After this assignment... 2125 rec.long = 'hello!' 2126 2127 ! The following is true: 2128 ! rec.w0 === 'he' 2129 ! rec.w1 === 'll' 2130 ! rec.w2 === 'o!' 2131 @end smallexample 2132 2133 The two maps share memory, and the size of the union is ultimately six bytes: 2134 2135 @example 2136 0 1 2 3 4 5 6 Byte offset 2137 ------------------------------- 2138 | | | | | | | 2139 ------------------------------- 2140 2141 ^ W0 ^ W1 ^ W2 ^ 2142 \-------/ \-------/ \-------/ 2143 2144 ^ LONG ^ 2145 \---------------------------/ 2146 @end example 2147 2148 Following is an example mirroring the layout of an Intel x86_64 register: 2149 2150 @example 2151 structure /reg/ 2152 union ! U0 ! rax 2153 map 2154 character(16) rx 2155 end map 2156 map 2157 character(8) rh ! rah 2158 union ! U1 2159 map 2160 character(8) rl ! ral 2161 end map 2162 map 2163 character(8) ex ! eax 2164 end map 2165 map 2166 character(4) eh ! eah 2167 union ! U2 2168 map 2169 character(4) el ! eal 2170 end map 2171 map 2172 character(4) x ! ax 2173 end map 2174 map 2175 character(2) h ! ah 2176 character(2) l ! al 2177 end map 2178 end union 2179 end map 2180 end union 2181 end map 2182 end union 2183 end structure 2184 record /reg/ a 2185 2186 ! After this assignment... 2187 a.rx = 'AAAAAAAA.BBB.C.D' 2188 2189 ! The following is true: 2190 a.rx === 'AAAAAAAA.BBB.C.D' 2191 a.rh === 'AAAAAAAA' 2192 a.rl === '.BBB.C.D' 2193 a.ex === '.BBB.C.D' 2194 a.eh === '.BBB' 2195 a.el === '.C.D' 2196 a.x === '.C.D' 2197 a.h === '.C' 2198 a.l === '.D' 2199 @end example 2200 2201 @node Type variants for integer intrinsics 2202 @subsection Type variants for integer intrinsics 2203 @cindex intrinsics, integer 2204 2205 Similar to the D/C prefixes to real functions to specify the input/output 2206 types, GNU Fortran offers B/I/J/K prefixes to integer functions for 2207 compatibility with DEC programs. The types implied by each are: 2208 2209 @example 2210 @code{B} - @code{INTEGER(kind=1)} 2211 @code{I} - @code{INTEGER(kind=2)} 2212 @code{J} - @code{INTEGER(kind=4)} 2213 @code{K} - @code{INTEGER(kind=8)} 2214 @end example 2215 2216 GNU Fortran supports these with the flag @option{-fdec-intrinsic-ints}. 2217 Intrinsics for which prefixed versions are available and in what form are noted 2218 in @ref{Intrinsic Procedures}. The complete list of supported intrinsics is 2219 here: 2220 2221 @multitable @columnfractions .2 .2 .2 .2 .2 2222 2223 @headitem Intrinsic @tab B @tab I @tab J @tab K 2224 2225 @item @code{@ref{ABS}} 2226 @tab @code{BABS} @tab @code{IIABS} @tab @code{JIABS} @tab @code{KIABS} 2227 @item @code{@ref{BTEST}} 2228 @tab @code{BBTEST} @tab @code{BITEST} @tab @code{BJTEST} @tab @code{BKTEST} 2229 @item @code{@ref{IAND}} 2230 @tab @code{BIAND} @tab @code{IIAND} @tab @code{JIAND} @tab @code{KIAND} 2231 @item @code{@ref{IBCLR}} 2232 @tab @code{BBCLR} @tab @code{IIBCLR} @tab @code{JIBCLR} @tab @code{KIBCLR} 2233 @item @code{@ref{IBITS}} 2234 @tab @code{BBITS} @tab @code{IIBITS} @tab @code{JIBITS} @tab @code{KIBITS} 2235 @item @code{@ref{IBSET}} 2236 @tab @code{BBSET} @tab @code{IIBSET} @tab @code{JIBSET} @tab @code{KIBSET} 2237 @item @code{@ref{IEOR}} 2238 @tab @code{BIEOR} @tab @code{IIEOR} @tab @code{JIEOR} @tab @code{KIEOR} 2239 @item @code{@ref{IOR}} 2240 @tab @code{BIOR} @tab @code{IIOR} @tab @code{JIOR} @tab @code{KIOR} 2241 @item @code{@ref{ISHFT}} 2242 @tab @code{BSHFT} @tab @code{IISHFT} @tab @code{JISHFT} @tab @code{KISHFT} 2243 @item @code{@ref{ISHFTC}} 2244 @tab @code{BSHFTC} @tab @code{IISHFTC} @tab @code{JISHFTC} @tab @code{KISHFTC} 2245 @item @code{@ref{MOD}} 2246 @tab @code{BMOD} @tab @code{IMOD} @tab @code{JMOD} @tab @code{KMOD} 2247 @item @code{@ref{NOT}} 2248 @tab @code{BNOT} @tab @code{INOT} @tab @code{JNOT} @tab @code{KNOT} 2249 @item @code{@ref{REAL}} 2250 @tab @code{--} @tab @code{FLOATI} @tab @code{FLOATJ} @tab @code{FLOATK} 2251 @end multitable 2252 2253 @node AUTOMATIC and STATIC attributes 2254 @subsection @code{AUTOMATIC} and @code{STATIC} attributes 2255 @cindex variable attributes 2256 @cindex @code{AUTOMATIC} 2257 @cindex @code{STATIC} 2258 2259 With @option{-fdec-static} GNU Fortran supports the DEC extended attributes 2260 @code{STATIC} and @code{AUTOMATIC} to provide explicit specification of entity 2261 storage. These follow the syntax of the Fortran standard @code{SAVE} attribute. 2262 2263 @code{STATIC} is exactly equivalent to @code{SAVE}, and specifies that 2264 an entity should be allocated in static memory. As an example, @code{STATIC} 2265 local variables will retain their values across multiple calls to a function. 2266 2267 Entities marked @code{AUTOMATIC} will be stack automatic whenever possible. 2268 @code{AUTOMATIC} is the default for local variables smaller than 2269 @option{-fmax-stack-var-size}, unless @option{-fno-automatic} is given. This 2270 attribute overrides @option{-fno-automatic}, @option{-fmax-stack-var-size}, and 2271 blanket @code{SAVE} statements. 2272 2273 2274 Examples: 2275 2276 @example 2277 subroutine f 2278 integer, automatic :: i ! automatic variable 2279 integer x, y ! static variables 2280 save 2281 ... 2282 endsubroutine 2283 @end example 2284 @example 2285 subroutine f 2286 integer a, b, c, x, y, z 2287 static :: x 2288 save y 2289 automatic z, c 2290 ! a, b, c, and z are automatic 2291 ! x and y are static 2292 endsubroutine 2293 @end example 2294 @example 2295 ! Compiled with -fno-automatic 2296 subroutine f 2297 integer a, b, c, d 2298 automatic :: a 2299 ! a is automatic; b, c, and d are static 2300 endsubroutine 2301 @end example 2302 2303 @node Extended math intrinsics 2304 @subsection Extended math intrinsics 2305 @cindex intrinsics, math 2306 @cindex intrinsics, trigonometric functions 2307 2308 GNU Fortran supports an extended list of mathematical intrinsics with the 2309 compile flag @option{-fdec-math} for compatability with legacy code. 2310 These intrinsics are described fully in @ref{Intrinsic Procedures} where it is 2311 noted that they are extensions and should be avoided whenever possible. 2312 2313 Specifically, @option{-fdec-math} enables the @ref{COTAN} intrinsic, and 2314 trigonometric intrinsics which accept or produce values in degrees instead of 2315 radians. Here is a summary of the new intrinsics: 2316 2317 @multitable @columnfractions .5 .5 2318 @headitem Radians @tab Degrees 2319 @item @code{@ref{ACOS}} @tab @code{@ref{ACOSD}}* 2320 @item @code{@ref{ASIN}} @tab @code{@ref{ASIND}}* 2321 @item @code{@ref{ATAN}} @tab @code{@ref{ATAND}}* 2322 @item @code{@ref{ATAN2}} @tab @code{@ref{ATAN2D}}* 2323 @item @code{@ref{COS}} @tab @code{@ref{COSD}}* 2324 @item @code{@ref{COTAN}}* @tab @code{@ref{COTAND}}* 2325 @item @code{@ref{SIN}} @tab @code{@ref{SIND}}* 2326 @item @code{@ref{TAN}} @tab @code{@ref{TAND}}* 2327 @end multitable 2328 2329 * Enabled with @option{-fdec-math}. 2330 2331 For advanced users, it may be important to know the implementation of these 2332 functions. They are simply wrappers around the standard radian functions, which 2333 have more accurate builtin versions. These functions convert their arguments 2334 (or results) to degrees (or radians) by taking the value modulus 360 (or 2*pi) 2335 and then multiplying it by a constant radian-to-degree (or degree-to-radian) 2336 factor, as appropriate. The factor is computed at compile-time as 180/pi (or 2337 pi/180). 2338 2339 @node Form feed as whitespace 2340 @subsection Form feed as whitespace 2341 @cindex form feed whitespace 2342 2343 Historically, legacy compilers allowed insertion of form feed characters ('\f', 2344 ASCII 0xC) at the beginning of lines for formatted output to line printers, 2345 though the Fortran standard does not mention this. GNU Fortran supports the 2346 interpretation of form feed characters in source as whitespace for 2347 compatibility. 2348 2349 @node TYPE as an alias for PRINT 2350 @subsection TYPE as an alias for PRINT 2351 @cindex type alias print 2352 For compatibility, GNU Fortran will interpret @code{TYPE} statements as 2353 @code{PRINT} statements with the flag @option{-fdec}. With this flag asserted, 2354 the following two examples are equivalent: 2355 2356 @smallexample 2357 TYPE *, 'hello world' 2358 @end smallexample 2359 2360 @smallexample 2361 PRINT *, 'hello world' 2362 @end smallexample 2363 2364 @node %LOC as an rvalue 2365 @subsection %LOC as an rvalue 2366 @cindex LOC 2367 Normally @code{%LOC} is allowed only in parameter lists. However the intrinsic 2368 function @code{LOC} does the same thing, and is usable as the right-hand-side of 2369 assignments. For compatibility, GNU Fortran supports the use of @code{%LOC} as 2370 an alias for the builtin @code{LOC} with @option{-std=legacy}. With this 2371 feature enabled the following two examples are equivalent: 2372 2373 @smallexample 2374 integer :: i, l 2375 l = %loc(i) 2376 call sub(l) 2377 @end smallexample 2378 2379 @smallexample 2380 integer :: i 2381 call sub(%loc(i)) 2382 @end smallexample 2383 2384 @node .XOR. operator 2385 @subsection .XOR. operator 2386 @cindex operators, xor 2387 2388 GNU Fortran supports @code{.XOR.} as a logical operator with @code{-std=legacy} 2389 for compatibility with legacy code. @code{.XOR.} is equivalent to 2390 @code{.NEQV.}. That is, the output is true if and only if the inputs differ. 2391 2392 @node Bitwise logical operators 2393 @subsection Bitwise logical operators 2394 @cindex logical, bitwise 2395 2396 With @option{-fdec}, GNU Fortran relaxes the type constraints on 2397 logical operators to allow integer operands, and performs the corresponding 2398 bitwise operation instead. This flag is for compatibility only, and should be 2399 avoided in new code. Consider: 2400 2401 @smallexample 2402 INTEGER :: i, j 2403 i = z'33' 2404 j = z'cc' 2405 print *, i .AND. j 2406 @end smallexample 2407 2408 In this example, compiled with @option{-fdec}, GNU Fortran will 2409 replace the @code{.AND.} operation with a call to the intrinsic 2410 @code{@ref{IAND}} function, yielding the bitwise-and of @code{i} and @code{j}. 2411 2412 Note that this conversion will occur if at least one operand is of integral 2413 type. As a result, a logical operand will be converted to an integer when the 2414 other operand is an integer in a logical operation. In this case, 2415 @code{.TRUE.} is converted to @code{1} and @code{.FALSE.} to @code{0}. 2416 2417 Here is the mapping of logical operator to bitwise intrinsic used with 2418 @option{-fdec}: 2419 2420 @multitable @columnfractions .25 .25 .5 2421 @headitem Operator @tab Intrinsic @tab Bitwise operation 2422 @item @code{.NOT.} @tab @code{@ref{NOT}} @tab complement 2423 @item @code{.AND.} @tab @code{@ref{IAND}} @tab intersection 2424 @item @code{.OR.} @tab @code{@ref{IOR}} @tab union 2425 @item @code{.NEQV.} @tab @code{@ref{IEOR}} @tab exclusive or 2426 @item @code{.EQV.} @tab @code{@ref{NOT}(@ref{IEOR})} @tab complement of exclusive or 2427 @end multitable 2428 2429 @node Extended I/O specifiers 2430 @subsection Extended I/O specifiers 2431 @cindex @code{CARRIAGECONTROL} 2432 @cindex @code{READONLY} 2433 @cindex @code{SHARE} 2434 @cindex @code{SHARED} 2435 @cindex @code{NOSHARED} 2436 @cindex I/O specifiers 2437 2438 GNU Fortran supports the additional legacy I/O specifiers 2439 @code{CARRIAGECONTROL}, @code{READONLY}, and @code{SHARE} with the 2440 compile flag @option{-fdec}, for compatibility. 2441 2442 @table @code 2443 @item CARRIAGECONTROL 2444 The @code{CARRIAGECONTROL} specifier allows a user to control line 2445 termination settings between output records for an I/O unit. The specifier has 2446 no meaning for readonly files. When @code{CARRAIGECONTROL} is specified upon 2447 opening a unit for formatted writing, the exact @code{CARRIAGECONTROL} setting 2448 determines what characters to write between output records. The syntax is: 2449 2450 @smallexample 2451 OPEN(..., CARRIAGECONTROL=cc) 2452 @end smallexample 2453 2454 Where @emph{cc} is a character expression that evaluates to one of the 2455 following values: 2456 2457 @multitable @columnfractions .2 .8 2458 @item @code{'LIST'} @tab One line feed between records (default) 2459 @item @code{'FORTRAN'} @tab Legacy interpretation of the first character (see below) 2460 @item @code{'NONE'} @tab No separator between records 2461 @end multitable 2462 2463 With @code{CARRIAGECONTROL='FORTRAN'}, when a record is written, the first 2464 character of the input record is not written, and instead determines the output 2465 record separator as follows: 2466 2467 @multitable @columnfractions .3 .3 .4 2468 @headitem Leading character @tab Meaning @tab Output separating character(s) 2469 @item @code{'+'} @tab Overprinting @tab Carriage return only 2470 @item @code{'-'} @tab New line @tab Line feed and carriage return 2471 @item @code{'0'} @tab Skip line @tab Two line feeds and carriage return 2472 @item @code{'1'} @tab New page @tab Form feed and carriage return 2473 @item @code{'$'} @tab Prompting @tab Line feed (no carriage return) 2474 @item @code{CHAR(0)} @tab Overprinting (no advance) @tab None 2475 @end multitable 2476 2477 @item READONLY 2478 The @code{READONLY} specifier may be given upon opening a unit, and is 2479 equivalent to specifying @code{ACTION='READ'}, except that the file may not be 2480 deleted on close (i.e. @code{CLOSE} with @code{STATUS="DELETE"}). The syntax 2481 is: 2482 2483 @smallexample 2484 @code{OPEN(..., READONLY)} 2485 @end smallexample 2486 2487 @item SHARE 2488 The @code{SHARE} specifier allows system-level locking on a unit upon opening 2489 it for controlled access from multiple processes/threads. The @code{SHARE} 2490 specifier has several forms: 2491 2492 @smallexample 2493 OPEN(..., SHARE=sh) 2494 OPEN(..., SHARED) 2495 OPEN(..., NOSHARED) 2496 @end smallexample 2497 2498 Where @emph{sh} in the first form is a character expression that evaluates to 2499 a value as seen in the table below. The latter two forms are aliases 2500 for particular values of @emph{sh}: 2501 2502 @multitable @columnfractions .3 .3 .4 2503 @headitem Explicit form @tab Short form @tab Meaning 2504 @item @code{SHARE='DENYRW'} @tab @code{NOSHARED} @tab Exclusive (write) lock 2505 @item @code{SHARE='DENYNONE'} @tab @code{SHARED} @tab Shared (read) lock 2506 @end multitable 2507 2508 In general only one process may hold an exclusive (write) lock for a given file 2509 at a time, whereas many processes may hold shared (read) locks for the same 2510 file. 2511 2512 The behavior of locking may vary with your operating system. On POSIX systems, 2513 locking is implemented with @code{fcntl}. Consult your corresponding operating 2514 system's manual pages for further details. Locking via @code{SHARE=} is not 2515 supported on other systems. 2516 2517 @end table 2518 2519 @node Legacy PARAMETER statements 2520 @subsection Legacy PARAMETER statements 2521 @cindex PARAMETER 2522 2523 For compatibility, GNU Fortran supports legacy PARAMETER statements without 2524 parentheses with @option{-std=legacy}. A warning is emitted if used with 2525 @option{-std=gnu}, and an error is acknowledged with a real Fortran standard 2526 flag (@option{-std=f95}, etc...). These statements take the following form: 2527 2528 @smallexample 2529 implicit real (E) 2530 parameter e = 2.718282 2531 real c 2532 parameter c = 3.0e8 2533 @end smallexample 2534 2535 @node Default exponents 2536 @subsection Default exponents 2537 @cindex exponent 2538 2539 For compatibility, GNU Fortran supports a default exponent of zero in real 2540 constants with @option{-fdec}. For example, @code{9e} would be 2541 interpreted as @code{9e0}, rather than an error. 2542 2543 2544 @node Extensions not implemented in GNU Fortran 2545 @section Extensions not implemented in GNU Fortran 2546 @cindex extensions, not implemented 2547 2548 The long history of the Fortran language, its wide use and broad 2549 userbase, the large number of different compiler vendors and the lack of 2550 some features crucial to users in the first standards have lead to the 2551 existence of a number of important extensions to the language. While 2552 some of the most useful or popular extensions are supported by the GNU 2553 Fortran compiler, not all existing extensions are supported. This section 2554 aims at listing these extensions and offering advice on how best make 2555 code that uses them running with the GNU Fortran compiler. 2556 2557 @c More can be found here: 2558 @c -- https://gcc.gnu.org/onlinedocs/gcc-3.4.6/g77/Missing-Features.html 2559 @c -- the list of Fortran and libgfortran bugs closed as WONTFIX: 2560 @c http://tinyurl.com/2u4h5y 2561 2562 @menu 2563 * ENCODE and DECODE statements:: 2564 * Variable FORMAT expressions:: 2565 @c * TYPE and ACCEPT I/O Statements:: 2566 @c * DEFAULTFILE, DISPOSE and RECORDTYPE I/O specifiers:: 2567 @c * Omitted arguments in procedure call:: 2568 * Alternate complex function syntax:: 2569 * Volatile COMMON blocks:: 2570 * OPEN( ... NAME=):: 2571 * Q edit descriptor:: 2572 @end menu 2573 2574 @node ENCODE and DECODE statements 2575 @subsection @code{ENCODE} and @code{DECODE} statements 2576 @cindex @code{ENCODE} 2577 @cindex @code{DECODE} 2578 2579 GNU Fortran does not support the @code{ENCODE} and @code{DECODE} 2580 statements. These statements are best replaced by @code{READ} and 2581 @code{WRITE} statements involving internal files (@code{CHARACTER} 2582 variables and arrays), which have been part of the Fortran standard since 2583 Fortran 77. For example, replace a code fragment like 2584 2585 @smallexample 2586 INTEGER*1 LINE(80) 2587 REAL A, B, C 2588 c ... Code that sets LINE 2589 DECODE (80, 9000, LINE) A, B, C 2590 9000 FORMAT (1X, 3(F10.5)) 2591 @end smallexample 2592 2593 @noindent 2594 with the following: 2595 2596 @smallexample 2597 CHARACTER(LEN=80) LINE 2598 REAL A, B, C 2599 c ... Code that sets LINE 2600 READ (UNIT=LINE, FMT=9000) A, B, C 2601 9000 FORMAT (1X, 3(F10.5)) 2602 @end smallexample 2603 2604 Similarly, replace a code fragment like 2605 2606 @smallexample 2607 INTEGER*1 LINE(80) 2608 REAL A, B, C 2609 c ... Code that sets A, B and C 2610 ENCODE (80, 9000, LINE) A, B, C 2611 9000 FORMAT (1X, 'OUTPUT IS ', 3(F10.5)) 2612 @end smallexample 2613 2614 @noindent 2615 with the following: 2616 2617 @smallexample 2618 CHARACTER(LEN=80) LINE 2619 REAL A, B, C 2620 c ... Code that sets A, B and C 2621 WRITE (UNIT=LINE, FMT=9000) A, B, C 2622 9000 FORMAT (1X, 'OUTPUT IS ', 3(F10.5)) 2623 @end smallexample 2624 2625 2626 @node Variable FORMAT expressions 2627 @subsection Variable @code{FORMAT} expressions 2628 @cindex @code{FORMAT} 2629 2630 A variable @code{FORMAT} expression is format statement which includes 2631 angle brackets enclosing a Fortran expression: @code{FORMAT(I<N>)}. GNU 2632 Fortran does not support this legacy extension. The effect of variable 2633 format expressions can be reproduced by using the more powerful (and 2634 standard) combination of internal output and string formats. For example, 2635 replace a code fragment like this: 2636 2637 @smallexample 2638 WRITE(6,20) INT1 2639 20 FORMAT(I<N+1>) 2640 @end smallexample 2641 2642 @noindent 2643 with the following: 2644 2645 @smallexample 2646 c Variable declaration 2647 CHARACTER(LEN=20) FMT 2648 c 2649 c Other code here... 2650 c 2651 WRITE(FMT,'("(I", I0, ")")') N+1 2652 WRITE(6,FMT) INT1 2653 @end smallexample 2654 2655 @noindent 2656 or with: 2657 2658 @smallexample 2659 c Variable declaration 2660 CHARACTER(LEN=20) FMT 2661 c 2662 c Other code here... 2663 c 2664 WRITE(FMT,*) N+1 2665 WRITE(6,"(I" // ADJUSTL(FMT) // ")") INT1 2666 @end smallexample 2667 2668 2669 @node Alternate complex function syntax 2670 @subsection Alternate complex function syntax 2671 @cindex Complex function 2672 2673 Some Fortran compilers, including @command{g77}, let the user declare 2674 complex functions with the syntax @code{COMPLEX FUNCTION name*16()}, as 2675 well as @code{COMPLEX*16 FUNCTION name()}. Both are non-standard, legacy 2676 extensions. @command{gfortran} accepts the latter form, which is more 2677 common, but not the former. 2678 2679 2680 @node Volatile COMMON blocks 2681 @subsection Volatile @code{COMMON} blocks 2682 @cindex @code{VOLATILE} 2683 @cindex @code{COMMON} 2684 2685 Some Fortran compilers, including @command{g77}, let the user declare 2686 @code{COMMON} with the @code{VOLATILE} attribute. This is 2687 invalid standard Fortran syntax and is not supported by 2688 @command{gfortran}. Note that @command{gfortran} accepts 2689 @code{VOLATILE} variables in @code{COMMON} blocks since revision 4.3. 2690 2691 2692 @node OPEN( ... NAME=) 2693 @subsection @code{OPEN( ... NAME=)} 2694 @cindex @code{NAME} 2695 2696 Some Fortran compilers, including @command{g77}, let the user declare 2697 @code{OPEN( ... NAME=)}. This is 2698 invalid standard Fortran syntax and is not supported by 2699 @command{gfortran}. @code{OPEN( ... NAME=)} should be replaced 2700 with @code{OPEN( ... FILE=)}. 2701 2702 @node Q edit descriptor 2703 @subsection @code{Q} edit descriptor 2704 @cindex @code{Q} edit descriptor 2705 2706 Some Fortran compilers provide the @code{Q} edit descriptor, which 2707 transfers the number of characters left within an input record into an 2708 integer variable. 2709 2710 A direct replacement of the @code{Q} edit descriptor is not available 2711 in @command{gfortran}. How to replicate its functionality using 2712 standard-conforming code depends on what the intent of the original 2713 code is. 2714 2715 Options to replace @code{Q} may be to read the whole line into a 2716 character variable and then counting the number of non-blank 2717 characters left using @code{LEN_TRIM}. Another method may be to use 2718 formatted stream, read the data up to the position where the @code{Q} 2719 descriptor occurred, use @code{INQUIRE} to get the file position, 2720 count the characters up to the next @code{NEW_LINE} and then start 2721 reading from the position marked previously. 2722 2723 2724 @c --------------------------------------------------------------------- 2725 @c --------------------------------------------------------------------- 2726 @c Mixed-Language Programming 2727 @c --------------------------------------------------------------------- 2728 2729 @node Mixed-Language Programming 2730 @chapter Mixed-Language Programming 2731 @cindex Interoperability 2732 @cindex Mixed-language programming 2733 2734 @menu 2735 * Interoperability with C:: 2736 * GNU Fortran Compiler Directives:: 2737 * Non-Fortran Main Program:: 2738 * Naming and argument-passing conventions:: 2739 @end menu 2740 2741 This chapter is about mixed-language interoperability, but also 2742 applies if you link Fortran code compiled by different compilers. In 2743 most cases, use of the C Binding features of the Fortran 2003 and 2744 later standards is sufficient. 2745 2746 For example, it is possible to mix Fortran code with C++ code as well 2747 as C, if you declare the interface functions as @code{extern "C"} on 2748 the C++ side and @code{BIND(C)} on the Fortran side, and follow the 2749 rules for interoperability with C. Note that you cannot manipulate 2750 C++ class objects in Fortran or vice versa except as opaque pointers. 2751 2752 You can use the @command{gfortran} command to link both Fortran and 2753 non-Fortran code into the same program, or you can use @command{gcc} 2754 or @command{g++} if you also add an explicit @option{-lgfortran} option 2755 to link with the Fortran library. If your main program is written in 2756 C or some other language instead of Fortran, see 2757 @ref{Non-Fortran Main Program}, below. 2758 2759 @node Interoperability with C 2760 @section Interoperability with C 2761 @cindex interoperability with C 2762 @cindex C interoperability 2763 2764 @menu 2765 * Intrinsic Types:: 2766 * Derived Types and struct:: 2767 * Interoperable Global Variables:: 2768 * Interoperable Subroutines and Functions:: 2769 * Working with C Pointers:: 2770 * Further Interoperability of Fortran with C:: 2771 @end menu 2772 2773 Since Fortran 2003 (ISO/IEC 1539-1:2004(E)) there is a 2774 standardized way to generate procedure and derived-type 2775 declarations and global variables that are interoperable with C 2776 (ISO/IEC 9899:1999). The @code{BIND(C)} attribute has been added 2777 to inform the compiler that a symbol shall be interoperable with C; 2778 also, some constraints are added. Note, however, that not 2779 all C features have a Fortran equivalent or vice versa. For instance, 2780 neither C's unsigned integers nor C's functions with variable number 2781 of arguments have an equivalent in Fortran. 2782 2783 Note that array dimensions are reversely ordered in C and that arrays in 2784 C always start with index 0 while in Fortran they start by default with 2785 1. Thus, an array declaration @code{A(n,m)} in Fortran matches 2786 @code{A[m][n]} in C and accessing the element @code{A(i,j)} matches 2787 @code{A[j-1][i-1]}. The element following @code{A(i,j)} (C: @code{A[j-1][i-1]}; 2788 assuming @math{i < n}) in memory is @code{A(i+1,j)} (C: @code{A[j-1][i]}). 2789 2790 @node Intrinsic Types 2791 @subsection Intrinsic Types 2792 @cindex C intrinsic type interoperability 2793 @cindex intrinsic type interoperability with C 2794 @cindex interoperability, intrinsic type 2795 2796 In order to ensure that exactly the same variable type and kind is used 2797 in C and Fortran, you should use the named constants for kind parameters 2798 that are defined in the @code{ISO_C_BINDING} intrinsic module. 2799 That module contains named constants of character type representing 2800 the escaped special characters in C, such as newline. 2801 For a list of the constants, see @ref{ISO_C_BINDING}. 2802 2803 For logical types, please note that the Fortran standard only guarantees 2804 interoperability between C99's @code{_Bool} and Fortran's @code{C_Bool}-kind 2805 logicals and C99 defines that @code{true} has the value 1 and @code{false} 2806 the value 0. Using any other integer value with GNU Fortran's @code{LOGICAL} 2807 (with any kind parameter) gives an undefined result. (Passing other integer 2808 values than 0 and 1 to GCC's @code{_Bool} is also undefined, unless the 2809 integer is explicitly or implicitly casted to @code{_Bool}.) 2810 2811 @node Derived Types and struct 2812 @subsection Derived Types and struct 2813 @cindex C derived type and struct interoperability 2814 @cindex derived type interoperability with C 2815 @cindex interoperability, derived type and struct 2816 2817 For compatibility of derived types with @code{struct}, use 2818 the @code{BIND(C)} attribute in the type declaration. For instance, the 2819 following type declaration 2820 2821 @smallexample 2822 USE ISO_C_BINDING 2823 TYPE, BIND(C) :: myType 2824 INTEGER(C_INT) :: i1, i2 2825 INTEGER(C_SIGNED_CHAR) :: i3 2826 REAL(C_DOUBLE) :: d1 2827 COMPLEX(C_FLOAT_COMPLEX) :: c1 2828 CHARACTER(KIND=C_CHAR) :: str(5) 2829 END TYPE 2830 @end smallexample 2831 2832 @noindent 2833 matches the following @code{struct} declaration in C 2834 2835 @smallexample 2836 struct @{ 2837 int i1, i2; 2838 /* Note: "char" might be signed or unsigned. */ 2839 signed char i3; 2840 double d1; 2841 float _Complex c1; 2842 char str[5]; 2843 @} myType; 2844 @end smallexample 2845 2846 Derived types with the C binding attribute shall not have the @code{sequence} 2847 attribute, type parameters, the @code{extends} attribute, nor type-bound 2848 procedures. Every component must be of interoperable type and kind and may not 2849 have the @code{pointer} or @code{allocatable} attribute. The names of the 2850 components are irrelevant for interoperability. 2851 2852 As there exist no direct Fortran equivalents, neither unions nor structs 2853 with bit field or variable-length array members are interoperable. 2854 2855 @node Interoperable Global Variables 2856 @subsection Interoperable Global Variables 2857 @cindex C variable interoperability 2858 @cindex variable interoperability with C 2859 @cindex interoperability, variable 2860 2861 Variables can be made accessible from C using the C binding attribute, 2862 optionally together with specifying a binding name. Those variables 2863 have to be declared in the declaration part of a @code{MODULE}, 2864 be of interoperable type, and have neither the @code{pointer} nor 2865 the @code{allocatable} attribute. 2866 2867 @smallexample 2868 MODULE m 2869 USE myType_module 2870 USE ISO_C_BINDING 2871 integer(C_INT), bind(C, name="_MyProject_flags") :: global_flag 2872 type(myType), bind(C) :: tp 2873 END MODULE 2874 @end smallexample 2875 2876 Here, @code{_MyProject_flags} is the case-sensitive name of the variable 2877 as seen from C programs while @code{global_flag} is the case-insensitive 2878 name as seen from Fortran. If no binding name is specified, as for 2879 @var{tp}, the C binding name is the (lowercase) Fortran binding name. 2880 If a binding name is specified, only a single variable may be after the 2881 double colon. Note of warning: You cannot use a global variable to 2882 access @var{errno} of the C library as the C standard allows it to be 2883 a macro. Use the @code{IERRNO} intrinsic (GNU extension) instead. 2884 2885 @node Interoperable Subroutines and Functions 2886 @subsection Interoperable Subroutines and Functions 2887 @cindex C procedure interoperability 2888 @cindex procedure interoperability with C 2889 @cindex function interoperability with C 2890 @cindex subroutine interoperability with C 2891 @cindex interoperability, subroutine and function 2892 2893 Subroutines and functions have to have the @code{BIND(C)} attribute to 2894 be compatible with C. The dummy argument declaration is relatively 2895 straightforward. However, one needs to be careful because C uses 2896 call-by-value by default while Fortran behaves usually similar to 2897 call-by-reference. Furthermore, strings and pointers are handled 2898 differently. 2899 2900 To pass a variable by value, use the @code{VALUE} attribute. 2901 Thus, the following C prototype 2902 2903 @smallexample 2904 @code{int func(int i, int *j)} 2905 @end smallexample 2906 2907 @noindent 2908 matches the Fortran declaration 2909 2910 @smallexample 2911 integer(c_int) function func(i,j) 2912 use iso_c_binding, only: c_int 2913 integer(c_int), VALUE :: i 2914 integer(c_int) :: j 2915 @end smallexample 2916 2917 Note that pointer arguments also frequently need the @code{VALUE} attribute, 2918 see @ref{Working with C Pointers}. 2919 2920 Strings are handled quite differently in C and Fortran. In C a string 2921 is a @code{NUL}-terminated array of characters while in Fortran each string 2922 has a length associated with it and is thus not terminated (by e.g. 2923 @code{NUL}). For example, if you want to use the following C function, 2924 2925 @smallexample 2926 #include <stdio.h> 2927 void print_C(char *string) /* equivalent: char string[] */ 2928 @{ 2929 printf("%s\n", string); 2930 @} 2931 @end smallexample 2932 2933 @noindent 2934 to print ``Hello World'' from Fortran, you can call it using 2935 2936 @smallexample 2937 use iso_c_binding, only: C_CHAR, C_NULL_CHAR 2938 interface 2939 subroutine print_c(string) bind(C, name="print_C") 2940 use iso_c_binding, only: c_char 2941 character(kind=c_char) :: string(*) 2942 end subroutine print_c 2943 end interface 2944 call print_c(C_CHAR_"Hello World"//C_NULL_CHAR) 2945 @end smallexample 2946 2947 As the example shows, you need to ensure that the 2948 string is @code{NUL} terminated. Additionally, the dummy argument 2949 @var{string} of @code{print_C} is a length-one assumed-size 2950 array; using @code{character(len=*)} is not allowed. The example 2951 above uses @code{c_char_"Hello World"} to ensure the string 2952 literal has the right type; typically the default character 2953 kind and @code{c_char} are the same and thus @code{"Hello World"} 2954 is equivalent. However, the standard does not guarantee this. 2955 2956 The use of strings is now further illustrated using the C library 2957 function @code{strncpy}, whose prototype is 2958 2959 @smallexample 2960 char *strncpy(char *restrict s1, const char *restrict s2, size_t n); 2961 @end smallexample 2962 2963 @noindent 2964 The function @code{strncpy} copies at most @var{n} characters from 2965 string @var{s2} to @var{s1} and returns @var{s1}. In the following 2966 example, we ignore the return value: 2967 2968 @smallexample 2969 use iso_c_binding 2970 implicit none 2971 character(len=30) :: str,str2 2972 interface 2973 ! Ignore the return value of strncpy -> subroutine 2974 ! "restrict" is always assumed if we do not pass a pointer 2975 subroutine strncpy(dest, src, n) bind(C) 2976 import 2977 character(kind=c_char), intent(out) :: dest(*) 2978 character(kind=c_char), intent(in) :: src(*) 2979 integer(c_size_t), value, intent(in) :: n 2980 end subroutine strncpy 2981 end interface 2982 str = repeat('X',30) ! Initialize whole string with 'X' 2983 call strncpy(str, c_char_"Hello World"//C_NULL_CHAR, & 2984 len(c_char_"Hello World",kind=c_size_t)) 2985 print '(a)', str ! prints: "Hello WorldXXXXXXXXXXXXXXXXXXX" 2986 end 2987 @end smallexample 2988 2989 The intrinsic procedures are described in @ref{Intrinsic Procedures}. 2990 2991 @node Working with C Pointers 2992 @subsection Working with C Pointers 2993 @cindex C pointers 2994 @cindex pointers, C 2995 2996 C pointers are represented in Fortran via the special opaque derived 2997 type @code{type(c_ptr)} (with private components). C pointers are distinct 2998 from Fortran objects with the @code{POINTER} attribute. Thus one needs to 2999 use intrinsic conversion procedures to convert from or to C pointers. 3000 For some applications, using an assumed type (@code{TYPE(*)}) can be 3001 an alternative to a C pointer, and you can also use library routines 3002 to access Fortran pointers from C. See @ref{Further Interoperability 3003 of Fortran with C}. 3004 3005 Here is an example of using C pointers in Fortran: 3006 3007 @smallexample 3008 use iso_c_binding 3009 type(c_ptr) :: cptr1, cptr2 3010 integer, target :: array(7), scalar 3011 integer, pointer :: pa(:), ps 3012 cptr1 = c_loc(array(1)) ! The programmer needs to ensure that the 3013 ! array is contiguous if required by the C 3014 ! procedure 3015 cptr2 = c_loc(scalar) 3016 call c_f_pointer(cptr2, ps) 3017 call c_f_pointer(cptr2, pa, shape=[7]) 3018 @end smallexample 3019 3020 When converting C to Fortran arrays, the one-dimensional @code{SHAPE} argument 3021 has to be passed. 3022 3023 If a pointer is a dummy argument of an interoperable procedure, it usually 3024 has to be declared using the @code{VALUE} attribute. @code{void*} 3025 matches @code{TYPE(C_PTR), VALUE}, while @code{TYPE(C_PTR)} alone 3026 matches @code{void**}. 3027 3028 Procedure pointers are handled analogously to pointers; the C type is 3029 @code{TYPE(C_FUNPTR)} and the intrinsic conversion procedures are 3030 @code{C_F_PROCPOINTER} and @code{C_FUNLOC}. 3031 3032 Let us consider two examples of actually passing a procedure pointer from 3033 C to Fortran and vice versa. Note that these examples are also very 3034 similar to passing ordinary pointers between both languages. First, 3035 consider this code in C: 3036 3037 @smallexample 3038 /* Procedure implemented in Fortran. */ 3039 void get_values (void (*)(double)); 3040 3041 /* Call-back routine we want called from Fortran. */ 3042 void 3043 print_it (double x) 3044 @{ 3045 printf ("Number is %f.\n", x); 3046 @} 3047 3048 /* Call Fortran routine and pass call-back to it. */ 3049 void 3050 foobar () 3051 @{ 3052 get_values (&print_it); 3053 @} 3054 @end smallexample 3055 3056 A matching implementation for @code{get_values} in Fortran, that correctly 3057 receives the procedure pointer from C and is able to call it, is given 3058 in the following @code{MODULE}: 3059 3060 @smallexample 3061 MODULE m 3062 IMPLICIT NONE 3063 3064 ! Define interface of call-back routine. 3065 ABSTRACT INTERFACE 3066 SUBROUTINE callback (x) 3067 USE, INTRINSIC :: ISO_C_BINDING 3068 REAL(KIND=C_DOUBLE), INTENT(IN), VALUE :: x 3069 END SUBROUTINE callback 3070 END INTERFACE 3071 3072 CONTAINS 3073 3074 ! Define C-bound procedure. 3075 SUBROUTINE get_values (cproc) BIND(C) 3076 USE, INTRINSIC :: ISO_C_BINDING 3077 TYPE(C_FUNPTR), INTENT(IN), VALUE :: cproc 3078 3079 PROCEDURE(callback), POINTER :: proc 3080 3081 ! Convert C to Fortran procedure pointer. 3082 CALL C_F_PROCPOINTER (cproc, proc) 3083 3084 ! Call it. 3085 CALL proc (1.0_C_DOUBLE) 3086 CALL proc (-42.0_C_DOUBLE) 3087 CALL proc (18.12_C_DOUBLE) 3088 END SUBROUTINE get_values 3089 3090 END MODULE m 3091 @end smallexample 3092 3093 Next, we want to call a C routine that expects a procedure pointer argument 3094 and pass it a Fortran procedure (which clearly must be interoperable!). 3095 Again, the C function may be: 3096 3097 @smallexample 3098 int 3099 call_it (int (*func)(int), int arg) 3100 @{ 3101 return func (arg); 3102 @} 3103 @end smallexample 3104 3105 It can be used as in the following Fortran code: 3106 3107 @smallexample 3108 MODULE m 3109 USE, INTRINSIC :: ISO_C_BINDING 3110 IMPLICIT NONE 3111 3112 ! Define interface of C function. 3113 INTERFACE 3114 INTEGER(KIND=C_INT) FUNCTION call_it (func, arg) BIND(C) 3115 USE, INTRINSIC :: ISO_C_BINDING 3116 TYPE(C_FUNPTR), INTENT(IN), VALUE :: func 3117 INTEGER(KIND=C_INT), INTENT(IN), VALUE :: arg 3118 END FUNCTION call_it 3119 END INTERFACE 3120 3121 CONTAINS 3122 3123 ! Define procedure passed to C function. 3124 ! It must be interoperable! 3125 INTEGER(KIND=C_INT) FUNCTION double_it (arg) BIND(C) 3126 INTEGER(KIND=C_INT), INTENT(IN), VALUE :: arg 3127 double_it = arg + arg 3128 END FUNCTION double_it 3129 3130 ! Call C function. 3131 SUBROUTINE foobar () 3132 TYPE(C_FUNPTR) :: cproc 3133 INTEGER(KIND=C_INT) :: i 3134 3135 ! Get C procedure pointer. 3136 cproc = C_FUNLOC (double_it) 3137 3138 ! Use it. 3139 DO i = 1_C_INT, 10_C_INT 3140 PRINT *, call_it (cproc, i) 3141 END DO 3142 END SUBROUTINE foobar 3143 3144 END MODULE m 3145 @end smallexample 3146 3147 @node Further Interoperability of Fortran with C 3148 @subsection Further Interoperability of Fortran with C 3149 @cindex Further Interoperability of Fortran with C 3150 @cindex TS 29113 3151 @cindex array descriptor 3152 @cindex dope vector 3153 @cindex assumed-type 3154 @cindex assumed-rank 3155 3156 GNU Fortran implements the Technical Specification ISO/IEC TS 3157 29113:2012, which extends the interoperability support of Fortran 2003 3158 and Fortran 2008 and is now part of the 2018 Fortran standard. 3159 Besides removing some restrictions and constraints, the Technical 3160 Specification adds assumed-type (@code{TYPE(*)}) and assumed-rank 3161 (@code{DIMENSION(..)}) variables and allows for interoperability of 3162 assumed-shape, assumed-rank, and deferred-shape arrays, as well as 3163 allocatables and pointers. Objects of these types are passed to 3164 @code{BIND(C)} functions as descriptors with a standard interface, 3165 declared in the header file @code{<ISO_Fortran_binding.h>}. 3166 3167 Note: Currently, GNU Fortran does not use internally the array descriptor 3168 (dope vector) as specified in the Technical Specification, but uses 3169 an array descriptor with different fields in functions without the 3170 @code{BIND(C)} attribute. Arguments to functions marked @code{BIND(C)} 3171 are converted to the specified form. If you need to access GNU Fortran's 3172 internal array descriptor, you can use the Chasm Language Interoperability 3173 Tools, @url{http://chasm-interop.sourceforge.net/}. 3174 3175 @node GNU Fortran Compiler Directives 3176 @section GNU Fortran Compiler Directives 3177 3178 @menu 3179 * ATTRIBUTES directive:: 3180 * UNROLL directive:: 3181 * BUILTIN directive:: 3182 * IVDEP directive:: 3183 * VECTOR directive:: 3184 * NOVECTOR directive:: 3185 @end menu 3186 3187 @node ATTRIBUTES directive 3188 @subsection ATTRIBUTES directive 3189 3190 The Fortran standard describes how a conforming program shall 3191 behave; however, the exact implementation is not standardized. In order 3192 to allow the user to choose specific implementation details, compiler 3193 directives can be used to set attributes of variables and procedures 3194 which are not part of the standard. Whether a given attribute is 3195 supported and its exact effects depend on both the operating system and 3196 on the processor; see 3197 @ref{Top,,C Extensions,gcc,Using the GNU Compiler Collection (GCC)} 3198 for details. 3199 3200 For procedures and procedure pointers, the following attributes can 3201 be used to change the calling convention: 3202 3203 @itemize 3204 @item @code{CDECL} -- standard C calling convention 3205 @item @code{STDCALL} -- convention where the called procedure pops the stack 3206 @item @code{FASTCALL} -- part of the arguments are passed via registers 3207 instead using the stack 3208 @end itemize 3209 3210 Besides changing the calling convention, the attributes also influence 3211 the decoration of the symbol name, e.g., by a leading underscore or by 3212 a trailing at-sign followed by the number of bytes on the stack. When 3213 assigning a procedure to a procedure pointer, both should use the same 3214 calling convention. 3215 3216 On some systems, procedures and global variables (module variables and 3217 @code{COMMON} blocks) need special handling to be accessible when they 3218 are in a shared library. The following attributes are available: 3219 3220 @itemize 3221 @item @code{DLLEXPORT} -- provide a global pointer to a pointer in the DLL 3222 @item @code{DLLIMPORT} -- reference the function or variable using a 3223 global pointer 3224 @end itemize 3225 3226 For dummy arguments, the @code{NO_ARG_CHECK} attribute can be used; in 3227 other compilers, it is also known as @code{IGNORE_TKR}. For dummy arguments 3228 with this attribute actual arguments of any type and kind (similar to 3229 @code{TYPE(*)}), scalars and arrays of any rank (no equivalent 3230 in Fortran standard) are accepted. As with @code{TYPE(*)}, the argument 3231 is unlimited polymorphic and no type information is available. 3232 Additionally, the argument may only be passed to dummy arguments 3233 with the @code{NO_ARG_CHECK} attribute and as argument to the 3234 @code{PRESENT} intrinsic function and to @code{C_LOC} of the 3235 @code{ISO_C_BINDING} module. 3236 3237 Variables with @code{NO_ARG_CHECK} attribute shall be of assumed-type 3238 (@code{TYPE(*)}; recommended) or of type @code{INTEGER}, @code{LOGICAL}, 3239 @code{REAL} or @code{COMPLEX}. They shall not have the @code{ALLOCATE}, 3240 @code{CODIMENSION}, @code{INTENT(OUT)}, @code{POINTER} or @code{VALUE} 3241 attribute; furthermore, they shall be either scalar or of assumed-size 3242 (@code{dimension(*)}). As @code{TYPE(*)}, the @code{NO_ARG_CHECK} attribute 3243 requires an explicit interface. 3244 3245 @itemize 3246 @item @code{NO_ARG_CHECK} -- disable the type, kind and rank checking 3247 @item @code{DEPRECATED} -- print a warning when using a such-tagged 3248 deprecated procedure, variable or parameter; the warning can be suppressed 3249 with @option{-Wno-deprecated-declarations}. 3250 @end itemize 3251 3252 3253 The attributes are specified using the syntax 3254 3255 @code{!GCC$ ATTRIBUTES} @var{attribute-list} @code{::} @var{variable-list} 3256 3257 where in free-form source code only whitespace is allowed before @code{!GCC$} 3258 and in fixed-form source code @code{!GCC$}, @code{cGCC$} or @code{*GCC$} shall 3259 start in the first column. 3260 3261 For procedures, the compiler directives shall be placed into the body 3262 of the procedure; for variables and procedure pointers, they shall be in 3263 the same declaration part as the variable or procedure pointer. 3264 3265 3266 @node UNROLL directive 3267 @subsection UNROLL directive 3268 3269 The syntax of the directive is 3270 3271 @code{!GCC$ unroll N} 3272 3273 You can use this directive to control how many times a loop should be unrolled. 3274 It must be placed immediately before a @code{DO} loop and applies only to the 3275 loop that follows. N is an integer constant specifying the unrolling factor. 3276 The values of 0 and 1 block any unrolling of the loop. 3277 3278 3279 @node BUILTIN directive 3280 @subsection BUILTIN directive 3281 3282 The syntax of the directive is 3283 3284 @code{!GCC$ BUILTIN (B) attributes simd FLAGS IF('target')} 3285 3286 You can use this directive to define which middle-end built-ins provide vector 3287 implementations. @code{B} is name of the middle-end built-in. @code{FLAGS} 3288 are optional and must be either "(inbranch)" or "(notinbranch)". 3289 @code{IF} statement is optional and is used to filter multilib ABIs 3290 for the built-in that should be vectorized. Example usage: 3291 3292 @smallexample 3293 !GCC$ builtin (sinf) attributes simd (notinbranch) if('x86_64') 3294 @end smallexample 3295 3296 The purpose of the directive is to provide an API among the GCC compiler and 3297 the GNU C Library which would define vector implementations of math routines. 3298 3299 3300 @node IVDEP directive 3301 @subsection IVDEP directive 3302 3303 The syntax of the directive is 3304 3305 @code{!GCC$ ivdep} 3306 3307 This directive tells the compiler to ignore vector dependencies in the 3308 following loop. It must be placed immediately before a @code{DO} loop 3309 and applies only to the loop that follows. 3310 3311 Sometimes the compiler may not have sufficient information to decide 3312 whether a particular loop is vectorizable due to potential 3313 dependencies between iterations. The purpose of the directive is to 3314 tell the compiler that vectorization is safe. 3315 3316 This directive is intended for annotation of existing code. For new 3317 code it is recommended to consider OpenMP SIMD directives as potential 3318 alternative. 3319 3320 3321 @node VECTOR directive 3322 @subsection VECTOR directive 3323 3324 The syntax of the directive is 3325 3326 @code{!GCC$ vector} 3327 3328 This directive tells the compiler to vectorize the following loop. It 3329 must be placed immediately before a @code{DO} loop and applies only to 3330 the loop that follows. 3331 3332 3333 @node NOVECTOR directive 3334 @subsection NOVECTOR directive 3335 3336 The syntax of the directive is 3337 3338 @code{!GCC$ novector} 3339 3340 This directive tells the compiler to not vectorize the following loop. 3341 It must be placed immediately before a @code{DO} loop and applies only 3342 to the loop that follows. 3343 3344 3345 @node Non-Fortran Main Program 3346 @section Non-Fortran Main Program 3347 3348 @menu 3349 * _gfortran_set_args:: Save command-line arguments 3350 * _gfortran_set_options:: Set library option flags 3351 * _gfortran_set_convert:: Set endian conversion 3352 * _gfortran_set_record_marker:: Set length of record markers 3353 * _gfortran_set_fpe:: Set when a Floating Point Exception should be raised 3354 * _gfortran_set_max_subrecord_length:: Set subrecord length 3355 @end menu 3356 3357 Even if you are doing mixed-language programming, it is very 3358 likely that you do not need to know or use the information in this 3359 section. Since it is about the internal structure of GNU Fortran, 3360 it may also change in GCC minor releases. 3361 3362 When you compile a @code{PROGRAM} with GNU Fortran, a function 3363 with the name @code{main} (in the symbol table of the object file) 3364 is generated, which initializes the libgfortran library and then 3365 calls the actual program which uses the name @code{MAIN__}, for 3366 historic reasons. If you link GNU Fortran compiled procedures 3367 to, e.g., a C or C++ program or to a Fortran program compiled by 3368 a different compiler, the libgfortran library is not initialized 3369 and thus a few intrinsic procedures do not work properly, e.g. 3370 those for obtaining the command-line arguments. 3371 3372 Therefore, if your @code{PROGRAM} is not compiled with 3373 GNU Fortran and the GNU Fortran compiled procedures require 3374 intrinsics relying on the library initialization, you need to 3375 initialize the library yourself. Using the default options, 3376 gfortran calls @code{_gfortran_set_args} and 3377 @code{_gfortran_set_options}. The initialization of the former 3378 is needed if the called procedures access the command line 3379 (and for backtracing); the latter sets some flags based on the 3380 standard chosen or to enable backtracing. In typical programs, 3381 it is not necessary to call any initialization function. 3382 3383 If your @code{PROGRAM} is compiled with GNU Fortran, you shall 3384 not call any of the following functions. The libgfortran 3385 initialization functions are shown in C syntax but using C 3386 bindings they are also accessible from Fortran. 3387 3388 3389 @node _gfortran_set_args 3390 @subsection @code{_gfortran_set_args} --- Save command-line arguments 3391 @fnindex _gfortran_set_args 3392 @cindex libgfortran initialization, set_args 3393 3394 @table @asis 3395 @item @emph{Description}: 3396 @code{_gfortran_set_args} saves the command-line arguments; this 3397 initialization is required if any of the command-line intrinsics 3398 is called. Additionally, it shall be called if backtracing is 3399 enabled (see @code{_gfortran_set_options}). 3400 3401 @item @emph{Syntax}: 3402 @code{void _gfortran_set_args (int argc, char *argv[])} 3403 3404 @item @emph{Arguments}: 3405 @multitable @columnfractions .15 .70 3406 @item @var{argc} @tab number of command line argument strings 3407 @item @var{argv} @tab the command-line argument strings; argv[0] 3408 is the pathname of the executable itself. 3409 @end multitable 3410 3411 @item @emph{Example}: 3412 @smallexample 3413 int main (int argc, char *argv[]) 3414 @{ 3415 /* Initialize libgfortran. */ 3416 _gfortran_set_args (argc, argv); 3417 return 0; 3418 @} 3419 @end smallexample 3420 @end table 3421 3422 3423 @node _gfortran_set_options 3424 @subsection @code{_gfortran_set_options} --- Set library option flags 3425 @fnindex _gfortran_set_options 3426 @cindex libgfortran initialization, set_options 3427 3428 @table @asis 3429 @item @emph{Description}: 3430 @code{_gfortran_set_options} sets several flags related to the Fortran 3431 standard to be used, whether backtracing should be enabled 3432 and whether range checks should be performed. The syntax allows for 3433 upward compatibility since the number of passed flags is specified; for 3434 non-passed flags, the default value is used. See also 3435 @pxref{Code Gen Options}. Please note that not all flags are actually 3436 used. 3437 3438 @item @emph{Syntax}: 3439 @code{void _gfortran_set_options (int num, int options[])} 3440 3441 @item @emph{Arguments}: 3442 @multitable @columnfractions .15 .70 3443 @item @var{num} @tab number of options passed 3444 @item @var{argv} @tab The list of flag values 3445 @end multitable 3446 3447 @item @emph{option flag list}: 3448 @multitable @columnfractions .15 .70 3449 @item @var{option}[0] @tab Allowed standard; can give run-time errors 3450 if e.g. an input-output edit descriptor is invalid in a given 3451 standard. Possible values are (bitwise or-ed) @code{GFC_STD_F77} (1), 3452 @code{GFC_STD_F95_OBS} (2), @code{GFC_STD_F95_DEL} (4), 3453 @code{GFC_STD_F95} (8), @code{GFC_STD_F2003} (16), @code{GFC_STD_GNU} 3454 (32), @code{GFC_STD_LEGACY} (64), @code{GFC_STD_F2008} (128), 3455 @code{GFC_STD_F2008_OBS} (256), @code{GFC_STD_F2008_TS} (512), 3456 @code{GFC_STD_F2018} (1024), @code{GFC_STD_F2018_OBS} (2048), and 3457 @code{GFC_STD=F2018_DEL} (4096). Default: @code{GFC_STD_F95_OBS | 3458 GFC_STD_F95_DEL | GFC_STD_F95 | GFC_STD_F2003 | GFC_STD_F2008 | 3459 GFC_STD_F2008_TS | GFC_STD_F2008_OBS | GFC_STD_F77 | GFC_STD_F2018 | 3460 GFC_STD_F2018_OBS | GFC_STD_F2018_DEL | GFC_STD_GNU | GFC_STD_LEGACY}. 3461 @item @var{option}[1] @tab Standard-warning flag; prints a warning to 3462 standard error. Default: @code{GFC_STD_F95_DEL | GFC_STD_LEGACY}. 3463 @item @var{option}[2] @tab If non zero, enable pedantic checking. 3464 Default: off. 3465 @item @var{option}[3] @tab Unused. 3466 @item @var{option}[4] @tab If non zero, enable backtracing on run-time 3467 errors. Default: off. (Default in the compiler: on.) 3468 Note: Installs a signal handler and requires command-line 3469 initialization using @code{_gfortran_set_args}. 3470 @item @var{option}[5] @tab If non zero, supports signed zeros. 3471 Default: enabled. 3472 @item @var{option}[6] @tab Enables run-time checking. Possible values 3473 are (bitwise or-ed): GFC_RTCHECK_BOUNDS (1), GFC_RTCHECK_ARRAY_TEMPS (2), 3474 GFC_RTCHECK_RECURSION (4), GFC_RTCHECK_DO (8), GFC_RTCHECK_POINTER (16), 3475 GFC_RTCHECK_MEM (32), GFC_RTCHECK_BITS (64). 3476 Default: disabled. 3477 @item @var{option}[7] @tab Unused. 3478 @item @var{option}[8] @tab Show a warning when invoking @code{STOP} and 3479 @code{ERROR STOP} if a floating-point exception occurred. Possible values 3480 are (bitwise or-ed) @code{GFC_FPE_INVALID} (1), @code{GFC_FPE_DENORMAL} (2), 3481 @code{GFC_FPE_ZERO} (4), @code{GFC_FPE_OVERFLOW} (8), 3482 @code{GFC_FPE_UNDERFLOW} (16), @code{GFC_FPE_INEXACT} (32). Default: None (0). 3483 (Default in the compiler: @code{GFC_FPE_INVALID | GFC_FPE_DENORMAL | 3484 GFC_FPE_ZERO | GFC_FPE_OVERFLOW | GFC_FPE_UNDERFLOW}.) 3485 @end multitable 3486 3487 @item @emph{Example}: 3488 @smallexample 3489 /* Use gfortran 4.9 default options. */ 3490 static int options[] = @{68, 511, 0, 0, 1, 1, 0, 0, 31@}; 3491 _gfortran_set_options (9, &options); 3492 @end smallexample 3493 @end table 3494 3495 3496 @node _gfortran_set_convert 3497 @subsection @code{_gfortran_set_convert} --- Set endian conversion 3498 @fnindex _gfortran_set_convert 3499 @cindex libgfortran initialization, set_convert 3500 3501 @table @asis 3502 @item @emph{Description}: 3503 @code{_gfortran_set_convert} set the representation of data for 3504 unformatted files. 3505 3506 @item @emph{Syntax}: 3507 @code{void _gfortran_set_convert (int conv)} 3508 3509 @item @emph{Arguments}: 3510 @multitable @columnfractions .15 .70 3511 @item @var{conv} @tab Endian conversion, possible values: 3512 GFC_CONVERT_NATIVE (0, default), GFC_CONVERT_SWAP (1), 3513 GFC_CONVERT_BIG (2), GFC_CONVERT_LITTLE (3). 3514 @end multitable 3515 3516 @item @emph{Example}: 3517 @smallexample 3518 int main (int argc, char *argv[]) 3519 @{ 3520 /* Initialize libgfortran. */ 3521 _gfortran_set_args (argc, argv); 3522 _gfortran_set_convert (1); 3523 return 0; 3524 @} 3525 @end smallexample 3526 @end table 3527 3528 3529 @node _gfortran_set_record_marker 3530 @subsection @code{_gfortran_set_record_marker} --- Set length of record markers 3531 @fnindex _gfortran_set_record_marker 3532 @cindex libgfortran initialization, set_record_marker 3533 3534 @table @asis 3535 @item @emph{Description}: 3536 @code{_gfortran_set_record_marker} sets the length of record markers 3537 for unformatted files. 3538 3539 @item @emph{Syntax}: 3540 @code{void _gfortran_set_record_marker (int val)} 3541 3542 @item @emph{Arguments}: 3543 @multitable @columnfractions .15 .70 3544 @item @var{val} @tab Length of the record marker; valid values 3545 are 4 and 8. Default is 4. 3546 @end multitable 3547 3548 @item @emph{Example}: 3549 @smallexample 3550 int main (int argc, char *argv[]) 3551 @{ 3552 /* Initialize libgfortran. */ 3553 _gfortran_set_args (argc, argv); 3554 _gfortran_set_record_marker (8); 3555 return 0; 3556 @} 3557 @end smallexample 3558 @end table 3559 3560 3561 @node _gfortran_set_fpe 3562 @subsection @code{_gfortran_set_fpe} --- Enable floating point exception traps 3563 @fnindex _gfortran_set_fpe 3564 @cindex libgfortran initialization, set_fpe 3565 3566 @table @asis 3567 @item @emph{Description}: 3568 @code{_gfortran_set_fpe} enables floating point exception traps for 3569 the specified exceptions. On most systems, this will result in a 3570 SIGFPE signal being sent and the program being aborted. 3571 3572 @item @emph{Syntax}: 3573 @code{void _gfortran_set_fpe (int val)} 3574 3575 @item @emph{Arguments}: 3576 @multitable @columnfractions .15 .70 3577 @item @var{option}[0] @tab IEEE exceptions. Possible values are 3578 (bitwise or-ed) zero (0, default) no trapping, 3579 @code{GFC_FPE_INVALID} (1), @code{GFC_FPE_DENORMAL} (2), 3580 @code{GFC_FPE_ZERO} (4), @code{GFC_FPE_OVERFLOW} (8), 3581 @code{GFC_FPE_UNDERFLOW} (16), and @code{GFC_FPE_INEXACT} (32). 3582 @end multitable 3583 3584 @item @emph{Example}: 3585 @smallexample 3586 int main (int argc, char *argv[]) 3587 @{ 3588 /* Initialize libgfortran. */ 3589 _gfortran_set_args (argc, argv); 3590 /* FPE for invalid operations such as SQRT(-1.0). */ 3591 _gfortran_set_fpe (1); 3592 return 0; 3593 @} 3594 @end smallexample 3595 @end table 3596 3597 3598 @node _gfortran_set_max_subrecord_length 3599 @subsection @code{_gfortran_set_max_subrecord_length} --- Set subrecord length 3600 @fnindex _gfortran_set_max_subrecord_length 3601 @cindex libgfortran initialization, set_max_subrecord_length 3602 3603 @table @asis 3604 @item @emph{Description}: 3605 @code{_gfortran_set_max_subrecord_length} set the maximum length 3606 for a subrecord. This option only makes sense for testing and 3607 debugging of unformatted I/O. 3608 3609 @item @emph{Syntax}: 3610 @code{void _gfortran_set_max_subrecord_length (int val)} 3611 3612 @item @emph{Arguments}: 3613 @multitable @columnfractions .15 .70 3614 @item @var{val} @tab the maximum length for a subrecord; 3615 the maximum permitted value is 2147483639, which is also 3616 the default. 3617 @end multitable 3618 3619 @item @emph{Example}: 3620 @smallexample 3621 int main (int argc, char *argv[]) 3622 @{ 3623 /* Initialize libgfortran. */ 3624 _gfortran_set_args (argc, argv); 3625 _gfortran_set_max_subrecord_length (8); 3626 return 0; 3627 @} 3628 @end smallexample 3629 @end table 3630 3631 3632 @node Naming and argument-passing conventions 3633 @section Naming and argument-passing conventions 3634 3635 This section gives an overview about the naming convention of procedures 3636 and global variables and about the argument passing conventions used by 3637 GNU Fortran. If a C binding has been specified, the naming convention 3638 and some of the argument-passing conventions change. If possible, 3639 mixed-language and mixed-compiler projects should use the better defined 3640 C binding for interoperability. See @pxref{Interoperability with C}. 3641 3642 @menu 3643 * Naming conventions:: 3644 * Argument passing conventions:: 3645 @end menu 3646 3647 3648 @node Naming conventions 3649 @subsection Naming conventions 3650 3651 According the Fortran standard, valid Fortran names consist of a letter 3652 between @code{A} to @code{Z}, @code{a} to @code{z}, digits @code{0}, 3653 @code{1} to @code{9} and underscores (@code{_}) with the restriction 3654 that names may only start with a letter. As vendor extension, the 3655 dollar sign (@code{$}) is additionally permitted with the option 3656 @option{-fdollar-ok}, but not as first character and only if the 3657 target system supports it. 3658 3659 By default, the procedure name is the lower-cased Fortran name with an 3660 appended underscore (@code{_}); using @option{-fno-underscoring} no 3661 underscore is appended while @code{-fsecond-underscore} appends two 3662 underscores. Depending on the target system and the calling convention, 3663 the procedure might be additionally dressed; for instance, on 32bit 3664 Windows with @code{stdcall}, an at-sign @code{@@} followed by an integer 3665 number is appended. For the changing the calling convention, see 3666 @pxref{GNU Fortran Compiler Directives}. 3667 3668 For common blocks, the same convention is used, i.e. by default an 3669 underscore is appended to the lower-cased Fortran name. Blank commons 3670 have the name @code{__BLNK__}. 3671 3672 For procedures and variables declared in the specification space of a 3673 module, the name is formed by @code{__}, followed by the lower-cased 3674 module name, @code{_MOD_}, and the lower-cased Fortran name. Note that 3675 no underscore is appended. 3676 3677 3678 @node Argument passing conventions 3679 @subsection Argument passing conventions 3680 3681 Subroutines do not return a value (matching C99's @code{void}) while 3682 functions either return a value as specified in the platform ABI or 3683 the result variable is passed as hidden argument to the function and 3684 no result is returned. A hidden result variable is used when the 3685 result variable is an array or of type @code{CHARACTER}. 3686 3687 Arguments are passed according to the platform ABI. In particular, 3688 complex arguments might not be compatible to a struct with two real 3689 components for the real and imaginary part. The argument passing 3690 matches the one of C99's @code{_Complex}. Functions with scalar 3691 complex result variables return their value and do not use a 3692 by-reference argument. Note that with the @option{-ff2c} option, 3693 the argument passing is modified and no longer completely matches 3694 the platform ABI. Some other Fortran compilers use @code{f2c} 3695 semantic by default; this might cause problems with 3696 interoperablility. 3697 3698 GNU Fortran passes most arguments by reference, i.e. by passing a 3699 pointer to the data. Note that the compiler might use a temporary 3700 variable into which the actual argument has been copied, if required 3701 semantically (copy-in/copy-out). 3702 3703 For arguments with @code{ALLOCATABLE} and @code{POINTER} 3704 attribute (including procedure pointers), a pointer to the pointer 3705 is passed such that the pointer address can be modified in the 3706 procedure. 3707 3708 For dummy arguments with the @code{VALUE} attribute: Scalar arguments 3709 of the type @code{INTEGER}, @code{LOGICAL}, @code{REAL} and 3710 @code{COMPLEX} are passed by value according to the platform ABI. 3711 (As vendor extension and not recommended, using @code{%VAL()} in the 3712 call to a procedure has the same effect.) For @code{TYPE(C_PTR)} and 3713 procedure pointers, the pointer itself is passed such that it can be 3714 modified without affecting the caller. 3715 @c FIXME: Document how VALUE is handled for CHARACTER, TYPE, 3716 @c CLASS and arrays, i.e. whether the copy-in is done in the caller 3717 @c or in the callee. 3718 3719 For Boolean (@code{LOGICAL}) arguments, please note that GCC expects 3720 only the integer value 0 and 1. If a GNU Fortran @code{LOGICAL} 3721 variable contains another integer value, the result is undefined. 3722 As some other Fortran compilers use @math{-1} for @code{.TRUE.}, 3723 extra care has to be taken -- such as passing the value as 3724 @code{INTEGER}. (The same value restriction also applies to other 3725 front ends of GCC, e.g. to GCC's C99 compiler for @code{_Bool} 3726 or GCC's Ada compiler for @code{Boolean}.) 3727 3728 For arguments of @code{CHARACTER} type, the character length is passed 3729 as a hidden argument at the end of the argument list. For 3730 deferred-length strings, the value is passed by reference, otherwise 3731 by value. The character length has the C type @code{size_t} (or 3732 @code{INTEGER(kind=C_SIZE_T)} in Fortran). Note that this is 3733 different to older versions of the GNU Fortran compiler, where the 3734 type of the hidden character length argument was a C @code{int}. In 3735 order to retain compatibility with older versions, one can e.g. for 3736 the following Fortran procedure 3737 3738 @smallexample 3739 subroutine fstrlen (s, a) 3740 character(len=*) :: s 3741 integer :: a 3742 print*, len(s) 3743 end subroutine fstrlen 3744 @end smallexample 3745 3746 define the corresponding C prototype as follows: 3747 3748 @smallexample 3749 #if __GNUC__ > 7 3750 typedef size_t fortran_charlen_t; 3751 #else 3752 typedef int fortran_charlen_t; 3753 #endif 3754 3755 void fstrlen_ (char*, int*, fortran_charlen_t); 3756 @end smallexample 3757 3758 In order to avoid such compiler-specific details, for new code it is 3759 instead recommended to use the ISO_C_BINDING feature. 3760 3761 Note with C binding, @code{CHARACTER(len=1)} result variables are 3762 returned according to the platform ABI and no hidden length argument 3763 is used for dummy arguments; with @code{VALUE}, those variables are 3764 passed by value. 3765 3766 For @code{OPTIONAL} dummy arguments, an absent argument is denoted 3767 by a NULL pointer, except for scalar dummy arguments of type 3768 @code{INTEGER}, @code{LOGICAL}, @code{REAL} and @code{COMPLEX} 3769 which have the @code{VALUE} attribute. For those, a hidden Boolean 3770 argument (@code{logical(kind=C_bool),value}) is used to indicate 3771 whether the argument is present. 3772 3773 Arguments which are assumed-shape, assumed-rank or deferred-rank 3774 arrays or, with @option{-fcoarray=lib}, allocatable scalar coarrays use 3775 an array descriptor. All other arrays pass the address of the 3776 first element of the array. With @option{-fcoarray=lib}, the token 3777 and the offset belonging to nonallocatable coarrays dummy arguments 3778 are passed as hidden argument along the character length hidden 3779 arguments. The token is an opaque pointer identifying the coarray 3780 and the offset is a passed-by-value integer of kind @code{C_PTRDIFF_T}, 3781 denoting the byte offset between the base address of the coarray and 3782 the passed scalar or first element of the passed array. 3783 3784 The arguments are passed in the following order 3785 @itemize @bullet 3786 @item Result variable, when the function result is passed by reference 3787 @item Character length of the function result, if it is a of type 3788 @code{CHARACTER} and no C binding is used 3789 @item The arguments in the order in which they appear in the Fortran 3790 declaration 3791 @item The the present status for optional arguments with value attribute, 3792 which are internally passed by value 3793 @item The character length and/or coarray token and offset for the first 3794 argument which is a @code{CHARACTER} or a nonallocatable coarray dummy 3795 argument, followed by the hidden arguments of the next dummy argument 3796 of such a type 3797 @end itemize 3798 3799 3800 @c --------------------------------------------------------------------- 3801 @c Coarray Programming 3802 @c --------------------------------------------------------------------- 3803 3804 @node Coarray Programming 3805 @chapter Coarray Programming 3806 @cindex Coarrays 3807 3808 @menu 3809 * Type and enum ABI Documentation:: 3810 * Function ABI Documentation:: 3811 @end menu 3812 3813 3814 @node Type and enum ABI Documentation 3815 @section Type and enum ABI Documentation 3816 3817 @menu 3818 * caf_token_t:: 3819 * caf_register_t:: 3820 * caf_deregister_t:: 3821 * caf_reference_t:: 3822 * caf_team_t:: 3823 @end menu 3824 3825 @node caf_token_t 3826 @subsection @code{caf_token_t} 3827 3828 Typedef of type @code{void *} on the compiler side. Can be any data 3829 type on the library side. 3830 3831 @node caf_register_t 3832 @subsection @code{caf_register_t} 3833 3834 Indicates which kind of coarray variable should be registered. 3835 3836 @verbatim 3837 typedef enum caf_register_t { 3838 CAF_REGTYPE_COARRAY_STATIC, 3839 CAF_REGTYPE_COARRAY_ALLOC, 3840 CAF_REGTYPE_LOCK_STATIC, 3841 CAF_REGTYPE_LOCK_ALLOC, 3842 CAF_REGTYPE_CRITICAL, 3843 CAF_REGTYPE_EVENT_STATIC, 3844 CAF_REGTYPE_EVENT_ALLOC, 3845 CAF_REGTYPE_COARRAY_ALLOC_REGISTER_ONLY, 3846 CAF_REGTYPE_COARRAY_ALLOC_ALLOCATE_ONLY 3847 } 3848 caf_register_t; 3849 @end verbatim 3850 3851 The values @code{CAF_REGTYPE_COARRAY_ALLOC_REGISTER_ONLY} and 3852 @code{CAF_REGTYPE_COARRAY_ALLOC_ALLOCATE_ONLY} are for allocatable components 3853 in derived type coarrays only. The first one sets up the token without 3854 allocating memory for allocatable component. The latter one only allocates the 3855 memory for an allocatable component in a derived type coarray. The token 3856 needs to be setup previously by the REGISTER_ONLY. This allows to have 3857 allocatable components un-allocated on some images. The status whether an 3858 allocatable component is allocated on a remote image can be queried by 3859 @code{_caf_is_present} which used internally by the @code{ALLOCATED} 3860 intrinsic. 3861 3862 @node caf_deregister_t 3863 @subsection @code{caf_deregister_t} 3864 3865 @verbatim 3866 typedef enum caf_deregister_t { 3867 CAF_DEREGTYPE_COARRAY_DEREGISTER, 3868 CAF_DEREGTYPE_COARRAY_DEALLOCATE_ONLY 3869 } 3870 caf_deregister_t; 3871 @end verbatim 3872 3873 Allows to specifiy the type of deregistration of a coarray object. The 3874 @code{CAF_DEREGTYPE_COARRAY_DEALLOCATE_ONLY} flag is only allowed for 3875 allocatable components in derived type coarrays. 3876 3877 @node caf_reference_t 3878 @subsection @code{caf_reference_t} 3879 3880 The structure used for implementing arbitrary reference chains. 3881 A @code{CAF_REFERENCE_T} allows to specify a component reference or any kind 3882 of array reference of any rank supported by gfortran. For array references all 3883 kinds as known by the compiler/Fortran standard are supported indicated by 3884 a @code{MODE}. 3885 3886 @verbatim 3887 typedef enum caf_ref_type_t { 3888 /* Reference a component of a derived type, either regular one or an 3889 allocatable or pointer type. For regular ones idx in caf_reference_t is 3890 set to -1. */ 3891 CAF_REF_COMPONENT, 3892 /* Reference an allocatable array. */ 3893 CAF_REF_ARRAY, 3894 /* Reference a non-allocatable/non-pointer array. I.e., the coarray object 3895 has no array descriptor associated and the addressing is done 3896 completely using the ref. */ 3897 CAF_REF_STATIC_ARRAY 3898 } caf_ref_type_t; 3899 @end verbatim 3900 3901 @verbatim 3902 typedef enum caf_array_ref_t { 3903 /* No array ref. This terminates the array ref. */ 3904 CAF_ARR_REF_NONE = 0, 3905 /* Reference array elements given by a vector. Only for this mode 3906 caf_reference_t.u.a.dim[i].v is valid. */ 3907 CAF_ARR_REF_VECTOR, 3908 /* A full array ref (:). */ 3909 CAF_ARR_REF_FULL, 3910 /* Reference a range on elements given by start, end and stride. */ 3911 CAF_ARR_REF_RANGE, 3912 /* Only a single item is referenced given in the start member. */ 3913 CAF_ARR_REF_SINGLE, 3914 /* An array ref of the kind (i:), where i is an arbitrary valid index in the 3915 array. The index i is given in the start member. */ 3916 CAF_ARR_REF_OPEN_END, 3917 /* An array ref of the kind (:i), where the lower bound of the array ref 3918 is given by the remote side. The index i is given in the end member. */ 3919 CAF_ARR_REF_OPEN_START 3920 } caf_array_ref_t; 3921 @end verbatim 3922 3923 @verbatim 3924 /* References to remote components of a derived type. */ 3925 typedef struct caf_reference_t { 3926 /* A pointer to the next ref or NULL. */ 3927 struct caf_reference_t *next; 3928 /* The type of the reference. */ 3929 /* caf_ref_type_t, replaced by int to allow specification in fortran FE. */ 3930 int type; 3931 /* The size of an item referenced in bytes. I.e. in an array ref this is 3932 the factor to advance the array pointer with to get to the next item. 3933 For component refs this gives just the size of the element referenced. */ 3934 size_t item_size; 3935 union { 3936 struct { 3937 /* The offset (in bytes) of the component in the derived type. 3938 Unused for allocatable or pointer components. */ 3939 ptrdiff_t offset; 3940 /* The offset (in bytes) to the caf_token associated with this 3941 component. NULL, when not allocatable/pointer ref. */ 3942 ptrdiff_t caf_token_offset; 3943 } c; 3944 struct { 3945 /* The mode of the array ref. See CAF_ARR_REF_*. */ 3946 /* caf_array_ref_t, replaced by unsigend char to allow specification in 3947 fortran FE. */ 3948 unsigned char mode[GFC_MAX_DIMENSIONS]; 3949 /* The type of a static array. Unset for array's with descriptors. */ 3950 int static_array_type; 3951 /* Subscript refs (s) or vector refs (v). */ 3952 union { 3953 struct { 3954 /* The start and end boundary of the ref and the stride. */ 3955 index_type start, end, stride; 3956 } s; 3957 struct { 3958 /* nvec entries of kind giving the elements to reference. */ 3959 void *vector; 3960 /* The number of entries in vector. */ 3961 size_t nvec; 3962 /* The integer kind used for the elements in vector. */ 3963 int kind; 3964 } v; 3965 } dim[GFC_MAX_DIMENSIONS]; 3966 } a; 3967 } u; 3968 } caf_reference_t; 3969 @end verbatim 3970 3971 The references make up a single linked list of reference operations. The 3972 @code{NEXT} member links to the next reference or NULL to indicate the end of 3973 the chain. Component and array refs can be arbitrarily mixed as long as they 3974 comply to the Fortran standard. 3975 3976 @emph{NOTES} 3977 The member @code{STATIC_ARRAY_TYPE} is used only when the @code{TYPE} is 3978 @code{CAF_REF_STATIC_ARRAY}. The member gives the type of the data referenced. 3979 Because no array descriptor is available for a descriptor-less array and 3980 type conversion still needs to take place the type is transported here. 3981 3982 At the moment @code{CAF_ARR_REF_VECTOR} is not implemented in the front end for 3983 descriptor-less arrays. The library caf_single has untested support for it. 3984 3985 @node caf_team_t 3986 @subsection @code{caf_team_t} 3987 3988 Opaque pointer to represent a team-handle. This type is a stand-in for the 3989 future implementation of teams. It is about to change without further notice. 3990 3991 @node Function ABI Documentation 3992 @section Function ABI Documentation 3993 3994 @menu 3995 * _gfortran_caf_init:: Initialiation function 3996 * _gfortran_caf_finish:: Finalization function 3997 * _gfortran_caf_this_image:: Querying the image number 3998 * _gfortran_caf_num_images:: Querying the maximal number of images 3999 * _gfortran_caf_image_status :: Query the status of an image 4000 * _gfortran_caf_failed_images :: Get an array of the indexes of the failed images 4001 * _gfortran_caf_stopped_images :: Get an array of the indexes of the stopped images 4002 * _gfortran_caf_register:: Registering coarrays 4003 * _gfortran_caf_deregister:: Deregistering coarrays 4004 * _gfortran_caf_is_present:: Query whether an allocatable or pointer component in a derived type coarray is allocated 4005 * _gfortran_caf_send:: Sending data from a local image to a remote image 4006 * _gfortran_caf_get:: Getting data from a remote image 4007 * _gfortran_caf_sendget:: Sending data between remote images 4008 * _gfortran_caf_send_by_ref:: Sending data from a local image to a remote image using enhanced references 4009 * _gfortran_caf_get_by_ref:: Getting data from a remote image using enhanced references 4010 * _gfortran_caf_sendget_by_ref:: Sending data between remote images using enhanced references 4011 * _gfortran_caf_lock:: Locking a lock variable 4012 * _gfortran_caf_unlock:: Unlocking a lock variable 4013 * _gfortran_caf_event_post:: Post an event 4014 * _gfortran_caf_event_wait:: Wait that an event occurred 4015 * _gfortran_caf_event_query:: Query event count 4016 * _gfortran_caf_sync_all:: All-image barrier 4017 * _gfortran_caf_sync_images:: Barrier for selected images 4018 * _gfortran_caf_sync_memory:: Wait for completion of segment-memory operations 4019 * _gfortran_caf_error_stop:: Error termination with exit code 4020 * _gfortran_caf_error_stop_str:: Error termination with string 4021 * _gfortran_caf_fail_image :: Mark the image failed and end its execution 4022 * _gfortran_caf_atomic_define:: Atomic variable assignment 4023 * _gfortran_caf_atomic_ref:: Atomic variable reference 4024 * _gfortran_caf_atomic_cas:: Atomic compare and swap 4025 * _gfortran_caf_atomic_op:: Atomic operation 4026 * _gfortran_caf_co_broadcast:: Sending data to all images 4027 * _gfortran_caf_co_max:: Collective maximum reduction 4028 * _gfortran_caf_co_min:: Collective minimum reduction 4029 * _gfortran_caf_co_sum:: Collective summing reduction 4030 * _gfortran_caf_co_reduce:: Generic collective reduction 4031 @end menu 4032 4033 4034 @node _gfortran_caf_init 4035 @subsection @code{_gfortran_caf_init} --- Initialiation function 4036 @cindex Coarray, _gfortran_caf_init 4037 4038 @table @asis 4039 @item @emph{Description}: 4040 This function is called at startup of the program before the Fortran main 4041 program, if the latter has been compiled with @option{-fcoarray=lib}. 4042 It takes as arguments the command-line arguments of the program. It is 4043 permitted to pass two @code{NULL} pointers as argument; if non-@code{NULL}, 4044 the library is permitted to modify the arguments. 4045 4046 @item @emph{Syntax}: 4047 @code{void _gfortran_caf_init (int *argc, char ***argv)} 4048 4049 @item @emph{Arguments}: 4050 @multitable @columnfractions .15 .70 4051 @item @var{argc} @tab intent(inout) An integer pointer with the number of 4052 arguments passed to the program or @code{NULL}. 4053 @item @var{argv} @tab intent(inout) A pointer to an array of strings with the 4054 command-line arguments or @code{NULL}. 4055 @end multitable 4056 4057 @item @emph{NOTES} 4058 The function is modelled after the initialization function of the Message 4059 Passing Interface (MPI) specification. Due to the way coarray registration 4060 works, it might not be the first call to the library. If the main program is 4061 not written in Fortran and only a library uses coarrays, it can happen that 4062 this function is never called. Therefore, it is recommended that the library 4063 does not rely on the passed arguments and whether the call has been done. 4064 @end table 4065 4066 4067 @node _gfortran_caf_finish 4068 @subsection @code{_gfortran_caf_finish} --- Finalization function 4069 @cindex Coarray, _gfortran_caf_finish 4070 4071 @table @asis 4072 @item @emph{Description}: 4073 This function is called at the end of the Fortran main program, if it has 4074 been compiled with the @option{-fcoarray=lib} option. 4075 4076 @item @emph{Syntax}: 4077 @code{void _gfortran_caf_finish (void)} 4078 4079 @item @emph{NOTES} 4080 For non-Fortran programs, it is recommended to call the function at the end 4081 of the main program. To ensure that the shutdown is also performed for 4082 programs where this function is not explicitly invoked, for instance 4083 non-Fortran programs or calls to the system's exit() function, the library 4084 can use a destructor function. Note that programs can also be terminated 4085 using the STOP and ERROR STOP statements; those use different library calls. 4086 @end table 4087 4088 4089 @node _gfortran_caf_this_image 4090 @subsection @code{_gfortran_caf_this_image} --- Querying the image number 4091 @cindex Coarray, _gfortran_caf_this_image 4092 4093 @table @asis 4094 @item @emph{Description}: 4095 This function returns the current image number, which is a positive number. 4096 4097 @item @emph{Syntax}: 4098 @code{int _gfortran_caf_this_image (int distance)} 4099 4100 @item @emph{Arguments}: 4101 @multitable @columnfractions .15 .70 4102 @item @var{distance} @tab As specified for the @code{this_image} intrinsic 4103 in TS18508. Shall be a non-negative number. 4104 @end multitable 4105 4106 @item @emph{NOTES} 4107 If the Fortran intrinsic @code{this_image} is invoked without an argument, which 4108 is the only permitted form in Fortran 2008, GCC passes @code{0} as 4109 first argument. 4110 @end table 4111 4112 4113 @node _gfortran_caf_num_images 4114 @subsection @code{_gfortran_caf_num_images} --- Querying the maximal number of images 4115 @cindex Coarray, _gfortran_caf_num_images 4116 4117 @table @asis 4118 @item @emph{Description}: 4119 This function returns the number of images in the current team, if 4120 @var{distance} is 0 or the number of images in the parent team at the specified 4121 distance. If failed is -1, the function returns the number of all images at 4122 the specified distance; if it is 0, the function returns the number of 4123 nonfailed images, and if it is 1, it returns the number of failed images. 4124 4125 @item @emph{Syntax}: 4126 @code{int _gfortran_caf_num_images(int distance, int failed)} 4127 4128 @item @emph{Arguments}: 4129 @multitable @columnfractions .15 .70 4130 @item @var{distance} @tab the distance from this image to the ancestor. 4131 Shall be positive. 4132 @item @var{failed} @tab shall be -1, 0, or 1 4133 @end multitable 4134 4135 @item @emph{NOTES} 4136 This function follows TS18508. If the num_image intrinsic has no arguments, 4137 then the compiler passes @code{distance=0} and @code{failed=-1} to the function. 4138 @end table 4139 4140 4141 @node _gfortran_caf_image_status 4142 @subsection @code{_gfortran_caf_image_status} --- Query the status of an image 4143 @cindex Coarray, _gfortran_caf_image_status 4144 4145 @table @asis 4146 @item @emph{Description}: 4147 Get the status of the image given by the id @var{image} of the team given by 4148 @var{team}. Valid results are zero, for image is ok, @code{STAT_STOPPED_IMAGE} 4149 from the ISO_FORTRAN_ENV module to indicate that the image has been stopped and 4150 @code{STAT_FAILED_IMAGE} also from ISO_FORTRAN_ENV to indicate that the image 4151 has executed a @code{FAIL IMAGE} statement. 4152 4153 @item @emph{Syntax}: 4154 @code{int _gfortran_caf_image_status (int image, caf_team_t * team)} 4155 4156 @item @emph{Arguments}: 4157 @multitable @columnfractions .15 .70 4158 @item @var{image} @tab the positive scalar id of the image in the current TEAM. 4159 @item @var{team} @tab optional; team on the which the inquiry is to be 4160 performed. 4161 @end multitable 4162 4163 @item @emph{NOTES} 4164 This function follows TS18508. Because team-functionality is not yet 4165 implemented a null-pointer is passed for the @var{team} argument at the moment. 4166 @end table 4167 4168 4169 @node _gfortran_caf_failed_images 4170 @subsection @code{_gfortran_caf_failed_images} --- Get an array of the indexes of the failed images 4171 @cindex Coarray, _gfortran_caf_failed_images 4172 4173 @table @asis 4174 @item @emph{Description}: 4175 Get an array of image indexes in the current @var{team} that have failed. The 4176 array is sorted ascendingly. When @var{team} is not provided the current team 4177 is to be used. When @var{kind} is provided then the resulting array is of that 4178 integer kind else it is of default integer kind. The returns an unallocated 4179 size zero array when no images have failed. 4180 4181 @item @emph{Syntax}: 4182 @code{int _gfortran_caf_failed_images (caf_team_t * team, int * kind)} 4183 4184 @item @emph{Arguments}: 4185 @multitable @columnfractions .15 .70 4186 @item @var{team} @tab optional; team on the which the inquiry is to be 4187 performed. 4188 @item @var{image} @tab optional; the kind of the resulting integer array. 4189 @end multitable 4190 4191 @item @emph{NOTES} 4192 This function follows TS18508. Because team-functionality is not yet 4193 implemented a null-pointer is passed for the @var{team} argument at the moment. 4194 @end table 4195 4196 4197 @node _gfortran_caf_stopped_images 4198 @subsection @code{_gfortran_caf_stopped_images} --- Get an array of the indexes of the stopped images 4199 @cindex Coarray, _gfortran_caf_stopped_images 4200 4201 @table @asis 4202 @item @emph{Description}: 4203 Get an array of image indexes in the current @var{team} that have stopped. The 4204 array is sorted ascendingly. When @var{team} is not provided the current team 4205 is to be used. When @var{kind} is provided then the resulting array is of that 4206 integer kind else it is of default integer kind. The returns an unallocated 4207 size zero array when no images have failed. 4208 4209 @item @emph{Syntax}: 4210 @code{int _gfortran_caf_stopped_images (caf_team_t * team, int * kind)} 4211 4212 @item @emph{Arguments}: 4213 @multitable @columnfractions .15 .70 4214 @item @var{team} @tab optional; team on the which the inquiry is to be 4215 performed. 4216 @item @var{image} @tab optional; the kind of the resulting integer array. 4217 @end multitable 4218 4219 @item @emph{NOTES} 4220 This function follows TS18508. Because team-functionality is not yet 4221 implemented a null-pointer is passed for the @var{team} argument at the moment. 4222 @end table 4223 4224 4225 @node _gfortran_caf_register 4226 @subsection @code{_gfortran_caf_register} --- Registering coarrays 4227 @cindex Coarray, _gfortran_caf_register 4228 4229 @table @asis 4230 @item @emph{Description}: 4231 Registers memory for a coarray and creates a token to identify the coarray. The 4232 routine is called for both coarrays with @code{SAVE} attribute and using an 4233 explicit @code{ALLOCATE} statement. If an error occurs and @var{STAT} is a 4234 @code{NULL} pointer, the function shall abort with printing an error message 4235 and starting the error termination. If no error occurs and @var{STAT} is 4236 present, it shall be set to zero. Otherwise, it shall be set to a positive 4237 value and, if not-@code{NULL}, @var{ERRMSG} shall be set to a string describing 4238 the failure. The routine shall register the memory provided in the 4239 @code{DATA}-component of the array descriptor @var{DESC}, when that component 4240 is non-@code{NULL}, else it shall allocate sufficient memory and provide a 4241 pointer to it in the @code{DATA}-component of @var{DESC}. The array descriptor 4242 has rank zero, when a scalar object is to be registered and the array 4243 descriptor may be invalid after the call to @code{_gfortran_caf_register}. 4244 When an array is to be allocated the descriptor persists. 4245 4246 For @code{CAF_REGTYPE_COARRAY_STATIC} and @code{CAF_REGTYPE_COARRAY_ALLOC}, 4247 the passed size is the byte size requested. For @code{CAF_REGTYPE_LOCK_STATIC}, 4248 @code{CAF_REGTYPE_LOCK_ALLOC} and @code{CAF_REGTYPE_CRITICAL} it is the array 4249 size or one for a scalar. 4250 4251 When @code{CAF_REGTYPE_COARRAY_ALLOC_REGISTER_ONLY} is used, then only a token 4252 for an allocatable or pointer component is created. The @code{SIZE} parameter 4253 is not used then. On the contrary when 4254 @code{CAF_REGTYPE_COARRAY_ALLOC_ALLOCATE_ONLY} is specified, then the 4255 @var{token} needs to be registered by a previous call with regtype 4256 @code{CAF_REGTYPE_COARRAY_ALLOC_REGISTER_ONLY} and either the memory specified 4257 in the @var{DESC}'s data-ptr is registered or allocate when the data-ptr is 4258 @code{NULL}. 4259 4260 @item @emph{Syntax}: 4261 @code{void caf_register (size_t size, caf_register_t type, caf_token_t *token, 4262 gfc_descriptor_t *desc, int *stat, char *errmsg, size_t errmsg_len)} 4263 4264 @item @emph{Arguments}: 4265 @multitable @columnfractions .15 .70 4266 @item @var{size} @tab For normal coarrays, the byte size of the coarray to be 4267 allocated; for lock types and event types, the number of elements. 4268 @item @var{type} @tab one of the caf_register_t types. 4269 @item @var{token} @tab intent(out) An opaque pointer identifying the coarray. 4270 @item @var{desc} @tab intent(inout) The (pseudo) array descriptor. 4271 @item @var{stat} @tab intent(out) For allocatable coarrays, stores the STAT=; 4272 may be @code{NULL} 4273 @item @var{errmsg} @tab intent(out) When an error occurs, this will be set to 4274 an error message; may be @code{NULL} 4275 @item @var{errmsg_len} @tab the buffer size of errmsg. 4276 @end multitable 4277 4278 @item @emph{NOTES} 4279 Nonallocatable coarrays have to be registered prior use from remote images. 4280 In order to guarantee this, they have to be registered before the main 4281 program. This can be achieved by creating constructor functions. That is what 4282 GCC does such that also for nonallocatable coarrays the memory is allocated and 4283 no static memory is used. The token permits to identify the coarray; to the 4284 processor, the token is a nonaliasing pointer. The library can, for instance, 4285 store the base address of the coarray in the token, some handle or a more 4286 complicated struct. The library may also store the array descriptor 4287 @var{DESC} when its rank is non-zero. 4288 4289 For lock types, the value shall only be used for checking the allocation 4290 status. Note that for critical blocks, the locking is only required on one 4291 image; in the locking statement, the processor shall always pass an 4292 image index of one for critical-block lock variables 4293 (@code{CAF_REGTYPE_CRITICAL}). For lock types and critical-block variables, 4294 the initial value shall be unlocked (or, respectively, not in critical 4295 section) such as the value false; for event types, the initial state should 4296 be no event, e.g. zero. 4297 @end table 4298 4299 4300 @node _gfortran_caf_deregister 4301 @subsection @code{_gfortran_caf_deregister} --- Deregistering coarrays 4302 @cindex Coarray, _gfortran_caf_deregister 4303 4304 @table @asis 4305 @item @emph{Description}: 4306 Called to free or deregister the memory of a coarray; the processor calls this 4307 function for automatic and explicit deallocation. In case of an error, this 4308 function shall fail with an error message, unless the @var{STAT} variable is 4309 not null. The library is only expected to free memory it allocated itself 4310 during a call to @code{_gfortran_caf_register}. 4311 4312 @item @emph{Syntax}: 4313 @code{void caf_deregister (caf_token_t *token, caf_deregister_t type, 4314 int *stat, char *errmsg, size_t errmsg_len)} 4315 4316 @item @emph{Arguments}: 4317 @multitable @columnfractions .15 .70 4318 @item @var{token} @tab the token to free. 4319 @item @var{type} @tab the type of action to take for the coarray. A 4320 @code{CAF_DEREGTYPE_COARRAY_DEALLOCATE_ONLY} is allowed only for allocatable or 4321 pointer components of derived type coarrays. The action only deallocates the 4322 local memory without deleting the token. 4323 @item @var{stat} @tab intent(out) Stores the STAT=; may be NULL 4324 @item @var{errmsg} @tab intent(out) When an error occurs, this will be set 4325 to an error message; may be NULL 4326 @item @var{errmsg_len} @tab the buffer size of errmsg. 4327 @end multitable 4328 4329 @item @emph{NOTES} 4330 For nonalloatable coarrays this function is never called. If a cleanup is 4331 required, it has to be handled via the finish, stop and error stop functions, 4332 and via destructors. 4333 @end table 4334 4335 4336 @node _gfortran_caf_is_present 4337 @subsection @code{_gfortran_caf_is_present} --- Query whether an allocatable or pointer component in a derived type coarray is allocated 4338 @cindex Coarray, _gfortran_caf_is_present 4339 4340 @table @asis 4341 @item @emph{Description}: 4342 Used to query the coarray library whether an allocatable component in a derived 4343 type coarray is allocated on a remote image. 4344 4345 @item @emph{Syntax}: 4346 @code{void _gfortran_caf_is_present (caf_token_t token, int image_index, 4347 gfc_reference_t *ref)} 4348 4349 @item @emph{Arguments}: 4350 @multitable @columnfractions .15 .70 4351 @item @var{token} @tab An opaque pointer identifying the coarray. 4352 @item @var{image_index} @tab The ID of the remote image; must be a positive 4353 number. 4354 @item @var{ref} @tab A chain of references to address the allocatable or 4355 pointer component in the derived type coarray. The object reference needs to be 4356 a scalar or a full array reference, respectively. 4357 @end multitable 4358 4359 @end table 4360 4361 @node _gfortran_caf_send 4362 @subsection @code{_gfortran_caf_send} --- Sending data from a local image to a remote image 4363 @cindex Coarray, _gfortran_caf_send 4364 4365 @table @asis 4366 @item @emph{Description}: 4367 Called to send a scalar, an array section or a whole array from a local 4368 to a remote image identified by the image_index. 4369 4370 @item @emph{Syntax}: 4371 @code{void _gfortran_caf_send (caf_token_t token, size_t offset, 4372 int image_index, gfc_descriptor_t *dest, caf_vector_t *dst_vector, 4373 gfc_descriptor_t *src, int dst_kind, int src_kind, bool may_require_tmp, 4374 int *stat)} 4375 4376 @item @emph{Arguments}: 4377 @multitable @columnfractions .15 .70 4378 @item @var{token} @tab intent(in) An opaque pointer identifying the coarray. 4379 @item @var{offset} @tab intent(in) By which amount of bytes the actual data is 4380 shifted compared to the base address of the coarray. 4381 @item @var{image_index} @tab intent(in) The ID of the remote image; must be a 4382 positive number. 4383 @item @var{dest} @tab intent(in) Array descriptor for the remote image for the 4384 bounds and the size. The @code{base_addr} shall not be accessed. 4385 @item @var{dst_vector} @tab intent(in) If not NULL, it contains the vector 4386 subscript of the destination array; the values are relative to the dimension 4387 triplet of the dest argument. 4388 @item @var{src} @tab intent(in) Array descriptor of the local array to be 4389 transferred to the remote image 4390 @item @var{dst_kind} @tab intent(in) Kind of the destination argument 4391 @item @var{src_kind} @tab intent(in) Kind of the source argument 4392 @item @var{may_require_tmp} @tab intent(in) The variable is @code{false} when 4393 it is known at compile time that the @var{dest} and @var{src} either cannot 4394 overlap or overlap (fully or partially) such that walking @var{src} and 4395 @var{dest} in element wise element order (honoring the stride value) will not 4396 lead to wrong results. Otherwise, the value is @code{true}. 4397 @item @var{stat} @tab intent(out) when non-NULL give the result of the 4398 operation, i.e., zero on success and non-zero on error. When NULL and an error 4399 occurs, then an error message is printed and the program is terminated. 4400 @end multitable 4401 4402 @item @emph{NOTES} 4403 It is permitted to have @var{image_index} equal the current image; the memory 4404 of the send-to and the send-from might (partially) overlap in that case. The 4405 implementation has to take care that it handles this case, e.g. using 4406 @code{memmove} which handles (partially) overlapping memory. If 4407 @var{may_require_tmp} is true, the library might additionally create a 4408 temporary variable, unless additional checks show that this is not required 4409 (e.g. because walking backward is possible or because both arrays are 4410 contiguous and @code{memmove} takes care of overlap issues). 4411 4412 Note that the assignment of a scalar to an array is permitted. In addition, 4413 the library has to handle numeric-type conversion and for strings, padding 4414 and different character kinds. 4415 @end table 4416 4417 4418 @node _gfortran_caf_get 4419 @subsection @code{_gfortran_caf_get} --- Getting data from a remote image 4420 @cindex Coarray, _gfortran_caf_get 4421 4422 @table @asis 4423 @item @emph{Description}: 4424 Called to get an array section or a whole array from a remote, 4425 image identified by the image_index. 4426 4427 @item @emph{Syntax}: 4428 @code{void _gfortran_caf_get (caf_token_t token, size_t offset, 4429 int image_index, gfc_descriptor_t *src, caf_vector_t *src_vector, 4430 gfc_descriptor_t *dest, int src_kind, int dst_kind, bool may_require_tmp, 4431 int *stat)} 4432 4433 @item @emph{Arguments}: 4434 @multitable @columnfractions .15 .70 4435 @item @var{token} @tab intent(in) An opaque pointer identifying the coarray. 4436 @item @var{offset} @tab intent(in) By which amount of bytes the actual data is 4437 shifted compared to the base address of the coarray. 4438 @item @var{image_index} @tab intent(in) The ID of the remote image; must be a 4439 positive number. 4440 @item @var{dest} @tab intent(out) Array descriptor of the local array to store 4441 the data retrieved from the remote image 4442 @item @var{src} @tab intent(in) Array descriptor for the remote image for the 4443 bounds and the size. The @code{base_addr} shall not be accessed. 4444 @item @var{src_vector} @tab intent(in) If not NULL, it contains the vector 4445 subscript of the source array; the values are relative to the dimension 4446 triplet of the @var{src} argument. 4447 @item @var{dst_kind} @tab intent(in) Kind of the destination argument 4448 @item @var{src_kind} @tab intent(in) Kind of the source argument 4449 @item @var{may_require_tmp} @tab intent(in) The variable is @code{false} when 4450 it is known at compile time that the @var{dest} and @var{src} either cannot 4451 overlap or overlap (fully or partially) such that walking @var{src} and 4452 @var{dest} in element wise element order (honoring the stride value) will not 4453 lead to wrong results. Otherwise, the value is @code{true}. 4454 @item @var{stat} @tab intent(out) When non-NULL give the result of the 4455 operation, i.e., zero on success and non-zero on error. When NULL and an error 4456 occurs, then an error message is printed and the program is terminated. 4457 @end multitable 4458 4459 @item @emph{NOTES} 4460 It is permitted to have @var{image_index} equal the current image; the memory of 4461 the send-to and the send-from might (partially) overlap in that case. The 4462 implementation has to take care that it handles this case, e.g. using 4463 @code{memmove} which handles (partially) overlapping memory. If 4464 @var{may_require_tmp} is true, the library might additionally create a 4465 temporary variable, unless additional checks show that this is not required 4466 (e.g. because walking backward is possible or because both arrays are 4467 contiguous and @code{memmove} takes care of overlap issues). 4468 4469 Note that the library has to handle numeric-type conversion and for strings, 4470 padding and different character kinds. 4471 @end table 4472 4473 4474 @node _gfortran_caf_sendget 4475 @subsection @code{_gfortran_caf_sendget} --- Sending data between remote images 4476 @cindex Coarray, _gfortran_caf_sendget 4477 4478 @table @asis 4479 @item @emph{Description}: 4480 Called to send a scalar, an array section or a whole array from a remote image 4481 identified by the @var{src_image_index} to a remote image identified by the 4482 @var{dst_image_index}. 4483 4484 @item @emph{Syntax}: 4485 @code{void _gfortran_caf_sendget (caf_token_t dst_token, size_t dst_offset, 4486 int dst_image_index, gfc_descriptor_t *dest, caf_vector_t *dst_vector, 4487 caf_token_t src_token, size_t src_offset, int src_image_index, 4488 gfc_descriptor_t *src, caf_vector_t *src_vector, int dst_kind, int src_kind, 4489 bool may_require_tmp, int *stat)} 4490 4491 @item @emph{Arguments}: 4492 @multitable @columnfractions .15 .70 4493 @item @var{dst_token} @tab intent(in) An opaque pointer identifying the 4494 destination coarray. 4495 @item @var{dst_offset} @tab intent(in) By which amount of bytes the actual data 4496 is shifted compared to the base address of the destination coarray. 4497 @item @var{dst_image_index} @tab intent(in) The ID of the destination remote 4498 image; must be a positive number. 4499 @item @var{dest} @tab intent(in) Array descriptor for the destination 4500 remote image for the bounds and the size. The @code{base_addr} shall not be 4501 accessed. 4502 @item @var{dst_vector} @tab intent(int) If not NULL, it contains the vector 4503 subscript of the destination array; the values are relative to the dimension 4504 triplet of the @var{dest} argument. 4505 @item @var{src_token} @tab intent(in) An opaque pointer identifying the source 4506 coarray. 4507 @item @var{src_offset} @tab intent(in) By which amount of bytes the actual data 4508 is shifted compared to the base address of the source coarray. 4509 @item @var{src_image_index} @tab intent(in) The ID of the source remote image; 4510 must be a positive number. 4511 @item @var{src} @tab intent(in) Array descriptor of the local array to be 4512 transferred to the remote image. 4513 @item @var{src_vector} @tab intent(in) Array descriptor of the local array to 4514 be transferred to the remote image 4515 @item @var{dst_kind} @tab intent(in) Kind of the destination argument 4516 @item @var{src_kind} @tab intent(in) Kind of the source argument 4517 @item @var{may_require_tmp} @tab intent(in) The variable is @code{false} when 4518 it is known at compile time that the @var{dest} and @var{src} either cannot 4519 overlap or overlap (fully or partially) such that walking @var{src} and 4520 @var{dest} in element wise element order (honoring the stride value) will not 4521 lead to wrong results. Otherwise, the value is @code{true}. 4522 @item @var{stat} @tab intent(out) when non-NULL give the result of the 4523 operation, i.e., zero on success and non-zero on error. When NULL and an error 4524 occurs, then an error message is printed and the program is terminated. 4525 @end multitable 4526 4527 @item @emph{NOTES} 4528 It is permitted to have the same image index for both @var{src_image_index} and 4529 @var{dst_image_index}; the memory of the send-to and the send-from might 4530 (partially) overlap in that case. The implementation has to take care that it 4531 handles this case, e.g. using @code{memmove} which handles (partially) 4532 overlapping memory. If @var{may_require_tmp} is true, the library 4533 might additionally create a temporary variable, unless additional checks show 4534 that this is not required (e.g. because walking backward is possible or because 4535 both arrays are contiguous and @code{memmove} takes care of overlap issues). 4536 4537 Note that the assignment of a scalar to an array is permitted. In addition, 4538 the library has to handle numeric-type conversion and for strings, padding and 4539 different character kinds. 4540 @end table 4541 4542 @node _gfortran_caf_send_by_ref 4543 @subsection @code{_gfortran_caf_send_by_ref} --- Sending data from a local image to a remote image with enhanced referencing options 4544 @cindex Coarray, _gfortran_caf_send_by_ref 4545 4546 @table @asis 4547 @item @emph{Description}: 4548 Called to send a scalar, an array section or a whole array from a local to a 4549 remote image identified by the @var{image_index}. 4550 4551 @item @emph{Syntax}: 4552 @code{void _gfortran_caf_send_by_ref (caf_token_t token, int image_index, 4553 gfc_descriptor_t *src, caf_reference_t *refs, int dst_kind, int src_kind, 4554 bool may_require_tmp, bool dst_reallocatable, int *stat, int dst_type)} 4555 4556 @item @emph{Arguments}: 4557 @multitable @columnfractions .15 .70 4558 @item @var{token} @tab intent(in) An opaque pointer identifying the coarray. 4559 @item @var{image_index} @tab intent(in) The ID of the remote image; must be a 4560 positive number. 4561 @item @var{src} @tab intent(in) Array descriptor of the local array to be 4562 transferred to the remote image 4563 @item @var{refs} @tab intent(in) The references on the remote array to store 4564 the data given by src. Guaranteed to have at least one entry. 4565 @item @var{dst_kind} @tab intent(in) Kind of the destination argument 4566 @item @var{src_kind} @tab intent(in) Kind of the source argument 4567 @item @var{may_require_tmp} @tab intent(in) The variable is @code{false} when 4568 it is known at compile time that the @var{dest} and @var{src} either cannot 4569 overlap or overlap (fully or partially) such that walking @var{src} and 4570 @var{dest} in element wise element order (honoring the stride value) will not 4571 lead to wrong results. Otherwise, the value is @code{true}. 4572 @item @var{dst_reallocatable} @tab intent(in) Set when the destination is of 4573 allocatable or pointer type and the refs will allow reallocation, i.e., the ref 4574 is a full array or component ref. 4575 @item @var{stat} @tab intent(out) When non-@code{NULL} give the result of the 4576 operation, i.e., zero on success and non-zero on error. When @code{NULL} and 4577 an error occurs, then an error message is printed and the program is terminated. 4578 @item @var{dst_type} @tab intent(in) Give the type of the destination. When 4579 the destination is not an array, than the precise type, e.g. of a component in 4580 a derived type, is not known, but provided here. 4581 @end multitable 4582 4583 @item @emph{NOTES} 4584 It is permitted to have @var{image_index} equal the current image; the memory of 4585 the send-to and the send-from might (partially) overlap in that case. The 4586 implementation has to take care that it handles this case, e.g. using 4587 @code{memmove} which handles (partially) overlapping memory. If 4588 @var{may_require_tmp} is true, the library might additionally create a 4589 temporary variable, unless additional checks show that this is not required 4590 (e.g. because walking backward is possible or because both arrays are 4591 contiguous and @code{memmove} takes care of overlap issues). 4592 4593 Note that the assignment of a scalar to an array is permitted. In addition, 4594 the library has to handle numeric-type conversion and for strings, padding 4595 and different character kinds. 4596 4597 Because of the more complicated references possible some operations may be 4598 unsupported by certain libraries. The library is expected to issue a precise 4599 error message why the operation is not permitted. 4600 @end table 4601 4602 4603 @node _gfortran_caf_get_by_ref 4604 @subsection @code{_gfortran_caf_get_by_ref} --- Getting data from a remote image using enhanced references 4605 @cindex Coarray, _gfortran_caf_get_by_ref 4606 4607 @table @asis 4608 @item @emph{Description}: 4609 Called to get a scalar, an array section or a whole array from a remote image 4610 identified by the @var{image_index}. 4611 4612 @item @emph{Syntax}: 4613 @code{void _gfortran_caf_get_by_ref (caf_token_t token, int image_index, 4614 caf_reference_t *refs, gfc_descriptor_t *dst, int dst_kind, int src_kind, 4615 bool may_require_tmp, bool dst_reallocatable, int *stat, int src_type)} 4616 4617 @item @emph{Arguments}: 4618 @multitable @columnfractions .15 .70 4619 @item @var{token} @tab intent(in) An opaque pointer identifying the coarray. 4620 @item @var{image_index} @tab intent(in) The ID of the remote image; must be a 4621 positive number. 4622 @item @var{refs} @tab intent(in) The references to apply to the remote structure 4623 to get the data. 4624 @item @var{dst} @tab intent(in) Array descriptor of the local array to store 4625 the data transferred from the remote image. May be reallocated where needed 4626 and when @var{DST_REALLOCATABLE} allows it. 4627 @item @var{dst_kind} @tab intent(in) Kind of the destination argument 4628 @item @var{src_kind} @tab intent(in) Kind of the source argument 4629 @item @var{may_require_tmp} @tab intent(in) The variable is @code{false} when 4630 it is known at compile time that the @var{dest} and @var{src} either cannot 4631 overlap or overlap (fully or partially) such that walking @var{src} and 4632 @var{dest} in element wise element order (honoring the stride value) will not 4633 lead to wrong results. Otherwise, the value is @code{true}. 4634 @item @var{dst_reallocatable} @tab intent(in) Set when @var{DST} is of 4635 allocatable or pointer type and its refs allow reallocation, i.e., the full 4636 array or a component is referenced. 4637 @item @var{stat} @tab intent(out) When non-@code{NULL} give the result of the 4638 operation, i.e., zero on success and non-zero on error. When @code{NULL} and an 4639 error occurs, then an error message is printed and the program is terminated. 4640 @item @var{src_type} @tab intent(in) Give the type of the source. When the 4641 source is not an array, than the precise type, e.g. of a component in a 4642 derived type, is not known, but provided here. 4643 @end multitable 4644 4645 @item @emph{NOTES} 4646 It is permitted to have @code{image_index} equal the current image; the memory 4647 of the send-to and the send-from might (partially) overlap in that case. The 4648 implementation has to take care that it handles this case, e.g. using 4649 @code{memmove} which handles (partially) overlapping memory. If 4650 @var{may_require_tmp} is true, the library might additionally create a 4651 temporary variable, unless additional checks show that this is not required 4652 (e.g. because walking backward is possible or because both arrays are 4653 contiguous and @code{memmove} takes care of overlap issues). 4654 4655 Note that the library has to handle numeric-type conversion and for strings, 4656 padding and different character kinds. 4657 4658 Because of the more complicated references possible some operations may be 4659 unsupported by certain libraries. The library is expected to issue a precise 4660 error message why the operation is not permitted. 4661 @end table 4662 4663 4664 @node _gfortran_caf_sendget_by_ref 4665 @subsection @code{_gfortran_caf_sendget_by_ref} --- Sending data between remote images using enhanced references on both sides 4666 @cindex Coarray, _gfortran_caf_sendget_by_ref 4667 4668 @table @asis 4669 @item @emph{Description}: 4670 Called to send a scalar, an array section or a whole array from a remote image 4671 identified by the @var{src_image_index} to a remote image identified by the 4672 @var{dst_image_index}. 4673 4674 @item @emph{Syntax}: 4675 @code{void _gfortran_caf_sendget_by_ref (caf_token_t dst_token, 4676 int dst_image_index, caf_reference_t *dst_refs, 4677 caf_token_t src_token, int src_image_index, caf_reference_t *src_refs, 4678 int dst_kind, int src_kind, bool may_require_tmp, int *dst_stat, 4679 int *src_stat, int dst_type, int src_type)} 4680 4681 @item @emph{Arguments}: 4682 @multitable @columnfractions .15 .70 4683 @item @var{dst_token} @tab intent(in) An opaque pointer identifying the 4684 destination coarray. 4685 @item @var{dst_image_index} @tab intent(in) The ID of the destination remote 4686 image; must be a positive number. 4687 @item @var{dst_refs} @tab intent(in) The references on the remote array to store 4688 the data given by the source. Guaranteed to have at least one entry. 4689 @item @var{src_token} @tab intent(in) An opaque pointer identifying the source 4690 coarray. 4691 @item @var{src_image_index} @tab intent(in) The ID of the source remote image; 4692 must be a positive number. 4693 @item @var{src_refs} @tab intent(in) The references to apply to the remote 4694 structure to get the data. 4695 @item @var{dst_kind} @tab intent(in) Kind of the destination argument 4696 @item @var{src_kind} @tab intent(in) Kind of the source argument 4697 @item @var{may_require_tmp} @tab intent(in) The variable is @code{false} when 4698 it is known at compile time that the @var{dest} and @var{src} either cannot 4699 overlap or overlap (fully or partially) such that walking @var{src} and 4700 @var{dest} in element wise element order (honoring the stride value) will not 4701 lead to wrong results. Otherwise, the value is @code{true}. 4702 @item @var{dst_stat} @tab intent(out) when non-@code{NULL} give the result of 4703 the send-operation, i.e., zero on success and non-zero on error. When 4704 @code{NULL} and an error occurs, then an error message is printed and the 4705 program is terminated. 4706 @item @var{src_stat} @tab intent(out) When non-@code{NULL} give the result of 4707 the get-operation, i.e., zero on success and non-zero on error. When 4708 @code{NULL} and an error occurs, then an error message is printed and the 4709 program is terminated. 4710 @item @var{dst_type} @tab intent(in) Give the type of the destination. When 4711 the destination is not an array, than the precise type, e.g. of a component in 4712 a derived type, is not known, but provided here. 4713 @item @var{src_type} @tab intent(in) Give the type of the source. When the 4714 source is not an array, than the precise type, e.g. of a component in a 4715 derived type, is not known, but provided here. 4716 @end multitable 4717 4718 @item @emph{NOTES} 4719 It is permitted to have the same image index for both @var{src_image_index} and 4720 @var{dst_image_index}; the memory of the send-to and the send-from might 4721 (partially) overlap in that case. The implementation has to take care that it 4722 handles this case, e.g. using @code{memmove} which handles (partially) 4723 overlapping memory. If @var{may_require_tmp} is true, the library 4724 might additionally create a temporary variable, unless additional checks show 4725 that this is not required (e.g. because walking backward is possible or because 4726 both arrays are contiguous and @code{memmove} takes care of overlap issues). 4727 4728 Note that the assignment of a scalar to an array is permitted. In addition, 4729 the library has to handle numeric-type conversion and for strings, padding and 4730 different character kinds. 4731 4732 Because of the more complicated references possible some operations may be 4733 unsupported by certain libraries. The library is expected to issue a precise 4734 error message why the operation is not permitted. 4735 @end table 4736 4737 4738 @node _gfortran_caf_lock 4739 @subsection @code{_gfortran_caf_lock} --- Locking a lock variable 4740 @cindex Coarray, _gfortran_caf_lock 4741 4742 @table @asis 4743 @item @emph{Description}: 4744 Acquire a lock on the given image on a scalar locking variable or for the 4745 given array element for an array-valued variable. If the @var{acquired_lock} 4746 is @code{NULL}, the function returns after having obtained the lock. If it is 4747 non-@code{NULL}, then @var{acquired_lock} is assigned the value true (one) when 4748 the lock could be obtained and false (zero) otherwise. Locking a lock variable 4749 which has already been locked by the same image is an error. 4750 4751 @item @emph{Syntax}: 4752 @code{void _gfortran_caf_lock (caf_token_t token, size_t index, int image_index, 4753 int *acquired_lock, int *stat, char *errmsg, size_t errmsg_len)} 4754 4755 @item @emph{Arguments}: 4756 @multitable @columnfractions .15 .70 4757 @item @var{token} @tab intent(in) An opaque pointer identifying the coarray. 4758 @item @var{index} @tab intent(in) Array index; first array index is 0. For 4759 scalars, it is always 0. 4760 @item @var{image_index} @tab intent(in) The ID of the remote image; must be a 4761 positive number. 4762 @item @var{acquired_lock} @tab intent(out) If not NULL, it returns whether lock 4763 could be obtained. 4764 @item @var{stat} @tab intent(out) Stores the STAT=; may be NULL. 4765 @item @var{errmsg} @tab intent(out) When an error occurs, this will be set to 4766 an error message; may be NULL. 4767 @item @var{errmsg_len} @tab intent(in) the buffer size of errmsg 4768 @end multitable 4769 4770 @item @emph{NOTES} 4771 This function is also called for critical blocks; for those, the array index 4772 is always zero and the image index is one. Libraries are permitted to use other 4773 images for critical-block locking variables. 4774 @end table 4775 4776 @node _gfortran_caf_unlock 4777 @subsection @code{_gfortran_caf_lock} --- Unlocking a lock variable 4778 @cindex Coarray, _gfortran_caf_unlock 4779 4780 @table @asis 4781 @item @emph{Description}: 4782 Release a lock on the given image on a scalar locking variable or for the 4783 given array element for an array-valued variable. Unlocking a lock variable 4784 which is unlocked or has been locked by a different image is an error. 4785 4786 @item @emph{Syntax}: 4787 @code{void _gfortran_caf_unlock (caf_token_t token, size_t index, int image_index, 4788 int *stat, char *errmsg, size_t errmsg_len)} 4789 4790 @item @emph{Arguments}: 4791 @multitable @columnfractions .15 .70 4792 @item @var{token} @tab intent(in) An opaque pointer identifying the coarray. 4793 @item @var{index} @tab intent(in) Array index; first array index is 0. For 4794 scalars, it is always 0. 4795 @item @var{image_index} @tab intent(in) The ID of the remote image; must be a 4796 positive number. 4797 @item @var{stat} @tab intent(out) For allocatable coarrays, stores the STAT=; 4798 may be NULL. 4799 @item @var{errmsg} @tab intent(out) When an error occurs, this will be set to 4800 an error message; may be NULL. 4801 @item @var{errmsg_len} @tab intent(in) the buffer size of errmsg 4802 @end multitable 4803 4804 @item @emph{NOTES} 4805 This function is also called for critical block; for those, the array index 4806 is always zero and the image index is one. Libraries are permitted to use other 4807 images for critical-block locking variables. 4808 @end table 4809 4810 @node _gfortran_caf_event_post 4811 @subsection @code{_gfortran_caf_event_post} --- Post an event 4812 @cindex Coarray, _gfortran_caf_event_post 4813 4814 @table @asis 4815 @item @emph{Description}: 4816 Increment the event count of the specified event variable. 4817 4818 @item @emph{Syntax}: 4819 @code{void _gfortran_caf_event_post (caf_token_t token, size_t index, 4820 int image_index, int *stat, char *errmsg, size_t errmsg_len)} 4821 4822 @item @emph{Arguments}: 4823 @multitable @columnfractions .15 .70 4824 @item @var{token} @tab intent(in) An opaque pointer identifying the coarray. 4825 @item @var{index} @tab intent(in) Array index; first array index is 0. For 4826 scalars, it is always 0. 4827 @item @var{image_index} @tab intent(in) The ID of the remote image; must be a 4828 positive number; zero indicates the current image, when accessed noncoindexed. 4829 @item @var{stat} @tab intent(out) Stores the STAT=; may be NULL. 4830 @item @var{errmsg} @tab intent(out) When an error occurs, this will be set to 4831 an error message; may be NULL. 4832 @item @var{errmsg_len} @tab intent(in) the buffer size of errmsg 4833 @end multitable 4834 4835 @item @emph{NOTES} 4836 This acts like an atomic add of one to the remote image's event variable. 4837 The statement is an image-control statement but does not imply sync memory. 4838 Still, all preceeding push communications of this image to the specified 4839 remote image have to be completed before @code{event_wait} on the remote 4840 image returns. 4841 @end table 4842 4843 4844 4845 @node _gfortran_caf_event_wait 4846 @subsection @code{_gfortran_caf_event_wait} --- Wait that an event occurred 4847 @cindex Coarray, _gfortran_caf_event_wait 4848 4849 @table @asis 4850 @item @emph{Description}: 4851 Wait until the event count has reached at least the specified 4852 @var{until_count}; if so, atomically decrement the event variable by this 4853 amount and return. 4854 4855 @item @emph{Syntax}: 4856 @code{void _gfortran_caf_event_wait (caf_token_t token, size_t index, 4857 int until_count, int *stat, char *errmsg, size_t errmsg_len)} 4858 4859 @item @emph{Arguments}: 4860 @multitable @columnfractions .15 .70 4861 @item @var{token} @tab intent(in) An opaque pointer identifying the coarray. 4862 @item @var{index} @tab intent(in) Array index; first array index is 0. For 4863 scalars, it is always 0. 4864 @item @var{until_count} @tab intent(in) The number of events which have to be 4865 available before the function returns. 4866 @item @var{stat} @tab intent(out) Stores the STAT=; may be NULL. 4867 @item @var{errmsg} @tab intent(out) When an error occurs, this will be set to 4868 an error message; may be NULL. 4869 @item @var{errmsg_len} @tab intent(in) the buffer size of errmsg 4870 @end multitable 4871 4872 @item @emph{NOTES} 4873 This function only operates on a local coarray. It acts like a loop checking 4874 atomically the value of the event variable, breaking if the value is greater 4875 or equal the requested number of counts. Before the function returns, the 4876 event variable has to be decremented by the requested @var{until_count} value. 4877 A possible implementation would be a busy loop for a certain number of spins 4878 (possibly depending on the number of threads relative to the number of available 4879 cores) followed by another waiting strategy such as a sleeping wait (possibly 4880 with an increasing number of sleep time) or, if possible, a futex wait. 4881 4882 The statement is an image-control statement but does not imply sync memory. 4883 Still, all preceeding push communications of this image to the specified 4884 remote image have to be completed before @code{event_wait} on the remote 4885 image returns. 4886 @end table 4887 4888 4889 4890 @node _gfortran_caf_event_query 4891 @subsection @code{_gfortran_caf_event_query} --- Query event count 4892 @cindex Coarray, _gfortran_caf_event_query 4893 4894 @table @asis 4895 @item @emph{Description}: 4896 Return the event count of the specified event variable. 4897 4898 @item @emph{Syntax}: 4899 @code{void _gfortran_caf_event_query (caf_token_t token, size_t index, 4900 int image_index, int *count, int *stat)} 4901 4902 @item @emph{Arguments}: 4903 @multitable @columnfractions .15 .70 4904 @item @var{token} @tab intent(in) An opaque pointer identifying the coarray. 4905 @item @var{index} @tab intent(in) Array index; first array index is 0. For 4906 scalars, it is always 0. 4907 @item @var{image_index} @tab intent(in) The ID of the remote image; must be a 4908 positive number; zero indicates the current image when accessed noncoindexed. 4909 @item @var{count} @tab intent(out) The number of events currently posted to 4910 the event variable. 4911 @item @var{stat} @tab intent(out) Stores the STAT=; may be NULL. 4912 @end multitable 4913 4914 @item @emph{NOTES} 4915 The typical use is to check the local event variable to only call 4916 @code{event_wait} when the data is available. However, a coindexed variable 4917 is permitted; there is no ordering or synchronization implied. It acts like 4918 an atomic fetch of the value of the event variable. 4919 @end table 4920 4921 4922 4923 @node _gfortran_caf_sync_all 4924 @subsection @code{_gfortran_caf_sync_all} --- All-image barrier 4925 @cindex Coarray, _gfortran_caf_sync_all 4926 4927 @table @asis 4928 @item @emph{Description}: 4929 Synchronization of all images in the current team; the program only continues 4930 on a given image after this function has been called on all images of the 4931 current team. Additionally, it ensures that all pending data transfers of 4932 previous segment have completed. 4933 4934 @item @emph{Syntax}: 4935 @code{void _gfortran_caf_sync_all (int *stat, char *errmsg, size_t errmsg_len)} 4936 4937 @item @emph{Arguments}: 4938 @multitable @columnfractions .15 .70 4939 @item @var{stat} @tab intent(out) Stores the status STAT= and may be NULL. 4940 @item @var{errmsg} @tab intent(out) When an error occurs, this will be set to 4941 an error message; may be NULL. 4942 @item @var{errmsg_len} @tab intent(in) the buffer size of errmsg 4943 @end multitable 4944 @end table 4945 4946 4947 4948 @node _gfortran_caf_sync_images 4949 @subsection @code{_gfortran_caf_sync_images} --- Barrier for selected images 4950 @cindex Coarray, _gfortran_caf_sync_images 4951 4952 @table @asis 4953 @item @emph{Description}: 4954 Synchronization between the specified images; the program only continues on a 4955 given image after this function has been called on all images specified for 4956 that image. Note that one image can wait for all other images in the current 4957 team (e.g. via @code{sync images(*)}) while those only wait for that specific 4958 image. Additionally, @code{sync images} ensures that all pending data 4959 transfers of previous segments have completed. 4960 4961 @item @emph{Syntax}: 4962 @code{void _gfortran_caf_sync_images (int count, int images[], int *stat, 4963 char *errmsg, size_t errmsg_len)} 4964 4965 @item @emph{Arguments}: 4966 @multitable @columnfractions .15 .70 4967 @item @var{count} @tab intent(in) The number of images which are provided in 4968 the next argument. For a zero-sized array, the value is zero. For 4969 @code{sync images (*)}, the value is @math{-1}. 4970 @item @var{images} @tab intent(in) An array with the images provided by the 4971 user. If @var{count} is zero, a NULL pointer is passed. 4972 @item @var{stat} @tab intent(out) Stores the status STAT= and may be NULL. 4973 @item @var{errmsg} @tab intent(out) When an error occurs, this will be set to 4974 an error message; may be NULL. 4975 @item @var{errmsg_len} @tab intent(in) the buffer size of errmsg 4976 @end multitable 4977 @end table 4978 4979 4980 4981 @node _gfortran_caf_sync_memory 4982 @subsection @code{_gfortran_caf_sync_memory} --- Wait for completion of segment-memory operations 4983 @cindex Coarray, _gfortran_caf_sync_memory 4984 4985 @table @asis 4986 @item @emph{Description}: 4987 Acts as optimization barrier between different segments. It also ensures that 4988 all pending memory operations of this image have been completed. 4989 4990 @item @emph{Syntax}: 4991 @code{void _gfortran_caf_sync_memory (int *stat, char *errmsg, size_t errmsg_len)} 4992 4993 @item @emph{Arguments}: 4994 @multitable @columnfractions .15 .70 4995 @item @var{stat} @tab intent(out) Stores the status STAT= and may be NULL. 4996 @item @var{errmsg} @tab intent(out) When an error occurs, this will be set to 4997 an error message; may be NULL. 4998 @item @var{errmsg_len} @tab intent(in) the buffer size of errmsg 4999 @end multitable 5000 5001 @item @emph{NOTE} A simple implementation could be 5002 @code{__asm__ __volatile__ ("":::"memory")} to prevent code movements. 5003 @end table 5004 5005 5006 5007 @node _gfortran_caf_error_stop 5008 @subsection @code{_gfortran_caf_error_stop} --- Error termination with exit code 5009 @cindex Coarray, _gfortran_caf_error_stop 5010 5011 @table @asis 5012 @item @emph{Description}: 5013 Invoked for an @code{ERROR STOP} statement which has an integer argument. The 5014 function should terminate the program with the specified exit code. 5015 5016 5017 @item @emph{Syntax}: 5018 @code{void _gfortran_caf_error_stop (int error)} 5019 5020 @item @emph{Arguments}: 5021 @multitable @columnfractions .15 .70 5022 @item @var{error} @tab intent(in) The exit status to be used. 5023 @end multitable 5024 @end table 5025 5026 5027 5028 @node _gfortran_caf_error_stop_str 5029 @subsection @code{_gfortran_caf_error_stop_str} --- Error termination with string 5030 @cindex Coarray, _gfortran_caf_error_stop_str 5031 5032 @table @asis 5033 @item @emph{Description}: 5034 Invoked for an @code{ERROR STOP} statement which has a string as argument. The 5035 function should terminate the program with a nonzero-exit code. 5036 5037 @item @emph{Syntax}: 5038 @code{void _gfortran_caf_error_stop (const char *string, size_t len)} 5039 5040 @item @emph{Arguments}: 5041 @multitable @columnfractions .15 .70 5042 @item @var{string} @tab intent(in) the error message (not zero terminated) 5043 @item @var{len} @tab intent(in) the length of the string 5044 @end multitable 5045 @end table 5046 5047 5048 5049 @node _gfortran_caf_fail_image 5050 @subsection @code{_gfortran_caf_fail_image} --- Mark the image failed and end its execution 5051 @cindex Coarray, _gfortran_caf_fail_image 5052 5053 @table @asis 5054 @item @emph{Description}: 5055 Invoked for an @code{FAIL IMAGE} statement. The function should terminate the 5056 current image. 5057 5058 @item @emph{Syntax}: 5059 @code{void _gfortran_caf_fail_image ()} 5060 5061 @item @emph{NOTES} 5062 This function follows TS18508. 5063 @end table 5064 5065 5066 5067 @node _gfortran_caf_atomic_define 5068 @subsection @code{_gfortran_caf_atomic_define} --- Atomic variable assignment 5069 @cindex Coarray, _gfortran_caf_atomic_define 5070 5071 @table @asis 5072 @item @emph{Description}: 5073 Assign atomically a value to an integer or logical variable. 5074 5075 @item @emph{Syntax}: 5076 @code{void _gfortran_caf_atomic_define (caf_token_t token, size_t offset, 5077 int image_index, void *value, int *stat, int type, int kind)} 5078 5079 @item @emph{Arguments}: 5080 @multitable @columnfractions .15 .70 5081 @item @var{token} @tab intent(in) An opaque pointer identifying the coarray. 5082 @item @var{offset} @tab intent(in) By which amount of bytes the actual data is 5083 shifted compared to the base address of the coarray. 5084 @item @var{image_index} @tab intent(in) The ID of the remote image; must be a 5085 positive number; zero indicates the current image when used noncoindexed. 5086 @item @var{value} @tab intent(in) the value to be assigned, passed by reference 5087 @item @var{stat} @tab intent(out) Stores the status STAT= and may be NULL. 5088 @item @var{type} @tab intent(in) The data type, i.e. @code{BT_INTEGER} (1) or 5089 @code{BT_LOGICAL} (2). 5090 @item @var{kind} @tab intent(in) The kind value (only 4; always @code{int}) 5091 @end multitable 5092 @end table 5093 5094 5095 5096 @node _gfortran_caf_atomic_ref 5097 @subsection @code{_gfortran_caf_atomic_ref} --- Atomic variable reference 5098 @cindex Coarray, _gfortran_caf_atomic_ref 5099 5100 @table @asis 5101 @item @emph{Description}: 5102 Reference atomically a value of a kind-4 integer or logical variable. 5103 5104 @item @emph{Syntax}: 5105 @code{void _gfortran_caf_atomic_ref (caf_token_t token, size_t offset, 5106 int image_index, void *value, int *stat, int type, int kind)} 5107 5108 @item @emph{Arguments}: 5109 @multitable @columnfractions .15 .70 5110 @item @var{token} @tab intent(in) An opaque pointer identifying the coarray. 5111 @item @var{offset} @tab intent(in) By which amount of bytes the actual data is 5112 shifted compared to the base address of the coarray. 5113 @item @var{image_index} @tab intent(in) The ID of the remote image; must be a 5114 positive number; zero indicates the current image when used noncoindexed. 5115 @item @var{value} @tab intent(out) The variable assigned the atomically 5116 referenced variable. 5117 @item @var{stat} @tab intent(out) Stores the status STAT= and may be NULL. 5118 @item @var{type} @tab the data type, i.e. @code{BT_INTEGER} (1) or 5119 @code{BT_LOGICAL} (2). 5120 @item @var{kind} @tab The kind value (only 4; always @code{int}) 5121 @end multitable 5122 @end table 5123 5124 5125 5126 @node _gfortran_caf_atomic_cas 5127 @subsection @code{_gfortran_caf_atomic_cas} --- Atomic compare and swap 5128 @cindex Coarray, _gfortran_caf_atomic_cas 5129 5130 @table @asis 5131 @item @emph{Description}: 5132 Atomic compare and swap of a kind-4 integer or logical variable. Assigns 5133 atomically the specified value to the atomic variable, if the latter has 5134 the value specified by the passed condition value. 5135 5136 @item @emph{Syntax}: 5137 @code{void _gfortran_caf_atomic_cas (caf_token_t token, size_t offset, 5138 int image_index, void *old, void *compare, void *new_val, int *stat, 5139 int type, int kind)} 5140 5141 @item @emph{Arguments}: 5142 @multitable @columnfractions .15 .70 5143 @item @var{token} @tab intent(in) An opaque pointer identifying the coarray. 5144 @item @var{offset} @tab intent(in) By which amount of bytes the actual data is 5145 shifted compared to the base address of the coarray. 5146 @item @var{image_index} @tab intent(in) The ID of the remote image; must be a 5147 positive number; zero indicates the current image when used noncoindexed. 5148 @item @var{old} @tab intent(out) The value which the atomic variable had 5149 just before the cas operation. 5150 @item @var{compare} @tab intent(in) The value used for comparision. 5151 @item @var{new_val} @tab intent(in) The new value for the atomic variable, 5152 assigned to the atomic variable, if @code{compare} equals the value of the 5153 atomic variable. 5154 @item @var{stat} @tab intent(out) Stores the status STAT= and may be NULL. 5155 @item @var{type} @tab intent(in) the data type, i.e. @code{BT_INTEGER} (1) or 5156 @code{BT_LOGICAL} (2). 5157 @item @var{kind} @tab intent(in) The kind value (only 4; always @code{int}) 5158 @end multitable 5159 @end table 5160 5161 5162 5163 @node _gfortran_caf_atomic_op 5164 @subsection @code{_gfortran_caf_atomic_op} --- Atomic operation 5165 @cindex Coarray, _gfortran_caf_atomic_op 5166 5167 @table @asis 5168 @item @emph{Description}: 5169 Apply an operation atomically to an atomic integer or logical variable. 5170 After the operation, @var{old} contains the value just before the operation, 5171 which, respectively, adds (GFC_CAF_ATOMIC_ADD) atomically the @code{value} to 5172 the atomic integer variable or does a bitwise AND, OR or exclusive OR 5173 between the atomic variable and @var{value}; the result is then stored in the 5174 atomic variable. 5175 5176 @item @emph{Syntax}: 5177 @code{void _gfortran_caf_atomic_op (int op, caf_token_t token, size_t offset, 5178 int image_index, void *value, void *old, int *stat, int type, int kind)} 5179 5180 @item @emph{Arguments}: 5181 @multitable @columnfractions .15 .70 5182 @item @var{op} @tab intent(in) the operation to be performed; possible values 5183 @code{GFC_CAF_ATOMIC_ADD} (1), @code{GFC_CAF_ATOMIC_AND} (2), 5184 @code{GFC_CAF_ATOMIC_OR} (3), @code{GFC_CAF_ATOMIC_XOR} (4). 5185 @item @var{token} @tab intent(in) An opaque pointer identifying the coarray. 5186 @item @var{offset} @tab intent(in) By which amount of bytes the actual data is 5187 shifted compared to the base address of the coarray. 5188 @item @var{image_index} @tab intent(in) The ID of the remote image; must be a 5189 positive number; zero indicates the current image when used noncoindexed. 5190 @item @var{old} @tab intent(out) The value which the atomic variable had 5191 just before the atomic operation. 5192 @item @var{val} @tab intent(in) The new value for the atomic variable, 5193 assigned to the atomic variable, if @code{compare} equals the value of the 5194 atomic variable. 5195 @item @var{stat} @tab intent(out) Stores the status STAT= and may be NULL. 5196 @item @var{type} @tab intent(in) the data type, i.e. @code{BT_INTEGER} (1) or 5197 @code{BT_LOGICAL} (2) 5198 @item @var{kind} @tab intent(in) the kind value (only 4; always @code{int}) 5199 @end multitable 5200 @end table 5201 5202 5203 5204 5205 @node _gfortran_caf_co_broadcast 5206 @subsection @code{_gfortran_caf_co_broadcast} --- Sending data to all images 5207 @cindex Coarray, _gfortran_caf_co_broadcast 5208 5209 @table @asis 5210 @item @emph{Description}: 5211 Distribute a value from a given image to all other images in the team. Has to 5212 be called collectively. 5213 5214 @item @emph{Syntax}: 5215 @code{void _gfortran_caf_co_broadcast (gfc_descriptor_t *a, 5216 int source_image, int *stat, char *errmsg, size_t errmsg_len)} 5217 5218 @item @emph{Arguments}: 5219 @multitable @columnfractions .15 .70 5220 @item @var{a} @tab intent(inout) An array descriptor with the data to be 5221 broadcasted (on @var{source_image}) or to be received (other images). 5222 @item @var{source_image} @tab intent(in) The ID of the image from which the 5223 data should be broadcasted. 5224 @item @var{stat} @tab intent(out) Stores the status STAT= and may be NULL. 5225 @item @var{errmsg} @tab intent(out) When an error occurs, this will be set to 5226 an error message; may be NULL. 5227 @item @var{errmsg_len} @tab intent(in) the buffer size of errmsg. 5228 @end multitable 5229 @end table 5230 5231 5232 5233 @node _gfortran_caf_co_max 5234 @subsection @code{_gfortran_caf_co_max} --- Collective maximum reduction 5235 @cindex Coarray, _gfortran_caf_co_max 5236 5237 @table @asis 5238 @item @emph{Description}: 5239 Calculates for each array element of the variable @var{a} the maximum 5240 value for that element in the current team; if @var{result_image} has the 5241 value 0, the result shall be stored on all images, otherwise, only on the 5242 specified image. This function operates on numeric values and character 5243 strings. 5244 5245 @item @emph{Syntax}: 5246 @code{void _gfortran_caf_co_max (gfc_descriptor_t *a, int result_image, 5247 int *stat, char *errmsg, int a_len, size_t errmsg_len)} 5248 5249 @item @emph{Arguments}: 5250 @multitable @columnfractions .15 .70 5251 @item @var{a} @tab intent(inout) An array descriptor for the data to be 5252 processed. On the destination image(s) the result overwrites the old content. 5253 @item @var{result_image} @tab intent(in) The ID of the image to which the 5254 reduced value should be copied to; if zero, it has to be copied to all images. 5255 @item @var{stat} @tab intent(out) Stores the status STAT= and may be NULL. 5256 @item @var{errmsg} @tab intent(out) When an error occurs, this will be set to 5257 an error message; may be NULL. 5258 @item @var{a_len} @tab intent(in) the string length of argument @var{a} 5259 @item @var{errmsg_len} @tab intent(in) the buffer size of errmsg 5260 @end multitable 5261 5262 @item @emph{NOTES} 5263 If @var{result_image} is nonzero, the data in the array descriptor @var{a} on 5264 all images except of the specified one become undefined; hence, the library may 5265 make use of this. 5266 @end table 5267 5268 5269 5270 @node _gfortran_caf_co_min 5271 @subsection @code{_gfortran_caf_co_min} --- Collective minimum reduction 5272 @cindex Coarray, _gfortran_caf_co_min 5273 5274 @table @asis 5275 @item @emph{Description}: 5276 Calculates for each array element of the variable @var{a} the minimum 5277 value for that element in the current team; if @var{result_image} has the 5278 value 0, the result shall be stored on all images, otherwise, only on the 5279 specified image. This function operates on numeric values and character 5280 strings. 5281 5282 @item @emph{Syntax}: 5283 @code{void _gfortran_caf_co_min (gfc_descriptor_t *a, int result_image, 5284 int *stat, char *errmsg, int a_len, size_t errmsg_len)} 5285 5286 @item @emph{Arguments}: 5287 @multitable @columnfractions .15 .70 5288 @item @var{a} @tab intent(inout) An array descriptor for the data to be 5289 processed. On the destination image(s) the result overwrites the old content. 5290 @item @var{result_image} @tab intent(in) The ID of the image to which the 5291 reduced value should be copied to; if zero, it has to be copied to all images. 5292 @item @var{stat} @tab intent(out) Stores the status STAT= and may be NULL. 5293 @item @var{errmsg} @tab intent(out) When an error occurs, this will be set to 5294 an error message; may be NULL. 5295 @item @var{a_len} @tab intent(in) the string length of argument @var{a} 5296 @item @var{errmsg_len} @tab intent(in) the buffer size of errmsg 5297 @end multitable 5298 5299 @item @emph{NOTES} 5300 If @var{result_image} is nonzero, the data in the array descriptor @var{a} on 5301 all images except of the specified one become undefined; hence, the library may 5302 make use of this. 5303 @end table 5304 5305 5306 5307 @node _gfortran_caf_co_sum 5308 @subsection @code{_gfortran_caf_co_sum} --- Collective summing reduction 5309 @cindex Coarray, _gfortran_caf_co_sum 5310 5311 @table @asis 5312 @item @emph{Description}: 5313 Calculates for each array element of the variable @var{a} the sum of all 5314 values for that element in the current team; if @var{result_image} has the 5315 value 0, the result shall be stored on all images, otherwise, only on the 5316 specified image. This function operates on numeric values only. 5317 5318 @item @emph{Syntax}: 5319 @code{void _gfortran_caf_co_sum (gfc_descriptor_t *a, int result_image, 5320 int *stat, char *errmsg, size_t errmsg_len)} 5321 5322 @item @emph{Arguments}: 5323 @multitable @columnfractions .15 .70 5324 @item @var{a} @tab intent(inout) An array descriptor with the data to be 5325 processed. On the destination image(s) the result overwrites the old content. 5326 @item @var{result_image} @tab intent(in) The ID of the image to which the 5327 reduced value should be copied to; if zero, it has to be copied to all images. 5328 @item @var{stat} @tab intent(out) Stores the status STAT= and may be NULL. 5329 @item @var{errmsg} @tab intent(out) When an error occurs, this will be set to 5330 an error message; may be NULL. 5331 @item @var{errmsg_len} @tab intent(in) the buffer size of errmsg 5332 @end multitable 5333 5334 @item @emph{NOTES} 5335 If @var{result_image} is nonzero, the data in the array descriptor @var{a} on 5336 all images except of the specified one become undefined; hence, the library may 5337 make use of this. 5338 @end table 5339 5340 5341 5342 @node _gfortran_caf_co_reduce 5343 @subsection @code{_gfortran_caf_co_reduce} --- Generic collective reduction 5344 @cindex Coarray, _gfortran_caf_co_reduce 5345 5346 @table @asis 5347 @item @emph{Description}: 5348 Calculates for each array element of the variable @var{a} the reduction 5349 value for that element in the current team; if @var{result_image} has the 5350 value 0, the result shall be stored on all images, otherwise, only on the 5351 specified image. The @var{opr} is a pure function doing a mathematically 5352 commutative and associative operation. 5353 5354 The @var{opr_flags} denote the following; the values are bitwise ored. 5355 @code{GFC_CAF_BYREF} (1) if the result should be returned 5356 by reference; @code{GFC_CAF_HIDDENLEN} (2) whether the result and argument 5357 string lengths shall be specified as hidden arguments; 5358 @code{GFC_CAF_ARG_VALUE} (4) whether the arguments shall be passed by value, 5359 @code{GFC_CAF_ARG_DESC} (8) whether the arguments shall be passed by descriptor. 5360 5361 5362 @item @emph{Syntax}: 5363 @code{void _gfortran_caf_co_reduce (gfc_descriptor_t *a, 5364 void * (*opr) (void *, void *), int opr_flags, int result_image, 5365 int *stat, char *errmsg, int a_len, size_t errmsg_len)} 5366 5367 @item @emph{Arguments}: 5368 @multitable @columnfractions .15 .70 5369 @item @var{a} @tab intent(inout) An array descriptor with the data to be 5370 processed. On the destination image(s) the result overwrites the old content. 5371 @item @var{opr} @tab intent(in) Function pointer to the reduction function 5372 @item @var{opr_flags} @tab intent(in) Flags regarding the reduction function 5373 @item @var{result_image} @tab intent(in) The ID of the image to which the 5374 reduced value should be copied to; if zero, it has to be copied to all images. 5375 @item @var{stat} @tab intent(out) Stores the status STAT= and may be NULL. 5376 @item @var{errmsg} @tab intent(out) When an error occurs, this will be set to 5377 an error message; may be NULL. 5378 @item @var{a_len} @tab intent(in) the string length of argument @var{a} 5379 @item @var{errmsg_len} @tab intent(in) the buffer size of errmsg 5380 @end multitable 5381 5382 @item @emph{NOTES} 5383 If @var{result_image} is nonzero, the data in the array descriptor @var{a} on 5384 all images except of the specified one become undefined; hence, the library may 5385 make use of this. 5386 5387 For character arguments, the result is passed as first argument, followed 5388 by the result string length, next come the two string arguments, followed 5389 by the two hidden string length arguments. With C binding, there are no hidden 5390 arguments and by-reference passing and either only a single character is passed 5391 or an array descriptor. 5392 @end table 5393 5394 5395 @c Intrinsic Procedures 5396 @c --------------------------------------------------------------------- 5397 5398 @include intrinsic.texi 5399 5400 5401 @tex 5402 \blankpart 5403 @end tex 5404 5405 @c --------------------------------------------------------------------- 5406 @c Contributing 5407 @c --------------------------------------------------------------------- 5408 5409 @node Contributing 5410 @unnumbered Contributing 5411 @cindex Contributing 5412 5413 Free software is only possible if people contribute to efforts 5414 to create it. 5415 We're always in need of more people helping out with ideas 5416 and comments, writing documentation and contributing code. 5417 5418 If you want to contribute to GNU Fortran, 5419 have a look at the long lists of projects you can take on. 5420 Some of these projects are small, 5421 some of them are large; 5422 some are completely orthogonal to the rest of what is 5423 happening on GNU Fortran, 5424 but others are ``mainstream'' projects in need of enthusiastic hackers. 5425 All of these projects are important! 5426 We will eventually get around to the things here, 5427 but they are also things doable by someone who is willing and able. 5428 5429 @menu 5430 * Contributors:: 5431 * Projects:: 5432 @end menu 5433 5434 5435 @node Contributors 5436 @section Contributors to GNU Fortran 5437 @cindex Contributors 5438 @cindex Credits 5439 @cindex Authors 5440 5441 Most of the parser was hand-crafted by @emph{Andy Vaught}, who is 5442 also the initiator of the whole project. Thanks Andy! 5443 Most of the interface with GCC was written by @emph{Paul Brook}. 5444 5445 The following individuals have contributed code and/or 5446 ideas and significant help to the GNU Fortran project 5447 (in alphabetical order): 5448 5449 @itemize @minus 5450 @item Janne Blomqvist 5451 @item Steven Bosscher 5452 @item Paul Brook 5453 @item Tobias Burnus 5454 @item Fran@,{c}ois-Xavier Coudert 5455 @item Bud Davis 5456 @item Jerry DeLisle 5457 @item Erik Edelmann 5458 @item Bernhard Fischer 5459 @item Daniel Franke 5460 @item Richard Guenther 5461 @item Richard Henderson 5462 @item Katherine Holcomb 5463 @item Jakub Jelinek 5464 @item Niels Kristian Bech Jensen 5465 @item Steven Johnson 5466 @item Steven G. Kargl 5467 @item Thomas Koenig 5468 @item Asher Langton 5469 @item H. J. Lu 5470 @item Toon Moene 5471 @item Brooks Moses 5472 @item Andrew Pinski 5473 @item Tim Prince 5474 @item Christopher D. Rickett 5475 @item Richard Sandiford 5476 @item Tobias Schl@"uter 5477 @item Roger Sayle 5478 @item Paul Thomas 5479 @item Andy Vaught 5480 @item Feng Wang 5481 @item Janus Weil 5482 @item Daniel Kraft 5483 @end itemize 5484 5485 The following people have contributed bug reports, 5486 smaller or larger patches, 5487 and much needed feedback and encouragement for the 5488 GNU Fortran project: 5489 5490 @itemize @minus 5491 @item Bill Clodius 5492 @item Dominique d'Humi@`eres 5493 @item Kate Hedstrom 5494 @item Erik Schnetter 5495 @item Gerhard Steinmetz 5496 @item Joost VandeVondele 5497 @end itemize 5498 5499 Many other individuals have helped debug, 5500 test and improve the GNU Fortran compiler over the past few years, 5501 and we welcome you to do the same! 5502 If you already have done so, 5503 and you would like to see your name listed in the 5504 list above, please contact us. 5505 5506 5507 @node Projects 5508 @section Projects 5509 5510 @table @emph 5511 5512 @item Help build the test suite 5513 Solicit more code for donation to the test suite: the more extensive the 5514 testsuite, the smaller the risk of breaking things in the future! We can 5515 keep code private on request. 5516 5517 @item Bug hunting/squishing 5518 Find bugs and write more test cases! Test cases are especially very 5519 welcome, because it allows us to concentrate on fixing bugs instead of 5520 isolating them. Going through the bugzilla database at 5521 @url{https://gcc.gnu.org/@/bugzilla/} to reduce testcases posted there and 5522 add more information (for example, for which version does the testcase 5523 work, for which versions does it fail?) is also very helpful. 5524 5525 @item Missing features 5526 For a larger project, consider working on the missing features required for 5527 Fortran language standards compliance (@pxref{Standards}), or contributing 5528 to the implementation of extensions such as OpenMP (@pxref{OpenMP}) or 5529 OpenACC (@pxref{OpenACC}) that are under active development. Again, 5530 contributing test cases for these features is useful too! 5531 5532 @end table 5533 5534 5535 @c --------------------------------------------------------------------- 5536 @c GNU General Public License 5537 @c --------------------------------------------------------------------- 5538 5539 @include gpl_v3.texi 5540 5541 5542 5543 @c --------------------------------------------------------------------- 5544 @c GNU Free Documentation License 5545 @c --------------------------------------------------------------------- 5546 5547 @include fdl.texi 5548 5549 5550 5551 @c --------------------------------------------------------------------- 5552 @c Funding Free Software 5553 @c --------------------------------------------------------------------- 5554 5555 @include funding.texi 5556 5557 @c --------------------------------------------------------------------- 5558 @c Indices 5559 @c --------------------------------------------------------------------- 5560 5561 @node Option Index 5562 @unnumbered Option Index 5563 @command{gfortran}'s command line options are indexed here without any 5564 initial @samp{-} or @samp{--}. Where an option has both positive and 5565 negative forms (such as -foption and -fno-option), relevant entries in 5566 the manual are indexed under the most appropriate form; it may sometimes 5567 be useful to look up both forms. 5568 @printindex op 5569 5570 @node Keyword Index 5571 @unnumbered Keyword Index 5572 @printindex cp 5573 5574 @bye 5575