     CUDABlocks = ceil((REAL)(qureg.numAmpsPerChunk>>1)/threadsPerCUDABlock);
     statevec_sigmaXKernel<<<CUDABlocks, threadsPerCUDABlock>>>(qureg, targetQubit);
 }
-
 
 __global__ void statevec_sigmaYKernel(QubitRegister qureg, const int targetQubit, const int conjFac){
 

    vecImag[stateInd] *= fac;
}

void densmatr_twoQubitDephase(QubitRegister qureg, int qubit1, int qubit2, REAL dephase) {
    if (dephase == 0)
        return;

        part1, part2, part3, part4, part5, colBit1, rowBit1, colBit2, rowBit2);
}

/** Works like oneQubitDephase but modifies every other element, and elements are averaged in pairs */
__global__ void densmatr_oneQubitDepolariseKernel(
    REAL depolLevel, REAL* vecReal, REAL *vecImag, long long int numAmpsToVisit,
    long long int part1, long long int part2, long long int part3, 
