Hi,
When running kspaceFirstOrder3DC() (or 3DG) with a sensor.record that includes 'I' or 'I_avg', the computation of acoustic intensities seems to be performed in private/kspaceFirstOrder_saveIntensity.m by commands of the form:
ux_sgt = fourierShift(sensor_data(sensor_mask_index).ux_non_staggered, 1/2);
sensor_data(sensor_mask_index).Ix = sensor_data(sensor_mask_index).p .* ux_sgt;
Does this mean that the --u_non_staggered_raw option passed to kspaceFirstOrder3D-OMP (or -CUDA) only affects spatial but not temporal staggering (which then needs to be fixed with fourierShift())?
Thank you very much!