1 files changed, 87 insertions, 0 deletions
diff --git a/src/armadillo/include/armadillo_bits/xtrans_mat_meat.hpp b/src/armadillo/include/armadillo_bits/xtrans_mat_meat.hpp
new file mode 100644
index 0000000..1872c30
--- /dev/null
+++ b/src/armadillo/include/armadillo_bits/xtrans_mat_meat.hpp
@@ -0,0 +1,87 @@
+// SPDX-License-Identifier: Apache-2.0
+// 
+// Copyright 2008-2016 Conrad Sanderson (http://conradsanderson.id.au)
+// Copyright 2008-2016 National ICT Australia (NICTA)
+// 
+// Licensed under the Apache License, Version 2.0 (the "License");
+// you may not use this file except in compliance with the License.
+// You may obtain a copy of the License at
+// http://www.apache.org/licenses/LICENSE-2.0
+// 
+// Unless required by applicable law or agreed to in writing, software
+// distributed under the License is distributed on an "AS IS" BASIS,
+// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+// See the License for the specific language governing permissions and
+// limitations under the License.
+// ------------------------------------------------------------------------
+
+
+//! \addtogroup xtrans_mat
+//! @{
+
+
+template<typename eT, bool do_conj>
+inline
+xtrans_mat<eT,do_conj>::xtrans_mat(const Mat<eT>& in_X)
+  : X     (in_X       )
+  , n_rows(in_X.n_cols)  // deliberately swapped
+  , n_cols(in_X.n_rows)
+  , n_elem(in_X.n_elem)
+  {
+  arma_extra_debug_sigprint();
+  }
+
+
+
+template<typename eT, bool do_conj>
+inline
+void
+xtrans_mat<eT,do_conj>::extract(Mat<eT>& out) const
+  {
+  arma_extra_debug_sigprint();
+  
+  really_do_conj ? op_htrans::apply_mat(out, X) : op_strans::apply_mat(out, X);
+  }
+
+
+
+template<typename eT, bool do_conj>
+inline
+eT
+xtrans_mat<eT,do_conj>::operator[](const uword ii) const
+  {
+  if(Y.n_elem > 0)
+    {
+    return Y[ii];
+    }
+  else
+    {
+    really_do_conj ? op_htrans::apply_mat(Y, X) : op_strans::apply_mat(Y, X);
+    return Y[ii];
+    }
+  }
+
+
+
+template<typename eT, bool do_conj>
+inline
+eT
+xtrans_mat<eT,do_conj>::at_alt(const uword ii) const
+  {
+  return (*this).operator[](ii);
+  }
+
+
+
+template<typename eT, bool do_conj>
+arma_inline
+eT
+xtrans_mat<eT,do_conj>::at(const uword in_row, const uword in_col) const
+  {
+  return really_do_conj ? eT(access::alt_conj(X.at(in_col, in_row))) : eT(X.at(in_col, in_row));
+  // in_row and in_col deliberately swapped above
+  }
+
+
+
+//! @}