|
|
@ -160,7 +160,7 @@ program cf_sample
|
|
|
|
|
|
|
|
|
|
|
|
b_col_glob => aux_b(:,1)
|
|
|
|
b_col_glob => aux_b(:,1)
|
|
|
|
do i=1, m_problem
|
|
|
|
do i=1, m_problem
|
|
|
|
b_col_glob(i) = (1.d0,1.d0)
|
|
|
|
b_col_glob(i) = (1.0,1.0)
|
|
|
|
enddo
|
|
|
|
enddo
|
|
|
|
endif
|
|
|
|
endif
|
|
|
|
call psb_bcast(ictxt,b_col_glob(1:m_problem))
|
|
|
|
call psb_bcast(ictxt,b_col_glob(1:m_problem))
|
|
|
@ -180,7 +180,7 @@ program cf_sample
|
|
|
|
! switch over different partition types
|
|
|
|
! switch over different partition types
|
|
|
|
if (ipart == 0) then
|
|
|
|
if (ipart == 0) then
|
|
|
|
call psb_barrier(ictxt)
|
|
|
|
call psb_barrier(ictxt)
|
|
|
|
if (iam == psb_root_) write(psb_out_unit,'("Partition type: block")')
|
|
|
|
if (iam == psb_root_) write(psb_out_unit,'("Partition type: block vector")')
|
|
|
|
allocate(ivg(m_problem),ipv(np))
|
|
|
|
allocate(ivg(m_problem),ipv(np))
|
|
|
|
do i=1,m_problem
|
|
|
|
do i=1,m_problem
|
|
|
|
call part_block(i,m_problem,np,ipv,nv)
|
|
|
|
call part_block(i,m_problem,np,ipv,nv)
|
|
|
@ -191,7 +191,7 @@ program cf_sample
|
|
|
|
|
|
|
|
|
|
|
|
else if (ipart == 2) then
|
|
|
|
else if (ipart == 2) then
|
|
|
|
if (iam == psb_root_) then
|
|
|
|
if (iam == psb_root_) then
|
|
|
|
write(psb_out_unit,'("Partition type: graph")')
|
|
|
|
write(psb_out_unit,'("Partition type: graph vector")')
|
|
|
|
write(psb_out_unit,'(" ")')
|
|
|
|
write(psb_out_unit,'(" ")')
|
|
|
|
! write(psb_err_unit,'("Build type: graph")')
|
|
|
|
! write(psb_err_unit,'("Build type: graph")')
|
|
|
|
call build_mtpart(aux_a,np)
|
|
|
|
call build_mtpart(aux_a,np)
|
|
|
@ -204,7 +204,7 @@ program cf_sample
|
|
|
|
& desc_a,b_col_glob,b_col,info,fmt=afmt,v=ivg)
|
|
|
|
& desc_a,b_col_glob,b_col,info,fmt=afmt,v=ivg)
|
|
|
|
|
|
|
|
|
|
|
|
else
|
|
|
|
else
|
|
|
|
if (iam == psb_root_) write(psb_out_unit,'("Partition type: block")')
|
|
|
|
if (iam == psb_root_) write(psb_out_unit,'("Partition type: block subroutine")')
|
|
|
|
call psb_matdist(aux_a, a, ictxt, &
|
|
|
|
call psb_matdist(aux_a, a, ictxt, &
|
|
|
|
& desc_a,b_col_glob,b_col,info,fmt=afmt,parts=part_block)
|
|
|
|
& desc_a,b_col_glob,b_col,info,fmt=afmt,parts=part_block)
|
|
|
|
end if
|
|
|
|
end if
|
|
|
|