def drot(x, y, c, s): """ This function applies a Givens rotation \M{(x', y') = (c x + s y, -s x + c y)} to the vectors x, y. """ xn = array_typed_copy(x, Float) yn = array_typed_copy(y, Float) _gslwrap.gsl_blas_drot(xn, yn, c, s) return (xn, yn)
def drot_cr(x_CR, y_CR, c, s): """ This function applies a Givens rotation \M{(x', y') = (c x + s y, -s x + c y)} to the vectors x, y. """ _gslwrap.gsl_blas_drot(x_CR, y_CR, c, s)