diff --git a/examples/fileread/mld_cexample_1lev.f90 b/examples/fileread/mld_cexample_1lev.f90
index f2522c4c..1b652efa 100644
--- a/examples/fileread/mld_cexample_1lev.f90
+++ b/examples/fileread/mld_cexample_1lev.f90
@@ -193,7 +193,7 @@ program mld_cexample_1lev
 
   ! set RAS
 
-  call P%init('AS',info)
+  call P%init(ictxt,'AS',info)
 
   ! set number of overlaps
 
diff --git a/examples/fileread/mld_cexample_ml.f90 b/examples/fileread/mld_cexample_ml.f90
index 3cbc42d9..b95c0617 100644
--- a/examples/fileread/mld_cexample_ml.f90
+++ b/examples/fileread/mld_cexample_ml.f90
@@ -212,7 +212,7 @@ program mld_cexample_ml
     ! GS sweep as pre/post-smoother and UMFPACK as coarsest-level
     ! solver
 
-    call P%init('ML',info)
+    call P%init(ictxt,'ML',info)
     kmethod = 'CG'
 
   case(2)
diff --git a/examples/fileread/mld_dexample_1lev.f90 b/examples/fileread/mld_dexample_1lev.f90
index 900d070e..126fa0ce 100644
--- a/examples/fileread/mld_dexample_1lev.f90
+++ b/examples/fileread/mld_dexample_1lev.f90
@@ -193,7 +193,7 @@ program mld_dexample_1lev
 
   ! set RAS
 
-  call P%init('AS',info)
+  call P%init(ictxt,'AS',info)
 
   ! set number of overlaps
 
diff --git a/examples/fileread/mld_dexample_ml.f90 b/examples/fileread/mld_dexample_ml.f90
index 7fe06f7e..6a3e6b1a 100644
--- a/examples/fileread/mld_dexample_ml.f90
+++ b/examples/fileread/mld_dexample_ml.f90
@@ -212,7 +212,7 @@ program mld_dexample_ml
     ! GS sweep as pre/post-smoother and UMFPACK as coarsest-level
     ! solver
 
-    call P%init('ML',info)
+    call P%init(ictxt,'ML',info)
     kmethod = 'CG'
 
   case(2)
diff --git a/examples/fileread/mld_sexample_1lev.f90 b/examples/fileread/mld_sexample_1lev.f90
index d4708e08..899c8734 100644
--- a/examples/fileread/mld_sexample_1lev.f90
+++ b/examples/fileread/mld_sexample_1lev.f90
@@ -193,7 +193,7 @@ program mld_sexample_1lev
 
   ! set RAS
 
-  call P%init('AS',info)
+  call P%init(ictxt,'AS',info)
 
   ! set number of overlaps
 
diff --git a/examples/fileread/mld_sexample_ml.f90 b/examples/fileread/mld_sexample_ml.f90
index 7b24bc31..88deffad 100644
--- a/examples/fileread/mld_sexample_ml.f90
+++ b/examples/fileread/mld_sexample_ml.f90
@@ -212,7 +212,7 @@ program mld_sexample_ml
     ! GS sweep as pre/post-smoother and UMFPACK as coarsest-level
     ! solver
 
-    call P%init('ML',info)
+    call P%init(ictxt,'ML',info)
     kmethod = 'CG'
 
   case(2)
diff --git a/examples/fileread/mld_zexample_1lev.f90 b/examples/fileread/mld_zexample_1lev.f90
index cc659515..0fc10b7a 100644
--- a/examples/fileread/mld_zexample_1lev.f90
+++ b/examples/fileread/mld_zexample_1lev.f90
@@ -193,7 +193,7 @@ program mld_zexample_1lev
 
   ! set RAS
 
-  call P%init('AS',info)
+  call P%init(ictxt,'AS',info)
 
   ! set number of overlaps
 
diff --git a/examples/fileread/mld_zexample_ml.f90 b/examples/fileread/mld_zexample_ml.f90
index 9c6cd54c..b711b9f4 100644
--- a/examples/fileread/mld_zexample_ml.f90
+++ b/examples/fileread/mld_zexample_ml.f90
@@ -212,7 +212,7 @@ program mld_zexample_ml
     ! GS sweep as pre/post-smoother and UMFPACK as coarsest-level
     ! solver
 
-    call P%init('ML',info)
+    call P%init(ictxt,'ML',info)
     kmethod = 'CG'
 
   case(2)
diff --git a/examples/pdegen/mld_dexample_1lev.f90 b/examples/pdegen/mld_dexample_1lev.f90
index 0b74503b..d26d4887 100644
--- a/examples/pdegen/mld_dexample_1lev.f90
+++ b/examples/pdegen/mld_dexample_1lev.f90
@@ -140,7 +140,7 @@ program mld_dexample_1lev
 
   ! set RAS
 
-  call P%init('AS',info)
+  call P%init(ictxt,'AS',info)
 
   ! set number of overlaps
 
diff --git a/examples/pdegen/mld_dexample_ml.f90 b/examples/pdegen/mld_dexample_ml.f90
index 4df65dd5..1895b6f1 100644
--- a/examples/pdegen/mld_dexample_ml.f90
+++ b/examples/pdegen/mld_dexample_ml.f90
@@ -168,7 +168,7 @@ program mld_dexample_ml
     ! GS sweep as pre/post-smoother and UMFPACK as coarsest-level
     ! solver
 
-    call P%init('ML',info)
+    call P%init(ictxt,'ML',info)
     kmethod = 'CG'
 
   case(2)
@@ -177,7 +177,7 @@ program mld_dexample_ml
     ! ILU(0) on the blocks) as pre- and post-smoother, and 8  block-Jacobi
     ! sweeps (with ILU(0) on the blocks) as coarsest-level solver
 
-    call P%init('ML',info)
+    call P%init(ictxt,'ML',info)
     call P%set('SMOOTHER_TYPE','BJAC',info)
     call P%set('COARSE_SOLVE','BJAC',info)
     call P%set('COARSE_SWEEPS',8,info)
@@ -189,7 +189,7 @@ program mld_dexample_ml
    ! GS sweeps as pre/post-smoother, a distributed coarsest
    ! matrix, and MUMPS as coarsest-level solver
 
-    call P%init('ML',info)
+    call P%init(ictxt,'ML',info)
     call P%set('ML_CYCLE','WCYCLE',info)
     call P%set('SMOOTHER_SWEEPS',2,info)
     call P%set('COARSE_SOLVE','MUMPS',info)
diff --git a/examples/pdegen/mld_sexample_1lev.f90 b/examples/pdegen/mld_sexample_1lev.f90
index 7c7966bc..40cc13e0 100644
--- a/examples/pdegen/mld_sexample_1lev.f90
+++ b/examples/pdegen/mld_sexample_1lev.f90
@@ -140,7 +140,7 @@ program mld_sexample_1lev
 
   ! set RAS
 
-  call P%init('AS',info)
+  call P%init(ictxt,'AS',info)
 
   ! set number of overlaps
 
diff --git a/examples/pdegen/mld_sexample_ml.f90 b/examples/pdegen/mld_sexample_ml.f90
index 5d5bc7f6..c65176d8 100644
--- a/examples/pdegen/mld_sexample_ml.f90
+++ b/examples/pdegen/mld_sexample_ml.f90
@@ -168,7 +168,7 @@ program mld_sexample_ml
     ! GS sweep as pre/post-smoother and UMFPACK as coarsest-level
     ! solver
 
-    call P%init('ML',info)
+    call P%init(ictxt,'ML',info)
     kmethod = 'CG'
 
   case(2)
@@ -177,7 +177,7 @@ program mld_sexample_ml
     ! ILU(0) on the blocks) as pre- and post-smoother, and 8  block-Jacobi
     ! sweeps (with ILU(0) on the blocks) as coarsest-level solver
 
-    call P%init('ML',info)
+    call P%init(ictxt,'ML',info)
     call P%set('SMOOTHER_TYPE','BJAC',info)
     call P%set('COARSE_SOLVE','BJAC',info)
     call P%set('COARSE_SWEEPS',8,info)
@@ -189,7 +189,7 @@ program mld_sexample_ml
    ! GS sweeps as pre/post-smoother, a distributed coarsest
    ! matrix, and MUMPS as coarsest-level solver
 
-    call P%init('ML',info)
+    call P%init(ictxt,'ML',info)
     call P%set('ML_CYCLE','WCYCLE',info)
     call P%set('SMOOTHER_SWEEPS',2,info)
     call P%set('COARSE_SOLVE','MUMPS',info)
diff --git a/tests/fileread/mld_cf_sample.f90 b/tests/fileread/mld_cf_sample.f90
index c494301b..490e460a 100644
--- a/tests/fileread/mld_cf_sample.f90
+++ b/tests/fileread/mld_cf_sample.f90
@@ -361,7 +361,7 @@ program mld_cf_sample
   !
   ! initialize the preconditioner
   !
-  call prec%init(p_choice%ptype,info)
+  call prec%init(ictxt,p_choice%ptype,info)
   select case(trim(psb_toupper(p_choice%ptype)))
   case ('NONE','NOPREC')
     ! Do nothing, keep defaults
diff --git a/tests/fileread/mld_df_sample.f90 b/tests/fileread/mld_df_sample.f90
index 7a84e045..c40f083b 100644
--- a/tests/fileread/mld_df_sample.f90
+++ b/tests/fileread/mld_df_sample.f90
@@ -361,7 +361,7 @@ program mld_df_sample
   !
   ! initialize the preconditioner
   !
-  call prec%init(p_choice%ptype,info)
+  call prec%init(ictxt,p_choice%ptype,info)
   select case(trim(psb_toupper(p_choice%ptype)))
   case ('NONE','NOPREC')
     ! Do nothing, keep defaults
diff --git a/tests/fileread/mld_sf_sample.f90 b/tests/fileread/mld_sf_sample.f90
index 58ace825..f2a9b7b6 100644
--- a/tests/fileread/mld_sf_sample.f90
+++ b/tests/fileread/mld_sf_sample.f90
@@ -361,7 +361,7 @@ program mld_sf_sample
   !
   ! initialize the preconditioner
   !
-  call prec%init(p_choice%ptype,info)
+  call prec%init(ictxt,p_choice%ptype,info)
   select case(trim(psb_toupper(p_choice%ptype)))
   case ('NONE','NOPREC')
     ! Do nothing, keep defaults
diff --git a/tests/fileread/mld_zf_sample.f90 b/tests/fileread/mld_zf_sample.f90
index 5474305f..9fdcb2dd 100644
--- a/tests/fileread/mld_zf_sample.f90
+++ b/tests/fileread/mld_zf_sample.f90
@@ -361,7 +361,7 @@ program mld_zf_sample
   !
   ! initialize the preconditioner
   !
-  call prec%init(p_choice%ptype,info)
+  call prec%init(ictxt,p_choice%ptype,info)
   select case(trim(psb_toupper(p_choice%ptype)))
   case ('NONE','NOPREC')
     ! Do nothing, keep defaults
diff --git a/tests/pdegen/mld_d_pde2d.f90 b/tests/pdegen/mld_d_pde2d.f90
index ce866028..d17c83c9 100644
--- a/tests/pdegen/mld_d_pde2d.f90
+++ b/tests/pdegen/mld_d_pde2d.f90
@@ -248,7 +248,7 @@ contains
     m   = idim*idim
     n   = m
     nnz = 7*((n+np-1)/np)
-    if(iam == psb_root_) write(psb_out_unit,'("Generating Matrix (size=",i0,")...")')n
+    if (iam == psb_root_) write(psb_out_unit,'("Generating Matrix (size=",i0,")...")')n
     t0 = psb_wtime()
     select case(partition_)
     case(1)
@@ -687,7 +687,7 @@ program mld_d_pde2d
   !
   ! initialize the preconditioner
   !
-  call prec%init(p_choice%ptype,info)
+  call prec%init(ictxt,p_choice%ptype,info)
   select case(trim(psb_toupper(p_choice%ptype)))
   case ('NONE','NOPREC')
     ! Do nothing, keep defaults
diff --git a/tests/pdegen/mld_d_pde3d.f90 b/tests/pdegen/mld_d_pde3d.f90
index 633d4fe0..6b4f491d 100644
--- a/tests/pdegen/mld_d_pde3d.f90
+++ b/tests/pdegen/mld_d_pde3d.f90
@@ -727,7 +727,7 @@ program mld_d_pde3d
   !
   ! initialize the preconditioner
   !
-  call prec%init(p_choice%ptype,info)
+  call prec%init(ictxt,p_choice%ptype,info)
   select case(trim(psb_toupper(p_choice%ptype)))
   case ('NONE','NOPREC')
     ! Do nothing, keep defaults
diff --git a/tests/pdegen/mld_s_pde2d.f90 b/tests/pdegen/mld_s_pde2d.f90
index a7fea3e5..e97e6210 100644
--- a/tests/pdegen/mld_s_pde2d.f90
+++ b/tests/pdegen/mld_s_pde2d.f90
@@ -248,7 +248,7 @@ contains
     m   = idim*idim
     n   = m
     nnz = 7*((n+np-1)/np)
-    if(iam == psb_root_) write(psb_out_unit,'("Generating Matrix (size=",i0,")...")')n
+    if (iam == psb_root_) write(psb_out_unit,'("Generating Matrix (size=",i0,")...")')n
     t0 = psb_wtime()
     select case(partition_)
     case(1)
@@ -687,7 +687,7 @@ program mld_s_pde2d
   !
   ! initialize the preconditioner
   !
-  call prec%init(p_choice%ptype,info)
+  call prec%init(ictxt,p_choice%ptype,info)
   select case(trim(psb_toupper(p_choice%ptype)))
   case ('NONE','NOPREC')
     ! Do nothing, keep defaults
diff --git a/tests/pdegen/mld_s_pde3d.f90 b/tests/pdegen/mld_s_pde3d.f90
index 1978fb81..3104cb02 100644
--- a/tests/pdegen/mld_s_pde3d.f90
+++ b/tests/pdegen/mld_s_pde3d.f90
@@ -727,7 +727,7 @@ program mld_s_pde3d
   !
   ! initialize the preconditioner
   !
-  call prec%init(p_choice%ptype,info)
+  call prec%init(ictxt,p_choice%ptype,info)
   select case(trim(psb_toupper(p_choice%ptype)))
   case ('NONE','NOPREC')
     ! Do nothing, keep defaults