psblas3-cbind:

cbind/prec/Makefile cbind/prec/psb_c_dprec.h cbind/prec/psb_cprec_cbind_mod.f90 cbind/prec/psb_prec_cbind_mod.f90 cbind/prec/psb_sprec_cbind_mod.f90 cbind/prec/psb_zprec_cbind_mod.f90 Start on PREC.
2026-10-06 22:55:08 +00:00 · 2017-05-04 20:23:40 +00:00
parent 39c90a8362
commit fbc41dace7
6 changed files with 368 additions and 6 deletions
@@ -7,11 +7,16 @@ HERE=..
 FINCLUDES=$(FMFLAG). $(FMFLAG)$(HERE) $(FMFLAG)$(INCLUDEDIR) 
 CINCLUDES=-I. -I$(HERE) -I$(INCLUDEDIR)

-OBJS=psb_prec_cbind_mod.o psb_dprec_cbind_mod.o psb_c_dprec.o
+OBJS=psb_prec_cbind_mod.o \
+psb_sprec_cbind_mod.o psb_dprec_cbind_mod.o psb_cprec_cbind_mod.o psb_zprec_cbind_mod.o \
+psb_c_dprec.o
 CMOD=psb_prec_cbind.h psb_c_dprec.h


-LIBMOD=psb_prec_cbind_mod$(.mod)  psb_prec_cbind_mod$(.mod) 
+LIBMOD=psb_prec_cbind_mod$(.mod)  \
+psb_sprec_cbind_mod$(.mod) psb_dprec_cbind_mod$(.mod) \
+psb_cprec_cbind_mod$(.mod) psb_zprec_cbind_mod$(.mod) 
+
 LOCAL_MODS=$(LIBMOD)
 LIBNAME=$(CPRECLIBNAME)	

@@ -22,7 +27,7 @@ lib:  $(OBJS) $(CMOD)
 	/bin/cp -p $(HERE)/$(LIBNAME) $(LIBDIR)
 	/bin/cp -p $(LIBMOD) $(CMOD) $(INCLUDEDIR)

-psb_prec_cbind_mod.o: psb_dprec_cbind_mod.o
+psb_prec_cbind_mod.o: psb_sprec_cbind_mod.o psb_dprec_cbind_mod.o psb_cprec_cbind_mod.o psb_zprec_cbind_mod.o
 veryclean: clean
 	/bin/rm -f $(HERE)/$(LIBNAME)

@@ -14,9 +14,9 @@ typedef struct PSB_C_DPREC {

 psb_c_dprec* psb_c_new_dprec();

-int  psb_c_dprecinit(psb_c_dprec *ph, const char *ptype);
-int  psb_c_dprecbld(psb_c_dspmat *ah, psb_c_descriptor *cdh, psb_c_dprec *ph);
-int  psb_c_dprecfree(psb_c_dprec *ph);
+psb_i_t  psb_c_dprecinit(psb_c_dprec *ph, const char *ptype);
+psb_i_t  psb_c_dprecbld(psb_c_dspmat *ah, psb_c_descriptor *cdh, psb_c_dprec *ph);
+psb_i_t  psb_c_dprecfree(psb_c_dprec *ph);
 #ifdef __cplusplus
 }
 #endif
@@ -0,0 +1,118 @@
+module psb_cprec_cbind_mod
+
+  use iso_c_binding
+  use psb_prec_mod, only : psb_cprec_type
+  use psb_objhandle_mod
+  use psb_base_string_cbind_mod
+
+  type, bind(c) :: psb_c_cprec
+    type(c_ptr) :: item = c_null_ptr
+  end type psb_c_cprec
+  
+  
+contains 
+
+
+  function  psb_c_cprecinit(ph,ptype) bind(c) result(res)
+    use psb_base_mod
+    use psb_prec_mod
+    use psb_base_string_cbind_mod
+    implicit none 
+    integer(c_int)          :: res
+    type(psb_c_cprec) :: ph
+    character(c_char)       :: ptype(*)
+    type(psb_cprec_type), pointer :: precp
+    integer                 :: info
+    character(len=80)       :: fptype
+
+    res = -1
+    if (c_associated(ph%item)) then 
+      return 
+    end if
+
+    allocate(precp,stat=info)
+    if (info /= 0) return
+    ph%item = c_loc(precp)
+
+    call stringc2f(ptype,fptype)
+    
+    call psb_precinit(precp,fptype,info) 
+    
+    res = min(0,info)
+    return
+  end function psb_c_cprecinit
+    
+
+
+  function  psb_c_cprecbld(ah,cdh,ph) bind(c) result(res)
+    use psb_base_mod
+    use psb_prec_mod
+    use psb_objhandle_mod
+    use psb_base_string_cbind_mod
+    implicit none 
+    
+    integer(c_int) :: res
+    type(psb_c_cspmat) :: ah
+    type(psb_c_cprec) :: ph
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer :: descp
+    type(psb_cspmat_type), pointer :: ap
+    type(psb_cprec_type), pointer :: precp
+
+    integer               :: info
+
+    res = -1
+!!$    write(*,*) 'Entry:   ', psb_c_cd_get_local_rows(cdh)
+    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
+    if (c_associated(ph%item)) then 
+      call c_f_pointer(ph%item,precp)
+    else
+      return 
+    end if
+
+    call psb_precbld(ap,descp, precp, info)
+
+    res = min(info,0)
+    
+  end function psb_c_cprecbld
+
+
+  function  psb_c_cprecfree(ph) bind(c) result(res)
+    use psb_base_mod
+    use psb_prec_mod
+    use psb_objhandle_mod
+    use psb_base_string_cbind_mod
+    implicit none 
+    
+    integer(c_int) :: res
+    type(psb_c_cprec) :: ph
+
+    type(psb_cprec_type), pointer :: precp
+
+    integer               :: info
+
+    res = -1
+    if (c_associated(ph%item)) then 
+      call c_f_pointer(ph%item,precp)
+    else
+      return 
+    end if
+    
+    call psb_precfree(precp, info)
+
+    res = min(info,0)
+    
+  end function psb_c_cprecfree
+
+
+end module psb_cprec_cbind_mod
@@ -1,3 +1,6 @@
 module psb_prec_cbind_mod
+  use psb_sprec_cbind_mod
  use psb_dprec_cbind_mod
+  use psb_cprec_cbind_mod
+  use psb_zprec_cbind_mod
 end module psb_prec_cbind_mod
@@ -0,0 +1,118 @@
+module psb_sprec_cbind_mod
+
+  use iso_c_binding
+  use psb_prec_mod, only : psb_sprec_type
+  use psb_objhandle_mod
+  use psb_base_string_cbind_mod
+
+  type, bind(c) :: psb_c_sprec
+    type(c_ptr) :: item = c_null_ptr
+  end type psb_c_sprec
+  
+  
+contains 
+
+
+  function  psb_c_sprecinit(ph,ptype) bind(c) result(res)
+    use psb_base_mod
+    use psb_prec_mod
+    use psb_base_string_cbind_mod
+    implicit none 
+    integer(c_int)          :: res
+    type(psb_c_sprec) :: ph
+    character(c_char)       :: ptype(*)
+    type(psb_sprec_type), pointer :: precp
+    integer                 :: info
+    character(len=80)       :: fptype
+
+    res = -1
+    if (c_associated(ph%item)) then 
+      return 
+    end if
+
+    allocate(precp,stat=info)
+    if (info /= 0) return
+    ph%item = c_loc(precp)
+
+    call stringc2f(ptype,fptype)
+    
+    call psb_precinit(precp,fptype,info) 
+    
+    res = min(0,info)
+    return
+  end function psb_c_sprecinit
+    
+
+
+  function  psb_c_sprecbld(ah,cdh,ph) bind(c) result(res)
+    use psb_base_mod
+    use psb_prec_mod
+    use psb_objhandle_mod
+    use psb_base_string_cbind_mod
+    implicit none 
+    
+    integer(c_int) :: res
+    type(psb_c_sspmat) :: ah
+    type(psb_c_sprec) :: ph
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer :: descp
+    type(psb_sspmat_type), pointer :: ap
+    type(psb_sprec_type), pointer :: precp
+
+    integer               :: info
+
+    res = -1
+!!$    write(*,*) 'Entry:   ', psb_c_cd_get_local_rows(cdh)
+    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
+    if (c_associated(ph%item)) then 
+      call c_f_pointer(ph%item,precp)
+    else
+      return 
+    end if
+
+    call psb_precbld(ap,descp, precp, info)
+
+    res = min(info,0)
+    
+  end function psb_c_sprecbld
+
+
+  function  psb_c_sprecfree(ph) bind(c) result(res)
+    use psb_base_mod
+    use psb_prec_mod
+    use psb_objhandle_mod
+    use psb_base_string_cbind_mod
+    implicit none 
+    
+    integer(c_int) :: res
+    type(psb_c_sprec) :: ph
+
+    type(psb_sprec_type), pointer :: precp
+
+    integer               :: info
+
+    res = -1
+    if (c_associated(ph%item)) then 
+      call c_f_pointer(ph%item,precp)
+    else
+      return 
+    end if
+    
+    call psb_precfree(precp, info)
+
+    res = min(info,0)
+    
+  end function psb_c_sprecfree
+
+
+end module psb_sprec_cbind_mod
@@ -0,0 +1,118 @@
+module psb_zprec_cbind_mod
+
+  use iso_c_binding
+  use psb_prec_mod, only : psb_zprec_type
+  use psb_objhandle_mod
+  use psb_base_string_cbind_mod
+
+  type, bind(c) :: psb_c_zprec
+    type(c_ptr) :: item = c_null_ptr
+  end type psb_c_zprec
+  
+  
+contains 
+
+
+  function  psb_c_zprecinit(ph,ptype) bind(c) result(res)
+    use psb_base_mod
+    use psb_prec_mod
+    use psb_base_string_cbind_mod
+    implicit none 
+    integer(c_int)          :: res
+    type(psb_c_zprec) :: ph
+    character(c_char)       :: ptype(*)
+    type(psb_zprec_type), pointer :: precp
+    integer                 :: info
+    character(len=80)       :: fptype
+
+    res = -1
+    if (c_associated(ph%item)) then 
+      return 
+    end if
+
+    allocate(precp,stat=info)
+    if (info /= 0) return
+    ph%item = c_loc(precp)
+
+    call stringc2f(ptype,fptype)
+    
+    call psb_precinit(precp,fptype,info) 
+    
+    res = min(0,info)
+    return
+  end function psb_c_zprecinit
+    
+
+
+  function  psb_c_zprecbld(ah,cdh,ph) bind(c) result(res)
+    use psb_base_mod
+    use psb_prec_mod
+    use psb_objhandle_mod
+    use psb_base_string_cbind_mod
+    implicit none 
+    
+    integer(c_int) :: res
+    type(psb_c_zspmat) :: ah
+    type(psb_c_zprec) :: ph
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer :: descp
+    type(psb_zspmat_type), pointer :: ap
+    type(psb_zprec_type), pointer :: precp
+
+    integer               :: info
+
+    res = -1
+!!$    write(*,*) 'Entry:   ', psb_c_cd_get_local_rows(cdh)
+    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
+    if (c_associated(ph%item)) then 
+      call c_f_pointer(ph%item,precp)
+    else
+      return 
+    end if
+
+    call psb_precbld(ap,descp, precp, info)
+
+    res = min(info,0)
+    
+  end function psb_c_zprecbld
+
+
+  function  psb_c_zprecfree(ph) bind(c) result(res)
+    use psb_base_mod
+    use psb_prec_mod
+    use psb_objhandle_mod
+    use psb_base_string_cbind_mod
+    implicit none 
+    
+    integer(c_int) :: res
+    type(psb_c_zprec) :: ph
+
+    type(psb_zprec_type), pointer :: precp
+
+    integer               :: info
+
+    res = -1
+    if (c_associated(ph%item)) then 
+      call c_f_pointer(ph%item,precp)
+    else
+      return 
+    end if
+    
+    call psb_precfree(precp, info)
+
+    res = min(info,0)
+    
+  end function psb_c_zprecfree
+
+
+end module psb_zprec_cbind_mod