ComputeMdot weirdness
A few weird things in ComputeMdot:
- Division by itself (I think that's
1
?)
Work.values() = Rho.values()/Rho.values();
- Int the loop:
for (int i = 0; i < nr-1; i++) {
for (int h =0; h < ni; h++) {
real surfr = (cosPhi[h]-cosPhi[h+1]); // works because cosPhi goes one past the end,
for (int j = 0; j < ns; j++) {
int l=j+h*ns+i*ni*ns;
// this :
mdotmap[l] = 0.;
mdotmap[l] = mdotmap[l]+rhos[l]*qrs[l]*vr[l]*Rinf2[i]*surfr;
// is equivalent to:
mdotmap[l] = rhos[l]*qrs[l]*vr[l]*Rinf2[i]*surfr;
}
}
}