The implicit statement (including implicit none) applies to a scoping unit. Such a thing is defined as BLOCK construct, derived-type definition, interface body, program unit, or subprogram, excluding all nested scoping units in it This “excluding all nested scoping units in it” suggests that it may be necessary/desirable to have implicit none in each function … Read more
ANY should actually be the right choice if ( ANY( lastNeighArray==”n” ) ) then there is also ALL if you wanted the whole array to contain that value.
As indicated in comments, real*8 isn’t standard Fortran. This FAQ entry has more details. That said, once we’re into the non-standard realm… The 8 refers to the number of bytes that the data type uses. So a 32-bit integer is integer*4 along the same lines. (But is also non-standard.) A quick search found this guide … Read more
Yeah, pointer arrays are funny in Fortran. The problem is that this: TYPE(domain),DIMENSION(:),POINTER :: dom does not define an array of pointers, as you might think, but a pointer to an array. There’s a number of cool things you can do with these things in Fortran – pointing to slices of large arrays, even with … Read more
According to The Fortran 2003 Handbook by Adams, et al., there is one difference between an intent(inout) argument and argument without specified intent. The actual argument (i.e., in the caller) in the intent(inout) case must always be definable. If the intent is not specified, the argument must be definable if execution of the subroutine attempts … Read more
From https://groups.google.com/forum/#!msg/android-ndk/QR1qiN0jIpE/g0MHkhTd4YMJ as selalerer suggested. I didn’t try this, so I’m posting as a community wiki for reference purposes. Fortran for x86 Android ================= The guide is based on this one, many thanks to Phil: Compiling Android NDK with Objective-C-enabled gcc errors 1) Download and unpack Android NDK ‘android-ndk-r8c’, (the older -r8b NDK won’t work … Read more
See also What are extern variables in C. This is mentioned in the C standard in informative Annex J as a common extension: J.5.11 Multiple external definitions There may be more than one external definition for the identifier of an object, with or without the explicit use of the keyword extern; if the definitions disagree, … Read more
I used to just do use modulename – then, as my application grew, I found it more and more difficult to find the source to functions (without turning to grep) – some of the other code floating around the office still uses a one-subroutine-per-file, which has its own set of problems, but it makes it … Read more
You can do that using reshape and shape intrinsics. Something like: INTEGER, DIMENSION(3, 3) :: array array = reshape((/ 1, 2, 3, 4, 5, 6, 7, 8, 9 /), shape(array)) But remember the column-major order. The array will be 1 4 7 2 5 8 3 6 9 after reshaping. So to get: 1 2 … Read more
I have followed the Julia project for a while now, and I have some comments to the code that might be relevant. It seems like you run a substantial amount of the code in global scope. The global environment is currently very slow in Julia, because the types all variables have to be checked on … Read more