Hi,
Thanks for the reply and interest.
I think what you are trying to say is that, in terms of above notation,
q(i,j) = p(i,j)*s(i,i)*s(j,j)
This is what I did in the simulation I mentioned in the post.
I randomly generated a symmetric positive definite matrix, and scaled (s(i,i)>=1...