fixed rand() bug in randfib=2

tractvolsx.h

@@ -128,26 +128,26 @@ namespace TRACTVOLSX{
 	
 	///////new xyz values from probabilistic interpolation
 	int newx,newy,newz; 
-	float tmp=rand(); tmp/=float(RAND_MAX);
+	float tmp=(float)rand()/float(RAND_MAX);
 	if(tmp>pcx)
 	  newx=fx;
 	else
 	  newx=cx;
 	
-	tmp=rand(); tmp/=float(RAND_MAX);
+	tmp=(float)rand()/float(RAND_MAX);
 	if(tmp>pcy)
 	  newy=fy;
 	else
 	  newy=cy;
 	
-	tmp=rand(); tmp/=float(RAND_MAX);
+	tmp=(float)rand()/float(RAND_MAX);
 	if(tmp>pcz)
 	  newz=fz;
 	else
 	  newz=cz;
  
 	ColumnVector th_ph_f(3);	
-	float samp=rand(); samp/=float(RAND_MAX);
+	float samp=(float)rand()/float(RAND_MAX);
 	samp=round(samp*((*thsamples[0]).tsize()-1));
 	float theta=0,phi=0;
 	float dotmax=0,dottmp=0;
@@ -167,7 +167,7 @@ namespace TRACTVOLSX{
 		fibst=0;
 	      }
 	      else{
-		float rtmp=rand()/float(RAND_MAX) * float(fibvec.size()-1);
+		float rtmp=(float)rand()/float(RAND_MAX) * float(fibvec.size()-1);
 		fibst = fibvec[ (int)round(rtmp) ];	      
 	      }
 	      
@@ -188,9 +188,9 @@ namespace TRACTVOLSX{
 	      else{
 		float fsumtmp2=0;
 		int fib=0;
-		float rtmp=rand()/float(RAND_MAX);
+		float rtmp=(float)rand()/float(RAND_MAX);
 		
-		while( fsumtmp2<rtmp){
+		while( fsumtmp2<=rtmp){
 		  float ft=(*fsamples[fib])(int(newx),int(newy),int(newz),int(samp));
 		  if(ft>opts.fibthresh.value()){
 		    fsumtmp2+=(ft/fsumtmp);