I currently searching for a subroutine from BLAS or LAPACK which realizes the following operation
A = alpha*A + beta * B
where A and B have different leading dimensions, i.e. A is stored as
and B is stored as
with LDB > N.
Is there a way to use existing BLAS or LAPACK operation to compute this or must I create an own subroutine for this kind of problems?