Commit dae0aea9 authored by Johan Richard's avatar Johan Richard
Browse files

Corrected bug in flux from cleanset 4

parent 3e7d4b30
......@@ -73,11 +73,11 @@ void do_itos(double **im, struct pixlist *pl, int npl, double dlsds, double z
if (is >= 0 && js >= 0 && is < ps.ny && js < ps.nx)
{
imult[is][js] += 1;
source[is][js] = ech*((imult[is][js] - 1) * source[is][js] +
pl[k].flux)
source[is][js] = ((imult[is][js] - 1) * source[is][js] +
pl[k].flux*ech)
/ ((double)imult[is][js]);
erreur[is][js] = ech*((imult[is][js] - 1) * erreur[is][js] +
pl[k].flux * pl[k].flux)
erreur[is][js] = ((imult[is][js] - 1) * erreur[is][js] +
pl[k].flux * pl[k].flux*ech*ech)
/ ((double)imult[is][js]);
}
......@@ -113,11 +113,11 @@ void do_itos(double **im, struct pixlist *pl, int npl, double dlsds, double z
if ((is >= 0) && (js >= 0) && (is < ps.ny) && (js < ps.nx))
{
imult[is][js] += 1;
source[is][js] = ech*((imult[is][js] - 1) * source[is][js] +
pl[k].flux)
source[is][js] = ((imult[is][js] - 1) * source[is][js] +
pl[k].flux*ech)
/ ((double)imult[is][js]);
erreur[is][js] = ech*((imult[is][js] - 1) * erreur[is][js] +
pl[k].flux * pl[k].flux)
erreur[is][js] = ((imult[is][js] - 1) * erreur[is][js] +
pl[k].flux * pl[k].flux*ech*ech)
/ ((double)imult[is][js]);
};
A.x += imFrame.pixelx / imFrame.ech * ((double) j);
......
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