View source for Mathematica Parallel Computing Configuration

This page contains information regarding how to setup nodes in the cluster as Remote Kernels in different methods to run parallel computation

===First Method : Launch the kernels using commands in the notebook===
Use the code given below in a Mathematica notebook. Change "kernels" variable to the number of kernels wanted for parallel computation
 Needs["SubKernels`RemoteKernels`"]
 kernels = 4; (*Number of kernels on each machine*)
 configureKernel[host_] := 
  SubKernels`RemoteKernels`RemoteMachine[host, kernels];
 LaunchKernels[
  Map[configureKernel, {"<nowiki>ssh://ramsey/</nowiki>", "<nowiki>ssh://fibonacci/</nowiki>", 
  "<nowiki>ssh://hopper/</nowiki>", "<nowiki>ssh://boole/</nowiki>", "<nowiki>ssh://heaviside/</nowiki>", 
  "<nowiki>ssh://squid1/</nowiki>", "<nowiki>ssh://squid2/</nowiki>"}]]; (*all the available machines, change to whatever is needed*)


Printing "Hello World" in parallel after setting up 4 kernels per machine with above code 

[[File:Image1.png|frameless|663x663px]]

[[File:Image.png|frameless|667x667px]]


===Second Method: Configure Remote Kernels in Mathematica Menus===
In this method, we will set up the nodes in the cluster as Remote Kernels from the Mathematica Parallelize menu. This setup is required only once, every subsequent time, simply running LaunchKernels

command should set up the remote kernels for parallel evaluation. 

====Step 1====
Open Mathematica and open Parallel Kernel Configuration from the Evaluation tab

[[File:Mathematica parallel.png|frameless|660x660px]]

====Step 2====
Under Parallel, choose Remote Kernels, click on add hosts and add as shown below. Tick enable next to the number of kernels. For each, you can configure the number of kernels you want running on the machine.

[[File:Image2.png|frameless|669x669px]]


Add all the machines and number of kernels you want to add. This is a one time setup, you won't need to do this again and again.


[[File:Image3.png|frameless|406x406px]]

====Step 3====
Disable RemoteKernel Object here (if it is enabled)

[[File:Image4.png|frameless|699x699px]]


If you don't disable this, subkernels will be created on the local kernel and not on the remote hosts.

====Step 4====
Once all these settings are done, running LaunchKernels[] from any notebook will launch the specification in the menu.

In this example 3 machines were configured in menu as shown above:

[[File:Image6.png|frameless|732x732px]]

[[File:Image5.png|frameless|732x732px]]


===Third Method: Launching the kernels in a loop===
For this method, the kernels will be launched by running the below code:

 Needs["SubKernels`
 
 RemoteKernels`"]
 
 kernels = 1;(*number of kernels on each machine*)
 
 configureKernel[host_] := 
 
 SubKernels`RemoteKernels`RemoteMachine[host, kernels];
 
 For[i = 0, i < 15, i++; 
 
 If[i < 10, 
 
  LaunchKernels[
 
   Map[configureKernel, {StringJoin[{"<nowiki>ssh://pnode0</nowiki>", ToString[i], 
 
       "/"}]}]], 
 
  LaunchKernels[
 
    Map[configureKernel, {StringJoin[{"<nowiki>ssh://pnode</nowiki>", ToString[i], 
 
       "/"}]}]]]]
We are launching 15 kernels here. Change the number in the loop condition to launch that number of kernels.


===Parallel Evaluation Example: Generating and Plotting Mandelbrot Set===


In this section we will take a look at how the time for Mandelbrot set generation is greatly improved when we run on parallel remote kernels.

Code for generating Mandelbrot set without parallel kernels:

[[File:Image7.png|frameless|527x527px]] 


Output and Timing:

[[File:Image8.png|frameless|547x547px]]


To parallelize this, we will launch 15 cores using the third method:

[[File:Image9.png|frameless|557x557px]]


Output and Timing for Parallel Evaluation:

[[File:Image10.png|frameless|566x566px]]



Code for running parallel evaluation of Mandelbrot set:

 f[c_, maxiter_Integer : 10^4] := 
 
 Abs[NestWhile[#^2 + c &, 0, Abs[#] < 2 &, 1, maxiter]] < 2
 
 f[0]
 
 (*True*)
 
 f[1]
 
 (*False*)
 
 f[I]
 
 (*True*)
 
 f[I, 10^6]
 
 (*True*)
 
 ArrayPlot[
 
 ParallelTable[
 
  Boole[f[N[x + I y], 10^3]], {y, -2, 2, 1/50}, {x, -2, 2, 1/50}]]