anchor.diff:      do i=2,it
         if (alpha(i)==0.or.beta(i)<0) then
           print *,'Unable to compute eigenvalue number',i,' and '
  -        do j=1,it
  +        do j=1,i
             print *,j,' alpha:',alpha(j),' beta:',beta
           enddo
           return
  
positive.code:       rho_prev = rho_curr
        ! check
        res_norm = Vect_dot_product(r,r)
        write (stream,*) "Norm is", res_norm
        ratio_norm = res_norm / init_norm
        if (ismaster()) &
             write(stream, '(i5,a,e22.15)') it,': res_norm=',dsqrt(res_norm)
  
negative.code:           endif
  
            if (.not.present(cdat)) then
              if (sctls%smoothers==-1) then
                allocate(tmpsol2(A%nrows))
                tmpsol2=0.0_rk
                call exact_sparse_multismoother(tmpsol2,A,rhs)
                call SpMtx_Ax(crhs,Restrict,tmpsol2,dozero=.true.) ! restriction
              else
                if (sctls%method>1.and.sctls%method/=5) then ! multiplicative Schwarz
                  call SpMtx_Ax(crhs,Restrict,res,dozero=.true.) ! restriction
                else
                  call SpMtx_Ax(crhs,Restrict,rhs,dozero=.true.) ! restriction
                endif
              endif
            endif
            !csol=0.0_rk
            write (stream,*) &
  