X.positiveDiagonal( Y, pvt );
PaceMatrix sol = (PaceMatrix) Y.clone();
X.rsolve( sol, pvt, pvt.size() );
DoubleVector r = Y.getColumn( pvt.size(), n-1, 0);
double sde = Math.sqrt(r.sum2() / r.size());
DoubleVector aHat = Y.getColumn( 0, pvt.size()-1, 0).times( 1./sde );
DoubleVector aTilde = null;
switch( paceEstimator) {