diff --git a/mlprec/impl/solver/mld_c_mumps_solver_apply.F90 b/mlprec/impl/solver/mld_c_mumps_solver_apply.F90
index 4dbfc292..ee00b9f5 100644
--- a/mlprec/impl/solver/mld_c_mumps_solver_apply.F90
+++ b/mlprec/impl/solver/mld_c_mumps_solver_apply.F90
@@ -107,7 +107,7 @@ subroutine c_mumps_solver_apply(alpha,sv,x,beta,y,desc_data,&
            & a_err='complex(psb_spk_)')
       goto 9999      
     end if
-    call psb_gather(gx, x, desc_data, info, root=0)
+    call psb_gather(gx, x, desc_data, info, root=izero)
   else
       info=psb_err_internal_error_
       call psb_errpush(info,name,& 
diff --git a/mlprec/impl/solver/mld_d_mumps_solver_apply.F90 b/mlprec/impl/solver/mld_d_mumps_solver_apply.F90
index bddc35df..5ee8763d 100644
--- a/mlprec/impl/solver/mld_d_mumps_solver_apply.F90
+++ b/mlprec/impl/solver/mld_d_mumps_solver_apply.F90
@@ -107,7 +107,7 @@ subroutine d_mumps_solver_apply(alpha,sv,x,beta,y,desc_data,&
            & a_err='real(psb_dpk_)')
       goto 9999      
     end if
-    call psb_gather(gx, x, desc_data, info, root=0)
+    call psb_gather(gx, x, desc_data, info, root=izero)
   else
       info=psb_err_internal_error_
       call psb_errpush(info,name,& 
diff --git a/mlprec/impl/solver/mld_s_mumps_solver_apply.F90 b/mlprec/impl/solver/mld_s_mumps_solver_apply.F90
index d84817e1..8870dbe4 100644
--- a/mlprec/impl/solver/mld_s_mumps_solver_apply.F90
+++ b/mlprec/impl/solver/mld_s_mumps_solver_apply.F90
@@ -107,7 +107,7 @@ subroutine s_mumps_solver_apply(alpha,sv,x,beta,y,desc_data,&
            & a_err='real(psb_spk_)')
       goto 9999      
     end if
-    call psb_gather(gx, x, desc_data, info, root=0)
+    call psb_gather(gx, x, desc_data, info, root=izero)
   else
       info=psb_err_internal_error_
       call psb_errpush(info,name,& 
diff --git a/mlprec/impl/solver/mld_z_mumps_solver_apply.F90 b/mlprec/impl/solver/mld_z_mumps_solver_apply.F90
index 2a0f6214..94e6758e 100644
--- a/mlprec/impl/solver/mld_z_mumps_solver_apply.F90
+++ b/mlprec/impl/solver/mld_z_mumps_solver_apply.F90
@@ -107,7 +107,7 @@ subroutine z_mumps_solver_apply(alpha,sv,x,beta,y,desc_data,&
            & a_err='complex(psb_dpk_)')
       goto 9999      
     end if
-    call psb_gather(gx, x, desc_data, info, root=0)
+    call psb_gather(gx, x, desc_data, info, root=izero)
   else
       info=psb_err_internal_error_
       call psb_errpush(info,name,&