#include <contrib/Eigen/src/GPU/DeviceBlasExpr.h>
C.rankUpdate(A), i.e. C += A * A^H, dispatched to cublasXsyrk (real) or cublasXherk (complex).