comp.lang.fortran
  Home FAQ Contact Sign in
Your Ad Here
comp.lang.fortran only
 
Advanced search
May 2008
motuwethfrsasuw
   1234 18
567891011 19
12131415161718 20
19202122232425 21
262728293031  22
2008
 Jan   Feb   Mar   Apr 
 May   Jun   Jul   Aug 
 Sep   Oct   Nov   Dec 
2008 2007 2006  
total
comp.lang.fortran Profile…
RELATED GROUPS

POPULAR GROUPS

more...


 Up
  Bad start on complex tests         


Author: James Van Buskirk
Date: May 16, 2008 02:24

C:\gfortran\james\intrinsics\func1>type abs.f90
program abstest
complex(10) c10

c10 = ( 1.7951958020513104220_10, 2.6457513110645905904_10 )
write(*,*) 'Alive to here'
write(*,*) abs(c10)
write(*,*) 'Survived the ordeal'
end program abstest

C:\gfortran\james\intrinsics\func1>gfortran abs.f90 -oabs

C:\gfortran\james\intrinsics\func1>abs
Alive to here

C:\gfortran\james\intrinsics\func1>

--
write(*,*) transfer((/17.392111325966148d0,6.5794487871554595D-85, &
6.0134700243160014d-154/),(/'x'/)); end
no comments
  KIND optional argument in transformational intrinsics         


Author: James Van Buskirk
Date: May 16, 2008 01:19

I found that gfortran is missing the KIND optional argument in some
of the transformational intrinsics. Also a couple of usages of KIND
were not documented, and gfortran does not yet accept f2008 usage for
the SIZE intrinsic.

C:\gfortran\james\intrinsics\func1>type kind.f90
program test
implicit none
integer, parameter :: ck = kind('A')
integer, parameter :: ik = selected_int_kind(18)
integer i
real x
integer array(1)

i = 42
x = 3.14
array = 13
Show full article (1.76Kb)
1 Comment
  Keyword arguments in intrinsics         


Author: James Van Buskirk
Date: May 16, 2008 01:15

I found some places where the gfortran manual had keyword arguments
inconsistent with standard. In all cases the actual compiler agreed
with the standard rather than gfortran documentation:

C:\gfortran\james\intrinsics\func1>type names.f90
program names
implicit none
integer i
real x
complex z
logical L
Show full article (1.13Kb)
no comments
  gfortran ---problem with random umber         


Author: rudra
Date: May 15, 2008 21:43

dear friends,
here is the full code that is creating the problem!!! while running
with ifort, the configuration is created as wished, but not in
gfortran....

!===========================================
! This is the main driver routine
! of the reverse monte carlo program
program rmonte
!===========================================
Show full article (4.37Kb)
7 Comments
  using gdb on windows         


Author: Ron Ford
Date: May 15, 2008 21:08

One reason that I keep silverfrost's f95 suite around is that if I can get
a program to compile, I can step through it with a debugger that I invoke
on one of the pull-downs. I'd like to emulate this capability with gcc's
gdb.

I googled for "gfortran debugger" and found a lot of information that was
often more relevant to the linux crowd. I *did* find this from a previous
exchange with FX:
> I read all the time about gfortran and gdb questions that seem to get
> solid answers. I believe they are available for linux.

OK then, a brief summary: gfortran is the GNU Fortran compiler, part of
GCC (the GNU Compiler Collection). GCC happens to be the system compiler
on Linux and most other free software operating systems, which means it's
widely used on these OS. gdb is the GNU debugger, also the default
debugger on Linux and other free software OS.

But that doesn't mean it's restricted to Linux. Both GCC (and gfortran)
and gdb are known to work on (free|net|open)bsd, Solaris, MacOS X,
Windows, AIX, Irix, HP-UX, Tru64 and probably others.

!end excerpt
Show full article (1.39Kb)
1 Comment
  IBM XLFortran 8.1 on MacOSX 10.5 (Leopard)         


Author: Al Greynolds
Date: May 15, 2008 07:52

I am currently running IBM's XLFortran 8.1 from a terminal shell on
MacOSX 10.4 (Tiger) without any issues (even though IBM has
discontinued it and never supported it on Tiger). I would like to
upgrade my system to MacOSX 10.5 (Leopard) and would like to know
ahead of time that XLF can be made to work from a terminal shell
(don't care about Xcode or even debugger support). Has anybody been
down this path already?

Al Greynolds
www.ruda.com
no comments
  fl32 vs df in cvf6.6         


Author: formulae translator
Date: May 15, 2008 07:47

Today I tried a very simple speed test using fl32 and df command
(compaq visual fortran 6.6)

program speed
integer :: i, j, k, n
real :: a, b, c
call cpu_time(a)
do i=1,1000
do j=1,100
do k=1,100
c=(sin(real(i))**2.0+cos(real(j))**2.0)+log(real(i))*log(real(j))
+real(k)
! write(*,*) c
end do
end do
end do
call cpu_time(b)
b=b-a
write(*,*) b,c
end program
Show full article (0.57Kb)
18 Comments
  fl32 vs df         


Author: formulae translator
Date: May 15, 2008 07:42

Today I tried a very simple speed test using fl32 and df command
(compaq visual fortran 6.6)

program speed
integer :: i, j, k, n
real :: a, b, c
call cpu_time(a)
do i=1,1000
do j=1,100
do k=1,100
c=(sin(real(i))**2.0+cos(real(j))**2.0)+log(real(i))*log(real(j))
+real(k)
!write(*,*) c
end do
end do
end do
call cpu_time(b)
b=b-a
write(*,*) b,c
end program
Show full article (0.50Kb)
no comments
  fl32 vs df         


Author: formulae translator
Date: May 15, 2008 07:42

Today I tried a very simple speed test using fl32 and df command
(compaq visual fortran 6.6)

program speed
integer :: i, j, k, n
real :: a, b, c
call cpu_time(a)
do i=1,1000
do j=1,100
do k=1,100
c=(sin(real(i))**2.0+cos(real(j))**2.0)+log(real(i))*log(real(j))
+real(k)
!write(*,*) c
end do
end do
end do
call cpu_time(b)
b=b-a
write(*,*) b,c
end program
Show full article (0.50Kb)
no comments
Your Ad Here
  Pointers question         


Author: tomguest
Date: May 15, 2008 07:32

This is a follow on from an earlier post of mine regarding sorting.

I've never used pointers before and I'm finding the lingo a bit heavy.

Here's my problem.

I have a potentially very large matrix that looks like

1 2 3 4 5 6
2 3 4 5 6 7
3 4 5 6 7 8
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0

For my largest problem the matrix will have 11 columns and about
10,000,000 rows.

I'm only interested in the section that is nonzero but am not aware of
how to change the matrix size quickly and efficiently as more rows are
needed. I was told that pointers maybe the way to go. The example I
was shown was to declare a single column vector which represents my
first column i.e.
Show full article (0.95Kb)
5 Comments
1 2 3 4 5 6 7 8 9