pr32512 - efficiency of RESHAPE and SPREAD

Paul Thomas paulthomas2@wanadoo.fr
Mon Aug 13 20:34:00 GMT 2007


Dear All,

In contemplating this excellent suggestion from Jaroslav Hajek, I have 
run into an interesting puzzle.

trans-array.c (gfc_conv_array_transpose) implements a transpose in which 
the operation is entirely performed inline and only acts on the 
descriptor.  Jaroslav suggests that we do the same for reshape and spread.

In investigating this, I found the rather wierd and wonderful property 
of gfortran; gfc_conv_array_transpose only works when it is the actual 
argument of something.  If I lift this constraint, it still works fine.

Does anybody know why this is so?

Cheers

Paul



More information about the Fortran mailing list