
Nicholas Andrews wrote:
void multi_state_kf::compute_gain(const matrix<double>& dev) { matrix<double> temp=prod(M,newP); matrix<double> transM=trans(M); matrix<double> A=(prod(temp,transM)+dev); permutation_matrix
pm(A.size2()); lu_factorize(A, pm); matrix<double> el = prod(newP,transM); lu_substitute(A, pm, el); K=el; }
This function was completely hacked together as I'm new to uBLAS (and rusty with linear algebra), but it got rid of most of my compile errors. Now only a few remain but I'm not sure what I'm doing wrong. I'm hoping that it will be obvious to someone here.
I tried to recreate your problem. Supposing that you defined somewhere
in the enclosing scope newP, K and M you can easily compile the above
sample by including
#include