Bugfix: enforce SPDness of cov matrices.

Conflicts:
	covcorr.cpp
Refactor to hoist logp_id computation where possible.

(speeds up all but geodesic and my distances.)

Conflicts:
	model.cpp
	train.cpp
	traintest.cpp
	warco.cpp
	warco.hpp