Fix CUDA fused transeq implementation for Neumann BCs

[semi-h]

Current fused transeq function doesn't input the _sym (npaire) variations of the derivative operations. This cannot cause an issue for periodic or Dirichlet BCs but we can't get the correct results with Neumann BCs. A simple fix but important before running simulations with Neumann BCs.