diff --git a/mlprec/mld_dslu_bld.f90 b/mlprec/mld_dslu_bld.f90
index 4fb0e7a0..ef8c500f 100644
--- a/mlprec/mld_dslu_bld.f90
+++ b/mlprec/mld_dslu_bld.f90
@@ -103,7 +103,7 @@ subroutine mld_dslu_bld(a,desc_a,p,info)
   !
   ! Compute the LU factorization
   !
-  call mld_dslu_factor(a%m,nzt,&
+  call mld_dslu_fact(a%m,nzt,&
        & a%aspk,a%ia2,a%ia1,p%iprcparm(mld_slu_ptr_),info)
 
   if (info /= 0) then
diff --git a/mlprec/mld_dslud_bld.f90 b/mlprec/mld_dslud_bld.f90
index 4d610842..de65c3f8 100644
--- a/mlprec/mld_dslud_bld.f90
+++ b/mlprec/mld_dslud_bld.f90
@@ -121,11 +121,11 @@ subroutine mld_dsludist_bld(a,desc_a,p,info)
   !
   ! Compute the LU factorization
   !
-  call mld_dsludist_factor(mglob,nrow,nzt,ifrst,&
+  call mld_dsludist_fact(mglob,nrow,nzt,ifrst,&
        & a%aspk,a%ia2,a%ia1,p%iprcparm(mld_slud_ptr_),&
        & npr, npc, info)
   if (info /= 0) then
-    ch_err='psb_slud_fact'
+    ch_err='psb_sludist_fact'
     call psb_errpush(4110,name,a_err=ch_err,i_err=(/info,0,0,0,0/))
     goto 9999
   end if
diff --git a/mlprec/mld_dumf_bld.f90 b/mlprec/mld_dumf_bld.f90
index d47a4b7c..0a2cd086 100644
--- a/mlprec/mld_dumf_bld.f90
+++ b/mlprec/mld_dumf_bld.f90
@@ -61,7 +61,7 @@
 !               to be factorized. Note that a is intent(inout), and not only
 !               intent(in), since the row and column indices of the      matrix
 !               stored in a are shifted by -1, and then again by +1, by the
-!               routine mld_dumf_factor, which is an interface to the UMFPACK
+!               routine mld_dumf_fact, which is an interface to the UMFPACK
 !               C code performing the factorization.
 !    desc_a  -  type(psb_desc_type), input.
 !               The communication descriptor associated to a.
@@ -108,7 +108,7 @@ subroutine mld_dumf_bld(a,desc_a,p,info)
   !
   ! Compute the LU factorization
   !
-  call mld_dumf_factor(a%m,nzt,&
+  call mld_dumf_fact(a%m,nzt,&
        & a%aspk,a%ia1,a%ia2,&
        & p%iprcparm(mld_umf_symptr_),p%iprcparm(mld_umf_numptr_),info)
 
diff --git a/mlprec/mld_slu_impl.c b/mlprec/mld_slu_impl.c
index 8d48379c..2c3ae623 100644
--- a/mlprec/mld_slu_impl.c
+++ b/mlprec/mld_slu_impl.c
@@ -36,7 +36,7 @@
  *
  * File: mld_slu_impl.c
  *
- * Functions: mld_dslu_factor_, mld_dslu_solve_, mld_dslu_free_.
+ * Functions: mld_dslu_fact_, mld_dslu_solve_, mld_dslu_free_.
  *
  * This file is an interface to the SuperLU routines for sparse factorization and
  * solve. It was obtained by modifying the c_fortran_dgssv.c file from the SuperLU
@@ -115,17 +115,17 @@ typedef struct {
 
 
 #ifdef Add_
-#define mld_dslu_factor_ mld_dslu_factor_
+#define mld_dslu_fact_   mld_dslu_fact_
 #define mld_dslu_solve_  mld_dslu_solve_
 #define mld_dslu_free_   mld_dslu_free_
 #endif
 #ifdef AddDouble_
-#define mld_dslu_factor_ mld_dslu_factor__
+#define mld_dslu_fact_   mld_dslu_fact__
 #define mld_dslu_solve_  mld_dslu_solve__
 #define mld_dslu_free_   mld_dslu_free__
 #endif
 #ifdef NoChange
-#define mld_dslu_factor_ mld_dslu_factor
+#define mld_dslu_fact_   mld_dslu_fact
 #define mld_dslu_solve_  mld_dslu_solve
 #define mld_dslu_free_   mld_dslu_free
 #endif
@@ -134,7 +134,7 @@ typedef struct {
 
 
 void
-mld_dslu_factor_(int *n, int *nnz,
+mld_dslu_fact_(int *n, int *nnz,
                  double *values, int *rowptr, int *colind,
 #ifdef Have_SLU_		 
 		 fptr *f_factors, /* a handle containing the address
diff --git a/mlprec/mld_slud_impl.c b/mlprec/mld_slud_impl.c
index 2fd77d18..2daa9b80 100644
--- a/mlprec/mld_slud_impl.c
+++ b/mlprec/mld_slud_impl.c
@@ -36,7 +36,7 @@
  *
  * File: mld_slud_impl.c
  *
- * Functions: mld_dsludist_factor_, mld_dsludist_solve_, mld_dsludist_free_.
+ * Functions: mld_dsludist_fact_, mld_dsludist_solve_, mld_dsludist_free_.
  *
  * This file is an interface to the SuperLU_dist routines for sparse factorization and
  * solve. It was obtained by modifying the c_fortran_dgssv.c file from the SuperLU_dist
@@ -114,17 +114,17 @@ typedef struct {
 
 
 #ifdef Add_
-#define mld_dsludist_factor_ mld_dsludist_factor_
+#define mld_dsludist_fact_   mld_dsludist_fact_
 #define mld_dsludist_solve_  mld_dsludist_solve_
 #define mld_dsludist_free_   mld_dsludist_free_
 #endif
 #ifdef AddDouble_
-#define mld_dsludist_factor_ mld_dsludist_factor__
+#define mld_dsludist_fact_   mld_dsludist_fact__
 #define mld_dsludist_solve_  mld_dsludist_solve__
 #define mld_dsludist_free_   mld_dsludist_free__
 #endif
 #ifdef NoChange
-#define mld_dsludist_factor_ mld_dsludist_factor
+#define mld_dsludist_fact_   mld_dsludist_fact
 #define mld_dsludist_solve_  mld_dsludist_solve
 #define mld_dsludist_free_   mld_dsludist_free
 #endif
@@ -133,7 +133,7 @@ typedef struct {
 
 
 void
-mld_dsludist_factor_(int *n, int *nl, int *nnzl, int *ffstr,
+mld_dsludist_fact_(int *n, int *nl, int *nnzl, int *ffstr,
 		     double *values, int *rowptr, int *colind,
 #ifdef Have_SLUDist_		 
 		     fptr *f_factors, /* a handle containing the address
@@ -174,7 +174,7 @@ mld_dsludist_factor_(int *n, int *nl, int *nnzl, int *ffstr,
     double *ival;
 
     trans = NOTRANS;
-/*     fprintf(stderr,"Entry to sludist_factor\n");     */
+/*     fprintf(stderr,"Entry to sludist_fact\n");     */
     grid = (gridinfo_t *) SUPERLU_MALLOC(sizeof(gridinfo_t));
     superlu_gridinit(MPI_COMM_WORLD, *nprow, *npcol, grid);
     /* Initialize the statistics variables. */
diff --git a/mlprec/mld_umf_impl.c b/mlprec/mld_umf_impl.c
index 82cdbcfd..ef5c8bf9 100644
--- a/mlprec/mld_umf_impl.c
+++ b/mlprec/mld_umf_impl.c
@@ -36,7 +36,7 @@
  *
  * File: mld_umf_impl.c
  *
- * Functions: mld_dumf_factor_, mld_dumf_solve_, mld_umf_free_.
+ * Functions: mld_dumf_fact_, mld_dumf_solve_, mld_umf_free_.
  *
  * This file is an interface to the UMFPACK routines for sparse factorization and
  * solve. It was obtained by adapting umfpack_di_demo under the original UMFPACK
@@ -76,17 +76,17 @@ Availability:
 
 
 #ifdef Add_
-#define mld_dumf_factor_ mld_dumf_factor_
+#define mld_dumf_fact_   mld_dumf_fact_
 #define mld_dumf_solve_  mld_dumf_solve_
 #define mld_dumf_free_   mld_dumf_free_
 #endif
 #ifdef AddDouble_
-#define mld_dumf_factor_ mld_dumf_factor__
+#define mld_dumf_fact_   mld_dumf_fact__
 #define mld_dumf_solve_  mld_dumf_solve__
 #define mld_dumf_free_   mld_dumf_free__
 #endif
 #ifdef NoChange
-#define mld_dumf_factor_ mld_dumf_factor
+#define mld_dumf_fact_   mld_dumf_fact
 #define mld_dumf_solve_  mld_dumf_solve
 #define mld_dumf_free_   mld_dumf_free
 #endif
@@ -104,7 +104,7 @@ typedef int fptr;  /* 32-bit by default */
 #endif
 
 void
-mld_dumf_factor_(int *n, int *nnz,
+mld_dumf_fact_(int *n, int *nnz,
                  double *values, int *rowind, int *colptr,
 #ifdef Have_UMF_		 
 		 fptr *symptr, 
diff --git a/mlprec/mld_zslu_bld.f90 b/mlprec/mld_zslu_bld.f90
index a6fe46e8..b5be9719 100644
--- a/mlprec/mld_zslu_bld.f90
+++ b/mlprec/mld_zslu_bld.f90
@@ -103,7 +103,7 @@ subroutine mld_zslu_bld(a,desc_a,p,info)
   !
   ! Compute the LU factorization
   !
-  call mld_zslu_factor(a%m,nzt,&
+  call mld_zslu_fact(a%m,nzt,&
        & a%aspk,a%ia2,a%ia1,p%iprcparm(mld_slu_ptr_),info)
 
   if (info /= 0) then
diff --git a/mlprec/mld_zslu_impl.c b/mlprec/mld_zslu_impl.c
index 298d1159..d4b7259e 100644
--- a/mlprec/mld_zslu_impl.c
+++ b/mlprec/mld_zslu_impl.c
@@ -36,7 +36,7 @@
  *
  * File: mld_zslu_impl.c
  *
- * Functions: mld_zslu_factor_, mld_zslu_solve_, mld_zslu_free_.
+ * Functions: mld_zslu_fact_, mld_zslu_solve_, mld_zslu_free_.
  *
  * This file is an interface to the SuperLU routines for sparse factorization and
  * solve. It was obtained by modifying the c_fortran_zgssv.c file from the SuperLU
@@ -114,17 +114,17 @@ typedef struct {
 
 
 #ifdef Add_
-#define mld_zslu_factor_ mld_zslu_factor_
+#define mld_zslu_fact_   mld_zslu_fact_
 #define mld_zslu_solve_  mld_zslu_solve_
 #define mld_zslu_free_   mld_zslu_free_
 #endif
 #ifdef AddDouble_
-#define mld_zslu_factor_ mld_zslu_factor__
+#define mld_zslu_fact_   mld_zslu_fact__
 #define mld_zslu_solve_  mld_zslu_solve__
 #define mld_zslu_free_   mld_zslu_free__
 #endif
 #ifdef NoChange
-#define mld_zslu_factor_ mld_zslu_factor
+#define mld_zslu_fact_   mld_zslu_fact
 #define mld_zslu_solve_  mld_zslu_solve
 #define mld_zslu_free_   mld_zslu_free
 #endif
@@ -133,7 +133,7 @@ typedef struct {
 
 
 void
-mld_zslu_factor_(int *n, int *nnz,
+mld_zslu_fact_(int *n, int *nnz,
 #ifdef Have_SLU_		 
                  doublecomplex *values, int *colind, int *rowptr,
 		 fptr *f_factors, /* a handle containing the address
diff --git a/mlprec/mld_zslud_bld.f90 b/mlprec/mld_zslud_bld.f90
index 330d98b8..b52c9d05 100644
--- a/mlprec/mld_zslud_bld.f90
+++ b/mlprec/mld_zslud_bld.f90
@@ -121,11 +121,11 @@ subroutine mld_zsludist_bld(a,desc_a,p,info)
   !
   ! Compute the LU factorization
   !
-  call mld_zsludist_factor(mglob,nrow,nzt,ifrst,&
+  call mld_zsludist_fact(mglob,nrow,nzt,ifrst,&
        & a%aspk,a%ia2,a%ia1,p%iprcparm(mld_slud_ptr_),&
        & npr, npc, info)
   if (info /= 0) then
-    ch_err='psb_slud_fact'
+    ch_err='psb_sludist_fact'
     call psb_errpush(4110,name,a_err=ch_err,i_err=(/info,0,0,0,0/))
     goto 9999
   end if
diff --git a/mlprec/mld_zslud_impl.c b/mlprec/mld_zslud_impl.c
index 66332eca..f5b13b84 100644
--- a/mlprec/mld_zslud_impl.c
+++ b/mlprec/mld_zslud_impl.c
@@ -36,7 +36,7 @@
  *
  * File: mld_zslud_impl.c
  *
- * Functions: mld_zsludist_factor_, mld_zsludist_solve_, mld_zsludist_free_.
+ * Functions: mld_zsludist_fact_, mld_zsludist_solve_, mld_zsludist_free_.
  *
  * This file is an interface to the SuperLU_dist routines for sparse factorization and
  * solve. It was obtained by modifying the c_fortran_zgssv.c file from the SuperLU_dist
@@ -114,17 +114,17 @@ typedef struct {
 
 
 #ifdef Add_
-#define mld_zsludist_factor_ mld_zsludist_factor_
+#define mld_zsludist_fact_   mld_zsludist_fact_
 #define mld_zsludist_solve_  mld_zsludist_solve_
 #define mld_zsludist_free_   mld_zsludist_free_
 #endif
 #ifdef AddDouble_
-#define mld_zsludist_factor_ mld_zsludist_factor__
+#define mld_zsludist_fact_   mld_zsludist_fact__
 #define mld_zsludist_solve_  mld_zsludist_solve__
 #define mld_zsludist_free_   mld_zsludist_free__
 #endif
 #ifdef NoChange
-#define mld_zsludist_factor_ mld_zsludist_factor
+#define mld_zsludist_fact_   mld_zsludist_fact
 #define mld_zsludist_solve_  mld_zsludist_solve
 #define mld_zsludist_free_   mld_zsludist_free
 #endif
@@ -133,7 +133,7 @@ typedef struct {
 
 
 void
-mld_zsludist_factor_(int *n, int *nl, int *nnzl, int *ffstr,
+mld_zsludist_fact_(int *n, int *nl, int *nnzl, int *ffstr,
 #ifdef Have_SLUDist_		 
 		     doublecomplex *values, int *rowptr, int *colind,
 		     fptr *f_factors, /* a handle containing the address
diff --git a/mlprec/mld_zumf_bld.f90 b/mlprec/mld_zumf_bld.f90
index 5012b076..6b2fb5d2 100644
--- a/mlprec/mld_zumf_bld.f90
+++ b/mlprec/mld_zumf_bld.f90
@@ -61,7 +61,7 @@
 !               to be factorized. Note that a is intent(inout), and not only
 !               intent(in), since the row and column indices of the      matrix
 !               stored in a are shifted by -1, and then again by +1, by the
-!               routine mld_zumf_factor, which is an interface to the UMFPACK
+!               routine mld_zumf_fact, which is an interface to the UMFPACK
 !               C code performing the factorization.
 !    desc_a  -  type(psb_desc_type), input.
 !               The communication descriptor associated to a.
@@ -108,7 +108,7 @@ subroutine mld_zumf_bld(a,desc_a,p,info)
   !
   ! Compute the LU factorization
   !
-  call mld_zumf_factor(a%m,nzt,&
+  call mld_zumf_fact(a%m,nzt,&
        & a%aspk,a%ia1,a%ia2,&
        & p%iprcparm(mld_umf_symptr_),p%iprcparm(mld_umf_numptr_),info)
 
diff --git a/mlprec/mld_zumf_impl.c b/mlprec/mld_zumf_impl.c
index 5eef3ca1..9b2733c9 100644
--- a/mlprec/mld_zumf_impl.c
+++ b/mlprec/mld_zumf_impl.c
@@ -36,7 +36,7 @@
  *
  * File: mld_zumf_impl.c
  *
- * Functions: mld_zumf_factor_, mld_zumf_solve_, mld_zumf_free_.
+ * Functions: mld_zumf_fact_, mld_zumf_solve_, mld_zumf_free_.
  *
  * This file is an interface to the UMFPACK routines for sparse factorization and
  * solve. It was obtained by adapting umfpack_zi_demo under the original UMFPACK
@@ -76,17 +76,17 @@ Availability:
 
 
 #ifdef Add_
-#define mld_zumf_factor_ mld_zumf_factor_
+#define mld_zumf_fact_   mld_zumf_fact_
 #define mld_zumf_solve_  mld_zumf_solve_
 #define mld_zumf_free_   mld_zumf_free_
 #endif
 #ifdef AddDouble_
-#define mld_zumf_factor_ mld_zumf_factor__
+#define mld_zumf_fact_   mld_zumf_fact__
 #define mld_zumf_solve_  mld_zumf_solve__
 #define mld_zumf_free_   mld_zumf_free__
 #endif
 #ifdef NoChange
-#define mld_zumf_factor_ mld_zumf_factor
+#define mld_zumf_fact_   mld_zumf_fact
 #define mld_zumf_solve_  mld_zumf_solve
 #define mld_zumf_free_   mld_zumf_free
 #endif
@@ -104,7 +104,7 @@ typedef int fptr;  /* 32-bit by default */
 #endif
 
 void
-mld_zumf_factor_(int *n, int *nnz,
+mld_zumf_fact_(int *n, int *nnz,
                  double *values, int *rowind, int *colptr,
 #ifdef Have_UMF_		 
 		 fptr *symptr,