Могу ли я безопасно использовать `spread` для копирования из среза массива в тот же массив?

Question

2D-данные хранятся в файле, который я хочу загрузить и распределить в 3D-массив по последнему индексу. Могу ли я безопасно загрузить данные в срез, а затем использовать спред из этого, например, так:

real, dimension(nx,ny,nz) :: data

call load_some_data(data(:,:,1))

data = spread(data(:,:,1), 3, nz)

Или мне нужен временный 2D-массив, чтобы спред соответствовал стандарту? Будет ли spread вызов работать также в предложении OpenMP workshare?

Answer 1

Да, вы можете сделать это. В зависимости от уровня оптимизации компилятора он будет использовать или не будет использовать временный массив, но это допустимо.

Нет, он не будет распараллелен в конструкции workshare. Простое присвоение массива будет, но я боюсь, что результат функции будет проблематичным. И звонок на spread, безусловно, будет последовательным.