4.10.9 Update

Some updates mostly in IOTOOLS and LINALG.
* introduced to_lower and to_upper procedures to change case of a given alphabetic string
* fixed some bugs, extended functionalities and exported operator .kx. for Kronecker product
of matrices (used in lattice real-space DMRG).

src/SF_IOTOOLS/IOFILE.f90

@@ -46,6 +46,8 @@ module IOFILE
   public :: str
   public :: txtfy !obsolete
   public :: reg
+  public :: to_lower
+  public :: to_upper
   !
   public :: file_size
   public :: file_length
@@ -70,7 +72,7 @@ module IOFILE
 
 
 
-  
+
   !+-----------------------------------------------------------------+
   !PURPOSE  : 
   !+-----------------------------------------------------------------+
@@ -397,6 +399,37 @@ end subroutine create_data_dir
 
 
 
+  function to_upper(StrIn) result(StrOut)
+    character(len=*), intent(in) :: strIn
+    character(len=len(strIn))    :: strOut
+    integer :: i
+    do i = 1,len(StrIn)
+       select case(StrIn(i:i))
+       case("a":"z")
+          StrOut(i:i) = achar(iachar(StrIn(i:i))-32)
+       case default
+          StrOut(i:i) = StrIn(i:i)
+       end select
+    end do
+  end function to_upper
+
+  function to_lower(StrIn) result(StrOut)
+    character(len=*), intent(in) :: strIn
+    character(len=len(strIn))    :: strOut
+    integer :: i
+    do i = 1,len(StrIn)
+       select case(StrIn(i:i))
+       case("A":"Z")
+          StrOut(i:i) = achar(iachar(StrIn(i:i))+32)
+       case default
+          StrOut(i:i) = StrIn(i:i)
+       end select
+    end do
+  end function to_lower
+
+
+
+
   ! function i_to_ch(i4) result(string)
   !   character(len=32) :: string
   !   integer           :: i4

src/SF_IOTOOLS/SF_IOTOOLS.f90

@@ -16,6 +16,8 @@ MODULE SF_IOTOOLS
   public :: str
   public :: txtfy !obsolete
   public :: reg
+  public :: to_lower
+  public :: to_upper
   !
   public :: file_size
   public :: file_length

src/SF_LINALG/SF_LINALG.f90

@@ -281,23 +281,35 @@ module SF_LINALG
   !Kroenecker product of matrices
   public :: kron
   public :: kronecker_product
-  public :: kroenecker_product
+  public :: operator(.kx.)
   !outer product of two 1d arrays to form a matrix
   public :: outerprod
   public :: cross_product
   public :: s3_product
   !
   interface kron
-     module procedure i_kronecker_product,d_kronecker_product,c_kronecker_product
+     module procedure :: i_kronecker_product
+     module procedure :: d_kronecker_product
+     module procedure :: dc_kronecker_product
+     module procedure :: cd_kronecker_product
+     module procedure :: c_kronecker_product
   end interface kron
   !
   interface kronecker_product
-     module procedure i_kronecker_product,d_kronecker_product,c_kronecker_product
+     module procedure :: i_kronecker_product
+     module procedure :: d_kronecker_product
+     module procedure :: dc_kronecker_product
+     module procedure :: cd_kronecker_product
+     module procedure :: c_kronecker_product
   end interface kronecker_product
   !
-  interface kroenecker_product
-     module procedure i_kronecker_product,d_kronecker_product,c_kronecker_product
-  end interface kroenecker_product
+  interface operator(.kx.)
+     module procedure :: i_kronecker_product
+     module procedure :: d_kronecker_product
+     module procedure :: dc_kronecker_product
+     module procedure :: cd_kronecker_product
+     module procedure :: c_kronecker_product
+  end interface operator(.kx.)
   !
   interface outerprod
      module procedure outerprod_d,outerprod_c

src/SF_LINALG/linalg_external_products.f90

@@ -23,21 +23,51 @@ function d_kronecker_product(A,B) result(AxB)
   integer            :: rowA,colA
   integer            :: rowB,colB
   real(8)            :: AxB(size(A,1)*size(B,1),size(A,2)*size(B,2))
-  AxB = 0
+  AxB = 0d0
   rowA=size(A,1) ; colA=size(A,2)
   rowB=size(B,1) ; colB=size(B,2)
   forall(i=1:rowA,j=1:colA)
      AxB(1+rowB*(i-1):rowB*i,1+colB*(j-1):colB*j)  =  A(i,j)*B(:,:)
   end forall
 end function d_kronecker_product
 !
+function dc_kronecker_product(A,B) result(AxB)
+  real(8),intent(in) :: A(:,:)
+  complex(8),intent(in) :: B(:,:)
+  integer            :: i,j
+  integer            :: rowA,colA
+  integer            :: rowB,colB
+  complex(8)            :: AxB(size(A,1)*size(B,1),size(A,2)*size(B,2))
+  AxB = zero
+  rowA=size(A,1) ; colA=size(A,2)
+  rowB=size(B,1) ; colB=size(B,2)
+  forall(i=1:rowA,j=1:colA)
+     AxB(1+rowB*(i-1):rowB*i,1+colB*(j-1):colB*j)  =  A(i,j)*B(:,:)
+  end forall
+end function dc_kronecker_product
+!
+function cd_kronecker_product(A,B) result(AxB)
+  complex(8),intent(in) :: A(:,:)
+  real(8),intent(in)    :: B(:,:)
+  integer               :: i,j
+  integer               :: rowA,colA
+  integer               :: rowB,colB
+  complex(8)            :: AxB(size(A,1)*size(B,1),size(A,2)*size(B,2))
+  AxB = zero
+  rowA=size(A,1) ; colA=size(A,2)
+  rowB=size(B,1) ; colB=size(B,2)
+  forall(i=1:rowA,j=1:colA)
+     AxB(1+rowB*(i-1):rowB*i,1+colB*(j-1):colB*j)  =  A(i,j)*B(:,:)
+  end forall
+end function cd_kronecker_product
+!
 function c_kronecker_product(A,B) result(AxB)
   complex(8),intent(in) :: A(:,:), B(:,:)
   integer               :: i,j
   integer               :: rowA,colA
   integer               :: rowB,colB
   complex(8)            :: AxB(size(A,1)*size(B,1),size(A,2)*size(B,2))
-  AxB = 0
+  AxB = zero
   rowA=size(A,1) ; colA=size(A,2)
   rowB=size(B,1) ; colB=size(B,2)
   forall(i=1:rowA,j=1:colA)