 output vector of the routine. z could be y, but not y + k (i.e. an overlapping area over y, but starting from a base index different from y). </td></tr>
    <tr><td class="paramname">y</td><td>The y input vector </td></tr>
    <tr><td class="paramname">alpha</td><td>The alpha scalar </td></tr>
    <tr><td class="paramname">dM</td><td>The stacked HDIA non zero values allocation pointer </td></tr>
    <tr><td class="paramname">offsets</td><td>The stacked HDIA diagonals offsets vector </td></tr>
    <tr><td c