code: def cublasChpr2(handle, uplo, n, alpha, x, inx, y, incy, AP):
      
  
      status = _libcublas.cublasChpr2_v2(handle, 
                                         _CUBLAS_FILL_MODE[uplo], 
                                         n, ctypes.byref(cuda.cuFloatComplex(alpha.real,
                                                                             alpha.imag)),
                                         int(x), incx, int(y), incy, int(AP))
      cublasCheckStatus(status)
docs: Rank-2 operation on Hermitian-packed matrix.