         apdq[mq] += smax[1]*wave[meqn + mq];
         apdq[mq] += smax[2]*wave[2*meqn + mq];
    }    
-#endif
+#else
     for(int mq = 0; mq < meqn; mq++)
     {
         amdq[mq] = 0;

        bpasdq[mq] += smax[1]*waveb[meqn + mq];
        bpasdq[mq] += smax[2]*waveb[2*meqn + mq];
    }
#endif

    for(int mq = 0; mq < meqn; mq++)
    {

    {        
        double smin = (s[mw] <  0) ? s[mw] : 0.0;
        double smax = (s[mw] >= 0) ? s[mw] : 0.0; 
        int idx = mw*meqn;
        
        amdq[0] += smin*wave[0 + idx];
        amdq[1] += smin*wave[1 + idx];
        amdq[2] += smin*wave[2 + idx];

        apdq[0] += smax*wave[0 + idx];
        apdq[1] += smax*wave[1 + idx];
        apdq[2] += smax*wave[2 + idx];
    }
#endif
}

__device__ cudaclaw_cuda_rpn2_t bump_rpn2 = bump_rpn2shallow;
