Inserting a Dirac column vector into a Dirac propagator matrix

Syntax	void QLA_P_eqop_diracvec_D( restrict QLA_DiracPropagator *r,
	QLA_DiracFermion *a, int jc, int js)
Meaning	$r_{i_c,i_s,j_c,j_s} = a_{i_c,i_s} \ \ \ \mbox{for $i_c = 0\ldots{}n_c-1, i_s = 0\ldots{}3$}$
eqop	eq