$$M=\sum a_i*x'_i*x_i^{T}$$
where $x_i$ is the $i^{th}$ column of the input matrix; $a_i$ is the i'th element of the weight vector. Note that $x$ and $x'$ might be different. It is for non-stationary kernel functions.
xixj(mat,mat.new=NULL,a=NULL)
cov.linear
,xixj_sta