anchor.diff:      deallocate(pa)
   
   
  -  end subroutine tt_kls
  -
  -  subroutine ztt_kls(d,n,m,ra,crA, crY0, ry, tau, rmax, kickrank, nswp, verb)
  +  end subroutine tt_ksl
  +
  +  subroutine ztt_ksl(d,n,m,ra,crA, crY0, ry, tau, rmax, kickrank, nswp, verb)
       use dispmodule
       use matrix_util ! dqr
       use ttals   ! als stuff
       use estnorm  !1-norm estimation
  -    use bfun_kls !bfun declarations
  -    use sfun_kls !We have to do this
  +    use bfun_ksl !bfun declarations
  +    use sfun_ksl !We have to do this
       use explib !Krylov exponential
       implicit none
       integer,intent(in) :: d,n(d),m(d),ra(d+1), rmax
  
positive.code:   end subroutine zinit_bfun_main
  
  
  end module bfun_kls
  
  module sfun_kls
  
    integer, private :: rx1, rx2, ra, ry1, ry2
    double precision, private, pointer :: phi1(:), phi2(:)
  
negative.code:     deallocate(pa)
  
  
    end subroutine ztt_ksl
    !! What we have: we have a starting vector + a matrix (no vector X!) 
    
  
  end module dyn_tt
  