Benchmark version.

cbind/mlprec/mld_c_dprec.h

@@ -29,7 +29,7 @@ extern "C" {
   psb_i_t mld_c_dprecfree(mld_c_dprec *ph);
   psb_i_t mld_c_dprecbld_opt(psb_c_dspmat *ah, psb_c_descriptor *cdh, 
 			  mld_c_dprec *ph, const char *afmt);
-
+  psb_i_t mld_c_ddescr(mld_c_dprec *ph);
 
   psb_i_t mld_c_dkrylov(const char *method, psb_c_dspmat *ah, mld_c_dprec *ph, 
 		  psb_c_dvector *bh, psb_c_dvector *xh,

cbind/mlprec/mld_c_zprec.h

@@ -30,6 +30,7 @@ extern "C" {
   psb_i_t mld_c_zprecbld_opt(psb_c_zspmat *ah, psb_c_descriptor *cdh, 
 			  mld_c_zprec *ph, const char *afmt);
 
+  psb_i_t mld_c_zdescr(mld_c_zprec *ph);
 
   psb_i_t mld_c_zkrylov(const char *method, psb_c_zspmat *ah, mld_c_zprec *ph, 
 		  psb_c_zvector *bh, psb_c_zvector *xh,

cbind/mlprec/mld_dprec_cbind_mod.F90

@@ -377,4 +377,32 @@ contains
     return
   end function mld_c_dprecfree
 
+  function mld_c_ddescr(ph) bind(c) result(res)
+   use psb_base_mod
+   use mld_prec_mod
+   implicit none
+ 
+   integer(psb_c_ipk_) :: res
+   type(psb_c_object_type) :: ph
+   integer               :: info
+   type(mld_dprec_type), pointer :: precp
+ 
+   res = -1
+   info = -1
+   if (c_associated(ph%item)) then
+     call c_f_pointer(ph%item,precp)
+   else
+     return
+   end if
+ 
+ 
+   call precp%descr()
+   call flush(output_unit)
+ 
+   info = 0
+   res = MLDC_ERR_FILTER(info)
+   MLDC_ERR_HANDLE(res)
+   return
+ end function mld_c_ddescr
+
 end module mld_dprec_cbind_mod

cbind/mlprec/mld_zprec_cbind_mod.F90

@@ -377,4 +377,32 @@ contains
     return
   end function mld_c_zprecfree
 
+  function mld_c_zdescr(ph) bind(c) result(res)
+   use psb_base_mod
+   use mld_prec_mod
+   implicit none
+ 
+   integer(psb_c_ipk_) :: res
+   type(psb_c_object_type) :: ph
+   integer               :: info
+   type(mld_zprec_type), pointer :: precp
+ 
+   res = -1
+   info = -1
+   if (c_associated(ph%item)) then
+     call c_f_pointer(ph%item,precp)
+   else
+     return
+   end if
+ 
+ 
+   call precp%descr()
+   call flush(output_unit)
+ 
+   info = 0
+   res = MLDC_ERR_FILTER(info)
+   MLDC_ERR_HANDLE(res)
+   return
+ end function mld_c_zdescr
+
 end module mld_zprec_cbind_mod

mlprec/impl/mld_c_hierarchy_bld.f90

@@ -443,7 +443,7 @@ subroutine mld_c_hierarchy_bld(a,desc_a,prec,info)
     iszv = newsz
     ! Fix the pointers, but the level 1 should
     ! be treated differently
-    if (.not.associated(prec%precv(i)%base_desc,desc_a)) then
+    if (.not.associated(prec%precv(1)%base_desc,desc_a)) then
       prec%precv(1)%base_desc => prec%precv(1)%desc_ac
     end if
     do i=2, iszv 

mlprec/impl/mld_d_hierarchy_bld.f90

@@ -443,7 +443,7 @@ subroutine mld_d_hierarchy_bld(a,desc_a,prec,info)
     iszv = newsz
     ! Fix the pointers, but the level 1 should
     ! be treated differently
-    if (.not.associated(prec%precv(i)%base_desc,desc_a)) then
+    if (.not.associated(prec%precv(1)%base_desc,desc_a)) then
       prec%precv(1)%base_desc => prec%precv(1)%desc_ac
     end if
     do i=2, iszv 

mlprec/impl/mld_s_hierarchy_bld.f90

@@ -443,7 +443,7 @@ subroutine mld_s_hierarchy_bld(a,desc_a,prec,info)
     iszv = newsz
     ! Fix the pointers, but the level 1 should
     ! be treated differently
-    if (.not.associated(prec%precv(i)%base_desc,desc_a)) then
+    if (.not.associated(prec%precv(1)%base_desc,desc_a)) then
       prec%precv(1)%base_desc => prec%precv(1)%desc_ac
     end if
     do i=2, iszv 

mlprec/impl/mld_z_hierarchy_bld.f90

@@ -443,7 +443,7 @@ subroutine mld_z_hierarchy_bld(a,desc_a,prec,info)
     iszv = newsz
     ! Fix the pointers, but the level 1 should
     ! be treated differently
-    if (.not.associated(prec%precv(i)%base_desc,desc_a)) then
+    if (.not.associated(prec%precv(1)%base_desc,desc_a)) then
       prec%precv(1)%base_desc => prec%precv(1)%desc_ac
     end if
     do i=2, iszv