Hello all,
In the following code two derived types, t1,t2 are defined and two functions whose arguments are of type t1, and output of type t2, an elemental and a "classic one". When the elemental is used wrong results are produced if t1 is defined as a class array. The code follows, along with correct and wrong results. Did something change about the class keyword from the previous version? Thank you!
module type_funs
type t1
integer :: i
end type t1
type t2
integer :: i
end type t2
type(t1) :: t1_start = t1(0)
contains
type(t2) elemental function bin_op(this,that) result(res)
type(t1), intent(in) :: this, that
res%i=this%i-that%i
end function bin_op
function bin_op_arr(this,that) result(res)
type(t1), dimension(:), intent(in) :: this
type(t1), intent(in) :: that
type(t2), dimension(size(this)) :: res
res%i=this%i-that%i
end function bin_op_arr
end module type_funs
program test_funs
use type_funs
implicit none
! this produces wrong results
class(t1), dimension(:), allocatable :: t1arr
! this produces correct results
!type(t1), dimension(:), allocatable :: t1arr
allocate(t1arr(10))
t1arr%i=[1:10]
print *, bin_op(t1arr,t1_start)
print *, bin_op_arr(t1arr,t1_start)
end program test_funs
correct result
1 2 3 4 5 6
7 8 9 10
1 2 3 4 5 6
7 8 9 10
wrong result
1 1 1 1 1 1
1 1 1 1
1 2 3 4 5 6
7 8 9 10
