Mario,
yeah, there was a typo actually this was just a
dummy created to reproduce my problem. I corrected
and had a re run:
//-----
Full Integral: 10000
Box1: 2655 Box2: 2556 Box3: 2535 Box4: 2439
> >> // 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%
>
> > 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:36:07 CEST
This archive was generated by hypermail 2.2.0 : Tue Jun 02 2009 - 17:50:02 CEST