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 acknowledge that you have read our updated terms of service, privacy policy and cookie policy, and that your continued use of the website is subject to these policies.

Browse other questions tagged or ask your own question.