Игорь Пашев - НОД на Фортране (2003)

Jun. 26th, 2011

02:23 pm - НОД на Фортране (2003)

Previous Entry Add to Memories Tell A Friend Next Entry

Язык в целом приятный, но доступ к командной строке,
стандартизированный только в Фортране 2003, убивает.

https://github.com/ip1981/GCD

! SYNOPSIS:
!
! # gfortran -o gcd-f gcd.f03
! # ./gcd-f 11 22 33 121
!

program GCD
implicit none

integer, allocatable :: ns(:)
integer :: i, n
character*20 :: tmpstr

n = command_argument_count()

allocate(ns(n)) ! allocate memory for numbers given in command line

do i = 1, n
    call get_command_argument(i, tmpstr)
    ns(i) = str2int(tmpstr)
end do

print *,  gcdn(ns)

deallocate(ns)


! If we declare functions first,
! we have to specify its types within
! the `program' section.
! See http://en.wikibooks.org/wiki/Fortran/Fortran_procedures_and_functions
contains

pure integer function str2int (s)
    character*(*), intent(in) :: s
    read (s, *) str2int
end function

pure recursive integer function gcd2 (a, b) result(GCD)
    integer, intent(in) :: a, b
    if (b == 0) then
        GCD = a
    else
        GCD = gcd2(b, mod(a, b))
    end if
end function gcd2

pure integer function gcdn(n)
    integer, intent(in) :: n(:) ! n is an array
    integer :: i
    gcdn = n(1)
    do i = 2, size(n)
        gcdn = gcd2(gcdn, n(i))
    end do
end function

end program

Tags: , ,
(4 комментария | Оставить комментарий)

Comments:

From:[info]phantom
Date:June 26th, 2011 - 02:21 pm
(Link)
Fortran многословен, однако.
(Reply to this) (Thread)
[User Picture]
From:[info]ketmar
Date:June 27th, 2011 - 04:42 pm
(Link)
это, в принципе, хорошо. си-подобный синтаксис довольно ублюдочен и нечитаем.
(Reply to this) (Parent) (Thread)
From:[info]phantom
Date:June 29th, 2011 - 10:51 pm
(Link)
Надо учиться видеть код сразу в AST-виде.
(Reply to this) (Parent) (Thread)
[User Picture]
From:[info]ketmar
Date:June 30th, 2011 - 12:31 am
(Link)
лисп уже придумали…
(Reply to this) (Parent)