Efficient matrix computation for tensor-product isogeometric analysis: The use of sum factorization