Пожалуйста, предложите эффективный и надежный метод (объем памяти, точность вычислений, скорость вычислений) для следующей матричной процедуры:
A + [ b * transpose ( b )] + transpose ([ b * transpose ( b *) 1012 *)])
где:
A является m -by- m плотной ранг-симметричной (ковариационной) матрицей
b - это m -by-1 vector
- обозначает умножение матрицы
1000 <<em> м <50000 </p>
Реализация будет в R