Avoid LX temporaries in X transpose

merge-paraggr-newops^2
Salvatore Filippone 5 years ago
parent f28e3a9ea9
commit 84a8b73416

@ -707,7 +707,7 @@ subroutine psb_c_simple_glob_transpose_ip(ain,desc_a,info)
! that the same DESC_A works for both A and A^T, which ! that the same DESC_A works for both A and A^T, which
! essentially means that A has a symmetric pattern. ! essentially means that A has a symmetric pattern.
! !
type(psb_lc_coo_sparse_mat) :: tmpc1, tmpc2 type(psb_c_coo_sparse_mat) :: tmpc1, tmpc2
integer(psb_ipk_) :: nz1, nz2, nzh, nz integer(psb_ipk_) :: nz1, nz2, nzh, nz
integer(psb_ipk_) :: ictxt, me, np integer(psb_ipk_) :: ictxt, me, np
integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz
@ -758,7 +758,7 @@ subroutine psb_c_simple_glob_transpose(ain,aout,desc_a,info)
! that the same DESC_A works for both A and A^T, which ! that the same DESC_A works for both A and A^T, which
! essentially means that A has a symmetric pattern. ! essentially means that A has a symmetric pattern.
! !
type(psb_lc_coo_sparse_mat) :: tmpc1, tmpc2 type(psb_c_coo_sparse_mat) :: tmpc1, tmpc2
integer(psb_ipk_) :: nz1, nz2, nzh, nz integer(psb_ipk_) :: nz1, nz2, nzh, nz
integer(psb_ipk_) :: ictxt, me, np integer(psb_ipk_) :: ictxt, me, np
integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz

@ -747,7 +747,7 @@ subroutine psb_d_simple_glob_transpose_ip(ain,desc_a,info)
! that the same DESC_A works for both A and A^T, which ! that the same DESC_A works for both A and A^T, which
! essentially means that A has a symmetric pattern. ! essentially means that A has a symmetric pattern.
! !
type(psb_ld_coo_sparse_mat) :: tmpc1, tmpc2 type(psb_d_coo_sparse_mat) :: tmpc1, tmpc2
integer(psb_ipk_) :: nz1, nz2, nzh, nz integer(psb_ipk_) :: nz1, nz2, nzh, nz
integer(psb_ipk_) :: ictxt, me, np integer(psb_ipk_) :: ictxt, me, np
integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz
@ -798,7 +798,7 @@ subroutine psb_d_simple_glob_transpose(ain,aout,desc_a,info)
! that the same DESC_A works for both A and A^T, which ! that the same DESC_A works for both A and A^T, which
! essentially means that A has a symmetric pattern. ! essentially means that A has a symmetric pattern.
! !
type(psb_ld_coo_sparse_mat) :: tmpc1, tmpc2 type(psb_d_coo_sparse_mat) :: tmpc1, tmpc2
integer(psb_ipk_) :: nz1, nz2, nzh, nz integer(psb_ipk_) :: nz1, nz2, nzh, nz
integer(psb_ipk_) :: ictxt, me, np integer(psb_ipk_) :: ictxt, me, np
integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz

@ -707,7 +707,7 @@ subroutine psb_s_simple_glob_transpose_ip(ain,desc_a,info)
! that the same DESC_A works for both A and A^T, which ! that the same DESC_A works for both A and A^T, which
! essentially means that A has a symmetric pattern. ! essentially means that A has a symmetric pattern.
! !
type(psb_ls_coo_sparse_mat) :: tmpc1, tmpc2 type(psb_s_coo_sparse_mat) :: tmpc1, tmpc2
integer(psb_ipk_) :: nz1, nz2, nzh, nz integer(psb_ipk_) :: nz1, nz2, nzh, nz
integer(psb_ipk_) :: ictxt, me, np integer(psb_ipk_) :: ictxt, me, np
integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz
@ -758,7 +758,7 @@ subroutine psb_s_simple_glob_transpose(ain,aout,desc_a,info)
! that the same DESC_A works for both A and A^T, which ! that the same DESC_A works for both A and A^T, which
! essentially means that A has a symmetric pattern. ! essentially means that A has a symmetric pattern.
! !
type(psb_ls_coo_sparse_mat) :: tmpc1, tmpc2 type(psb_s_coo_sparse_mat) :: tmpc1, tmpc2
integer(psb_ipk_) :: nz1, nz2, nzh, nz integer(psb_ipk_) :: nz1, nz2, nzh, nz
integer(psb_ipk_) :: ictxt, me, np integer(psb_ipk_) :: ictxt, me, np
integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz

@ -707,7 +707,7 @@ subroutine psb_z_simple_glob_transpose_ip(ain,desc_a,info)
! that the same DESC_A works for both A and A^T, which ! that the same DESC_A works for both A and A^T, which
! essentially means that A has a symmetric pattern. ! essentially means that A has a symmetric pattern.
! !
type(psb_lz_coo_sparse_mat) :: tmpc1, tmpc2 type(psb_z_coo_sparse_mat) :: tmpc1, tmpc2
integer(psb_ipk_) :: nz1, nz2, nzh, nz integer(psb_ipk_) :: nz1, nz2, nzh, nz
integer(psb_ipk_) :: ictxt, me, np integer(psb_ipk_) :: ictxt, me, np
integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz
@ -758,7 +758,7 @@ subroutine psb_z_simple_glob_transpose(ain,aout,desc_a,info)
! that the same DESC_A works for both A and A^T, which ! that the same DESC_A works for both A and A^T, which
! essentially means that A has a symmetric pattern. ! essentially means that A has a symmetric pattern.
! !
type(psb_lz_coo_sparse_mat) :: tmpc1, tmpc2 type(psb_z_coo_sparse_mat) :: tmpc1, tmpc2
integer(psb_ipk_) :: nz1, nz2, nzh, nz integer(psb_ipk_) :: nz1, nz2, nzh, nz
integer(psb_ipk_) :: ictxt, me, np integer(psb_ipk_) :: ictxt, me, np
integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz integer(psb_lpk_) :: i, j, k, nrow, ncol, nlz

Loading…
Cancel
Save