bugfixes is random

la_random.cc

@@ -5,7 +5,7 @@ namespace LA {
 
 WEAK_SYMBOL  double randdouble() {return random()/(1.+RAND_MAX);}
 WEAK_SYMBOL  double randdoublesigned() {return 2.*random()/(1.+RAND_MAX)-1.;}
-WEAK_SYMBOL  int randint32() {return random();}
-WEAK_SYMBOL  uint64_t randint64() {uint64_t r = random(); r<<=32; r|= random(); return r; }
+WEAK_SYMBOL  unsigned int randint32() {unsigned int r = random(); r<<=1; r^= random(); return r;}
+WEAK_SYMBOL  uint64_t randint64() {uint64_t r = random(); r<<=31; r|= random(); r<<=2; r|= random()&3; return r; } //random() gives 31 bits
 }//namespace
 

la_random.h

@@ -7,7 +7,7 @@ namespace LA {
 
 extern double randdouble();
 extern double randdoublesigned();
-extern int randint32();
+extern unsigned int randint32();
 extern uint64_t randint64();
 
 //RANDOM numbers defaulting to standard library but switchable to user's functions