[Patch, fortran] PR29422 and PR29428 - allocatable component wrinkles
François-Xavier Coudert
fxcoudert@gmail.com
Thu Oct 12 07:43:00 GMT 2006
Hi Paul,
> In discussing PR29391 with FX, he quoted:
>
> "That's because of F95 13.14.53:
>
> Case (i): For an array section or for an array expression other than a whole
> array or array structure component, LBOUND(ARRAY, DIM) has the value 1.
> For a
> whole array or array structure component, LBOUND(ARRAY, DIM) has the value:
> (a) equal to the lower bound for subscript DIM of ARRAY if dimension DIM of
> ARRAY does not have extent zero or if ARRAY is an assumed-size array of rank
> DIM, or
> (b) 1 otherwise."
About that: I've thought again, and I think the best we can do is to
actually add the extra code in the {L,U}BOUND intrinsic. gfortran uses
0-based temporaries and incorrect bounds for many array functions
return values, but it doesn't matter at all. It's all internal code,
the only place where it shows up is when you call {L,U}BOUND, because
it gets the value directly from the descriptor.
I'm now working on a very different solution to PR29391, where LBOUND
actually checks whether its argument is "a whole array or array
structure component", whether it is "an assumed-size array of rank
DIM, or has extent zero in dimension DIM" and fetch the descriptor
bounds only in this case (and in this case, I believe that they
descriptor bound values are always "right"); in other cases, it
returns 1 (or, for UBOUND, returns the number of elements in this
dimension).
What do you think about that?
FX
More information about the Fortran
mailing list