CODEWITHSUNDEEP
arrayspointersfortran
ClimateUnboxed
ClimateUnboxed

Reputation: 8067

How to use of an array of pointers to assign non-adjacent array entries in FORTRAN

I want to be able to reset elements in an array which are not in a consecutive block of memory. I thought to use an array of pointers for this rather than a pointer array, as my understanding of pointer arrays is that they must point to a coherent block of memory (e.g pointer(1:10,1:10) => target(1:100) )

My simple test program is as follows:

program test
  implicit none

  integer, target :: arr(4,4)

  type ptr
    integer, pointer :: p
  end type ptr

  type(ptr), dimension(2) :: idx 

  arr=0

  idx(1)%p=>arr(2,2) 
  idx(2)%p=>arr(4,2)

  idx(1)%p=5 ! this is okay 
  idx(2)%p=5 ! this is okay  
  idx(1:2)%p=5 ! this gives an error

  print *,arr

end program test

The first two statements idx(n)%p=5 are okay, but I want to be able to set a chunk of the array in one statement using the spanning method idx(1:n)%p=5 but when I do this I get the following compile error:

Error: Component to the right of a part reference with nonzero rank must not have the POINTER attribute at (1)

Can I set the value of a chunk of array entries using pointer somehow? Perhaps it is in fact possible with a pointer array rather than an array of pointers...

I think this is related to Fortran: using a vector to index a multidimensional array

but I can't see how to use that answer here.

Upvotes: 0

Views: 287

Answers (1)

High Performance Mark
High Performance Mark

Reputation: 78316

Probably an extended comment rather than an answer, but I need a little formatting ...

It's best not to think of Fortran supporting arrays of pointers. You've already grasped that (I think) and built the usual workaround of an array of derived type each instance of the type having an element which is a pointer.

What's not entirely clear is why you don't use multiple vector subscripts such as 

arr([2,4],[2]) = 5

(OK, the second subscript is something of a degenerate vector but the statement is intended to have the same effect as your attempt to use pointers.) If you would rather, you could use array subscript triplets instead

arr(2:4:2,2) = 5

Perhaps you've simplified your needs too far to express the necessity of using pointers in which case my feeble suggestions will not meet your undeclared needs.

Upvotes: 1

Related Questions