> the integral for each box, and add those values. Contrary to my
> expectations, the value does not turn out to be same as that
> obtained by integrating full histo.
<snip>
> // Box1
> float box1=BoxIntegral(h2,0,x_cut,0,y_cut);
>
> //Box2
> float box2=BoxIntegral(h2,0,x_cut,y_cut,1);
>
> //Box3
> float box3=BoxIntegral(h2,x_cut,1,0,y_cut);
>
> //Box4
> float box4=BoxIntegral(h2,x_cut,1,0,y_cut);
Aren't you taking a small overlap in the region where x_cut and y_cut lines overlap? I'd guess that if you exclude the overlap region you could well end up with the right result. Though the 25% increase is a bit more than I'd have expected from the overlap in case of 100 x 100 matrix.
Mario Received on Tue Jun 02 2009 - 15:27:45 CEST
This archive was generated by hypermail 2.2.0 : Tue Jun 02 2009 - 17:50:02 CEST