anchor.diff:      integer :: typ, order
       integer, optional :: typ0, order0
       real(8) :: ermax, min_res
  -    real(8) anorm, anorm_t
  +    real(8) anorm
       typ = 2 ! 1 - KSL, 2 - KSL-symm
       if ( present(typ0) ) then
           typ = typ0 
  
positive.code:           call init_bfun_sizes(ry(i),n(i),ry(i+1),ry(i),n(i),ry(i+1),ra(i),ra(i+1),ry(i)*n(i)*ry(i+1),ry(i)*n(i)*ry(i+1))
            call dinit_bfun_main(phinew(i)%p,crA(pa(i):pa(i+1)-1),phinew(i+1)%p)
            anorm = normest(ry(i)*n(i)*ry(i+1),4, dmatvec, dmatvec_transp)
            anorm_t = norm_true(ry(i)*n(i)*ry(i+1), dmatvec, dmatvec_transp)
            print *, 'anorm, anorm_t = ', anorm, anorm_t
            call dexp_mv(ry(i)*n(i)*ry(i+1),order,tau/2,crnew(i)%p,curcr,eps,anorm,dmatvec)
            if ( i < d ) then
               !In this case, we have to put S(i+1) backwards in time (heh)
  
negative.code:     integer :: max_phi_size, max_core_size, max_R_size
      integer :: order
      real(8) :: ermax, min_res
      real(8) anorm
      complex(8) ZERO, ONE
      parameter( ZERO=(0.0d0,0.0d0), ONE=(1.0d0,0.0d0) )
  
  