anchor.diff:          end do
         end if
   
  -      result = result * arbitrary_basis_element(tt)
  +      result = result * arbitrary_basis_element(t)
   
       end subroutine Lmult1
   
  
positive.code:     ! declarations
      real(dp), dimension(:,:), allocatable  :: Ham
      complex(dpc), allocatable:: rho1(:,:), prhox1(:,:), prhodx1(:,:)
      !complex(dpc), allocatable:: rho2(:,:), prhox2(:,:), prhodx2(:,:)
      complex(dpc), dimension(:,:,:,:,:), allocatable      :: Ueg, Uee, Ugg
      complex(dpc), dimension(:,:), private, allocatable   :: A, U, VT
      complex(dpc), dimension(:), private, allocatable     :: B
  
negative.code:     subroutine clean_lindblad_fit()
      end subroutine clean_lindblad_fit
  
      subroutine Lmult1 (t, rhoin, result)
        complex(dpc), intent(in)   :: rhoin(:,:)
        real(dp), intent(in)       :: t
        complex(dpc), intent(out)  :: result(:,:)
        integer(i4b) :: i
        complex(dpc), parameter:: iconst = dcmplx(0.0, 1.0)
  