add smilei_omp_threads in namelist

benchmarks/gpu/tst3d_gpu_o2_thermal_plasma_medium.py

@@ -66,7 +66,7 @@ def InitialChargeDensity(x, y, z):
      number_of_patches = kPatchPerGridDimension,
      EM_boundary_conditions = [ ["periodic"] ],
      print_every = 10,
-     random_seed = smilei_mpi_rank)
+     )
 
 Vectorization(mode = "off")
 

benchmarks/gpu/tst3d_gpu_o2_thermal_plasma_short.py

@@ -61,7 +61,6 @@
     gpu_computing = True,
 
     # random_seed = 0xDEADBEEF,
-    random_seed = smilei_mpi_rank,
 )
 
 Vectorization(

benchmarks/gpu/tst3d_v_o2_thermal_plasma_medium.py

@@ -67,7 +67,7 @@ def InitialChargeDensity(x, y, z):
      number_of_patches = kPatchPerGridDimension,
      EM_boundary_conditions = [ ["periodic"] ],
      print_every = 10,
-     random_seed = smilei_mpi_rank)
+     )
 
 Vectorization(mode = "on")
 

benchmarks/gpu/tst3d_v_o2_thermal_plasma_short.py

@@ -61,7 +61,6 @@
     gpu_computing = False,
 
     # random_seed = 0xDEADBEEF,
-    random_seed = smilei_mpi_rank,
 )
 
 Vectorization(

benchmarks/tst2d_18_em_pml.py

@@ -25,7 +25,6 @@
                               ['PML','PML'],
                              ],
     number_of_pml_cells             = [[10,10],[10,10]],
-    random_seed = smilei_mpi_rank
 )
 
 Antenna(

benchmarks/tst2d_s_o4_laser_wake_vay.py

@@ -28,8 +28,6 @@
 
     solve_poisson = False,
     print_every = 100,
-
-    random_seed = smilei_mpi_rank
 )
 
 MovingWindow(

benchmarks/tst2d_s_o4_radiation_pressure_acc.py

@@ -28,8 +28,6 @@
         ['silver-muller'],
         ['periodic'],
     ],
-
-    random_seed = smilei_mpi_rank
 )
 
 Vectorization(

benchmarks/tst2d_tasks_01_radiation_pressure_acc.py

@@ -39,7 +39,6 @@
         ['periodic'],
     ],
     cluster_width = 16,
-    random_seed = smilei_mpi_rank
 )
 
 

benchmarks/tst2d_v_o2_em_propagation.py

@@ -32,8 +32,6 @@
     ],
 
     EM_boundary_conditions_k = [[cos(ang), sin(ang)],[-1.,0.],[0.,1.],[0.,-1.]],
-
-    random_seed = smilei_mpi_rank
 )
 
 Vectorization(

benchmarks/tst2d_v_o4_em_propagation.py

@@ -32,8 +32,6 @@
     ],
 
     EM_boundary_conditions_k = [[cos(ang), sin(ang)],[-1.,0.],[0.,1.],[0.,-1.]],
-
-    random_seed = smilei_mpi_rank
 )
 
 Vectorization(

benchmarks/tst2d_v_o4_laser_wake_vay.py

@@ -28,8 +28,6 @@
 
     solve_poisson = False,
     print_every = 100,
-
-    random_seed = smilei_mpi_rank
 )
 
 Vectorization(

benchmarks/tst2d_v_o4_multiphoton_Breit_Wheeler.py

@@ -95,9 +95,6 @@ def n0_positron(x,y):
     simulation_time = Tsim,
 
     EM_boundary_conditions = [field_cond, field_cond],
-
-    random_seed = smilei_mpi_rank,
-
     reference_angular_frequency_SI = wr
 )
 

benchmarks/tst2d_v_o4_radiation_pressure_acc.py

@@ -28,8 +28,6 @@
         ['silver-muller'],
         ['periodic'],
     ],
-
-    random_seed = smilei_mpi_rank
 )
 
 

benchmarks/tst3d_s_o4_em_propagation.py

@@ -21,8 +21,6 @@
     simulation_time = Tsim,
 
     EM_boundary_conditions = [ ['silver-muller'] ],
-
-    random_seed = smilei_mpi_rank
 )
 
 LaserGaussian3D(

doc/Sphinx/Overview/releases.rst

@@ -30,6 +30,12 @@ Changes made in the repository (not released)
 * Features:
 
   * Relativistic field initialization now supports multiple species and both direction propagations.
+  * Added the argument ``phase_offset`` in laser definitions such as ``LaserGaussian2D``.
+  * The ``LaserGaussianAM`` definition will only use one coordinate for its ``focus`` argument 
+    (the transverse coordinate of the focus in this geometry is zero).
+  * Small improvements in PML for envelope model (AM and 2D).
+  * Deprecated ``smilei_rand_max``.
+  * New namelist variables ``smilei_omp_threads`` and ``smilei_total_cores``.
 
 * Happi:
 
@@ -44,11 +50,6 @@ Changes made in the repository (not released)
 
   * Dark theme (click the switch on the bottom left, or set browser preferences).
 
-* Added the argument ``phase_offset`` in laser definitions such as ``LaserGaussian2D``.
-* The ``LaserGaussianAM`` definition will only use one coordinate for its ``focus`` argument 
-  (the transverse coordinate of the focus in this geometry is zero).
-* Small improvements in PML for envelope model (AM and 2D).
-
 * Bug fixes:
 
   * ``dump_minutes`` often failed to write some checkpoint files.

doc/Sphinx/Use/namelist.rst

@@ -60,7 +60,8 @@ for each MPI process). The following steps are executed:
 
    * The rank of the current MPI process as :py:data:`smilei_mpi_rank`.
    * The total number of MPI processes as :py:data:`smilei_mpi_size`.
-   * The maximum random integer as :py:data:`smilei_rand_max`.
+   * The number of OpenMP threads per MPI :py:data:`smilei_omp_threads`.
+   * The total number of cores :py:data:`smilei_total_cores`.
 
 #. The namelist(s) is executed.
 
@@ -3619,9 +3620,15 @@ namelist. They should not be re-defined by the user!
 
   The total number of MPI processes.
 
-..
-  <<Not showing this anymore because of new rand system>>
-  .. py:data:: smilei_rand_max
+.. py:data:: smilei_omp_threads
+
+  The number of OpenMP threads per MPI.
+
+.. py:data:: smilei_total_cores
 
-    The largest random integer.
+  The total number of cores.
 
+.. note::
+
+  These variables can be access during ``happi`` post-processing, e.g.
+  ``S.namelist.smilei_mpi_size``.

src/Params/Params.cpp

@@ -129,16 +129,20 @@ Params::Params( SmileiMPI *smpi, std::vector<std::string> namelistsFiles ) :
     PyObject_SetAttrString( Py_main, "_test_mode", Py_False );
     PyTools::checkPyError();
 
-    // here we add the rank, in case some script need it
+    // we add the rank, in case some script needs it
     PyModule_AddIntConstant( Py_main, "smilei_mpi_rank", smpi->getRank() );
 
-    // here we add the MPI size, in case some script need it
+    // we add the MPI size, in case some script needs it
     PyModule_AddIntConstant( Py_main, "smilei_mpi_size", smpi->getSize() );
     namelist += string( "smilei_mpi_size = " ) + to_string( smpi->getSize() ) + "\n";
 
-    // here we add the larget int, important to get a valid seed for randomization
-    PyModule_AddIntConstant( Py_main, "smilei_rand_max", RAND_MAX );
-    namelist += string( "smilei_rand_max = " ) + to_string( RAND_MAX ) + "\n\n";
+    // we add the openMP size, in case some script needs it
+    PyModule_AddIntConstant( Py_main, "smilei_omp_threads", smpi->getOMPMaxThreads() );
+    namelist += string( "smilei_omp_threads = " ) + to_string( smpi->getOMPMaxThreads() ) + "\n";
+
+    // we add the total number of cores, in case some script needs it
+    PyModule_AddIntConstant( Py_main, "smilei_total_cores", smpi->getGlobalNumCores() );
+    namelist += string( "smilei_total_cores = " ) + to_string( smpi->getGlobalNumCores() ) + "\n";
 
     // Running pyprofiles.py
     runScript( string( reinterpret_cast<const char *>( pyprofiles_py ), pyprofiles_py_len ), "pyprofiles.py", globals );

src/Python/pyinit.py

@@ -645,7 +645,8 @@ class MultiphotonBreitWheeler(SmileiComponent):
 # Smilei-defined
 smilei_mpi_rank = 0
 smilei_mpi_size = 1
-smilei_rand_max = 2**31-1
+smilei_omp_threads = 1
+smilei_total_cores = 1
 
 # Variable to set to False for the actual run (useful for the test mode)
 _test_mode = True