From 58198de6faa85fadf1b6189f403ec0c5b8b81272 Mon Sep 17 00:00:00 2001
From: Cirdans-Home <f.durastante@na.iac.cnr.it>
Date: Mon, 6 Apr 2020 15:50:49 +0200
Subject: [PATCH] Added c interfaces to set_asb(), set_upd(), set_bld() for
 sparse matrices

---
 cbind/base/psb_c_cbase.h              |  3 +
 cbind/base/psb_c_dbase.h              |  3 +
 cbind/base/psb_c_psblas_cbind_mod.f90 | 85 +++++++++++++++++++++++++++
 cbind/base/psb_c_sbase.h              |  3 +
 cbind/base/psb_c_zbase.h              |  3 +
 cbind/base/psb_d_psblas_cbind_mod.f90 | 85 +++++++++++++++++++++++++++
 cbind/base/psb_s_psblas_cbind_mod.f90 | 85 +++++++++++++++++++++++++++
 cbind/base/psb_z_psblas_cbind_mod.f90 | 85 +++++++++++++++++++++++++++
 8 files changed, 352 insertions(+)

diff --git a/cbind/base/psb_c_cbase.h b/cbind/base/psb_c_cbase.h
index dbc81b42..452b74ef 100644
--- a/cbind/base/psb_c_cbase.h
+++ b/cbind/base/psb_c_cbase.h
@@ -44,6 +44,9 @@ psb_l_t    psb_c_cnnz(psb_c_cspmat *mh,psb_c_descriptor *cdh);
 bool    	 psb_c_cis_matupd(psb_c_cspmat *mh,psb_c_descriptor *cdh);
 bool    	 psb_c_cis_matasb(psb_c_cspmat *mh,psb_c_descriptor *cdh);
 bool    	 psb_c_cis_matbld(psb_c_cspmat *mh,psb_c_descriptor *cdh);
+psb_i_t    psb_c_cset_matupd(psb_c_cspmat *mh,psb_c_descriptor *cdh);
+psb_i_t    psb_c_cset_matasb(psb_c_cspmat *mh,psb_c_descriptor *cdh);
+psb_i_t    psb_c_cset_matbld(psb_c_cspmat *mh,psb_c_descriptor *cdh);
 
 /* psb_i_t    psb_c_cspasb_opt(psb_c_cspmat *mh, psb_c_descriptor *cdh,  */
 /* 			const char *afmt, psb_i_t upd, psb_i_t dupl); */
diff --git a/cbind/base/psb_c_dbase.h b/cbind/base/psb_c_dbase.h
index d510f0fa..0705b6ba 100644
--- a/cbind/base/psb_c_dbase.h
+++ b/cbind/base/psb_c_dbase.h
@@ -44,6 +44,9 @@ psb_l_t    psb_c_dnnz(psb_c_dspmat *mh,psb_c_descriptor *cdh);
 bool    	 psb_c_dis_matupd(psb_c_dspmat *mh,psb_c_descriptor *cdh);
 bool    	 psb_c_dis_matasb(psb_c_dspmat *mh,psb_c_descriptor *cdh);
 bool    	 psb_c_dis_matbld(psb_c_dspmat *mh,psb_c_descriptor *cdh);
+psb_i_t    psb_c_dset_matupd(psb_c_dspmat *mh,psb_c_descriptor *cdh);
+psb_i_t    psb_c_dset_matasb(psb_c_dspmat *mh,psb_c_descriptor *cdh);
+psb_i_t    psb_c_dset_matbld(psb_c_dspmat *mh,psb_c_descriptor *cdh);
 
 /* psb_i_t    psb_c_dspasb_opt(psb_c_dspmat *mh, psb_c_descriptor *cdh,  */
 /* 			const char *afmt, psb_i_t upd, psb_i_t dupl); */
diff --git a/cbind/base/psb_c_psblas_cbind_mod.f90 b/cbind/base/psb_c_psblas_cbind_mod.f90
index 6bfc9d26..c68626f5 100644
--- a/cbind/base/psb_c_psblas_cbind_mod.f90
+++ b/cbind/base/psb_c_psblas_cbind_mod.f90
@@ -972,4 +972,89 @@ contains
     res = psb_is_matbld(ap,descp,info)
   end function
 
+  function psb_c_cset_matupd(ah,cdh) bind(c) result(res)
+    implicit none
+    integer(psb_c_ipk_) :: res
+
+    type(psb_c_cspmat)   :: ah
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer      :: descp
+    type(psb_cspmat_type), pointer  :: ap
+
+    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
+
+    call ap%set_upd()
+
+    res = psb_success_
+  end function
+
+  function psb_c_cset_matasb(ah,cdh) bind(c) result(res)
+    implicit none
+    integer(psb_c_ipk_)    :: res
+
+    type(psb_c_cspmat)   :: ah
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer      :: descp
+    type(psb_cspmat_type), pointer  :: ap
+
+    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
+
+    call ap%set_asb()
+
+    res = psb_success_
+
+  end function
+
+  function psb_c_cset_matbld(ah,cdh) bind(c) result(res)
+    implicit none
+    integer(psb_c_ipk_)    :: res
+
+    type(psb_c_cspmat)   :: ah
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer      :: descp
+    type(psb_cspmat_type), pointer  :: ap
+
+    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
+
+    call ap%set_bld()
+
+    res = psb_success_
+  end function
+
 end module psb_c_psblas_cbind_mod
diff --git a/cbind/base/psb_c_sbase.h b/cbind/base/psb_c_sbase.h
index 5e89d14e..cc80bee3 100644
--- a/cbind/base/psb_c_sbase.h
+++ b/cbind/base/psb_c_sbase.h
@@ -44,6 +44,9 @@ psb_l_t    psb_c_snnz(psb_c_sspmat *mh,psb_c_descriptor *cdh);
 bool    	 psb_c_sis_matupd(psb_c_sspmat *mh,psb_c_descriptor *cdh);
 bool    	 psb_c_sis_matasb(psb_c_sspmat *mh,psb_c_descriptor *cdh);
 bool    	 psb_c_sis_matbld(psb_c_sspmat *mh,psb_c_descriptor *cdh);
+psb_i_t    psb_c_sset_matupd(psb_c_sspmat *mh,psb_c_descriptor *cdh);
+psb_i_t    psb_c_sset_matasb(psb_c_sspmat *mh,psb_c_descriptor *cdh);
+psb_i_t    psb_c_sset_matbld(psb_c_sspmat *mh,psb_c_descriptor *cdh);
 
 /* psb_i_t    psb_c_sspasb_opt(psb_c_sspmat *mh, psb_c_descriptor *cdh,  */
 /* 			const char *afmt, psb_i_t upd, psb_i_t dupl); */
diff --git a/cbind/base/psb_c_zbase.h b/cbind/base/psb_c_zbase.h
index c245fd83..38056aff 100644
--- a/cbind/base/psb_c_zbase.h
+++ b/cbind/base/psb_c_zbase.h
@@ -44,6 +44,9 @@ psb_l_t    psb_c_znnz(psb_c_zspmat *mh,psb_c_descriptor *cdh);
 bool    	 psb_c_zis_matupd(psb_c_zspmat *mh,psb_c_descriptor *cdh);
 bool    	 psb_c_zis_matasb(psb_c_zspmat *mh,psb_c_descriptor *cdh);
 bool    	 psb_c_zis_matbld(psb_c_zspmat *mh,psb_c_descriptor *cdh);
+psb_i_t    psb_c_zset_matupd(psb_c_zspmat *mh,psb_c_descriptor *cdh);
+psb_i_t    psb_c_zset_matasb(psb_c_zspmat *mh,psb_c_descriptor *cdh);
+psb_i_t    psb_c_zset_matbld(psb_c_zspmat *mh,psb_c_descriptor *cdh);
 
 /* psb_i_t    psb_c_zspasb_opt(psb_c_zspmat *mh, psb_c_descriptor *cdh,  */
 /* 			const char *afmt, psb_i_t upd, psb_i_t dupl); */
diff --git a/cbind/base/psb_d_psblas_cbind_mod.f90 b/cbind/base/psb_d_psblas_cbind_mod.f90
index 41257e29..50eca563 100644
--- a/cbind/base/psb_d_psblas_cbind_mod.f90
+++ b/cbind/base/psb_d_psblas_cbind_mod.f90
@@ -1073,4 +1073,89 @@ contains
     res = psb_is_matbld(ap,descp,info)
   end function
 
+  function psb_c_dset_matupd(ah,cdh) bind(c) result(res)
+    implicit none
+    integer(psb_c_ipk_) :: res
+
+    type(psb_c_dspmat)   :: ah
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer      :: descp
+    type(psb_dspmat_type), pointer  :: ap
+
+    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
+
+    call ap%set_upd()
+
+    res = psb_success_
+  end function
+
+  function psb_c_dset_matasb(ah,cdh) bind(c) result(res)
+    implicit none
+    integer(psb_c_ipk_)    :: res
+
+    type(psb_c_dspmat)   :: ah
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer      :: descp
+    type(psb_dspmat_type), pointer  :: ap
+
+    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
+
+    call ap%set_asb()
+
+    res = psb_success_
+
+  end function
+
+  function psb_c_dset_matbld(ah,cdh) bind(c) result(res)
+    implicit none
+    integer(psb_c_ipk_)    :: res
+
+    type(psb_c_dspmat)   :: ah
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer      :: descp
+    type(psb_dspmat_type), pointer  :: ap
+
+    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
+
+    call ap%set_bld()
+
+    res = psb_success_
+  end function
+
 end module psb_d_psblas_cbind_mod
diff --git a/cbind/base/psb_s_psblas_cbind_mod.f90 b/cbind/base/psb_s_psblas_cbind_mod.f90
index 5b25f77a..dbc5de0f 100644
--- a/cbind/base/psb_s_psblas_cbind_mod.f90
+++ b/cbind/base/psb_s_psblas_cbind_mod.f90
@@ -1073,4 +1073,89 @@ contains
     res = psb_is_matbld(ap,descp,info)
   end function
 
+  function psb_c_sset_matupd(ah,cdh) bind(c) result(res)
+    implicit none
+    integer(psb_c_ipk_) :: res
+
+    type(psb_c_sspmat)   :: ah
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer      :: descp
+    type(psb_sspmat_type), pointer  :: ap
+
+    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
+
+    call ap%set_upd()
+
+    res = psb_success_
+  end function
+
+  function psb_c_sset_matasb(ah,cdh) bind(c) result(res)
+    implicit none
+    integer(psb_c_ipk_)    :: res
+
+    type(psb_c_sspmat)   :: ah
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer      :: descp
+    type(psb_sspmat_type), pointer  :: ap
+
+    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
+
+    call ap%set_asb()
+
+    res = psb_success_
+
+  end function
+
+  function psb_c_sset_matbld(ah,cdh) bind(c) result(res)
+    implicit none
+    integer(psb_c_ipk_)    :: res
+
+    type(psb_c_sspmat)   :: ah
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer      :: descp
+    type(psb_sspmat_type), pointer  :: ap
+
+    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
+
+    call ap%set_bld()
+
+    res = psb_success_
+  end function
+
 end module psb_s_psblas_cbind_mod
diff --git a/cbind/base/psb_z_psblas_cbind_mod.f90 b/cbind/base/psb_z_psblas_cbind_mod.f90
index 7986b459..73a55a52 100644
--- a/cbind/base/psb_z_psblas_cbind_mod.f90
+++ b/cbind/base/psb_z_psblas_cbind_mod.f90
@@ -972,4 +972,89 @@ contains
     res = psb_is_matbld(ap,descp,info)
   end function
 
+  function psb_c_zset_matupd(ah,cdh) bind(c) result(res)
+    implicit none
+    integer(psb_c_ipk_) :: res
+
+    type(psb_c_zspmat)   :: ah
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer      :: descp
+    type(psb_zspmat_type), pointer  :: ap
+
+    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
+
+    call ap%set_upd()
+
+    res = psb_success_
+  end function
+
+  function psb_c_zset_matasb(ah,cdh) bind(c) result(res)
+    implicit none
+    integer(psb_c_ipk_)    :: res
+
+    type(psb_c_zspmat)   :: ah
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer      :: descp
+    type(psb_zspmat_type), pointer  :: ap
+
+    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
+
+    call ap%set_asb()
+
+    res = psb_success_
+
+  end function
+
+  function psb_c_zset_matbld(ah,cdh) bind(c) result(res)
+    implicit none
+    integer(psb_c_ipk_)    :: res
+
+    type(psb_c_zspmat)   :: ah
+    type(psb_c_descriptor) :: cdh
+
+    type(psb_desc_type), pointer      :: descp
+    type(psb_zspmat_type), pointer  :: ap
+
+    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
+
+    call ap%set_bld()
+
+    res = psb_success_
+  end function
+
 end module psb_z_psblas_cbind_mod