anchor.diff:        if (sctls%method/=0) then
           call Add_common_interf(z,A,Msh)
         endif
  -!write(stream,*)'localz==:',z
  +!call Print_Glob_Vect(z,Msh,'global z===')
   !call Print_Glob_Vect(z,Msh,'global z===')
         ! compute current rho
         rho_curr = Vect_dot_product(r,z)
  
positive.code:           !  deallocate(A%M_bound)
            !endif
            call SpMtx_arrange(A,D_SpMtx_ARRNG_ROWS,sort=.true.)
     write(stream,*)'AAAAAAAAAAAAAAAA is:'
     call SpMtx_printRaw(A)
     call SpMtx_printRaw(A_interf_)
     write(stream,*)'AAAAAAAAAAAAAAAA :'
     call flush(stream)
            call sparse_multisolve(sol=sol,A=A,rhs=rhs, &
                                   A_interf_=A_interf_, &
                                    refactor=refactor_) !fine solves 
            ! put the original structure and orders back:
            A%indi=A_tmp%indi
            A%indj=A_tmp%indj
  
negative.code:     real(kind=rk),dimension(:),pointer :: sol
      type(SpMtx)                        :: A
      real(kind=rk),dimension(:),pointer :: rhs
      type(Mesh),intent(in)              :: M ! Mesh
      real(kind=rk),dimension(:),pointer :: res ! -- residual vector, allocated
                                                ! here for multiplicative Schwarz
      type(SpMtx),optional               :: A_interf_  ! matr@interf.
      type(SpMtx),optional               :: CoarseMtx_ ! Coarse matrix
      type(SpMtx),optional               :: Restrict ! Restriction matrix
      logical,intent(inout),optional :: refactor_
      type(CoarseData),intent(inout),optional :: cdat
      ! ----- local: ------
      real(kind=rk),dimension(:),pointer,save :: csol,crhs,tmpsol,tmpsol2,clrhs,ctmp
  