One question: the cuda docs say that ``nrm2``, ``amax`` and related functions which return scalars need to have a synchronization call before the host accesses the data.  Is this done somewhere under the hood, or should a call to ``pycuda.autoinit.context.synchronize()`` be included here?