Имеется ли матричная операция, которая при поэлементном делении матрицы на матрицу не будет трогать главную диагональ.
В этом примере в делителе нули по диагонали дают Inf, а нужно чтобы было ноль:

octave:32> a = [1 2 3; 4 5 6; 7 8 9]
a =

   1   2   3
   4   5   6
   7   8   9

octave:33> r = [0 1 2 ; 1 0 2; 1 2 0]
r =

   0   1   2
   1   0   2
   1   2   0

octave:34> a./r
ans =

      Inf   2.0000   1.5000
   4.0000      Inf   3.0000
   7.0000   4.0000      Inf

P.S.1 Стараюсь обойти циклы.

P.S.2 Может в Питоне что-то есть?