===== Mathematica =====

==== Mathematica serial job ====

=== Execution of a Mathematica serial program ===

Script ''math-vrt.m'':

<code mathematica math-vrt.m>
Print[Integrate[1/(x^4 - a^4), x]];

(*Prints all Mersenne PRime numbers less than 10000*)
Print[Select[Range[10000],PrimeQ[2^#-1]&]];

A = Sum[i, {i,1,100}];
B = Mean[{25, 36, 22, 16, 8, 42}];
Print[A + B];

Quit[];
</code>

Script ''slurm-math-vrt.sh'':

<code bash slurm-math-vrt.sh>
#!/bin/bash
#SBATCH --job-name=MathVRT
#SBATCH --output=%x.o%j
#SBATCH --error=%x.e%j
#SBATCH --ntasks=1
#SBATCH --partition=vrt
#SBATCH --qos=vrt
#SBATCH --mem=8G
#SBATCH --time=0-00:30:00
#SBATCH --account=<account>

module load mathematica

math -script math-vrt.m
</code>

Submitting ''slurm-math-vrt.sh'':

<code>
sbatch slurm-math-vrt.sh
</code>

==== Mathematica parallel job ====

=== Execution of a parallel job with Mathematica ===

Script ''math-cpu.m'':

<code mathematica math-cpu.m>
GetEnvironment["SLURM_NTASKS_PER_NODE"]

(*Limits Mathematica to requested resources*)
Unprotect[$ProcessorCount];
$ProcessorCount = ToExpression[Environment["SLURM_NTASKS_PER_NODE"]];

(*Prints the machine name that each kernel is running on*)
Print[ParallelEvaluate[$MachineName]];

(*Prints all Mersenne PRime numbers less than 10000*)
Print[Parallelize[Select[Range[10000],PrimeQ[2^#-1]&]]];

Quit[];
</code>

Script ''slurm-math-cpu.sh'': 

<code bash slurm-math-cpu.sh>
#!/bin/bash
#SBATCH --job-name=MathCPU
#SBATCH --output=%x.o%j
#SBATCH --error=%x.e%j
#SBATCH --nodes=1
#SBATCH --ntasks-per-node=8
#SBATCH --partition=cpu
#SBATCH --qos=cpu
#SBATCH --mem=16G
#SBATCH --time=0-00:30:00
#SBATCH --account=<account>

module load mathematica

math -script math-cpu.m
</code>

Submitting ''slurm-math-cpu.sh'':

<code>
sbatch slurm-math-cpu.sh
</code>