base/tools/psb_c_map.f90
 base/tools/psb_d_map.f90
 base/tools/psb_s_map.f90
 base/tools/psb_z_map.f90

Silly bug with iam/np
trunk
Salvatore Filippone 8 years ago
parent 85da2dbd27
commit 854c3d8418

@ -306,7 +306,7 @@ subroutine psb_c_map_Y2X_vect(alpha,x,beta,y,map,info,work,vtx,vty)
type(psb_c_vect_type),pointer :: ptx, pty
complex(psb_spk_), allocatable :: xta(:), yta(:)
integer(psb_ipk_) :: i, j, nr1, nc1,nr2, nc2,&
& map_kind, nr, ictxt
& map_kind, nr, ictxt, iam, np
character(len=20), parameter :: name='psb_map_Y2Xv'
info = psb_success_
@ -322,6 +322,7 @@ subroutine psb_c_map_Y2X_vect(alpha,x,beta,y,map,info,work,vtx,vty)
case(psb_map_aggr_)
ictxt = map%p_desc_X%get_context()
call psb_info(ictxt,iam,np)
nr2 = map%p_desc_X%get_global_rows()
nc2 = map%p_desc_X%get_local_cols()
if (present(vty)) then
@ -348,6 +349,7 @@ subroutine psb_c_map_Y2X_vect(alpha,x,beta,y,map,info,work,vtx,vty)
case(psb_map_gen_linear_)
ictxt = map%desc_X%get_context()
call psb_info(ictxt,iam,np)
nr1 = map%desc_Y%get_local_rows()
nc1 = map%desc_Y%get_local_cols()
nr2 = map%desc_X%get_global_rows()

@ -306,7 +306,7 @@ subroutine psb_d_map_Y2X_vect(alpha,x,beta,y,map,info,work,vtx,vty)
type(psb_d_vect_type),pointer :: ptx, pty
real(psb_dpk_), allocatable :: xta(:), yta(:)
integer(psb_ipk_) :: i, j, nr1, nc1,nr2, nc2,&
& map_kind, nr, ictxt
& map_kind, nr, ictxt, iam, np
character(len=20), parameter :: name='psb_map_Y2Xv'
info = psb_success_
@ -322,6 +322,7 @@ subroutine psb_d_map_Y2X_vect(alpha,x,beta,y,map,info,work,vtx,vty)
case(psb_map_aggr_)
ictxt = map%p_desc_X%get_context()
call psb_info(ictxt,iam,np)
nr2 = map%p_desc_X%get_global_rows()
nc2 = map%p_desc_X%get_local_cols()
if (present(vty)) then
@ -348,6 +349,7 @@ subroutine psb_d_map_Y2X_vect(alpha,x,beta,y,map,info,work,vtx,vty)
case(psb_map_gen_linear_)
ictxt = map%desc_X%get_context()
call psb_info(ictxt,iam,np)
nr1 = map%desc_Y%get_local_rows()
nc1 = map%desc_Y%get_local_cols()
nr2 = map%desc_X%get_global_rows()

@ -306,7 +306,7 @@ subroutine psb_s_map_Y2X_vect(alpha,x,beta,y,map,info,work,vtx,vty)
type(psb_s_vect_type),pointer :: ptx, pty
real(psb_spk_), allocatable :: xta(:), yta(:)
integer(psb_ipk_) :: i, j, nr1, nc1,nr2, nc2,&
& map_kind, nr, ictxt
& map_kind, nr, ictxt, iam, np
character(len=20), parameter :: name='psb_map_Y2Xv'
info = psb_success_
@ -322,6 +322,7 @@ subroutine psb_s_map_Y2X_vect(alpha,x,beta,y,map,info,work,vtx,vty)
case(psb_map_aggr_)
ictxt = map%p_desc_X%get_context()
call psb_info(ictxt,iam,np)
nr2 = map%p_desc_X%get_global_rows()
nc2 = map%p_desc_X%get_local_cols()
if (present(vty)) then
@ -348,6 +349,7 @@ subroutine psb_s_map_Y2X_vect(alpha,x,beta,y,map,info,work,vtx,vty)
case(psb_map_gen_linear_)
ictxt = map%desc_X%get_context()
call psb_info(ictxt,iam,np)
nr1 = map%desc_Y%get_local_rows()
nc1 = map%desc_Y%get_local_cols()
nr2 = map%desc_X%get_global_rows()

@ -306,7 +306,7 @@ subroutine psb_z_map_Y2X_vect(alpha,x,beta,y,map,info,work,vtx,vty)
type(psb_z_vect_type),pointer :: ptx, pty
complex(psb_dpk_), allocatable :: xta(:), yta(:)
integer(psb_ipk_) :: i, j, nr1, nc1,nr2, nc2,&
& map_kind, nr, ictxt
& map_kind, nr, ictxt, iam, np
character(len=20), parameter :: name='psb_map_Y2Xv'
info = psb_success_
@ -322,6 +322,7 @@ subroutine psb_z_map_Y2X_vect(alpha,x,beta,y,map,info,work,vtx,vty)
case(psb_map_aggr_)
ictxt = map%p_desc_X%get_context()
call psb_info(ictxt,iam,np)
nr2 = map%p_desc_X%get_global_rows()
nc2 = map%p_desc_X%get_local_cols()
if (present(vty)) then
@ -348,6 +349,7 @@ subroutine psb_z_map_Y2X_vect(alpha,x,beta,y,map,info,work,vtx,vty)
case(psb_map_gen_linear_)
ictxt = map%desc_X%get_context()
call psb_info(ictxt,iam,np)
nr1 = map%desc_Y%get_local_rows()
nc1 = map%desc_Y%get_local_cols()
nr2 = map%desc_X%get_global_rows()

Loading…
Cancel
Save