===== Mathematica Jobs ===== 

==== Mathematica serial job ====

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

Script ''math-vrt.m'':

<code>
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}]
Answer = A + B

Quit[];
</code>

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

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

module load mathematica

math 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-bdw.m'':

<code>
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-bdw.sh'': 

<code bash>
#!/bin/sh
#SBATCH --job-name=MathBDW
#SBATCH --output=%x.o%j
#SBATCH --error=%x.e%j
#SBATCH --nodes=1
#SBATCH --tasks-per-node=8
#SBATCH --partition=bdw
#SBATCH --mem=16G
#SBATCH --time=0-00:30:00
#SBATCH --account=<accountname>

module load mathematica

math math-bdw.m
</code>

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

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