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