>> // 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.
Ehh ... it's even simpler :) You have a typo. Your Box3 and Box4 are identical regions. You should be integrating from y_cut to 1 in Box4 to cover the whole range.
Mario Received on Tue Jun 02 2009 - 15:29:21 CEST
This archive was generated by hypermail 2.2.0 : Tue Jun 02 2009 - 17:50:02 CEST