test cbrt func

BloubWorld/essai.f90

@@ -1,17 +1,11 @@
 program essai
 
-  use bloubspace
+!  use bloubspace
   use mathstuff
   implicit none
 
-  type(t_boundingbox)    :: bbox
 
-  call load_boundingbox("WS/boundinboxes.dat", bbox, "cube    ")
-
-  print *, bbox
-
-
-  ! call test_random(20)
+  call test_random(10)
 
 
   STOP ': BECAUSE JOB IS DONE'
@@ -21,17 +15,17 @@ program essai
 
 subroutine test_random(nbre)
   integer, intent(in)         :: nbre
-  integer                     :: foo, bar
-  real                        :: quux
+  integer                     :: foo
+  real                        :: quux, bar
   double precision            :: somme
 
   call init_random_seed()           ! in module 'mathstuff'
   somme = 0.0
   do foo=1, nbre
-    quux = rand()
+    quux = 10.0 * rand()
     somme = somme + quux
-    bar = mod(irand(), 7)
-    print *, foo, quux, somme/foo, bar
+    bar = quux ** (.1/.3)
+    print *, quux, bar, somme/foo
   enddo
 end subroutine test_random
 ! --------------------------------------------------------------