Dirac spin projection

Syntax	void QLA_H_eqop_spproj_D( restrict QLA_HalfFermion *r,
	QLA_DiracFermion *a, int d, int p)
Meaning	$r = (1 + p\gamma_d)a$
eqop	eq