本文整理汇总了C++中MATRIX::find_element方法的典型用法代码示例。如果您正苦于以下问题：C++ MATRIX::find_element方法的具体用法？C++ MATRIX::find_element怎么用？C++ MATRIX::find_element使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类MATRIX的用法示例。
在下文中一共展示了MATRIX::find_element方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: incomplete_cholesky_decompose
size_t incomplete_cholesky_decompose(MATRIX& A)
{
using namespace ublas;
typedef typename MATRIX::value_type T;
// read access to a const matrix is faster
const MATRIX& A_c(A);
const size_t n = A.size1();
for (size_t k=0 ; k < n; k++) {
double qL_kk = A_c(k,k) - inner_prod( project( row( A_c, k ), range(0, k) ),
project( row( A_c, k ), range(0, k) ) );
if (qL_kk <= 0) {
return 1 + k;
} else {
double L_kk = sqrt( qL_kk );
// aktualisieren
for (size_t i = k+1; i < A.size1(); ++i) {
T* Aik = A.find_element(i, k);
if (Aik != 0) {
*Aik = ( *Aik - inner_prod( project( row( A_c, k ), range(0, k) ),
project( row( A_c, i ), range(0, k) ) ) ) / L_kk;
}
}
A(k,k) = L_kk;
}
}
return 0;
}