Factorization in VanLeerTheta
In Transport::VanLeerTheta, I came accros:
for (int i = 0; i < nr; i++) {
real surfr = 0.5*(Rsup2[i]-Rinf2[i])*dt;
for (int h = 0; h < ni; h++) {
real dphi = phi[h+1]-phi[h];
for (int j = 0; j < ns-1; j++) {
int l = j+h*ns+i*ni*ns;
int ljp = l+1;
qb[l] += invVolume[l]*surfr*dphi*(qrsXrhosXvt[l]-qrsXrhosXvt[ljp]);
....
By unrolling invVolume[l]*surfr*dphi
(and droppin dt) we get something like:
1.5*(radSup[i]^2-radInf[i]^2)*(phi[h+1]-phi[h])
/ (radSup[i]^3-radInf[i]^3)*(cosPhi[h]-cosPhi[h+1])*(theta[j+1]-theta[j])
Which I'm sure has some meaning. What is it ?
Thanks