Fortran

Floating Point Multiplication - 32 and 64 bit differences
Hello, I need some help understanding floating point operations. I am seeing a small difference in a program when it is compiled for 32 and 64 bit machines. Here is the program: -------------------------------------------------------- PROGRAM DPStudy IMPLICIT NONE REAL(8) :: FVF REAL(8) :: THR REAL(8) ... 27 May 2010 10:20

negative zeros and dot_product matmul
m_b_metcalf <michaelmetcalf(a)compuserve.com> wrote: On May 16, 9:07�pm, Thomas Koenig <tkoe...(a)netcologne.de> wrote: Consider the following program: program main � implicit none � real, dimension(1) :: a, b � real, dimension(1,1) :: aa, bb � a(1) = -1.0 � b(1) = 0.0 � print *,a(1)... 17 May 2010 19:06

meaning of #if system/port/comp
Hi all, I would like to ask if somebody could explain to me the meaning of such a code: #if system #if comp #if port As long as I understand, with such a ifs I can differentiate between windows and linux/unix, but that's all. Thanks & regards peter ... 17 May 2010 09:05

How to define a user data type containing an array with dimensions decided in running time
My particular situation is that n, the dimension of vector, remains the same throughout the program and there are many variables (scalars, vectors, matrices) of type User_Data. And I do not perform "allocate" statement for each such variable. On May 15, 7:46 pm, nos...(a)see.signature (Richard Maine) wrote: ... 16 May 2010 22:20

Pointer for... pointing and array reallocation issue
Richard Maine wrote: Hifi-Comp <wenbinyu.heaven(a)gmail.com> wrote: I am thinking to create a user data type containing an array with dimensions decided in running time. Currently I have the following: INTEGER,PARAMETER:: n=3 TYPE,PUBLIC:: User_Data REAL(DBL_AD)::scal... 18 May 2010 13:39

negative zeros and dot_product / matmul
Consider the following program: program main implicit none real, dimension(1) :: a, b real, dimension(1,1) :: aa, bb a(1) = -1.0 b(1) = 0.0 print *,a(1)*b(1) print *,dot_product(a,b) aa(1,1) = -1.0 bb(1,1) = 0.0 print *,matmul(aa,bb) end program main Assuming that the processor supp... 17 May 2010 13:30

How to define a user data type containing an array with dimensions decided in running time
Hifi-Comp <wenbinyu.heaven(a)gmail.com> wrote: I am thinking to create a user data type containing an array with dimensions decided in running time. Currently I have the following: INTEGER,PARAMETER:: n=3 TYPE,PUBLIC:: User_Data REAL(DBL_AD)::scale REAL(DBL_AD)::vector(n) ... 15 May 2010 20:15

How to define a user data type containing an array with dimensions decided in running time
Dear All, I am thinking to create a user data type containing an array with dimensions decided in running time. Currently I have the following: INTEGER,PARAMETER:: n=3 TYPE,PUBLIC:: User_Data REAL(DBL_AD)::scale REAL(DBL_AD)::vector(n) END TYPE User_Data In ideal situation, I want the user to i... 15 May 2010 18:04

Need efficient Fortran equivalent to Matlab accumarray
Hello, I'm looking for a efficient, vectorized Fortran90/95 code that works like Matlab's accumarray function, which accumulates elements of an input (data) array, based on subscript arrays (of the same shape as the data array) which point to the location in a new target array that the input array are to acc... 15 May 2010 06:08

deferred length character as subroutine argument
I am struggling to get something like the following to work - is it possible at all? cheers -ralf program main character( : ), allocatable :: theStr call GetStr( theStr ) contains subroutine GetStr( theStr ) character( : ), allocatable, intent( out ) :: theStr ! open a file etc. ... 15 May 2010 16:58