rotm.hpp

Go to the documentation of this file.

 
//
// Copyright (c) 2017-2021, University of Tennessee. All rights reserved.
// Copyright (c) 2021-2023, University of Colorado Denver. All rights reserved.
//
// This file is part of <T>LAPACK.
// <T>LAPACK is free software: you can redistribute it and/or modify it under
// the terms of the BSD 3-Clause license. See the accompanying LICENSE file.
 
#ifndef TLAPACK_LEGACY_ROTM_HH
#define TLAPACK_LEGACY_ROTM_HH
 
#include "tlapack/blas/rotm.hpp"
#include "tlapack/legacy_api/base/types.hpp"
#include "tlapack/legacy_api/base/utils.hpp"
 
namespace tlapack {
namespace legacy {
 
    template <typename TX, typename TY>
    void rotm(idx_t n,
              TX* x,
              int_t incx,
              TY* y,
              int_t incy,
              scalar_type<TX, TY> const param[5])
    {
        // constants
        const int flag = (int)param[0];
        const scalar_type<TX, TY>* h = &param[1];
 
        // check arguments
        tlapack_check_false(incx == 0);
        tlapack_check_false(incy == 0);
        tlapack_check_false(flag < -2 && flag > 1);
 
        // quick return
        if (n <= 0) return;
 
        tlapack_expr_with_2vectors(
            x_, TX, n, x, incx, y_, TY, n, y, incy,
            if (flag == -2) return rotm<-2>(x_, y_, h);
            else if (flag == -1) return rotm<-1>(x_, y_, h);
            else if (flag == 0) return rotm<0>(x_, y_, h);
            else return rotm<1>(x_, y_, h););
    }
 
}  // namespace legacy
}  // namespace tlapack
 
#endif  //  #ifndef TLAPACK_LEGACY_ROTM_HH