  We consider a versatile matrix model of the form ${\bf A}+i {\bf B}$, where
${\bf A}$ and ${\bf B}$ are real random circulant matrices with independent
but, in general, non-identically distributed Gaussian entries. For this model,
we derive exact result for the joint probability density function and find that
it is a multivariate Gaussian. Consequently, exact expression for arbitrary
order marginal density also ensues. It is demonstrated that by adjusting the
averages and variances of the Gaussian elements of ${\bf A}$ and ${\bf B}$, we
can interpolate between a remarkably wide range of eigenvalue distributions in
the complex plane. In particular, we can examine the crossover between a random
real circulant matrix and a random complex circulant matrix. We also extend our
study to include Wigner-like and Wishart-like matrices constructed from our
general random circulant matrix. To validate our analytical findings, Monte
Carlo simulations are conducted, which confirm the accuracy of our results.
