Commit e17c0254 authored by Florian Goth's avatar Florian Goth
Browse files

remove unnecessary checks

parent f523f5aa
......@@ -118,17 +118,6 @@
endif
else
write (*,*) "MSCB, method:", method
if (op%N - 2* (op%N/2) /= 0) then
write (*,*) "operator dimension not divisible by two."
error stop 2
endif
! check for chemical potential
! do i = 1, op%N
! if (Op%O(i,i) /= 0) then
! write (*,*) "chemical potential not supported."
! error stop 2
! endif
! enddo
if (method == 1) then
allocate(eulerexp)
call eulerexp%init(op)
......@@ -138,7 +127,7 @@
call mscbexp%init(op, method)
call ExpOpT_vec%pushback(mscbexp)
endif
! Let's apply a quick check to see, if the resulting matrix gives a symmetric time evolution.
ndim = size(Op%O, 1)
allocate(mat(ndim, ndim))
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment