mirror of
https://github.com/sfilippone/psblas3.git
synced 2026-10-07 07:04:57 +00:00
psblas3:
base/serial/impl/psb_c_base_mat_impl.F90 base/serial/impl/psb_c_mat_impl.F90 base/serial/impl/psb_d_base_mat_impl.F90 base/serial/impl/psb_d_mat_impl.F90 base/serial/impl/psb_s_base_mat_impl.F90 base/serial/impl/psb_s_mat_impl.F90 base/serial/impl/psb_z_base_mat_impl.F90 base/serial/impl/psb_z_mat_impl.F90 Fixed silly bug
This commit is contained in:
@@ -636,6 +636,7 @@ subroutine psb_c_base_tril(a,b,info,&
|
||||
& nzin=nzin)
|
||||
if (info /= psb_success_) goto 9999
|
||||
call b%set_nzeros(nzin+nzout)
|
||||
nzin = nzin+nzout
|
||||
end do
|
||||
call b%fix(info)
|
||||
nzout = b%get_nzeros()
|
||||
@@ -745,6 +746,7 @@ subroutine psb_c_base_triu(a,b,info,&
|
||||
& nzin=nzin)
|
||||
if (info /= psb_success_) goto 9999
|
||||
call b%set_nzeros(nzin+nzout)
|
||||
nzin = nzin+nzout
|
||||
end do
|
||||
call b%fix(info)
|
||||
nzout = b%get_nzeros()
|
||||
|
||||
@@ -1025,7 +1025,7 @@ subroutine psb_c_triu(a,b,info,diag,imin,imax,&
|
||||
allocate(acoo,stat=info)
|
||||
|
||||
if (info == psb_success_) then
|
||||
call a%a%tril(acoo,info,diag,imin,imax,&
|
||||
call a%a%triu(acoo,info,diag,imin,imax,&
|
||||
& jmin,jmax,rscale,cscale)
|
||||
else
|
||||
info = psb_err_alloc_dealloc_
|
||||
|
||||
@@ -636,6 +636,7 @@ subroutine psb_d_base_tril(a,b,info,&
|
||||
& nzin=nzin)
|
||||
if (info /= psb_success_) goto 9999
|
||||
call b%set_nzeros(nzin+nzout)
|
||||
nzin = nzin+nzout
|
||||
end do
|
||||
call b%fix(info)
|
||||
nzout = b%get_nzeros()
|
||||
@@ -745,6 +746,7 @@ subroutine psb_d_base_triu(a,b,info,&
|
||||
& nzin=nzin)
|
||||
if (info /= psb_success_) goto 9999
|
||||
call b%set_nzeros(nzin+nzout)
|
||||
nzin = nzin+nzout
|
||||
end do
|
||||
call b%fix(info)
|
||||
nzout = b%get_nzeros()
|
||||
|
||||
@@ -1025,7 +1025,7 @@ subroutine psb_d_triu(a,b,info,diag,imin,imax,&
|
||||
allocate(acoo,stat=info)
|
||||
|
||||
if (info == psb_success_) then
|
||||
call a%a%tril(acoo,info,diag,imin,imax,&
|
||||
call a%a%triu(acoo,info,diag,imin,imax,&
|
||||
& jmin,jmax,rscale,cscale)
|
||||
else
|
||||
info = psb_err_alloc_dealloc_
|
||||
|
||||
@@ -636,6 +636,7 @@ subroutine psb_s_base_tril(a,b,info,&
|
||||
& nzin=nzin)
|
||||
if (info /= psb_success_) goto 9999
|
||||
call b%set_nzeros(nzin+nzout)
|
||||
nzin = nzin+nzout
|
||||
end do
|
||||
call b%fix(info)
|
||||
nzout = b%get_nzeros()
|
||||
@@ -745,6 +746,7 @@ subroutine psb_s_base_triu(a,b,info,&
|
||||
& nzin=nzin)
|
||||
if (info /= psb_success_) goto 9999
|
||||
call b%set_nzeros(nzin+nzout)
|
||||
nzin = nzin+nzout
|
||||
end do
|
||||
call b%fix(info)
|
||||
nzout = b%get_nzeros()
|
||||
|
||||
@@ -1025,7 +1025,7 @@ subroutine psb_s_triu(a,b,info,diag,imin,imax,&
|
||||
allocate(acoo,stat=info)
|
||||
|
||||
if (info == psb_success_) then
|
||||
call a%a%tril(acoo,info,diag,imin,imax,&
|
||||
call a%a%triu(acoo,info,diag,imin,imax,&
|
||||
& jmin,jmax,rscale,cscale)
|
||||
else
|
||||
info = psb_err_alloc_dealloc_
|
||||
|
||||
@@ -636,6 +636,7 @@ subroutine psb_z_base_tril(a,b,info,&
|
||||
& nzin=nzin)
|
||||
if (info /= psb_success_) goto 9999
|
||||
call b%set_nzeros(nzin+nzout)
|
||||
nzin = nzin+nzout
|
||||
end do
|
||||
call b%fix(info)
|
||||
nzout = b%get_nzeros()
|
||||
@@ -745,6 +746,7 @@ subroutine psb_z_base_triu(a,b,info,&
|
||||
& nzin=nzin)
|
||||
if (info /= psb_success_) goto 9999
|
||||
call b%set_nzeros(nzin+nzout)
|
||||
nzin = nzin+nzout
|
||||
end do
|
||||
call b%fix(info)
|
||||
nzout = b%get_nzeros()
|
||||
|
||||
@@ -1025,7 +1025,7 @@ subroutine psb_z_triu(a,b,info,diag,imin,imax,&
|
||||
allocate(acoo,stat=info)
|
||||
|
||||
if (info == psb_success_) then
|
||||
call a%a%tril(acoo,info,diag,imin,imax,&
|
||||
call a%a%triu(acoo,info,diag,imin,imax,&
|
||||
& jmin,jmax,rscale,cscale)
|
||||
else
|
||||
info = psb_err_alloc_dealloc_
|
||||
|
||||
Reference in New Issue
Block a user