"Fossies" - the Fresh Open Source Software Archive  

Source code changes of the file "include/armadillo_bits/op_det_meat.hpp" between
armadillo-10.8.2.tar.xz and armadillo-11.0.0.tar.xz

About: Armadillo is a C++ linear algebra library (matrix maths) aiming towards a good balance between speed and ease of use.

[ To the main armadillo source changes report ]

op_det_meat.hpp  (armadillo-10.8.2.tar.xz):op_det_meat.hpp  (armadillo-11.0.0.tar.xz)
skipping to change at line 53 skipping to change at line 53
out_val = op_det::apply_trimat(strip.M); out_val = op_det::apply_trimat(strip.M);
return true; return true;
} }
Mat<eT> A(expr.get_ref()); Mat<eT> A(expr.get_ref());
arma_debug_check( (A.is_square() == false), "det(): given matrix must be squar e sized" ); arma_debug_check( (A.is_square() == false), "det(): given matrix must be squar e sized" );
if((A.n_rows <= 4) && is_cx<eT>::no) const uword N = A.n_rows;
if(N == 0) { out_val = eT(1); return true; }
if(N == 1) { out_val = A[0]; return true; }
if((is_cx<eT>::no) && (N <= 4))
{ {
constexpr T det_min = std::numeric_limits<T>::epsilon(); constexpr T det_min = std::numeric_limits<T>::epsilon();
constexpr T det_max = T(1) / std::numeric_limits<T>::epsilon(); constexpr T det_max = T(1) / std::numeric_limits<T>::epsilon();
const eT det_val = op_det::apply_tiny(A); eT det_val = eT(0);
const T abs_det_val = std::abs(det_val);
if(N == 2) { det_val = op_det::apply_tiny_2x2(A); }
if(N == 3) { det_val = op_det::apply_tiny_3x3(A); }
if(N == 4) { det_val = op_det::apply_tiny_4x4(A); }
const T abs_det_val = std::abs(det_val);
if((abs_det_val > det_min) && (abs_det_val < det_max)) { out_val = det_val; return true; } if((abs_det_val > det_min) && (abs_det_val < det_max)) { out_val = det_val; return true; }
// fallthrough if det_val is suspect // fallthrough if det_val is suspect
} }
if(A.is_diagmat()) { out_val = op_det::apply_diagmat(A); return true; } if(A.is_diagmat()) { out_val = op_det::apply_diagmat(A); return true; }
const bool is_triu = trimat_helper::is_triu(A); const bool is_triu = trimat_helper::is_triu(A);
const bool is_tril = is_triu ? false : trimat_helper::is_tril(A); const bool is_tril = is_triu ? false : trimat_helper::is_tril(A);
skipping to change at line 124 skipping to change at line 134
for(uword i=0; i<N; ++i) { val *= P.at(i,i); } for(uword i=0; i<N; ++i) { val *= P.at(i,i); }
return val; return val;
} }
template<typename eT> template<typename eT>
arma_cold arma_cold
inline inline
eT eT
op_det::apply_tiny(const Mat<eT>& X) op_det::apply_tiny_2x2(const Mat<eT>& X)
{ {
arma_extra_debug_sigprint(); arma_extra_debug_sigprint();
// NOTE: assuming matrix X is square sized const eT* Xm = X.memptr();
const uword N = X.n_rows; return ( Xm[pos<0,0>::n2]*Xm[pos<1,1>::n2] - Xm[pos<0,1>::n2]*Xm[pos<1,0>::n2]
const eT* Xm = X.memptr(); );
}
if(N == 0) { return eT(1); } template<typename eT>
arma_cold
inline
eT
op_det::apply_tiny_3x3(const Mat<eT>& X)
{
arma_extra_debug_sigprint();
if(N == 1) { return Xm[0]; } const eT* Xm = X.memptr();
if(N == 2) // const double tmp1 = X.at(0,0) * X.at(1,1) * X.at(2,2);
{ // const double tmp2 = X.at(0,1) * X.at(1,2) * X.at(2,0);
return ( Xm[pos<0,0>::n2]*Xm[pos<1,1>::n2] - Xm[pos<0,1>::n2]*Xm[pos<1,0>::n // const double tmp3 = X.at(0,2) * X.at(1,0) * X.at(2,1);
2] ); // const double tmp4 = X.at(2,0) * X.at(1,1) * X.at(0,2);
} // const double tmp5 = X.at(2,1) * X.at(1,2) * X.at(0,0);
// const double tmp6 = X.at(2,2) * X.at(1,0) * X.at(0,1);
// return (tmp1+tmp2+tmp3) - (tmp4+tmp5+tmp6);
const eT val1 = Xm[pos<0,0>::n3]*(Xm[pos<2,2>::n3]*Xm[pos<1,1>::n3] - Xm[pos<2
,1>::n3]*Xm[pos<1,2>::n3]);
const eT val2 = Xm[pos<1,0>::n3]*(Xm[pos<2,2>::n3]*Xm[pos<0,1>::n3] - Xm[pos<2
,1>::n3]*Xm[pos<0,2>::n3]);
const eT val3 = Xm[pos<2,0>::n3]*(Xm[pos<1,2>::n3]*Xm[pos<0,1>::n3] - Xm[pos<1
,1>::n3]*Xm[pos<0,2>::n3]);
if(N == 3) return ( val1 - val2 + val3 );
{ }
// const double tmp1 = X.at(0,0) * X.at(1,1) * X.at(2,2);
// const double tmp2 = X.at(0,1) * X.at(1,2) * X.at(2,0);
// const double tmp3 = X.at(0,2) * X.at(1,0) * X.at(2,1);
// const double tmp4 = X.at(2,0) * X.at(1,1) * X.at(0,2);
// const double tmp5 = X.at(2,1) * X.at(1,2) * X.at(0,0);
// const double tmp6 = X.at(2,2) * X.at(1,0) * X.at(0,1);
// return (tmp1+tmp2+tmp3) - (tmp4+tmp5+tmp6);
const eT val1 = Xm[pos<0,0>::n3]*(Xm[pos<2,2>::n3]*Xm[pos<1,1>::n3] - Xm[pos
<2,1>::n3]*Xm[pos<1,2>::n3]);
const eT val2 = Xm[pos<1,0>::n3]*(Xm[pos<2,2>::n3]*Xm[pos<0,1>::n3] - Xm[pos
<2,1>::n3]*Xm[pos<0,2>::n3]);
const eT val3 = Xm[pos<2,0>::n3]*(Xm[pos<1,2>::n3]*Xm[pos<0,1>::n3] - Xm[pos
<1,1>::n3]*Xm[pos<0,2>::n3]);
return ( val1 - val2 + val3 ); template<typename eT>
} arma_cold
inline
eT
op_det::apply_tiny_4x4(const Mat<eT>& X)
{
arma_extra_debug_sigprint();
if(N == 4) const eT* Xm = X.memptr();
{
const eT val_03_12 = Xm[pos<0,3>::n4] * Xm[pos<1,2>::n4];
const eT val_02_13 = Xm[pos<0,2>::n4] * Xm[pos<1,3>::n4];
const eT val_03_11 = Xm[pos<0,3>::n4] * Xm[pos<1,1>::n4];
const eT val_01_13 = Xm[pos<0,1>::n4] * Xm[pos<1,3>::n4];
const eT val_02_11 = Xm[pos<0,2>::n4] * Xm[pos<1,1>::n4];
const eT val_01_12 = Xm[pos<0,1>::n4] * Xm[pos<1,2>::n4];
const eT val_03_10 = Xm[pos<0,3>::n4] * Xm[pos<1,0>::n4];
const eT val_00_13 = Xm[pos<0,0>::n4] * Xm[pos<1,3>::n4];
const eT val_02_10 = Xm[pos<0,2>::n4] * Xm[pos<1,0>::n4];
const eT val_00_12 = Xm[pos<0,0>::n4] * Xm[pos<1,2>::n4];
const eT val_01_10 = Xm[pos<0,1>::n4] * Xm[pos<1,0>::n4];
const eT val_00_11 = Xm[pos<0,0>::n4] * Xm[pos<1,1>::n4];
const eT val_21_30 = Xm[pos<2,1>::n4] * Xm[pos<3,0>::n4];
const eT val_22_30 = Xm[pos<2,2>::n4] * Xm[pos<3,0>::n4];
const eT val_23_30 = Xm[pos<2,3>::n4] * Xm[pos<3,0>::n4];
const eT val_20_31 = Xm[pos<2,0>::n4] * Xm[pos<3,1>::n4];
const eT val_22_31 = Xm[pos<2,2>::n4] * Xm[pos<3,1>::n4];
const eT val_23_31 = Xm[pos<2,3>::n4] * Xm[pos<3,1>::n4];
const eT val_20_32 = Xm[pos<2,0>::n4] * Xm[pos<3,2>::n4];
const eT val_21_32 = Xm[pos<2,1>::n4] * Xm[pos<3,2>::n4];
const eT val_23_32 = Xm[pos<2,3>::n4] * Xm[pos<3,2>::n4];
const eT val_20_33 = Xm[pos<2,0>::n4] * Xm[pos<3,3>::n4];
const eT val_21_33 = Xm[pos<2,1>::n4] * Xm[pos<3,3>::n4];
const eT val_22_33 = Xm[pos<2,2>::n4] * Xm[pos<3,3>::n4];
const eT val = \
val_03_12 * val_21_30 \
- val_02_13 * val_21_30 \
- val_03_11 * val_22_30 \
+ val_01_13 * val_22_30 \
+ val_02_11 * val_23_30 \
- val_01_12 * val_23_30 \
- val_03_12 * val_20_31 \
+ val_02_13 * val_20_31 \
+ val_03_10 * val_22_31 \
- val_00_13 * val_22_31 \
- val_02_10 * val_23_31 \
+ val_00_12 * val_23_31 \
+ val_03_11 * val_20_32 \
- val_01_13 * val_20_32 \
- val_03_10 * val_21_32 \
+ val_00_13 * val_21_32 \
+ val_01_10 * val_23_32 \
- val_00_11 * val_23_32 \
- val_02_11 * val_20_33 \
+ val_01_12 * val_20_33 \
+ val_02_10 * val_21_33 \
- val_00_12 * val_21_33 \
- val_01_10 * val_22_33 \
+ val_00_11 * val_22_33 \
;
return val; const eT val_03_12 = Xm[pos<0,3>::n4] * Xm[pos<1,2>::n4];
} const eT val_02_13 = Xm[pos<0,2>::n4] * Xm[pos<1,3>::n4];
const eT val_03_11 = Xm[pos<0,3>::n4] * Xm[pos<1,1>::n4];
const eT val_01_13 = Xm[pos<0,1>::n4] * Xm[pos<1,3>::n4];
const eT val_02_11 = Xm[pos<0,2>::n4] * Xm[pos<1,1>::n4];
const eT val_01_12 = Xm[pos<0,1>::n4] * Xm[pos<1,2>::n4];
const eT val_03_10 = Xm[pos<0,3>::n4] * Xm[pos<1,0>::n4];
const eT val_00_13 = Xm[pos<0,0>::n4] * Xm[pos<1,3>::n4];
const eT val_02_10 = Xm[pos<0,2>::n4] * Xm[pos<1,0>::n4];
const eT val_00_12 = Xm[pos<0,0>::n4] * Xm[pos<1,2>::n4];
const eT val_01_10 = Xm[pos<0,1>::n4] * Xm[pos<1,0>::n4];
const eT val_00_11 = Xm[pos<0,0>::n4] * Xm[pos<1,1>::n4];
const eT val_21_30 = Xm[pos<2,1>::n4] * Xm[pos<3,0>::n4];
const eT val_22_30 = Xm[pos<2,2>::n4] * Xm[pos<3,0>::n4];
const eT val_23_30 = Xm[pos<2,3>::n4] * Xm[pos<3,0>::n4];
const eT val_20_31 = Xm[pos<2,0>::n4] * Xm[pos<3,1>::n4];
const eT val_22_31 = Xm[pos<2,2>::n4] * Xm[pos<3,1>::n4];
const eT val_23_31 = Xm[pos<2,3>::n4] * Xm[pos<3,1>::n4];
const eT val_20_32 = Xm[pos<2,0>::n4] * Xm[pos<3,2>::n4];
const eT val_21_32 = Xm[pos<2,1>::n4] * Xm[pos<3,2>::n4];
const eT val_23_32 = Xm[pos<2,3>::n4] * Xm[pos<3,2>::n4];
const eT val_20_33 = Xm[pos<2,0>::n4] * Xm[pos<3,3>::n4];
const eT val_21_33 = Xm[pos<2,1>::n4] * Xm[pos<3,3>::n4];
const eT val_22_33 = Xm[pos<2,2>::n4] * Xm[pos<3,3>::n4];
const eT val = \
val_03_12 * val_21_30 \
- val_02_13 * val_21_30 \
- val_03_11 * val_22_30 \
+ val_01_13 * val_22_30 \
+ val_02_11 * val_23_30 \
- val_01_12 * val_23_30 \
- val_03_12 * val_20_31 \
+ val_02_13 * val_20_31 \
+ val_03_10 * val_22_31 \
- val_00_13 * val_22_31 \
- val_02_10 * val_23_31 \
+ val_00_12 * val_23_31 \
+ val_03_11 * val_20_32 \
- val_01_13 * val_20_32 \
- val_03_10 * val_21_32 \
+ val_00_13 * val_21_32 \
+ val_01_10 * val_23_32 \
- val_00_11 * val_23_32 \
- val_02_11 * val_20_33 \
+ val_01_12 * val_20_33 \
+ val_02_10 * val_21_33 \
- val_00_12 * val_21_33 \
- val_01_10 * val_22_33 \
+ val_00_11 * val_22_33 \
;
return eT(0); return val;
} }
//! @} //! @}
 End of changes. 13 change blocks. 
95 lines changed or deleted 109 lines changed or added
To the Header Files section of the armadillo source changes report or the top of this page
Mainly generated by pkgdiff using rfcdiff
Home  |  About  |  Features  |  All  |  Newest  |  Dox  |  Diffs  |  RSS Feeds  |  Screenshots  |  Comments  |  Imprint  |  Privacy  |  HTTP(S)