rotate.h

Go to the documentation of this file.

// ========================================================================= //
// This file is part of MyraMath, copyright (c) 2014-2019 by Ryan A Chilton  //
// and distributed by MyraCore, LLC. See LICENSE.txt for license terms.      //
// ========================================================================= //

#ifndef MYRAMATH_DENSE_ROTATE_H
#define MYRAMATH_DENSE_ROTATE_H

// For number types.
#include <myramath/MYRAMATH_EXPORT.h>
#include <myramath/utility/Number.h>

// For std::pair
#include <utility>

namespace myra {

// Forward declarations.
template<class Number> class Matrix22;

MYRAMATH_EXPORT std::pair<NumberS, NumberS> rotate1(NumberS a0, NumberS a1);
MYRAMATH_EXPORT std::pair<NumberD, NumberD> rotate1(NumberD a0, NumberD a1);
MYRAMATH_EXPORT std::pair<NumberC, NumberC> rotate1(NumberC a0, NumberC a1);
MYRAMATH_EXPORT std::pair<NumberZ, NumberZ> rotate1(NumberZ a0, NumberZ a1);

MYRAMATH_EXPORT Matrix22<NumberS> make_Matrix22(NumberS a0, NumberS a1);
MYRAMATH_EXPORT Matrix22<NumberD> make_Matrix22(NumberD a0, NumberD a1);
MYRAMATH_EXPORT Matrix22<NumberC> make_Matrix22(NumberC a0, NumberC a1);
MYRAMATH_EXPORT Matrix22<NumberZ> make_Matrix22(NumberZ a0, NumberZ a1);

MYRAMATH_EXPORT Matrix22<NumberS> sy_rotate2(NumberS a00, NumberS a11, NumberS a10);
MYRAMATH_EXPORT Matrix22<NumberD> sy_rotate2(NumberD a00, NumberD a11, NumberD a10);
MYRAMATH_EXPORT Matrix22<NumberC> sy_rotate2(NumberC a00, NumberC a11, NumberC a10);
MYRAMATH_EXPORT Matrix22<NumberZ> sy_rotate2(NumberZ a00, NumberZ a11, NumberZ a10);

MYRAMATH_EXPORT Matrix22<NumberS> sy_rotate2(const Matrix22<NumberS>& A);
MYRAMATH_EXPORT Matrix22<NumberD> sy_rotate2(const Matrix22<NumberD>& A);
MYRAMATH_EXPORT Matrix22<NumberC> sy_rotate2(const Matrix22<NumberC>& A);
MYRAMATH_EXPORT Matrix22<NumberZ> sy_rotate2(const Matrix22<NumberZ>& A);

MYRAMATH_EXPORT Matrix22<NumberC> he_rotate2(NumberS a00, NumberS a11, NumberC a10);
MYRAMATH_EXPORT Matrix22<NumberZ> he_rotate2(NumberD a00, NumberD a11, NumberZ a10);

MYRAMATH_EXPORT Matrix22<NumberC> he_rotate2(const Matrix22<NumberC>& A);
MYRAMATH_EXPORT Matrix22<NumberZ> he_rotate2(const Matrix22<NumberZ>& A);

template<class Number> class gesvd2_output
  {
  public:
  
    typedef typename ReflectPrecision<Number>::type Precision;

    Matrix22<Number> U;    // Unitary U, U'U = I
    Matrix22<Precision> S; // Real diagonal S, note S(1,0) = S(0,1) = 0
    Matrix22<Number> V;    // Unitary V, VV' = I
  };

MYRAMATH_EXPORT gesvd2_output<NumberS> gesvd2(const Matrix22<NumberS>& A);
MYRAMATH_EXPORT gesvd2_output<NumberD> gesvd2(const Matrix22<NumberD>& A);
MYRAMATH_EXPORT gesvd2_output<NumberC> gesvd2(const Matrix22<NumberC>& A);
MYRAMATH_EXPORT gesvd2_output<NumberZ> gesvd2(const Matrix22<NumberZ>& A);

} // namespace

#endif