View source for Third year CMP compulsory experiment/Accelerating the code

'''<big><span style="color:blue; ">This is the fourth section of the third year CMP experiment. You can return to the previous page, [[Third_year_CMP_compulsory_experiment/Introduction to Monte Carlo simulation|Introduction to Monte Carlo simulation]], or jump ahead to the next section, [[Third year CMP compulsory experiment/The effect of temperature|The effect of temperature]].</span></big>'''

So far, your energy() and magnetisation() functions probably use Python loops to go over every spin in your system (and a double loop for the energy). One of the drawbacks of Python's flexibility is a performance penalty &mdash; though developing code in Python tends to be faster than in C++, the resulting program is much, much slower.

Luckily, the NumPy module provides an interface to a number of functions written to give high performance. In this section, we are going to use these functions to improve the performance of our simulation.

'''<big>TASK 4a</big>: Use the script ILtimetrial.py to record how long your ''current'' version of IsingLattice.py takes to perform 2000 Monte Carlo steps. This will vary, depending on what else the computer happens to be doing, so perform repeats and report the error in your average!''' 

'''<big>TASK 4b</big>: Look at the documentation for the [http://docs.scipy.org/doc/numpy/reference/generated/numpy.sum.html NumPy sum function]. You should be able to modify your magnetisation() function so that it uses this to evaluate M. The energy is a little trickier. Familiarise yourself with the NumPy [http://docs.scipy.org/doc/numpy/reference/generated/numpy.roll.html roll] and [http://docs.scipy.org/doc/numpy/reference/generated/numpy.multiply.html multiply] functions, and use these to replace your energy double loop (you will need to call roll and multiply twice!).'''

'''<big>TASK 4c</big>: Use the script ILtimetrial.py to record how long your ''new'' version of IsingLattice.py takes to perform 2000 Monte Carlo steps. This will vary, depending on what else the computer happens to be doing, so perform repeats and report the error in your average!''' 

'''<big><span style="color:blue; ">This is the fourth section of the third year CMP experiment. You can return to the previous page, [[Third_year_CMP_compulsory_experiment/Introduction to Monte Carlo simulation|Introduction to Monte Carlo simulation]], or jump ahead to the next section, [[Third year CMP compulsory experiment/The effect of temperature|The effect of temperature]].</span></big>'''