Quantitative DCE-MRI analysis requires fitting a model to the acquired enhancing time-course data. These models involve a convolution operation, which has a variety of possible implementations. The goal of this work was to evaluate three different implementations of the convolution operation: (i) Fourier, (ii) Summation, (iii) Iterative. The accuracy and execution time of each implementation was compared on a virtual phantom. It was found that the iterative technique was the fastest while also having the best overall accuracy.