Accelerating Band Linear Algebra Operations on GPUs with Application in Model Reduction