utils.hpp

Go to the documentation of this file.

 
//
// 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_UTILS_HH
#define TLAPACK_LEGACY_UTILS_HH
 
#include "tlapack/legacy_api/base/legacyArray.hpp"
 
#define tlapack_expr_with_2vectors(x, TX, n, X, incx, y, TY, m, Y, incy, expr) \
    do {                                                                       \
        using tlapack::legacy::internal::create_vector;                        \
        using tlapack::legacy::internal::create_backward_vector;               \
        if (incx == 1) {                                                       \
            auto x = create_vector((TX*)X, n);                                 \
            tlapack_expr_with_vector(y, TY, m, Y, incy, expr);                 \
        }                                                                      \
        else if (incx == -1) {                                                 \
            auto x = create_backward_vector((TX*)X, n);                        \
            tlapack_expr_with_vector(y, TY, m, Y, incy, expr);                 \
        }                                                                      \
        else if (incx > 1) {                                                   \
            auto x = create_vector((TX*)X, n, incx);                           \
            tlapack_expr_with_vector(y, TY, m, Y, incy, expr);                 \
        }                                                                      \
        else {                                                                 \
            auto x = create_backward_vector((TX*)X, n, -incx);                 \
            tlapack_expr_with_vector(y, TY, m, Y, incy, expr);                 \
        }                                                                      \
    } while (false)
 
#define tlapack_expr_with_vector(x, TX, n, X, incx, expr)        \
    do {                                                         \
        using tlapack::legacy::internal::create_vector;          \
        using tlapack::legacy::internal::create_backward_vector; \
        if (incx == 1) {                                         \
            auto x = create_vector((TX*)X, n);                   \
            expr;                                                \
        }                                                        \
        else if (incx == -1) {                                   \
            auto x = create_backward_vector((TX*)X, n);          \
            expr;                                                \
        }                                                        \
        else if (incx > 1) {                                     \
            auto x = create_vector((TX*)X, n, incx);             \
            expr;                                                \
        }                                                        \
        else {                                                   \
            auto x = create_backward_vector((TX*)X, n, -incx);   \
            expr;                                                \
        }                                                        \
    } while (false)
 
#define tlapack_expr_with_vector_positiveInc(x, TX, n, X, incx, expr) \
    do {                                                              \
        using tlapack::legacy::internal::create_vector;               \
        if (incx == 1) {                                              \
            auto x = create_vector((TX*)X, n);                        \
            expr;                                                     \
        }                                                             \
        else {                                                        \
            auto x = create_vector((TX*)X, n, incx);                  \
            expr;                                                     \
        }                                                             \
    } while (false)
 
#endif  // TLAPACK_LEGACY_UTILS_HH