Chris@49
|
1 // Copyright (C) 2011-2012 NICTA (www.nicta.com.au)
|
Chris@49
|
2 // Copyright (C) 2011-2012 Conrad Sanderson
|
Chris@49
|
3 //
|
Chris@49
|
4 // This Source Code Form is subject to the terms of the Mozilla Public
|
Chris@49
|
5 // License, v. 2.0. If a copy of the MPL was not distributed with this
|
Chris@49
|
6 // file, You can obtain one at http://mozilla.org/MPL/2.0/.
|
Chris@49
|
7
|
Chris@49
|
8
|
Chris@49
|
9 //! \addtogroup fn_strans
|
Chris@49
|
10 //! @{
|
Chris@49
|
11
|
Chris@49
|
12
|
Chris@49
|
13
|
Chris@49
|
14 template<typename T1>
|
Chris@49
|
15 arma_inline
|
Chris@49
|
16 const Op<T1, op_strans>
|
Chris@49
|
17 strans
|
Chris@49
|
18 (
|
Chris@49
|
19 const T1& X,
|
Chris@49
|
20 const typename enable_if< is_arma_type<T1>::value == true >::result* junk1 = 0,
|
Chris@49
|
21 const typename arma_cx_only<typename T1::elem_type>::result* junk2 = 0
|
Chris@49
|
22 )
|
Chris@49
|
23 {
|
Chris@49
|
24 arma_extra_debug_sigprint();
|
Chris@49
|
25 arma_ignore(junk1);
|
Chris@49
|
26 arma_ignore(junk2);
|
Chris@49
|
27
|
Chris@49
|
28 return Op<T1, op_strans>(X);
|
Chris@49
|
29 }
|
Chris@49
|
30
|
Chris@49
|
31
|
Chris@49
|
32
|
Chris@49
|
33 // NOTE: for non-complex objects, deliberately returning op_htrans instead of op_strans,
|
Chris@49
|
34 // NOTE: due to currently more optimisations available when using op_htrans, especially by glue_times
|
Chris@49
|
35 template<typename T1>
|
Chris@49
|
36 arma_inline
|
Chris@49
|
37 const Op<T1, op_htrans>
|
Chris@49
|
38 strans
|
Chris@49
|
39 (
|
Chris@49
|
40 const T1& X,
|
Chris@49
|
41 const typename enable_if< is_arma_type<T1>::value == true >::result* junk1 = 0,
|
Chris@49
|
42 const typename arma_not_cx<typename T1::elem_type>::result* junk2 = 0
|
Chris@49
|
43 )
|
Chris@49
|
44 {
|
Chris@49
|
45 arma_extra_debug_sigprint();
|
Chris@49
|
46 arma_ignore(junk1);
|
Chris@49
|
47 arma_ignore(junk2);
|
Chris@49
|
48
|
Chris@49
|
49 return Op<T1, op_htrans>(X);
|
Chris@49
|
50 }
|
Chris@49
|
51
|
Chris@49
|
52
|
Chris@49
|
53
|
Chris@49
|
54 //! two consecutive transpose operations cancel each other
|
Chris@49
|
55 template<typename T1>
|
Chris@49
|
56 arma_inline
|
Chris@49
|
57 const T1&
|
Chris@49
|
58 strans(const Op<T1, op_strans>& X)
|
Chris@49
|
59 {
|
Chris@49
|
60 arma_extra_debug_sigprint();
|
Chris@49
|
61 arma_extra_debug_print("strans(): removing op_strans");
|
Chris@49
|
62
|
Chris@49
|
63 return X.m;
|
Chris@49
|
64 }
|
Chris@49
|
65
|
Chris@49
|
66
|
Chris@49
|
67
|
Chris@49
|
68 //
|
Chris@49
|
69 // handling of sparse matrices
|
Chris@49
|
70
|
Chris@49
|
71
|
Chris@49
|
72 template<typename T1>
|
Chris@49
|
73 inline
|
Chris@49
|
74 typename
|
Chris@49
|
75 enable_if2
|
Chris@49
|
76 <
|
Chris@49
|
77 is_arma_sparse_type<T1>::value,
|
Chris@49
|
78 const SpOp<T1,spop_strans>
|
Chris@49
|
79 >::result
|
Chris@49
|
80 strans(const T1& x)
|
Chris@49
|
81 {
|
Chris@49
|
82 arma_extra_debug_sigprint();
|
Chris@49
|
83
|
Chris@49
|
84 return SpOp<T1,spop_strans>(x);
|
Chris@49
|
85 }
|
Chris@49
|
86
|
Chris@49
|
87
|
Chris@49
|
88
|
Chris@49
|
89 //! @}
|