Profiling¶
The Profile
module provides tools to help developers improve the
performance of their code. When used, it takes measurements on running
code, and produces output that helps you understand how much time is
spent on individual line(s). The most common usage is to identify
“bottlenecks” as targets for optimization.
Profile
implements what is known as a “sampling” or statistical
profiler. It
works by periodically taking a backtrace during the execution of any
task. Each backtrace captures the currently-running function and
line number, plus the complete chain of function calls that led to
this line, and hence is a “snapshot” of the current state of
execution.
If much of your run time is spent executing a particular line of code, this line will show up frequently in the set of all backtraces. In other words, the “cost” of a given line—or really, the cost of the sequence of function calls up to and including this line—is proportional to how often it appears in the set of all backtraces.
A sampling profiler does not provide complete line-by-line coverage, because the backtraces occur at intervals (by default, 1 ms on Unix systems and 10 ms on Windows, although the actual scheduling is subject to operating system load). Moreover, as discussed further below, because samples are collected at a sparse subset of all execution points, the data collected by a sampling profiler is subject to statistical noise.
Despite these limitations, sampling profilers have substantial strengths:
- You do not have to make any modifications to your code to take timing measurements (in contrast to the alternative instrumenting profiler).
- It can profile into Julia’s core code and even (optionally) into C and Fortran libraries.
- By running “infrequently” there is very little performance overhead; while profiling, your code can run at nearly native speed.
For these reasons, it’s recommended that you try using the built-in sampling profiler before considering any alternatives.
Basic usage¶
Let’s work with a simple test case:
function myfunc()A=rand(100,100,200)maximum(A)end
It’s a good idea to first run the code you intend to profile at least once (unless you want to profile Julia’s JIT-compiler):
julia>myfunc()# run once to force compilation
Now we’re ready to profile this function:
julia>@profilemyfunc()
To see the profiling results, there is a graphical browser available, but here we’ll use the text-based display that comes with the standard library:
julia>Profile.print()23client.jl;_start;line:37323client.jl;run_repl;line:16623client.jl;eval_user_input;line:9123profile.jl;anonymous;line:148none;myfunc;line:28dSFMT.jl;dsfmt_gv_fill_array_close_open!;line:12815none;myfunc;line:32reduce.jl;max;line:352reduce.jl;max;line:3611reduce.jl;max;line:37
Each line of this display represents a particular spot (line number) in the code. Indentation is used to indicate the nested sequence of function calls, with more-indented lines being deeper in the sequence of calls. In each line, the first “field” indicates the number of backtraces (samples) taken at this line or in any functions executed by this line. The second field is the file name, followed by a semicolon; the third is the function name followed by a semicolon, and the fourth is the line number. Note that the specific line numbers may change as Julia’s code changes; if you want to follow along, it’s best to run this example yourself.
In this example, we can see that the top level is client.jl
‘s
_start
function. This is the first Julia function that gets called
when you launch Julia. If you examine line 373 of client.jl
,
you’ll see that (at the time of this writing) it calls run_repl()
,
mentioned on the second line. This in turn calls eval_user_input()
.
These are the functions in client.jl
that interpret what you type
at the REPL, and since we’re working interactively these functions
were invoked when we entered @profilemyfunc()
. The next line
reflects actions taken in the @profile
macro.
The first line shows that 23 backtraces were taken at line 373 of
client.jl
, but it’s not that this line was “expensive” on its own:
the second line reveals that all 23 of these backtraces were actually
triggered inside its call to run_repl
, and so on. To find out
which operations are actually taking the time, we need to look deeper
in the call chain.
The first “important” line in this output is this one:
8none;myfunc;line:2
none
refers to the fact that we defined myfunc
in the REPL,
rather than putting it in a file; if we had used a file, this would
show the file name. Line 2 of myfunc()
contains the call to
rand
, and there were 8 (out of 23) backtraces that occurred at
this line. Below that, you can see a call to
dsfmt_gv_fill_array_close_open!()
inside dSFMT.jl
. You might be surprised not to see the
rand
function listed explicitly: that’s because rand
is inlined,
and hence doesn’t appear in the backtraces.
A little further down, you see:
15none;myfunc;line:3
Line 3 of myfunc
contains the call to max
, and there were 15
(out of 23) backtraces taken here. Below that, you can see the
specific places in base/reduce.jl
that carry out the
time-consuming operations in the max
function for this type of
input data.
Overall, we can tentatively conclude that finding the maximum element is approximately twice as expensive as generating the random numbers. We could increase our confidence in this result by collecting more samples:
julia>@profile(fori=1:100;myfunc();end)julia>Profile.print()3121client.jl;_start;line:3733121client.jl;run_repl;line:1663121client.jl;eval_user_input;line:913121profile.jl;anonymous;line:1848none;myfunc;line:2842dSFMT.jl;dsfmt_gv_fill_array_close_open!;line:1281510none;myfunc;line:374reduce.jl;max;line:35122reduce.jl;max;line:361314reduce.jl;max;line:37
In general, if you have N
samples collected at a line, you can
expect an uncertainty on the order of sqrt(N)
(barring other
sources of noise, like how busy the computer is with other tasks). The
major exception to this rule is garbage collection, which runs
infrequently but tends to be quite expensive. (Since Julia’s garbage
collector is written in C, such events can be detected using the
C=true
output mode described below, or by using ProfileView.jl.)
This illustrates the default “tree” dump; an alternative is the “flat” dump, which accumulates counts independent of their nesting:
julia>Profile.print(format=:flat)CountFileFunctionLine3121client.jl_start3733121client.jleval_user_input913121client.jlrun_repl166842dSFMT.jldsfmt_gv_fill_array_close_open!128848nonemyfunc21510nonemyfunc33121profile.jlanonymous174reduce.jlmax35122reduce.jlmax361314reduce.jlmax37
If your code has recursion, one potentially-confusing point is that a line in a “child” function can accumulate more counts than there are total backtraces. Consider the following function definitions:
dumbsum(n::Integer)=n==1?1:1+dumbsum(n-1)dumbsum3()=dumbsum(3)
If you were to profile dumbsum3
, and a backtrace was taken while it was executing dumbsum(1)
, the backtrace would look like this:
dumbsum3dumbsum(3)dumbsum(2)dumbsum(1)
Consequently, this child function gets 3 counts, even though the parent only gets one. The “tree” representation makes this much clearer, and for this reason (among others) is probably the most useful way to view the results.
Accumulation and clearing¶
Results from @profile
accumulate in a buffer; if you run multiple
pieces of code under @profile
, then Profile.print()
will show
you the combined results. This can be very useful, but sometimes you
want to start fresh; you can do so with Profile.clear()
.
Options for controlling the display of profile results¶
Profile.print()
has more options than we’ve described so far.
Let’s see the full declaration:
function print(io::IO=STDOUT,data=fetch();format=:tree,C=false,combine=true,cols=tty_cols(),maxdepth=typemax(Int),sortedby=:filefuncline)
Let’s discuss these arguments in order:
The first argument allows you to save the results to a file, but the default is to print to
STDOUT
(the console).The second argument contains the data you want to analyze; by default that is obtained from
Profile.fetch()
, which pulls out the backtraces from a pre-allocated buffer. For example, if you want to profile the profiler, you could say:data=copy(Profile.fetch())Profile.clear()@profileProfile.print(STDOUT,data)# Prints the previous resultsProfile.print()# Prints results from Profile.print()
The first keyword argument,
format
, was introduced above. The possible choices are:tree
and:flat
.C
, if set totrue
, allows you to see even the calls to C code. Try running the introductory example withProfile.print(C=true)
. This can be extremely helpful in deciding whether it’s Julia code or C code that is causing a bottleneck; settingC=true
also improves the interpretability of the nesting, at the cost of longer profile dumps.Some lines of code contain multiple operations; for example,
s+=A[i]
contains both an array reference (A[i]
) and a sum operation. These correspond to different lines in the generated machine code, and hence there may be two or more different addresses captured during backtraces on this line.combine=true
lumps them together, and is probably what you typically want, but you can generate an output separately for each unique instruction pointer withcombine=false
.cols
allows you to control the number of columns that you are willing to use for display. When the text would be wider than the display, you might see output like this:33inference.jl;abstract_call;line:64533inference.jl;abstract_call;line:64533...rence.jl;abstract_call_gf;line:56733...nce.jl;typeinf;line:1201+15...nce.jl;...t_interpret;line:900+35...ence.jl;abstract_eval;line:758+45...ence.jl;...ct_eval_call;line:733+65...ence.jl;abstract_call;line:645
File/function names are sometimes truncated (with
...
), and indentation is truncated with a+n
at the beginning, wheren
is the number of extra spaces that would have been inserted, had there been room. If you want a complete profile of deeply-nested code, often a good idea is to save to a file and use a very widecols
setting:s=open("/tmp/prof.txt","w")Profile.print(s,cols=500)close(s)
maxdepth
can be used to limit the size of the output in:tree
format (it nests only up to levelmaxdepth
)sortedby=:count
sorts the:flat
format in order of increasing counts
Configuration¶
@profile
just accumulates backtraces, and the analysis happens
when you call Profile.print()
. For a long-running computation,
it’s entirely possible that the pre-allocated buffer for storing
backtraces will be filled. If that happens, the backtraces stop but
your computation continues. As a consequence, you may miss some
important profiling data (you will get a warning when that happens).
You can obtain and configure the relevant parameters this way:
Profile.init()# returns the current settingsProfile.init(n,delay)Profile.init(delay=0.01)
n
is the total number of instruction pointers you can store, with
a default value of 10^6
. If your typical backtrace is 20
instruction pointers, then you can collect 50000 backtraces, which
suggests a statistical uncertainty of less than 1%. This may be good
enough for most applications.
Consequently, you are more likely to need to modify delay
,
expressed in seconds, which sets the amount of time that Julia gets
between snapshots to perform the requested computations. A very
long-running job might not need frequent backtraces. The default
setting is delay=0.001
. Of course, you can decrease the delay
as well as increase it; however, the overhead of profiling grows once
the delay becomes similar to the amount of time needed to take a
backtrace (~30 microseconds on the author’s laptop).
Memory allocation analysis¶
One of the most common techniques to improve performance is to reduce
memory allocation. The total amount of allocation can be measured
with @time
and @allocated
, and specific lines triggering
allocation can often be inferred from profiling via the cost of
garbage collection that these lines incur. However, sometimes it is
more efficient to directly measure the amount of memory allocated by
each line of code.
To measure allocation line-by-line, start Julia with the
--track-allocation=<setting>
command-line option, for which you
can choose none
(the default, do not measure allocation), user
(measure memory allocation everywhere except Julia’s core code), or
all
(measure memory allocation at each line of Julia code).
Allocation gets measured for each line of compiled code. When you quit
Julia, the cumulative results are written to text files with .mem
appended after the file name, residing in the same directory as the
source file. Each line lists the total number of bytes allocated.
The Coverage
package contains some elementary analysis tools, for
example to sort the lines in order of number of bytes allocated.
In interpreting the results, there are a few important details. Under
the user
setting, the first line of any function directly called
from the REPL will exhibit allocation due to events that happen in the
REPL code itself. More significantly, JIT-compilation also adds to
allocation counts, because much of Julia’s compiler is written in
Julia (and compilation usually requires memory allocation). The
recommended procedure is to force compilation by executing all the
commands you want to analyze, then call Profile.clear_malloc_data()
to
reset all allocation counters. Finally, execute the desired commands
and quit Julia to trigger the generation of the .mem
files.