Updated parallel Mathematica examples

Author: pkrastev

Parallel_Computing/Mathematica/README.md

@@ -4,25 +4,57 @@ Program computes PI with parallel Monte-Carlo method.
 
 ### Contents:
 
-* parallel_pi.m: Mathematica source code
-* run.sbatch: Btach-job submission script to send the job to the queue.
+* `parallel_pi.m`: Mathematica source code
+* `run.sbatch`: Btach-job submission script to send the job to the queue.
 
-#### Example Batch-Job Submission Script:
+### Source code:
+
+```mathematica
+(* Compute PI in via parallel Monte-Carlo method *)
+tstart = AbsoluteTime[];
+Print [ "Parallel calculation of PI via Monte-Carlo method." ];
+
+nproc = 8;
+LaunchKernels[nproc];
+
+Print[ " Number of kernels: ", nproc];
+
+n = 10^8;
+m = n / nproc;
+Print [ " Total number of hits: ", n ];
+Print [ " Number of hits per core: ", m ];
+
+acceptpoint[j_] := Total[Table[ 1 - Floor[ Sqrt[ (Random[])^2 + (Random[])^2 ] ], {i,1,m} ] ];
+DistributeDefinitions[n,m,acceptpoint];
+t1 = AbsoluteTime[];
+hits = ParallelTable[acceptpoint[j], {j,1,nproc} ];
+t2 = AbsoluteTime[];
+tt = t2 - t1;
+hits = Total[ hits ];
+pi = hits / ( nproc * m ) * 4.0;
+Print [ " Computed PI = ", pi ];
+Print [ " Time in parallel calculation: ", tt ] ;
+tend = AbsoluteTime[];
+ttotal = tend - tstart;
+Print [ " Total time: ", ttotal  ];
+Quit [ ]
+```
+
+### Example Batch-Job Submission Script:
 
 ```bash
 #!/bin/bash
 #SBATCH -J parallel_pi
 #SBATCH -o parallel_pi.out
 #SBATCH -e parallel_pi.err
-#SBATCH -p shared
+#SBATCH -p test
 #SBATCH -N 1
 #SBATCH -c 8
 #SBATCH -t 0-00:30
 #SBATCH --mem=4000
 
 # Load required software modules
-source new-modules.sh
-module load mathematica/11.1.1-fasrc01
+module load mathematica/13.3.0-fasrc01
 
 # Run program
 srun -n 1 -c 8 math -script parallel_pi.m
@@ -31,8 +63,6 @@ srun -n 1 -c 8 math -script parallel_pi.m
 ### Example Usage:
 
 ```bash
-source new-modules.sh
-module load mathematica/11.1.1-fasrc01
 sbatch run.sbatch
 ```
 

Parallel_Computing/Mathematica/run.sbtach

@@ -2,15 +2,14 @@
 #SBATCH -J parallel_pi
 #SBATCH -o parallel_pi.out
 #SBATCH -e parallel_pi.err
-#SBATCH -p shared
+#SBATCH -p test
 #SBATCH -N 1
 #SBATCH -c 8
 #SBATCH -t 0-00:30
 #SBATCH --mem=4000
 
 # Load required software modules
-module load centos6/0.0.1-fasrc01
-module load mathematica/11.1.1-fasrc01
+module load mathematica/13.3.0-fasrc01
 
 # Run program
 srun -n 1 -c 8 math -script parallel_pi.m