I am in need to optimize the performance of the 3D solver that integrates a system of equations for seismic wave propagation. Unsurprisingly, the function that implements the finite-difference stencil computations, is the bottleneck of the computations.

I read this article, where in Section 3.2 the authors say that the register blocking allowed them to obtain substantial speedup of the 3D stencil computations.

I do not understand how to implement this technique. Is it possible to find somewhere the description of register blocking in clear form?

  • $\begingroup$ Is it any different from regular blocking, except instead of fitting a block into cache, you fit it into some SSE registers? $\endgroup$ – Kirill Dec 23 '18 at 14:56
  • $\begingroup$ This is what I try to understand. $\endgroup$ – Dmitry Kabanov Dec 26 '18 at 9:17

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.