diff options
| -rw-r--r-- | fortran/navy_four_levels/FourLevelPulseProp_Double.f95 | 70 |
1 files changed, 0 insertions, 70 deletions
diff --git a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 index 21249de..5eb6af8 100644 --- a/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 +++ b/fortran/navy_four_levels/FourLevelPulseProp_Double.f95 @@ -508,76 +508,6 @@ stop end -! - subroutine IMatConstruct(Om,Imat) -! -! NOTE: This subroutine actually calculates -Imat because we need to solve LPsi=-Imat. - implicit none - integer i,nmat - parameter (nmat=3) - real*8 Om - complex*8 Imat(nmat) - - Imat(1)=-cmplx(0.,-0.5*Om) - Imat(2)=-cmplx(0.,0.5*Om) - Imat(3)=cmplx(0.,0.) - return - end - - -! - subroutine LMatConstruct(Om,del,L) -! - - implicit none - integer nmat - parameter (nmat=3) - real*8 del,ga12,Om,W21 - complex*8 L(nmat,nmat) - common/para/ga12,W21 - - - L(1,1)=cmplx (-ga12,-del) - L(1,3)=cmplx(0.,Om) - L(2,2)=cmplx(-ga12,del) - L(2,3)=cmplx(0.,-Om) - L(3,1)=cmplx(0.,0.5*Om) - L(3,2)=cmplx(0.,-0.5*Om) - L(3,3)=cmplx(-W21,0.) - return - end - -! - subroutine MatCheck(L,Linv) -! - - implicit none - integer i,j,k,nmat - parameter (nmat=3) - complex*8 L(nmat,nmat),Linv(nmat,nmat),Res(nmat,nmat) - - write (*,*)'L = ' - do 10 i=1,nmat - write (*,120)(L(i,j),j=1,nmat) -10 continue - write (*,*)'Linv = ' - do 20 i=1,nmat - write (*,120)(Linv(i,j),j=1,nmat) -20 continue - write (*,*)'Res = ' - do 50 i=1,nmat - do 40 j=1,nmat - Res(i,j)=cmplx(0.,0.) - do 30 k=1,nmat - Res(i,j)=Res(i,j)+Linv(i,k)*L(k,j) -30 continue -40 continue - write (*,120)(Res(i,j),j=1,nmat) -50 continue -120 format(1x,3(f8.4,'+i',f8.4,' ')) - - return - end |