base/tools/psb_ccdbldext.F90
 base/tools/psb_dcdbldext.F90
 base/tools/psb_scdbldext.F90
 base/tools/psb_zcdbldext.F90

Revert: in tight loops on small sizes it's better to call directly
into %indxmap%
psblas-testmv
Salvatore Filippone 11 years ago
parent fb8e8ae927
commit 81b9285e85

@ -254,7 +254,7 @@ Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info, extype)
Do j=0,n_elem_recv-1 Do j=0,n_elem_recv-1
idx = desc_a%ovrlap_index(counter+psb_elem_recv_+j) idx = desc_a%ovrlap_index(counter+psb_elem_recv_+j)
call desc_ov%l2g(idx,gidx,info) call desc_ov%indxmap%l2g(idx,gidx,info)
If (gidx < 0) then If (gidx < 0) then
info=-3 info=-3
call psb_errpush(info,name) call psb_errpush(info,name)
@ -611,7 +611,6 @@ Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info, extype)
End If End If
End Do End Do
select case(extype_) select case(extype_)
case(psb_ovt_xhal_) case(psb_ovt_xhal_)
! !

@ -254,7 +254,7 @@ Subroutine psb_dcdbldext(a,desc_a,novr,desc_ov,info, extype)
Do j=0,n_elem_recv-1 Do j=0,n_elem_recv-1
idx = desc_a%ovrlap_index(counter+psb_elem_recv_+j) idx = desc_a%ovrlap_index(counter+psb_elem_recv_+j)
call desc_ov%l2g(idx,gidx,info) call desc_ov%indxmap%l2g(idx,gidx,info)
If (gidx < 0) then If (gidx < 0) then
info=-3 info=-3
call psb_errpush(info,name) call psb_errpush(info,name)
@ -611,7 +611,6 @@ Subroutine psb_dcdbldext(a,desc_a,novr,desc_ov,info, extype)
End If End If
End Do End Do
select case(extype_) select case(extype_)
case(psb_ovt_xhal_) case(psb_ovt_xhal_)
! !

@ -254,7 +254,7 @@ Subroutine psb_scdbldext(a,desc_a,novr,desc_ov,info, extype)
Do j=0,n_elem_recv-1 Do j=0,n_elem_recv-1
idx = desc_a%ovrlap_index(counter+psb_elem_recv_+j) idx = desc_a%ovrlap_index(counter+psb_elem_recv_+j)
call desc_ov%l2g(idx,gidx,info) call desc_ov%indxmap%l2g(idx,gidx,info)
If (gidx < 0) then If (gidx < 0) then
info=-3 info=-3
call psb_errpush(info,name) call psb_errpush(info,name)
@ -611,7 +611,6 @@ Subroutine psb_scdbldext(a,desc_a,novr,desc_ov,info, extype)
End If End If
End Do End Do
select case(extype_) select case(extype_)
case(psb_ovt_xhal_) case(psb_ovt_xhal_)
! !

@ -254,7 +254,7 @@ Subroutine psb_zcdbldext(a,desc_a,novr,desc_ov,info, extype)
Do j=0,n_elem_recv-1 Do j=0,n_elem_recv-1
idx = desc_a%ovrlap_index(counter+psb_elem_recv_+j) idx = desc_a%ovrlap_index(counter+psb_elem_recv_+j)
call desc_ov%l2g(idx,gidx,info) call desc_ov%indxmap%l2g(idx,gidx,info)
If (gidx < 0) then If (gidx < 0) then
info=-3 info=-3
call psb_errpush(info,name) call psb_errpush(info,name)
@ -611,7 +611,6 @@ Subroutine psb_zcdbldext(a,desc_a,novr,desc_ov,info, extype)
End If End If
End Do End Do
select case(extype_) select case(extype_)
case(psb_ovt_xhal_) case(psb_ovt_xhal_)
! !

Loading…
Cancel
Save