psblas-3.3-maint:

Changelog
 base/tools/psb_ccdbldext.F90
 base/tools/psb_dcdbldext.F90
 base/tools/psb_scdbldext.F90
 base/tools/psb_zcdbldext.F90

Fix bug in cdbldext: must call indxmap%reinit.

Changelog

@@ -1,5 +1,6 @@
 Changelog. A lot less detailed than usual, at least for past
 	history.
+2015/07/17: Fix bug in cdbldext. 	
 2015/04/18: Fix bug in spmv transpose case.	
 2015/03/28: Fix bug in transpose.
 2015/03/20: Introduce some optimizations for conversions to/from.	

base/tools/psb_ccdbldext.F90

@@ -177,8 +177,9 @@ Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info, extype)
     ! in a separate method. 
     call psb_cd_switch_ovl_indxmap(desc_ov,info) 
   end if
-
-  call psb_cd_set_ovl_bld(desc_ov,info)
+  if (info == 0) call desc_ov%indxmap%reinit(info)
+  if (info == 0) call psb_cd_set_ovl_bld(desc_ov,info)
+  if (info /= 0) goto 9999
 
   If (debug_level >= psb_debug_outer_)then 
     Write(debug_unit,*) me,' ',trim(name),&

base/tools/psb_dcdbldext.F90

@@ -177,8 +177,9 @@ Subroutine psb_dcdbldext(a,desc_a,novr,desc_ov,info, extype)
     ! in a separate method. 
     call psb_cd_switch_ovl_indxmap(desc_ov,info) 
   end if
-
-  call psb_cd_set_ovl_bld(desc_ov,info)
+  if (info == 0) call desc_ov%indxmap%reinit(info)
+  if (info == 0) call psb_cd_set_ovl_bld(desc_ov,info)
+  if (info /= 0) goto 9999
 
   If (debug_level >= psb_debug_outer_)then 
     Write(debug_unit,*) me,' ',trim(name),&

base/tools/psb_scdbldext.F90

@@ -177,8 +177,9 @@ Subroutine psb_scdbldext(a,desc_a,novr,desc_ov,info, extype)
     ! in a separate method. 
     call psb_cd_switch_ovl_indxmap(desc_ov,info) 
   end if
-
-  call psb_cd_set_ovl_bld(desc_ov,info)
+  if (info == 0) call desc_ov%indxmap%reinit(info)
+  if (info == 0) call psb_cd_set_ovl_bld(desc_ov,info)
+  if (info /= 0) goto 9999
 
   If (debug_level >= psb_debug_outer_)then 
     Write(debug_unit,*) me,' ',trim(name),&

base/tools/psb_zcdbldext.F90

@@ -177,8 +177,9 @@ Subroutine psb_zcdbldext(a,desc_a,novr,desc_ov,info, extype)
     ! in a separate method. 
     call psb_cd_switch_ovl_indxmap(desc_ov,info) 
   end if
-
-  call psb_cd_set_ovl_bld(desc_ov,info)
+  if (info == 0) call desc_ov%indxmap%reinit(info)
+  if (info == 0) call psb_cd_set_ovl_bld(desc_ov,info)
+  if (info /= 0) goto 9999
 
   If (debug_level >= psb_debug_outer_)then 
     Write(debug_unit,*) me,' ',trim(name),&