Commit f523f5e2 authored by Jean-Eric Campagne's avatar Jean-Eric Campagne
Browse files

(JEC/MR) 11/4/15 transpose the matrix LnkMtx in the LaguerreTransform::MultiAnalysis function

parent d7584663
......@@ -96,7 +96,8 @@ void LaguerreTransform::MultiAnalysis(const vector< complex<r_8> >& fi,
r_8 rk = nodes_[k];
lag.Values(rk,LnAll);
for (int n = 0; n<N_; n++ ){
LnkMtx[n*N_+k] = LnAll[n]*facts[n];
// LnkMtx[n*N_+k] = LnAll[n]*facts[n];
LnkMtx[n+N_*k] = LnAll[n]*facts[n];
}
}
......@@ -110,7 +111,8 @@ void LaguerreTransform::MultiAnalysis(const vector< complex<r_8> >& fi,
complex<r_8> fli = fi[l+ i*stride];
r_8 wi = weights_[i];
for (int n=0; n<N_; n++){
vtmp[n] += fli * wi * LnkMtx[n*N_+i];
// vtmp[n] += fli * wi * LnkMtx[n*N_+i];
vtmp[n] += fli * wi * LnkMtx[n+N_*i];
}//loop on k
}//loop on n
for (int n=0; n<N_; n++) {
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment