Non-Gaussian positive-definite matrix-valued random fields for elliptic stochastic partial differential operators