Я хочу разреженную матрицу только с диагональными блоками большой матрицы D
В настоящее время я делаю это следующим образом:
hadamard=bdiag(matrix(1,500,500),matrix(1,300,300),...))
D*hadamard
Проблема в том, что D очень большой иэто умножение элементов требует много времени.Есть ли лучший способ сделать это.