function x = CGPre(A,b,P)
  x = zeros(size(A,1),1);
  x = CGxPre(A,b,x,P);
