S1 = matrix(c(1, 0.9, 0.8, 0.5,
0.9, 1.1, 0.7, 0.4,
0.8, 0.7, 1.5, 0.2,
0.5, 0.4, 0.2, 1), ncol = 4)
S2 = t(S1)%*% S1
S2 = S2/2
eval_objective(PC = c(1,0,0,0,sqrt(2),0,0,-sqrt(2)),
eta = 1, gamma = 0.5,
COVS = list(S1, S2))
Run the code above in your browser using DataLab