#include <contrib/Eigen/src/GPU/DeviceSolverExpr.h>
d_A.lu().solve(d_B), dispatched to cusolverDnXgetrf + cusolverDnXgetrs.