You cannot select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
psblas3/cbind/util/psb_z_util_cbind_mod.f90

89 lines
2.1 KiB
Fortran

module psb_zutil_cbind_mod
use iso_c_binding
use psb_util_mod
use psb_base_mod
use psb_objhandle_mod
use psb_base_string_cbind_mod
contains
function psb_c_zmm_mat_write(ah,matrixtitle,filename) bind(c) result(res)
use psb_base_mod
use psb_util_mod
use psb_base_string_cbind_mod
implicit none
integer(psb_c_ipk_) :: res
type(psb_c_zspmat) :: ah
character(c_char) :: matrixtitle(*)
character(c_char) :: filename(*)
type(psb_zspmat_type), pointer :: ap
character(len=1024) :: mtitle
character(len=1024) :: fname
integer(psb_c_ipk_) :: info
res = -1
if (c_associated(ah%item)) then
call c_f_pointer(ah%item,ap)
else
return
end if
call stringc2f(matrixtitle,mtitle)
call stringc2f(filename,fname)
call mm_mat_write(ap,mtitle,info,filename=fname)
res = info
end function psb_c_zmm_mat_write
function psb_c_zglobal_mat_write(ah,cdh) bind(c) result(res)
use psb_base_mod
use psb_util_mod
use psb_base_string_cbind_mod
implicit none
integer(psb_c_ipk_) :: res
type(psb_c_zspmat) :: ah
type(psb_c_descriptor) :: cdh
type(psb_zspmat_type), pointer :: ap
type(psb_desc_type), pointer :: descp
! Local variables
type(psb_zspmat_type) :: aglobal
integer(psb_ipk_) :: info, iam, np
type(psb_ctxt_type) :: ctxt
character(len=40) :: matrixname
res = -1
if (c_associated(cdh%item)) then
call c_f_pointer(cdh%item,descp)
else
return
end if
if (c_associated(ah%item)) then
call c_f_pointer(ah%item,ap)
else
return
end if
ctxt = descp%get_ctxt()
call psb_info(ctxt,iam,np)
call psb_gather(aglobal,ap,descp,info)
if (iam == psb_root_) then
write(matrixname,'("A-np-",I1,".mtx")') np
call mm_mat_write(aglobal,"Global matrix",info,filename=trim(matrixname))
end if
call psb_spfree(aglobal,descp,info)
res = info
end function psb_c_zglobal_mat_write
end module psb_zutil_cbind_mod