base/serial/impl/psb_c_csr_impl.f90
 base/serial/impl/psb_d_csr_impl.f90
 base/serial/impl/psb_s_csr_impl.f90
 base/serial/impl/psb_z_csr_impl.f90

Fix silly bug in cp_from_coo
psblas3-accel
Salvatore Filippone 10 years ago
parent ca4a919959
commit e1cd09b457

@ -2702,9 +2702,9 @@ subroutine psb_c_cp_csr_from_coo(a,b,info)
! Dirty trick: call move_alloc to have the new data allocated just once.
call psb_safe_ab_cpy(b%ia,itemp,info)
if (info /= psb_success_) call psb_safe_ab_cpy(b%ja,a%ja,info)
if (info /= psb_success_) call psb_safe_ab_cpy(b%val,a%val,info)
if (info /= psb_success_) call psb_realloc(max(nr+1,nc+1),a%irp,info)
if (info == psb_success_) call psb_safe_ab_cpy(b%ja,a%ja,info)
if (info == psb_success_) call psb_safe_ab_cpy(b%val,a%val,info)
if (info == psb_success_) call psb_realloc(max(nr+1,nc+1),a%irp,info)
endif

@ -2702,9 +2702,9 @@ subroutine psb_d_cp_csr_from_coo(a,b,info)
! Dirty trick: call move_alloc to have the new data allocated just once.
call psb_safe_ab_cpy(b%ia,itemp,info)
if (info /= psb_success_) call psb_safe_ab_cpy(b%ja,a%ja,info)
if (info /= psb_success_) call psb_safe_ab_cpy(b%val,a%val,info)
if (info /= psb_success_) call psb_realloc(max(nr+1,nc+1),a%irp,info)
if (info == psb_success_) call psb_safe_ab_cpy(b%ja,a%ja,info)
if (info == psb_success_) call psb_safe_ab_cpy(b%val,a%val,info)
if (info == psb_success_) call psb_realloc(max(nr+1,nc+1),a%irp,info)
endif

@ -2702,9 +2702,9 @@ subroutine psb_s_cp_csr_from_coo(a,b,info)
! Dirty trick: call move_alloc to have the new data allocated just once.
call psb_safe_ab_cpy(b%ia,itemp,info)
if (info /= psb_success_) call psb_safe_ab_cpy(b%ja,a%ja,info)
if (info /= psb_success_) call psb_safe_ab_cpy(b%val,a%val,info)
if (info /= psb_success_) call psb_realloc(max(nr+1,nc+1),a%irp,info)
if (info == psb_success_) call psb_safe_ab_cpy(b%ja,a%ja,info)
if (info == psb_success_) call psb_safe_ab_cpy(b%val,a%val,info)
if (info == psb_success_) call psb_realloc(max(nr+1,nc+1),a%irp,info)
endif

@ -2702,9 +2702,9 @@ subroutine psb_z_cp_csr_from_coo(a,b,info)
! Dirty trick: call move_alloc to have the new data allocated just once.
call psb_safe_ab_cpy(b%ia,itemp,info)
if (info /= psb_success_) call psb_safe_ab_cpy(b%ja,a%ja,info)
if (info /= psb_success_) call psb_safe_ab_cpy(b%val,a%val,info)
if (info /= psb_success_) call psb_realloc(max(nr+1,nc+1),a%irp,info)
if (info == psb_success_) call psb_safe_ab_cpy(b%ja,a%ja,info)
if (info == psb_success_) call psb_safe_ab_cpy(b%val,a%val,info)
if (info == psb_success_) call psb_realloc(max(nr+1,nc+1),a%irp,info)
endif

Loading…
Cancel
Save