Recursive procedures in subroutines

Help
1

Dear friends, I have to perform a recursive procedure within a subroutine, something, for example, like this:

module xxxx
…
subroutine rrrr ! (for simplicity I don’t state the intent of the variables)

A = A + b*50

end subroutine rrrr
end module xxxx

program yyyy
use xxxx
…
do i = 1,100
call rrrr()
end do
…
end program yyyy

My doubt is that I do not know where to put the initial value of A = 0. For example, if I place it within the module before the subroutine, every time the subroutine is called, will A be set = 0? In this case, it won’t work. Therefore, could you please help me understand where A should be placed in such a code?

Many thanks for your kind help.

2

I don’t see recursion in your example.
A rather useless recursive program is given below

program xxxx

  implicit none

  print*, rrrr(100.)

contains

recursive function rrrr(u) result(A)
  real, intent(in) :: u
  real :: A

  print*, 'u', u
  if (u < .5) then
    A = u
  else
    A = rrrr(u/2.)
  endif

end function rrrr

end program xxxx
1 Like
3

Placing it within the module and initialising it to 0 would work.
Placing it inside the recursive procedure and initialing it to 0 would work.

Initialization can only be done

a) as part of a declaration statement, e.g.
REAL:: A=0.

b)separately, in a DATA statement.e.g.
REAL A
DATA A/0./

c) for a derived type, as default initialization, e.g.
TYPE widget
REAL ::A=0.
END TYPE

TYPE(widget) ::x

an assignment statement (A=0.) does NOT count as initialization.

1 Like
4

unfortunately not Fortran, but the concepts are similar:

Currently, the course is free of charge.

5

@themos, @MarDie, many thanks for your help.