Merge branch 'bugfixes_current_release'

Makefile.am

@@ -196,6 +196,10 @@ dist_files_DATA = \
   test/Fortran/test_real.F90 \
   test/Fortran/test_real_with_c.F90 \
   test/Fortran/test_toeplitz.F90 \
+  test/Fortran/test_transpose_multiply_real.F90 \
+  test/Fortran/test_transpose_multiply_complex.F90 \
+  test/Fortran/test_cholesky_real.F90 \
+  test/Fortran/test_cholesky_complex.F90 \
   src/elpa2_print_kernels.F90
 
 dist_doc_DATA = README.md USERS_GUIDE.md INSTALL.md CONTRIBUTING.md LICENSE Changelog COPYING/COPYING COPYING/gpl.txt COPYING/lgpl.txt
@@ -219,6 +223,10 @@ noinst_PROGRAMS = \
   elpa2_test_real_api@SUFFIX@ \
   elpa2_test_complex_api@SUFFIX@ \
   elpa1_real_toeplitz@SUFFIX@ \
+  elpa1_real_transpose_multiply@SUFFIX@ \
+  elpa1_complex_transpose_multiply@SUFFIX@ \
+  elpa1_real_cholesky@SUFFIX@ \
+  elpa1_complex_cholesky@SUFFIX@ \
   elpa1_test_real_with_c@SUFFIX@ \
   elpa1_test_real_c_version@SUFFIX@ \
   elpa1_test_complex_c_version@SUFFIX@ \
@@ -278,6 +286,26 @@ elpa1_real_toeplitz@SUFFIX@_LDADD = $(build_lib)
 elpa1_real_toeplitz@SUFFIX@_FCFLAGS = $(AM_FCFLAGS) @FC_MODOUT@private_modules @FC_MODINC@private_modules
 EXTRA_elpa1_real_toeplitz@SUFFIX@_DEPENDENCIES = test/Fortran/elpa_print_headers.X90
 
+elpa1_real_transpose_multiply@SUFFIX@_SOURCES = test/Fortran/test_transpose_multiply_real.F90
+elpa1_real_transpose_multiply@SUFFIX@_LDADD = $(build_lib)
+elpa1_real_transpose_multiply@SUFFIX@_FCFLAGS = $(AM_FCFLAGS) @FC_MODOUT@private_modules @FC_MODINC@private_modules
+EXTRA_elpa1_real_transpose_multiply@SUFFIX@_DEPENDENCIES = test/Fortran/elpa_print_headers.X90
+
+elpa1_complex_transpose_multiply@SUFFIX@_SOURCES = test/Fortran/test_transpose_multiply_complex.F90
+elpa1_complex_transpose_multiply@SUFFIX@_LDADD = $(build_lib)
+elpa1_complex_transpose_multiply@SUFFIX@_FCFLAGS = $(AM_FCFLAGS) @FC_MODOUT@private_modules @FC_MODINC@private_modules
+EXTRA_elpa1_complex_transpose_multiply@SUFFIX@_DEPENDENCIES = test/Fortran/elpa_print_headers.X90
+
+elpa1_real_cholesky@SUFFIX@_SOURCES = test/Fortran/test_cholesky_real.F90
+elpa1_real_cholesky@SUFFIX@_LDADD = $(build_lib)
+elpa1_real_cholesky@SUFFIX@_FCFLAGS = $(AM_FCFLAGS) @FC_MODOUT@private_modules @FC_MODINC@private_modules
+EXTRA_elpa1_real_cholesky@SUFFIX@_DEPENDENCIES = test/Fortran/elpa_print_headers.X90
+
+elpa1_complex_cholesky@SUFFIX@_SOURCES = test/Fortran/test_cholesky_complex.F90
+elpa1_complex_cholesky@SUFFIX@_LDADD = $(build_lib)
+elpa1_complex_cholesky@SUFFIX@_FCFLAGS = $(AM_FCFLAGS) @FC_MODOUT@private_modules @FC_MODINC@private_modules
+EXTRA_elpa1_complex_cholesky@SUFFIX@_DEPENDENCIES = test/Fortran/elpa_print_headers.X90
+
 elpa1_test_real_with_c@SUFFIX@_SOURCES = test/Fortran/test_real_with_c.F90
 elpa1_test_real_with_c@SUFFIX@_LDADD = $(build_lib)
 elpa1_test_real_with_c@SUFFIX@_FCFLAGS = $(AM_FCFLAGS) @FC_MODOUT@private_modules @FC_MODINC@private_modules
@@ -341,6 +369,10 @@ check_SCRIPTS = \
   elpa2_test_real_api@SUFFIX@.sh \
   elpa2_test_complex_api@SUFFIX@.sh \
   elpa1_real_toeplitz@SUFFIX@.sh \
+  elpa1_real_transpose_multiply@SUFFIX@.sh \
+  elpa1_complex_transpose_multiply@SUFFIX@.sh \
+  elpa1_real_cholesky@SUFFIX@.sh \
+  elpa1_complex_cholesky@SUFFIX@.sh \
   elpa2_print_kernels@SUFFIX@ \
   elpa1_test_real_c_version@SUFFIX@.sh \
   elpa1_test_complex_c_version@SUFFIX@.sh \
@@ -388,6 +420,7 @@ CLEANFILES = \
   elpa2_test*\
   elpa2_real* \
   elpa1_real* \
+  *.sh \
   *.i
 
 clean-local:

README.md

@@ -30,7 +30,7 @@ consisting of the following organizations:
 - Technische Universität München, Lehrstuhl für Informatik mit
   Schwerpunkt Wissenschaftliches Rechnen ,
 - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-- Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+- Max-Plack-Institut für Mathematik in den Naturwissenschaften,
   Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
   and
 - IBM Deutschland GmbH

src/aligned_mem.F90

@@ -10,7 +10,7 @@
 !    - Technische Universität München, Lehrstuhl für Informatik mit
 !      Schwerpunkt Wissenschaftliches Rechnen ,
 !    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-!    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+!    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 !      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 !      and
 !    - IBM Deutschland GmbH

src/elpa1.F90

@@ -9,7 +9,7 @@
 !    - Technische Universität München, Lehrstuhl für Informatik mit
 !      Schwerpunkt Wissenschaftliches Rechnen ,
 !    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-!    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+!    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 !      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 !      and
 !    - IBM Deutschland GmbH
@@ -65,7 +65,7 @@
 !>    - Technische Universität München, Lehrstuhl für Informatik mit
 !>      Schwerpunkt Wissenschaftliches Rechnen ,
 !>    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-!>    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+!>    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 !>      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 !>      and
 !>    - IBM Deutschland GmbH

src/elpa1_auxiliary.F90

@@ -10,7 +10,7 @@
 !    - Technische Universität München, Lehrstuhl für Informatik mit
 !      Schwerpunkt Wissenschaftliches Rechnen ,
 !    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-!    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+!    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 !      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 !      and
 !    - IBM Deutschland GmbH
@@ -105,12 +105,12 @@ module elpa1_auxiliary
 !>                              Only upper triangle is needs to be set.
 !>                              The lower triangle is not referenced.
 !> \param  lda                  Leading dimension of a
-!> \param                       matrixCols  local columns of matrix a
 !> \param  nblk                 blocksize of cyclic distribution, must be the same in both directions!
+!> \param  matrixCols           local columns of matrix a
 !> \param  mpi_comm_rows        MPI communicator for rows
 !> \param  mpi_comm_cols        MPI communicator for columns
 !> \param wantDebug             logical, more debug information on failure
-!> \result succes                logical, reports success or failure
+!> \result succes               logical, reports success or failure
     interface invert_trm_real
       module procedure elpa_invert_trm_real
     end interface
@@ -124,8 +124,8 @@ module elpa1_auxiliary
 !>                              On return, the upper triangle contains the Cholesky factor
 !>                              and the lower triangle is set to 0.
 !> \param  lda                  Leading dimension of a
-!> \param                       matrixCols  local columns of matrix a
 !> \param  nblk                 blocksize of cyclic distribution, must be the same in both directions!
+!> \param  matrixCols           local columns of matrix a
 !> \param  mpi_comm_rows        MPI communicator for rows
 !> \param  mpi_comm_cols        MPI communicator for columns
 !> \param wantDebug             logical, more debug information on failure
@@ -142,8 +142,8 @@ module elpa1_auxiliary
 !>                              Only upper triangle is needs to be set.
 !>                              The lower triangle is not referenced.
 !> \param  lda                  Leading dimension of a
-!> \param                       matrixCols  local columns of matrix a
 !> \param  nblk                 blocksize of cyclic distribution, must be the same in both directions!
+!> \param  matrixCols           local columns of matrix a
 !> \param  mpi_comm_rows        MPI communicator for rows
 !> \param  mpi_comm_cols        MPI communicator for columns
 !> \param wantDebug             logical, more debug information on failure
@@ -236,8 +236,8 @@ module elpa1_auxiliary
 !>                              On return, the upper triangle contains the Cholesky factor
 !>                              and the lower triangle is set to 0.
 !> \param  lda                  Leading dimension of a
-!> \param                       matrixCols  local columns of matrix a
 !> \param  nblk                 blocksize of cyclic distribution, must be the same in both directions!
+!> \param  matrixCols           local columns of matrix a
 !> \param  mpi_comm_rows        MPI communicator for rows
 !> \param  mpi_comm_cols        MPI communicator for columns
 !> \param wantDebug             logical, more debug information on failure
@@ -322,7 +322,6 @@ module elpa1_auxiliary
       tmatc = 0
 
       do n = 1, na, nblk
-
         ! Calculate first local row and column of the still remaining matrix
         ! on the local processor
 
@@ -341,7 +340,7 @@ module elpa1_auxiliary
 
             call dpotrf('U',na-n+1,a(l_row1,l_col1),lda,info)
             if (info/=0) then
-              if (wantDebug) write(error_unit,*) "elpa_cholesky_real: Error in dpotrf"
+              if (wantDebug) write(error_unit,*) "elpa_cholesky_real: Error in dpotrf 1: ",info
               success = .false.
               return
             endif
@@ -361,7 +360,7 @@ module elpa1_auxiliary
 
             call dpotrf('U',nblk,a(l_row1,l_col1),lda,info)
             if (info/=0) then
-              if (wantDebug) write(error_unit,*) "elpa_cholesky_real: Error in dpotrf"
+              if (wantDebug) write(error_unit,*) "elpa_cholesky_real: Error in dpotrf 2: ",info
               success = .false.
               return
             endif
@@ -442,12 +441,12 @@ module elpa1_auxiliary
 !>                              Only upper triangle is needs to be set.
 !>                              The lower triangle is not referenced.
 !> \param  lda                  Leading dimension of a
-!> \param                       matrixCols  local columns of matrix a
 !> \param  nblk                 blocksize of cyclic distribution, must be the same in both directions!
+!> \param  matrixCols           local columns of matrix a
 !> \param  mpi_comm_rows        MPI communicator for rows
 !> \param  mpi_comm_cols        MPI communicator for columns
 !> \param wantDebug             logical, more debug information on failure
-!> \result succes                logical, reports success or failure
+!> \result succes               logical, reports success or failure
      function elpa_invert_trm_real(na, a, lda, nblk, matrixCols, mpi_comm_rows, mpi_comm_cols, wantDebug) result(success)
        use precision
        use elpa1_compute
@@ -600,8 +599,8 @@ module elpa1_auxiliary
 !>                              On return, the upper triangle contains the Cholesky factor
 !>                              and the lower triangle is set to 0.
 !> \param  lda                  Leading dimension of a
-!> \param                       matrixCols  local columns of matrix a
 !> \param  nblk                 blocksize of cyclic distribution, must be the same in both directions!
+!> \param  matrixCols           local columns of matrix a
 !> \param  mpi_comm_rows        MPI communicator for rows
 !> \param  mpi_comm_cols        MPI communicator for columns
 !> \param wantDebug             logical, more debug information on failure
@@ -805,12 +804,12 @@ module elpa1_auxiliary
 !>                              Only upper triangle is needs to be set.
 !>                              The lower triangle is not referenced.
 !> \param  lda                  Leading dimension of a
-!> \param                       matrixCols  local columns of matrix a
 !> \param  nblk                 blocksize of cyclic distribution, must be the same in both directions!
+!> \param  matrixCols           local columns of matrix a
 !> \param  mpi_comm_rows        MPI communicator for rows
 !> \param  mpi_comm_cols        MPI communicator for columns
 !> \param wantDebug             logical, more debug information on failure
-!> \result succes                logical, reports success or failure
+!> \result succes               logical, reports success or failure
     function elpa_invert_trm_complex(na, a, lda, nblk, matrixCols, mpi_comm_rows, mpi_comm_cols, wantDebug) result(success)
 
        use precision

src/elpa1_compute_private.F90

@@ -10,7 +10,7 @@
 !    - Technische Universität München, Lehrstuhl für Informatik mit
 !      Schwerpunkt Wissenschaftliches Rechnen ,
 !    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-!    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+!    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 !      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 !      and
 !    - IBM Deutschland GmbH

src/elpa2.F90

@@ -10,7 +10,7 @@
 !    - Technische Universität München, Lehrstuhl für Informatik mit
 !      Schwerpunkt Wissenschaftliches Rechnen ,
 !    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-!    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+!    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 !      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 !      and
 !    - IBM Deutschland GmbH

src/elpa2_compute.F90

@@ -10,7 +10,7 @@
 !    - Technische Universität München, Lehrstuhl für Informatik mit
 !      Schwerpunkt Wissenschaftliches Rechnen ,
 !    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-!    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+!    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 !      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 !      and
 !    - IBM Deutschland GmbH

src/elpa2_kernels/elpa2_kernels_asm_x86_64.s

@@ -10,7 +10,7 @@
 #    - Technische Universität München, Lehrstuhl für Informatik mit
 #      Schwerpunkt Wissenschaftliches Rechnen ,
 #    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-#    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+#    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 #      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 #      and
 #    - IBM Deutschland GmbH

src/elpa2_kernels/elpa2_kernels_complex.F90

@@ -10,7 +10,7 @@
 !    - Technische Universität München, Lehrstuhl für Informatik mit
 !      Schwerpunkt Wissenschaftliches Rechnen ,
 !    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-!    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+!    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 !      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 !      and
 !    - IBM Deutschland GmbH

src/elpa2_kernels/elpa2_kernels_complex_avx-avx2_1hv.c

@@ -10,7 +10,7 @@
 //    - Technische Universität München, Lehrstuhl für Informatik mit
 //      Schwerpunkt Wissenschaftliches Rechnen ,
 //    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-//    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+//    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 //      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 //      and
 //    - IBM Deutschland GmbH

src/elpa2_kernels/elpa2_kernels_complex_avx-avx2_2hv.c

@@ -10,7 +10,7 @@
 //    - Technische Universität München, Lehrstuhl für Informatik mit
 //      Schwerpunkt Wissenschaftliches Rechnen ,
 //    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-//    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+//    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 //      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 //      and
 //    - IBM Deutschland GmbH

src/elpa2_kernels/elpa2_kernels_complex_simple.F90

@@ -10,7 +10,7 @@
 !    - Technische Universität München, Lehrstuhl für Informatik mit
 !      Schwerpunkt Wissenschaftliches Rechnen ,
 !    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-!    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+!    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 !      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 !      and
 !    - IBM Deutschland GmbH

src/elpa2_kernels/elpa2_kernels_complex_sse_1hv.c

@@ -10,7 +10,7 @@
 //    - Technische Universität München, Lehrstuhl für Informatik mit
 //      Schwerpunkt Wissenschaftliches Rechnen ,
 //    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-//    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+//    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 //      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 //      and
 //    - IBM Deutschland GmbH

src/elpa2_kernels/elpa2_kernels_complex_sse_2hv.c

@@ -10,7 +10,7 @@
 //    - Technische Universität München, Lehrstuhl für Informatik mit
 //      Schwerpunkt Wissenschaftliches Rechnen ,
 //    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-//    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+//    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 //      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 //      and
 //    - IBM Deutschland GmbH

src/elpa2_kernels/elpa2_kernels_real.F90

@@ -10,7 +10,7 @@
 !    - Technische Universität München, Lehrstuhl für Informatik mit
 !      Schwerpunkt Wissenschaftliches Rechnen ,
 !    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-!    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+!    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 !      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 !      and
 !    - IBM Deutschland GmbH

src/elpa2_kernels/elpa2_kernels_real_avx-avx2_2hv.c

@@ -10,7 +10,7 @@
 //    - Technische Universität München, Lehrstuhl für Informatik mit
 //      Schwerpunkt Wissenschaftliches Rechnen ,
 //    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-//    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+//    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 //      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 //      and
 //    - IBM Deutschland GmbH

src/elpa2_kernels/elpa2_kernels_real_avx-avx2_4hv.c

@@ -10,7 +10,7 @@
 //    - Technische Universität München, Lehrstuhl für Informatik mit
 //      Schwerpunkt Wissenschaftliches Rechnen ,
 //    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-//    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+//    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 //      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 //      and
 //    - IBM Deutschland GmbH

src/elpa2_kernels/elpa2_kernels_real_avx-avx2_6hv.c

@@ -10,7 +10,7 @@
 //    - Technische Universität München, Lehrstuhl für Informatik mit
 //      Schwerpunkt Wissenschaftliches Rechnen ,
 //    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-//    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+//    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 //      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 //      and
 //    - IBM Deutschland GmbH

src/elpa2_kernels/elpa2_kernels_real_bgp.f90

@@ -10,7 +10,7 @@
 !    - Technische Universität München, Lehrstuhl für Informatik mit
 !      Schwerpunkt Wissenschaftliches Rechnen ,
 !    - Fritz-Haber-Institut, Berlin, Abt. Theorie,
-!    - Max-Plack-Institut für Mathematik in den Naturwissenschaftrn,
+!    - Max-Plack-Institut für Mathematik in den Naturwissenschaften,
 !      Leipzig, Abt. Komplexe Strukutren in Biologie und Kognition,
 !      and
 !    - IBM Deutschland GmbH