Description
When running a simulation using an initial plane wave, the flags to add Bielak absorbing conditions are ignored by the GPU solver. The simulation results using CPU and GPU are different.
Affected SPECFEM2D version
8.1.0
Your software and hardware environment
ubuntu22.04, RTX 5090, cuda 12.8
Reproduction steps
# use an existing initial wave field as source or start from zero
initialfield = .true.
add_Bielak_conditions_bottom = .true. # add Bielak conditions or not if initial plane wave
add_Bielak_conditions_right = .true.
add_Bielak_conditions_top = .true.
add_Bielak_conditions_left = .true.
# set to true to use GPUs
GPU_MODE = .true. / .false.