  The fractional Laplacian $\Delta^{\beta/2}$ is the generator of
$\beta$-stable L\'evy process, which is the scaling limit of the L\'evy fight.
Due to the divergence of the second moment of the jump length of the L\'evy
fight it is not appropriate as a physical model in many practical applications.
However, using a parameter $\lambda$ to exponentially temper the isotropic
power law measure of the jump length leads to the tempered L\'evy fight, which
has finite second moment. For short time the tempered L\'evy fight exhibits the
dynamics of L\'evy fight while after sufficiently long time it turns to normal
diffusion. The generator of tempered $\beta$-stable L\'evy process is the
tempered fractional Laplacian $(\Delta+\lambda)^{\beta/2}$ [W.H. Deng, B.Y. Li,
W.Y. Tian, and P.W. Zhang, Multiscale Model. Simul., in press, 2017]. In the
current work, we present new computational methods for the tempered fractional
Laplacian equation, including the cases with the homogeneous and nonhomogeneous
generalized Dirichlet type boundary conditions. We prove the well-posedness of
the Galerkin weak formulation and provide convergence analysis of the single
scaling B-spline and multiscale Riesz bases finite element methods. We propose
a technique for efficiently generating the entries of the dense stiffness
matrix and for solving the resulting algebraic equation by preconditioning. We
also present several numerical experiments to verify the theoretical results.
