max@0: // Copyright (C) 2010-2011 NICTA (www.nicta.com.au)
max@0: // Copyright (C) 2010-2011 Conrad Sanderson
max@0: // 
max@0: // This file is part of the Armadillo C++ library.
max@0: // It is provided without any warranty of fitness
max@0: // for any purpose. You can redistribute this file
max@0: // and/or modify it under the terms of the GNU
max@0: // Lesser General Public License (LGPL) as published
max@0: // by the Free Software Foundation, either version 3
max@0: // of the License or (at your option) any later version.
max@0: // (see http://www.opensource.org/licenses for more info)
max@0: 
max@0: 
max@0: 
max@0: //! \addtogroup glue_cross
max@0: //! @{
max@0: 
max@0: 
max@0: 
max@0: template<typename T1, typename T2>
max@0: inline
max@0: void
max@0: glue_cross::apply(Mat<typename T1::elem_type>& out, const Glue<T1, T2, glue_cross>& X)
max@0:   {
max@0:   arma_extra_debug_sigprint();
max@0:   
max@0:   typedef typename T1::elem_type      eT;
max@0:   typedef typename Proxy<T1>::ea_type ea_type1;
max@0:   typedef typename Proxy<T2>::ea_type ea_type2;
max@0:   
max@0:   const Proxy<T1> A(X.A);
max@0:   const Proxy<T2> B(X.B);
max@0:   
max@0:   arma_debug_check( ((A.get_n_elem() != 3) || (B.get_n_elem() != 3)), "cross(): input vectors must have 3 elements" );
max@0:   
max@0:   out.set_size(A.get_n_rows(), A.get_n_cols());
max@0:   
max@0:   eT*      out_mem = out.memptr();
max@0:   ea_type1 PA      = A.get_ea();
max@0:   ea_type2 PB      = B.get_ea();
max@0:   
max@0:   const eT ax = PA[0];
max@0:   const eT ay = PA[1];
max@0:   const eT az = PA[2];
max@0:   
max@0:   const eT bx = PB[0];
max@0:   const eT by = PB[1];
max@0:   const eT bz = PB[2];
max@0:   
max@0:   out_mem[0] = ay*bz - az*by;
max@0:   out_mem[1] = az*bx - ax*bz;
max@0:   out_mem[2] = ax*by - ay*bx;
max@0:   }
max@0: 
max@0: 
max@0: 
max@0: //! @}