A novel and comprehensive method is proposed for calculating the dislocation Love numbers (DLNs), Green's functions (GFs), and the corresponding deformation in a transversely isotropic and layered elastic half-space. It is based on the newly introduced Fourier-Bessel series system of vector functions, along with the dual variable and position method. Two important features associated with this new system are: (1) it is much faster than the conventional cylindrical system of vector functions; (2) we can even pre-calculate the DLNs which are only possible in terms of this new system. This is due to the fact that the variables to be solved in the new system are functions of the simple discrete zero points of the Bessel functions, instead of the numerical integration of continuous Bessel functions between the neighboring zero points as in the conventional system. The introduced dual variable and position method is unconditionally stable as compared to the traditional propagator matrix method in dealing with layering. Exact asymptotic expressions of the DLNs for large wavenumber are further derived, which makes the Kummer's transformation applicable in accelerating the convergence of the corresponding GFs. For the reduced case of homogeneous and isotropic half-space, the present solutions amazingly reduce to the existing exact closed-form solutions. These new features are further seamlessly combined for calculating the deformation due to a finite dislocation (or a finite fault in geophysics) in the layered structure, which are demonstrated to be accurate and efficient.