#include <contrib/Eigen/src/GPU/DeviceSolverExpr.h>
d_A.llt().solve(d_B), dispatched to cusolverDnXpotrf + cusolverDnXpotrs.