2020 Running batch jobs at TACC
Reservations
Use our summer school reservation (intro_NGS) when submitting batch jobs to get higher priority on the Lonestar5 normal queue during this course:
sbatch --reservation=intro_NGS <batch_file>.slurm
Note that a reservation is different from the TACC allocation/project for this class, which is UT-2015-05-18.
- 1 Compute cluster overview
- 2 Software at TACC
- 3 Job Execution
- 3.1 SLURM at a glance
- 3.2 Simple example
- 3.3 Job parameters
- 3.3.1 launcher_creator.py
- 3.3.2 job name and commands file
- 3.3.3 queues and runtime
- 3.3.4 allocation and SUs
- 3.3.4.1 ALLOCATION setting in .bashrc
- 3.3.5 wayness (tasks per node)
- 3.4 Wayness example
- 4 Some best practices
- 5 Interactive sessions (idev)
Compute cluster overview
When you SSH into ls5, your session is assigned to one of a small set of login nodes (also called head nodes). These are not the compute nodes that will run your jobs.
Think of a node as a computer, like your laptop, but probably with more cores and memory. Now multiply that computer a thousand or more, and you have a cluster.
The small set of login nodes are a shared resource (type the users command to see everyone currently logged in) and are not meant for running interactive programs – for that you submit a description of what you want done to a batch system, which farms the work out to one or more compute nodes.
On the other hand, the login nodes are intended for copying files to and from TACC, so they have a lot of network bandwidth while compute nodes have limited network bandwidth.
So follow these guidelines:
Do not perform substantial computation on the login nodes.
They are closely monitored, and you will get warnings from the TACC admin folks!
Code is usually developed and tested somewhere other than TACC, and only moved over when pretty solid.
Do not perform significant network access from your batch jobs.
Instead, stage your data onto $SCRATCH from a login node before submitting your job.
Lonestar5 and Stampede2 overview and comparison
Here is a comparison of the configurations and ls5 and stampede2. As you can see, stampede2 is the larger cluster, launched in 2017.
| ls5 | stampede2 |
|---|---|---|
login nodes | 6 20 cores each | 6 28 cores each |
standard compute nodes | 1,252 24 cores per node (48 virtual) | 4,200 KNL (Knights Landing)
1,736 SKX (Skylake)
|
large memory nodes | 10 total 2 w/1 TB memory, 48 cores | -- |
batch system | SLURM | SLURM |
maximum job run time | 48 hours | 96 hours on KNL nodes 48 hours on SKX nodes |
Note the use of the term virtual core above. Compute cores are standalone processors – mini CPUs, each of which can execute separate sets of instructions. However modern cores may also have hyper-threading enabled, where a single core can appear as more than one virtual processor to the Operating system (see https://en.wikipedia.org/wiki/Hyper-threading for more on hyper-threading). For example, Lonestar5 nodes have 2 hyperthreads (HTs) per core. So with 2 HTs for each of the 24 physical cores, each node has a total of 48 virtual cores.
User guides for ls5 and stampede2 can be found at:
Unfortunately, the TACC user guides are aimed towards a different user community – the weather modelers and aerodynamic flow simulators who need very fast matrix manipulation and other high performance computing (HPC) features. The usage patterns for bioinformatics – generally running 3rd party tools on many different datasets – is rather a special case for HPC. TACC calls our type of processing "parameter sweep jobs" and has a special process for running them, using ther launcher module.
Software at TACC
Programs and your $PATH
When you type in the name of an arbitrary program (ls for example), how does the shell know where to find that program? The answer is your $PATH. $PATH is a pre-defined environment variable whose value is a list of directories.The shell looks for program names in that list, in the order the directories appear.
To determine where the shell will find a particular program, use the which command:
Using which to search $PATH
which rsync
which catThe module system
The module system is an incredibly powerful way to have literally thousands of software packages available, some of which are incompatible with each other, without causing complete havoc. The TACC staff builds the desired package from source code in well-known locations that are NOT on your $PATH. Then, when a module is loaded, its binaries are added to your $PATH.
For example, the following module load command makes the bwa aligner available to you:
How module load affects $PATH
# first type "bwa" to show that it is not present in your environment:
bwa
# it's not on your $PATH either:
which bwa
# now add bwa to your environment and try again:
module load bwa
bwa
# and see how it's now on your $PATH:
which bwa
# you can see the new directory at the front of $PATH
echo $PATH
# to remove it, use "unload"
module unload bwa
bwa
# gone from $PATH again...
which bwamodule spider
These days the TACC module system includes hundreds of useful bioinformatics programs. To see if your favorite software package has been installed at TACC, use module spider:
module spider fastqc
module spider samtools
module spider bedtools
TACC BioContainers modules
It is quite a large systems administration task to install software at TACC and configure it for the module system. As a result, TACC was always behind in making important bioinformatics software availble. To address this problem, TACC moved to providing bioinformatics software via containers, which are virtual machines like VMware and Virtual Box, but are lighter weight: they require less disk space because they rely more on the host's base Linux environment. Specifically, TACC (and many other High Performance Computing clusters) use Singularity containers, which are similar to Docker containers but are more suited to the HPC environment (in fact one can build a Docker container then easily convert it to Singularity for use at TACC).
TACC obtains its containers from BioContainers (https://biocontainers.pro/ and https://github.com/BioContainers/containers), a large public repository of bioinformatics tool Singularity containers. This has allowed TACC to easily provision thousands of such tools!
These BioContainers are not visible in TACC's "standard" module system, but only after the master biocontainers module is loaded:
# Make sure the non-biocontainers version of bwa is not loaded
module unload bwa
# Verify that bwa is not available
bwa
# Verify that these programs are not in the standard TACC module system
module spider kallisto
module spider bowtie2
module spider minimap2
module spider multiqc
module spider GATK
module spider velvet
# Load the Biocontainers master module (this takes a while)
module load biocontainers
# Now look for those programs
module spider kallisto
module spider bowtie2
module spider minimap2
module spider multiqc
module spider GATK
module spider velvetNotice how the BioContainers module names have "ctr" in their names, version numbers, and other identifying information.
loading a biocontainer module
Once the biocontainers module has been loaded, you can just load the desired tool module, as with the kallisto pseudo-aligner program below.
# Load the Biocontainers master module
module load biocontainers
# Load the default kallisto biocontainer
module load kallisto
# Verify you can now execute kallisto
kallistoNote that loading a BioContainer does not add anything to your $PATH. Instead, it defines an alias, which is just a shortcut for executing the command. You can see the alias definition using the type command. And you can ensure the program is available using the command -v utility.
# Note that kallisto has not been added to your $PATH
which kallisto
# Instead, an alias has been defined. Use type to see its definition
type kallisto
# Ensure kallisto is available with command -v
command -v kallistoinstalling custom software
Even with all the tools available at TACC, inevitably you'll need something they don't have. In this case you can build the tool yourself and install it in a local TACC directory. While building 3rd party tools is beyond the scope of this course, it's really not that hard. The trick is keeping it all organized.
For one thing, remember that your $HOME directory quota is fairly small (10 GB on ls5), and that can fill up quickly if you install many programs. We recommend creating an installation area in your $WORK directory and installing programs there. You can then make symbolic links to the binaries you need in your $HOME/local/bin directory (which was added to your $PATH in your .bashrc).
See how we used a similar trick to make the launcher_creator.py program available to you. Using the ls -l option shows you where symbolic links point to:
Real location of launcher_creator.py
ls -l ~/local/bin$PATH caveat
Remember that the order of locations in the $PATH environment variable is the order in which the locations will be searched. In particular, the module load command adds to the front of your path. This can mask similarly-named programs, for example, in your $HOME/local/bin directory.
Job Execution
Job execution is controlled by the SLURM batch system on both ls5 and stampede2.
To run a job you prepare 2 files:
a commands file file containing the commands to run, one command per line (<job_name>.cmds)
a job control file that describes how to run the job (<job_name>.slurm)
The process of running the job involves these steps:
Create a commands file containing exactly one command per line.
Prepare a job control file for the commands file that describes how the job should be run.
You submit the job control file to the batch system. The job is then said to be queued to run.
The batch system prioritizes the job based on the number of compute nodes needed and the job run time requested.
When compute nodes become available, the job tasks (command lines in the <job_name>.cmds file) are assigned to one or more compute nodes and begin to run in parallel.
The job completes when either:
you cancel the job manually
all tasks in the job complete (successfully or not!)
the requested job run time has expired
SLURM at a glance
Here are the main components of the SLURM batch system.
| ls5, stampede2 |
|---|---|
batch system | SLURM |
batch control file name | < |
job submission command |
|
job monitoring command |
|
job stop command |
|
Simple example
Let's go through a simple example. Execute the following commands to copy a pre-made simple.cmds commands file:
Copy simple commands
mkdir -p $SCRATCH/core_ngs/slurm/simple
cd $SCRATCH/core_ngs/slurm/simple
cp $CORENGS/tacc/simple.cmds .What are the tasks we want to do? Each task corresponds to one line in the simple.cmds file, so let's take a look at it using the cat (concatenate) command that simply reads a file and writes each line of content to standard output (here, your Terminal):
View simple commands
cat simple.cmdsThe tasks we want to perform look like this:
echo "Command 1 on `hostname` - `date`" > cmd1.log 2>&1
echo "Command 2 on `hostname` - `date`" > cmd2.log 2>&1
echo "Command 3 on `hostname` - `date`" > cmd3.log 2>&1
echo "Command 4 on `hostname` - `date`" > cmd4.log 2>&1
echo "Command 5 on `hostname` - `date`" > cmd5.log 2>&1
echo "Command 6 on `hostname` - `date`" > cmd6.log 2>&1
echo "Command 7 on `hostname` - `date`" > cmd7.log 2>&1
echo "Command 8 on `hostname` - `date`" > cmd8.log 2>&1There are 8 tasks. Each is a simple echo command that just outputs string containing the task number and date to a different file.
Use the handy launcher_creator.py program to create the job submission script.
Create batch submission script for simple commands
launcher_creator.py -j simple.cmds -n simple -t 00:00:05 -w 8 -a UT-2015-05-18 -q developmentYou should see output something like the following, and you should see a simple.slurm batch submission file in the current directory.
Project simple.
Using job file simple.cmds.
Using development queue.
For 00:00:05 time.
Using UT-2015-05-18 allocation.
Not sending start/stop email.
Launcher successfully created. Type "sbatch simple.slurm" to queue your job.Submit your batch job like this, then check the batch queue to see the job's status.
Submit simple job to batch queue
sbatch --reservation=intro_NGS simple.slurm
showq -uIf you're quick, you'll see a queue status something like this:
SUMMARY OF JOBS FOR USER: <abattenh>
ACTIVE JOBS--------------------
JOBID JOBNAME USERNAME STATE NODES REMAINING STARTTIME
================================================================================
2916562 simple abattenh Running 1 0:00:57 Sat Jun 13 13:52:32
WAITING JOBS------------------------
JOBID JOBNAME USERNAME STATE NODES WCLIMIT QUEUETIME
================================================================================
Total Jobs: 1 Active Jobs: 1 Idle Jobs: 0 Blocked Jobs: 0If you don't see your simple job in either the ACTIVE or WAITING sections of your queue, it probably already finished – it should only run for a second or two!
Notice in my queue status, where the STATE is Running, there is only one node assigned. Why is this, since there were 8 tasks?
Every job, no matter how few tasks requested, will be assigned at least one node. Each ls5 node has 24 physical cores, so each of the 8 tasks can be assigned to a core.
Exercise: What files were created by your job?
filename wildcarding
Here's a cute trick for viewing the contents all your output files at once, using the cat command and filename wildcarding.
Multi-character filename wildcarding
cat cmd*.logThe cat command actually takes a list of one or more files (if you're giving it files rather than standard input – more on this shortly) and outputs the concatenation of them to standard output. The asterisk ( * ) in cmd*.log is a multi-character wildcard that matches any filename starting with cmd then ending with .log. So it would match cmd_hello_world.log.
You can also specify single-character matches inside brackets ( [ ] ) in either of the ways below, this time using the ls command so you can better see what is matching:
Single character filename wildcarding
ls cmd[12345678].log
ls cmd[1-8].logThis technique is sometimes called filename globbing, and the pattern a glob. Don't ask me why – it's a Unix thing. Globbing – translating a glob pattern into a list of files – is one of the handy thing the bash shell does for you. (Read more about Wildcards and special filenames.)
Exercise: How would you list all files starting with simple?
Here's what my cat output looks like. Notice the times are all nearly the same because all the tasks ran in parallel. That's the power of cluster computing!
Command 1 on nid00022 - Sat Jun 13 13:52:52 CDT 2020
Command 2 on nid00022 - Sat Jun 13 13:52:47 CDT 2020
Command 3 on nid00022 - Sat Jun 13 13:52:58 CDT 2020
Command 4 on nid00022 - Sat Jun 13 13:52:53 CDT 2020
Command 5 on nid00022 - Sat Jun 13 13:52:41 CDT 2020
Command 6 on nid00022 - Sat Jun 13 13:52:49 CDT 2020
Command 7 on nid00022 - Sat Jun 13 13:52:40 CDT 2020
Command 8 on nid00022 - Sat Jun 13 13:52:48 CDT 2020echo
Lets take a closer look at a typical task in the simple.cmds file.
An echo command
echo "Command 3 `date`" > cmd3.log 2>&1The echo command is like a print statement in the bash shell. Echo takes its arguments and writes them to one line of standard output. While not always required, it is a good idea to put the output string in double quotes.
backtick evaluation
So what is this funny looking `date` bit doing? Well, date is just another Linux command (try just typing it in). Here we don't want the shell to put the string "date" in the output, we want it to execute the date command and put the result text into the output. The backquotes ( ` ` also called backticks) around the date command tell the shell we want that command executed and its output substituted into the string. (Read more about Quoting in the shell.)
Backtick evaluation
# These are equivalent:
date
echo `date`
# But different from this:
echo dateoutput redirection
There's still more to learn from one of our simple tasks, something called output redirection:
echo "Command 3 `date`" > cmd3.log 2>&1Normally echo writes its string to standard output. If you invoke echo in an interactive shell like Terminal, standard output is displayed to the Terminal window.
Usually we want to separate the outputs of all our commands. Why is this important? Suppose we run a job with 100 commands, each one a whole pipeline (alignment, for example). 88 finish fine but 12 do not. Just try figuring out which ones had the errors, and where the errors occurred, if all the normal output is in one intermingled file and all the error in another intermingled file!
So in the above example the first '>' says to redirect the standard output of the echo command to the cmd3.log file. The '2>&1' part says to redirect standard error to the same place. Technically, it says to redirect standard error (built-in Linux stream 2) to the same place as standard output (built-in Linux stream 1); and since standard output is going to cmd3.log, any standard error will go there also. (Read more about Standard I/O streams.)
So what happens when output is generated by tasks in a batch job? Well, you may have noticed the files with names like simple.e2916562 and simple.o2916562 were created by your job. simple.o2916562 contains all standard output and simple.o2916562 contains all standard error generated by your tasks that was not redirected elsewhere, as well as information relating to running your job and its tasks.
Job parameters
Now that we've executed a really simple job, let's take a look at some important job submission parameters. These correspond to arguments to the launcher_creator.py script.
A bit of background. Historically, TACC was set up to cater to researchers writing their own C or Fortran codes highly optimized to exploit parallelism (the HPC crowd). Much of TACC's documentation is aimed at this audience, which makes it difficult to pick out the important parts for us.
The kind of jobs we biologists generally run are relatively new to TACC. They even have a special name for them: "parametric sweeps", by which they mean the same program running on different data sets.
In fact there is a special software module required to run our jobs, called the launcher module. You don't need to worry about activating the launcher module – that's done by the <job_name>.slurm script created by launcher_creator.py like this:
module load launcherThe launcher module knows how to interpret various job parameters in the <job_name>.slurm batch SLURM submission script and use them to create your job and assign its tasks to compute nodes. Our launcher_creator.py program is a simple Python script that lets you specify job parameters and writes out a valid <job_name>.slurm submission script.
launcher_creator.py
If you call launcher_creator.py with no arguments it gives you its usage description:
launcher_creator.py usage
usage: launcher_creator.py [-h] -n NAME -t TIME_REQUEST [-j JOB_FILE]
[-b SHELL_COMMANDS] [-B SHELL_COMMANDS_FILE]
[-q QUEUE] [-a [ALLOCATION]] [-m MODULES]
[-M MODULES_FILE] [-w WAYNESS] [-N NUM_NODES]
[-e [EMAIL]] [-l LAUNCHER] [-s]
Create launchers for TACC clusters. Report problems to rt-other@ccbb.utexas.edu
optional arguments:
-h, --help show this help message and exit
Required:
-n NAME, --name NAME The name of your job.
-t TIME_REQUEST, --time TIME_REQUEST
The time you want to give to your job. Format:
hh:mm:ss
Commands:
You must use at least one of these options to submit your commands for
TACC.
-j JOB_FILE, --jobs JOB_FILE
The name of the job file containing your commands.
-b SHELL_COMMANDS, --bash SHELL_COMMANDS
A string of shell (Bash, zsh, etc) commands that are
executed before any parametric jobs are launched.
-B SHELL_COMMANDS_FILE, --bash_file SHELL_COMMANDS_FILE
A file containing shell (Bash, zsh, etc) commands that
are executed before any parametric jobs are launched.
Optional:
-q QUEUE, --queue QUEUE
The TACC allocation for job submission.
Default="development"
-a [ALLOCATION], -A [ALLOCATION], --allocation [ALLOCATION]
The TACC allocation for job submission. You can set a
default ALLOCATION environment variable.
-m MODULES, --modules MODULES
A list of module commands. The "launcher" module is
always automatically included. Example: -m "module
swap intel gcc; module load bedtools"
-M MODULES_FILE, --modules_file MODULES_FILE
A file containing module commands.
-w WAYNESS, --wayness WAYNESS
Wayness: the number of commands you want to give each
node. The default is the number of cores per node.
-N NUM_NODES, --num_nodes NUM_NODES
Number of nodes to request. You probably don't need
this option. Use wayness instead. You ONLY need it if
you want to run a job list that isn't defined at the
time you submit the launcher.
-e [EMAIL], --email [EMAIL]
Your email address if you want to receive an email
from Lonestar when your job starts and ends. Without
an argument, it will use a default EMAIL_ADDRESS
environment variable.
-l LAUNCHER, --launcher_name LAUNCHER
The name of the launcher script that will be created.
Default="<name>.slurm"
-s Echoes the launcher filename to stdout.Because it is a long help message, we may want to pipe the output to more, a pager that displays one screen of text at a time. Type the spacebar to advance to the next page, and Ctrl-c to exit from more.
Get usage information for launcher_creator.py
# Use spacebar to page forward; Ctrl-c to exit
launcher_creator.py -h | morejob name and commands file
Recall how the simple.slurm batch file was created:
Create batch submission script for simple commands
launcher_creator.py -j simple.cmds -n simple -t 00:00:05 -w 8 -a UT-2015-05-18 -q developmenThe name of your commands file is given with the -j simple.cmds argument.
Your desired job name is given with the -n <job_name> argument.
The <job_name> (here simple) is the job name you will see in your queue.
By default a corresponding <job_name>.slurm batch file is created for you.
It contains the name of the commands file that the batch system will execute.
queues and runtime
TACC resources are partitioned into queues: a named set of compute nodes with different characteristics. The major ones on ls5 are listed below. Generally you use development (-q development) when you are writing and testing your code, then normal once you're sure your commands will execute properly.
queue name | maximum runtime | purpose |
|---|---|---|
development | 2 hrs | development (short queue wait times) |
normal | 48 hrs | normal priority (queue waits can sometimes be long) |
largemem | 48 hrs | large memory jobs |
In launcher_creator.py, the queue is specified by the -q argument.