Sum over an array component of a multi-dimensional derived type (without do loop)? [duplicate]

Question

I have a derived type array, where the derived types can also be arrays. Is there any way to sum over elements without resorting to a do loop?

program main

   type t_couple
      real :: income(2) = 5
   end type t_couple

   type(t_couple) :: couple(2)

   print *, " sum 1 ", sum(couple(1)%income)  ! works fine
   print *, " sum 2 ", sum(couple%income)     ! doesn't work

end program main

%F90-E-ERROR, A component cannot be an array if the encompassing structure is an array.

Of course we have all seen this error message many times and can write a do loop to solve this, but I'm wondering if there is a simpler and more concise way?

Alexander Vogt · Accepted Answer · 2016-04-25 19:37:51Z

1

You could use an implicit loop:

print *, " sum 2 ", sum( [( couple(i)%income, i=1,size(couple))] )

answered Apr 25, 2016 at 19:37

Alexander Vogt

18.2k13 gold badges56 silver badges72 bronze badges

Sign up to request clarification or add additional context in comments.

2 Comments

JohnE Over a year ago

Well, +1 for doing it on one line but it's still a loop... But probably there is no way to do it without a loop.

Alexander Vogt Over a year ago

@JohnE At least I don't know of any other way...

Collectives™ on Stack Overflow

Sum over an array component of a multi-dimensional derived type (without do loop)? [duplicate]

1 Answer 1

2 Comments

Linked

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

2 Comments

Linked

Related