SUMMARYIn this paper, we investigate the normwise, mixed, and componentwise condition numbers and their upper bounds for the Moore-Penrose inverse of the Kronecker product and more general matrix function compositions involving Kronecker products. We also present the condition numbers and their upper bounds for the associated Kronecker product linear least squares solution with full column rank. In practice, the derived upper bounds for the mixed and componentwise condition numbers for Kronecker product linear least squares solution can be efficiently estimated using the Hager-Higham Algorithm.