fortran90

问题 I'm enough of a novice to not know the terminology, so I can't search the Web for the answer to this. More than once, in programming, I've wanted to do something like this. A and B are subroutines, c and d are functions. A and B each call a function multiple times inside them. call A(c(x)) call A(d(x)) call B(c(x)) call B(d(x)) This structure doesn't work. I'm told that Fortran doesn't support aliasing of functions, at least in this context. (Most search results involving "aliasing" refer to

问题 I am trying to collect and process all source files in the top directory, to deal with the fact that variables are not passed to sub directories. I have cobbled together the following CMakeLists.txt files: cmake_minimum_required(VERSION 2.8) project(main) enable_language(Fortran) enable_testing() set (CMAKE_Fortran_COMPILER ifort) set (CMAKE_Fortran_FLAGS " -g -C -fixed") set (CMAKE_Fortran95_FLAGS " -openmp ") # function to collect all the sources from sub-directories # into a single list

问题 I am trying to collect and process all source files in the top directory, to deal with the fact that variables are not passed to sub directories. I have cobbled together the following CMakeLists.txt files: cmake_minimum_required(VERSION 2.8) project(main) enable_language(Fortran) enable_testing() set (CMAKE_Fortran_COMPILER ifort) set (CMAKE_Fortran_FLAGS " -g -C -fixed") set (CMAKE_Fortran95_FLAGS " -openmp ") # function to collect all the sources from sub-directories # into a single list

问题 I am using MPI for solving PDE. For this, I breakdown the 2D domain into different cells (size of each of these cells is " xcell,ycell " with xcell = size_x_domain/(number of X subdomains) and ycell = size_y_domain/(number of Y subdomains) . So, I am running the code with number of processes = (number of X subdomains)*(number of Y subdomains) The gain relatively to sequential version is that I communicate between each process representing the sub-domains. Here a figure illustrating my

问题 I have an array with multiple dimensions (the goal is to allow for about 100) and each dimension has a size of about 2^10 and I only need to store in it about 1000 double precision coefficients. I don't need to do any operation with this array aside from reading and writing into it. The code is written in Fortran 90. I assume that if I a library like one of the ones mentioned in this answer I would be able to store the do this, but would this be optimized for the simple reading and writing

问题 Edit: Gfortran 6 now supports these extensions :) I have some old f77 code that extensively uses UNIONs and MAPs. I need to compile this using gfortran, which does not support these extensions. I have figured out how to convert all non-supported extensions except for these and I am at a loss. I have had several thoughts on possible approaches, but haven't been able to successfully implement anything. I need for the existing UDTs to be accessed in the same way that they currently are; I can

问题 This question already has answers here : Array declaration in Fortran (2 answers) Closed last year . I'm using gfortran version 7.2.0. I'm quite new to Fortran. I know there are different versions of Fortran. In the code below, I'm declaring arrays (or actually tensors) using different syntaxes program arrays implicit none integer :: m(3, 4) integer, dimension(3, 4) :: n print *, "m = ", m print *, "n = ", n end program arrays In one case, I'm using the dimension statement, in the other I am

问题 I'm learning Fortran90 after a brief introduction to Fortran77 a few years ago. When printing integers in Fortran, you must specify how many spaces you want to reserve for printing the integer. Consider this program... implicit none integer :: i i = 123 write(*, '(A, I3, A)') "'", i, "'" !3 spaces for output = no padding write(*, '(A, I5, A)') "'", i, "'" !5 is too many, so output is padded write(*, '(A, I2, A)') "'", i, "'" !2 is too few, so output is jibberish write(*, '(A, I:, A)') "'", i,

问题 I'm learning Fortran90 after a brief introduction to Fortran77 a few years ago. When printing integers in Fortran, you must specify how many spaces you want to reserve for printing the integer. Consider this program... implicit none integer :: i i = 123 write(*, '(A, I3, A)') "'", i, "'" !3 spaces for output = no padding write(*, '(A, I5, A)') "'", i, "'" !5 is too many, so output is padded write(*, '(A, I2, A)') "'", i, "'" !2 is too few, so output is jibberish write(*, '(A, I:, A)') "'", i,

问题 I have had problems calling successive subroutines with assumed-shape arrays in Fortran 90. More specifically, I call two levels of subroutines, passing an assumed-shape array as a parameter, but in the end the array is lost. To demonstrate it, one can follow the code below. program main INTERFACE subroutine sub1(x) real, dimension(:):: x real C end subroutine sub1 subroutine sub2(x) real, dimension(:):: x real C end subroutine sub2 END INTERFACE real, dimension(:), allocatable:: x allocate(x