jiri
/
LA_library
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

*** empty log message ***

This commit is contained in:
jiri 2005-09-07 20:40:54 +00:00
parent 19ad78d6aa
commit 6be91ad771
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
bisection.h
View File

@ -1,6 +1,6 @@
#ifndef _BISECTION_H
#define _BISECTION_H
//general bisection search
//general bisection search between dm and hm
//returns dm-1 on failure, otherwise number between dm and hm
//cmp returns 0 on equal, >0 if first > second argument
@ -39,15 +39,15 @@ int interpolation_find(INDEX dm0, INDEX high, const SUBJECT *x, const SUBJECT *b
{
d1=(*dist)(x,base+low*lead_dimension_base);
d2=(*dist)(base+high*lead_dimension_base,x);
cout << "intlp "<<low<<" "<<high<<" "<<d1<<" "<<d2<,endl
if(d1<=0 || d2 >=0) break;
cout << "intlp "<<low<<" "<<high<<" "<<d1<<" "<<d2<<endl;
if(d1<=0 || d2 <=0) break;
d3=(*dist)(base+high*lead_dimension_base,base+low*lead_dimension_base);
mid = low + (int)((d1 * (high-low)) / d3+.5);
mid = low + (int)((d1 * (high-low)) / d3);
if(mid<low||mid>high) laerror("interpolation_find: error in distance function");
d4=(*dist)(x,base+mid*lead_dimension_base);
cout << "intlp2 "<<mid<<" "<<d3<<" "<<d4<<endl;
if(d4>0)
low = mid + 1;
else if(d3<0)
else if(d4<0)
high = mid - 1;
else
return mid;