I need to create random vectors of real numbers a_i satisfying the following constraints:
abs(a_i) < c_i;
sum(a_i)< A; # sum of elements smaller than A
sum(b_i * a_i) < B; # weighted sum is smaller than B
aT*A*a < D # quadratic multiplication with A smaller than D
where c_i, b_i, A, B, D are constants.
What would be the typical algorithm to generate efficiently this kind of vector?