Есть ли в Octave способ вычислить и сохранить только диагональ матричного произведения?

В основном, как делать:vector = diag(A*B);

Меня не волнует ни одно из значенийA*B кроме тех, что по диагонали. Размеры матрицы около80k x 12 а также12 x 80k, так что даже если бы я не заботился о скорости / дополнительной памяти, он просто не помещался в RAM.

Странно, поскольку Octave представляет собой пакет для огромных наборов данных, и диагонали очень важны, поэтому это должно быть возможно.