anchor.diff:  	real, intent(out),optional :: lpc(size(ll,2),nkeep), &
   		& leof(size(ll,1),nkeep), rpc(size(rr,2),nkeep), &
   		& reof(size(rr,1),nkeep),ev(nkeep)
  -	logical, intent(in),  optional :: bLargeMatrix, usecorr
  +	logical, intent(in),  optional :: bigmat, usecorr
   	real, intent(out), optional :: ev_sum
   	integer, intent(out),  optional :: info
   
  
positive.code: 	!	- lw:    Left weights
  	!	- rw:    Right weights
  	!	- usecorr:  Use correlations instead of covariances
  	!	- bLargeMatrix: Use gesdd instead of gesvd (faster for large matrices, but uses more workspace) [default:.false.]
  	!
  	! Dependencies:
  	!	[sd]gemm(BLAS) gesvd(LAPACK95) gesdd(LAPACK95)
  
  
  
negative.code: 
  	! Use divide/conquer algo?
  	! ------------------------
  	if(.not.present(bigmat))then
  		zbigmat = .false.
  	else
  		zbigmat = bigmat
  	end if
  
  	! Weights
  