12345qiupeng
/
gtsam
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

18% performance boost by using pointer to address destination matrix A

release/4.3a0
Chris Beall 2009-12-11 23:59:37 +00:00
parent 631248f902
commit be10fc2a90
1 changed files with 6 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
cpp/Matrix.cpp
View File

@ -345,8 +345,12 @@ weighted_eliminate(Matrix& A, Vector& b, const Vector& sigmas) {
for (int i=0;i<m;++i) { // update all rows for (int i=0;i<m;++i) { // update all rows
double ai = a(i); double ai = a(i);
b(i) -= ai*d; b(i) -= ai*d;
for (int j2=j+1;j2<n;++j2) // limit to only columns in separator double *Aptr = A.data().begin()+i*n+j+1;
A(i,j2) -= ai*r(j2); for (int j2=j+1;j2<n;++j2){ // limit to only columns in separator
//A(i,j2) -= ai*r(j2);
*Aptr -= ai*r(j2);
Aptr++;
}
} }
if (verbose) print(sub(A,0,m,j+1,n), "updated A"); if (verbose) print(sub(A,0,m,j+1,n), "updated A");
if (verbose) print(b, "updated b "); if (verbose) print(b, "updated b ");