Alternativa a diag (X '* C * X) en matlab
Considere un N x 1 vector x y una N x N matriz C. Me gustaría evaluar
s = x'*C*x;
en matlab para muchas muestras del vector x, p. considere M muestras de x como una matriz N x M X; esto se puede hacer usando
S = diag(X'*C*X);
pero esta es una solución pobre ya que se asigna una matriz M x M en el proceso, y esto se rompe para M> 1e5. ¿Existe alguna funcionalidad de matlab que pueda sugerir una alternativa?