                         #endif
 
                         #if use_shared == 1
-                        //if (offset >= 0 && offset < shared_memory_size-1)
                             sh_col_idx[offset++] = i+j;
                         #else
                         col_idx[offset++] = i+j;

        if (k-block_start >= 0 && k-block_start < shared_memory_size-1)
            col_idx[k] = sh_col_idx[k-block_start];
    }

    #endif




}


        #endif


            int laneid = tx & (32-1);
            if (output) {
            for (int j=0; j < 32-laneid && output; j++) {
                if (i+j >= 0 && i+j<N) {

                float diffct = ct_i - LDG(ct,i+j);
