Step 4: Poisson solve — divide by waves x and y are cosine transformed, so neither has a Nyquist mode and the only singular entry is the mean, which the eps_wave test handles. The x Nyquist zeroing of the 010 case must not be copied here: in this layout x index nx/2 + 1 is the regular cosine mode N/2 (#465).
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| complex(kind=dp), | intent(inout), | device, dimension(:, :, :) | :: | div_u | ||
| complex(kind=dp), | intent(in), | device, dimension(:, :, :) | :: | waves | ||
| integer, | intent(in), | value | :: | nz_h | ||
| integer, | intent(in), | value | :: | nx | ||
| integer, | intent(in), | value | :: | ny |