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

small performance improvement by using a pointer to address r

release/4.3a0
Chris Beall 2009-12-14 16:56:50 +00:00
parent 71bc9ca0f7
commit a7a5e5e816
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
cpp/Matrix.cpp
View File

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