library("CVXR");
A <- Variable(3,3);
D <- matrix(c(A[1,1],0,0,0,A[2,2],0,0,0,A[3,3]),nrow=3);
A-D
Здесь у нас есть матрица переменных, A (я не уверен, правильно ли это вызывать объект), и D как матрица переменных, в которой диагональ состоит из диагонали A, а все остальное - 0.00.
Я хотел бы иметь возможность устанавливать ограничения, такие как (A-D)%*%x == y,
как обычно, например, A%*%x == y будет работать нормально.
Однако R не похоже на A-D - "объект s4 не вычитается". Как это сделать?