diff --git a/util/Makefile b/util/Makefile
index f19650899..355b69e8e 100644
--- a/util/Makefile
+++ b/util/Makefile
@@ -8,10 +8,12 @@ MODDIR=../modules
 HERE=.
 
 BASEOBJS= psb_blockpart_mod.o psb_metispart_mod.o psb_partidx_mod.o \
-      psb_hbio_mod.o psb_mmio_mod.o psb_mat_dist_mod.o \
-      psb_s_mat_dist_mod.o psb_d_mat_dist_mod.o psb_c_mat_dist_mod.o psb_z_mat_dist_mod.o \
-      psb_renum_mod.o psb_gps_mod.o \
-      psb_s_renum_mod.o psb_d_renum_mod.o psb_c_renum_mod.o psb_z_renum_mod.o
+	psb_hbio_mod.o psb_mmio_mod.o \
+	psb_i_mmio_mod.o psb_s_mmio_mod.o psb_d_mmio_mod.o psb_c_mmio_mod.o psb_z_mmio_mod.o \
+	psb_mat_dist_mod.o \
+        psb_s_mat_dist_mod.o psb_d_mat_dist_mod.o psb_c_mat_dist_mod.o psb_z_mat_dist_mod.o \
+        psb_renum_mod.o psb_gps_mod.o \
+        psb_s_renum_mod.o psb_d_renum_mod.o psb_c_renum_mod.o psb_z_renum_mod.o
 IMPLOBJS= psb_s_hbio_impl.o psb_d_hbio_impl.o \
 	psb_c_hbio_impl.o psb_z_hbio_impl.o \
 	psb_s_mmio_impl.o psb_d_mmio_impl.o \
@@ -46,7 +48,7 @@ psb_util_mod.o: $(BASEOBJS)
 psb_metispart_mod.o: psb_metis_int.o
 psb_mat_dist_mod.o:  psb_s_mat_dist_mod.o psb_d_mat_dist_mod.o psb_c_mat_dist_mod.o psb_z_mat_dist_mod.o 
 psb_renum_mod.o: psb_s_renum_mod.o psb_d_renum_mod.o psb_c_renum_mod.o psb_z_renum_mod.o
-
+psb_mmio_mod.o:	psb_i_mmio_mod.o psb_s_mmio_mod.o psb_d_mmio_mod.o psb_c_mmio_mod.o psb_z_mmio_mod.o
 $(IMPLOBJS): $(BASEOBJS)
 
 veryclean: clean
diff --git a/util/psb_c_mmio_mod.F90 b/util/psb_c_mmio_mod.F90
new file mode 100644
index 000000000..f1bb58977
--- /dev/null
+++ b/util/psb_c_mmio_mod.F90
@@ -0,0 +1,155 @@
+!   
+!                Parallel Sparse BLAS  version 3.5
+!      (C) Copyright 2006-2018
+!        Salvatore Filippone    
+!        Alfredo Buttari      
+!   
+!    Redistribution and use in source and binary forms, with or without
+!    modification, are permitted provided that the following conditions
+!    are met:
+!      1. Redistributions of source code must retain the above copyright
+!         notice, this list of conditions and the following disclaimer.
+!      2. Redistributions in binary form must reproduce the above copyright
+!         notice, this list of conditions, and the following disclaimer in the
+!         documentation and/or other materials provided with the distribution.
+!      3. The name of the PSBLAS group or the names of its contributors may
+!         not be used to endorse or promote products derived from this
+!         software without specific prior written permission.
+!   
+!    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+!    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
+!    TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+!    PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
+!    BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+!    CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+!    SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+!    INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+!    CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+!    ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+!    POSSIBILITY OF SUCH DAMAGE.
+!   
+!    
+module psb_c_mmio_mod
+
+
+  use psb_base_mod, only :  psb_ipk_, psb_lpk_, psb_spk_,&
+       & psb_c_vect_type, psb_cspmat_type, psb_lcspmat_type
+
+  !public mm_mat_read, mm_mat_write, mm_array_read, mm_array_write
+  interface mm_array_read
+    subroutine mm_cvet_read(b, info, iunit, filename)   
+      import :: psb_spk_, psb_ipk_
+      implicit none
+      complex(psb_spk_), allocatable, intent(out)  :: b(:)
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_cvet_read
+    subroutine mm_cvet2_read(b, info, iunit, filename)   
+      import :: psb_spk_, psb_ipk_
+      implicit none
+      complex(psb_spk_), allocatable, intent(out)  :: b(:,:)
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_cvet2_read
+    subroutine mm_cvect_read(b, info, iunit, filename)   
+      import :: psb_ipk_,psb_c_vect_type
+      implicit none
+      type(psb_c_vect_type), intent(inout)  :: b
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in) :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_cvect_read
+  end interface mm_array_read
+
+
+#if 0
+#if ! defined(PSB_HAVE_BUGGY_GENERICS)
+  interface mm_vet_read
+    procedure mm_cvet_read, mm_cvet2_read
+  end interface
+#endif
+#endif
+  
+  interface mm_array_write
+    subroutine mm_cvet2_write(b, header, info, iunit, filename)   
+      import :: psb_spk_, psb_ipk_
+      implicit none
+      complex(psb_spk_), intent(in)  :: b(:,:)
+      character(len=*), intent(in) :: header
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_cvet2_write
+    subroutine mm_cvet1_write(b, header, info, iunit, filename)   
+      import :: psb_spk_, psb_ipk_
+      implicit none
+      complex(psb_spk_), intent(in)  :: b(:)
+      character(len=*), intent(in) :: header
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_cvet1_write
+    subroutine mm_cvect_write(b, header, info, iunit, filename)   
+      import :: psb_ipk_, psb_c_vect_type
+      implicit none
+      type(psb_c_vect_type), intent(inout)    :: b
+      character(len=*), intent(in)            :: header
+      integer(psb_ipk_), intent(out)          :: info
+      integer(psb_ipk_), optional, intent(in) :: iunit
+      character(len=*), optional, intent(in)  :: filename
+    end subroutine mm_cvect_write
+  end interface
+
+#if ! defined(PSB_HAVE_BUGGY_GENERICS)
+  interface mm_vet_write
+    procedure mm_cvet1_write, mm_cvet2_write
+  end interface
+#endif
+
+  interface mm_mat_read
+    subroutine cmm_mat_read(a, info, iunit, filename)   
+      import :: psb_cspmat_type, psb_ipk_
+      implicit none
+      type(psb_cspmat_type), intent(out)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine cmm_mat_read
+    subroutine lcmm_mat_read(a, info, iunit, filename)   
+      import :: psb_lcspmat_type, psb_ipk_
+      implicit none
+      type(psb_lcspmat_type), intent(out)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine lcmm_mat_read
+  end interface
+
+  interface mm_mat_write
+    subroutine cmm_mat_write(a,mtitle,info,iunit,filename)
+      import :: psb_cspmat_type, psb_ipk_ 
+      implicit none
+      type(psb_cspmat_type), intent(in)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      character(len=*), intent(in) :: mtitle
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine cmm_mat_write
+    subroutine lcmm_mat_write(a,mtitle,info,iunit,filename)
+      import :: psb_lcspmat_type, psb_ipk_ 
+      implicit none
+      type(psb_lcspmat_type), intent(in)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      character(len=*), intent(in) :: mtitle
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine lcmm_mat_write
+  end interface
+#if 0
+#if ! defined(PSB_HAVE_BUGGY_GENERICS)
+  public mm_vet_read, mm_vet_write
+#endif
+#endif
+end module psb_c_mmio_mod
diff --git a/util/psb_d_mmio_mod.F90 b/util/psb_d_mmio_mod.F90
new file mode 100644
index 000000000..98ec9dcc6
--- /dev/null
+++ b/util/psb_d_mmio_mod.F90
@@ -0,0 +1,155 @@
+!   
+!                Parallel Sparse BLAS  version 3.5
+!      (C) Copyright 2006-2018
+!        Salvatore Filippone    
+!        Alfredo Buttari      
+!   
+!    Redistribution and use in source and binary forms, with or without
+!    modification, are permitted provided that the following conditions
+!    are met:
+!      1. Redistributions of source code must retain the above copyright
+!         notice, this list of conditions and the following disclaimer.
+!      2. Redistributions in binary form must reproduce the above copyright
+!         notice, this list of conditions, and the following disclaimer in the
+!         documentation and/or other materials provided with the distribution.
+!      3. The name of the PSBLAS group or the names of its contributors may
+!         not be used to endorse or promote products derived from this
+!         software without specific prior written permission.
+!   
+!    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+!    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
+!    TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+!    PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
+!    BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+!    CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+!    SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+!    INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+!    CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+!    ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+!    POSSIBILITY OF SUCH DAMAGE.
+!   
+!    
+module psb_d_mmio_mod
+
+
+  use psb_base_mod, only :  psb_ipk_, psb_lpk_, psb_dpk_,&
+       & psb_d_vect_type, psb_dspmat_type, psb_ldspmat_type
+
+  !public mm_mat_read, mm_mat_write, mm_array_read, mm_array_write
+  interface mm_array_read
+    subroutine mm_dvet_read(b, info, iunit, filename)   
+      import :: psb_dpk_, psb_ipk_
+      implicit none
+      real(psb_dpk_), allocatable, intent(out)  :: b(:)
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_dvet_read
+    subroutine mm_dvet2_read(b, info, iunit, filename)   
+      import :: psb_dpk_, psb_ipk_
+      implicit none
+      real(psb_dpk_), allocatable, intent(out)  :: b(:,:)
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_dvet2_read
+    subroutine mm_dvect_read(b, info, iunit, filename)   
+      import :: psb_ipk_,psb_d_vect_type
+      implicit none
+      type(psb_d_vect_type), intent(inout)  :: b
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in) :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_dvect_read
+  end interface mm_array_read
+
+
+#if 0
+#if ! defined(PSB_HAVE_BUGGY_GENERICS)
+  interface mm_vet_read
+    procedure mm_dvet_read, mm_dvet2_read
+  end interface
+#endif
+#endif
+  
+  interface mm_array_write
+    subroutine mm_dvet2_write(b, header, info, iunit, filename)   
+      import :: psb_dpk_, psb_ipk_
+      implicit none
+      real(psb_dpk_), intent(in)  :: b(:,:)
+      character(len=*), intent(in) :: header
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_dvet2_write
+    subroutine mm_dvet1_write(b, header, info, iunit, filename)   
+      import :: psb_dpk_, psb_ipk_
+      implicit none
+      real(psb_dpk_), intent(in)  :: b(:)
+      character(len=*), intent(in) :: header
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_dvet1_write
+    subroutine mm_dvect_write(b, header, info, iunit, filename)   
+      import :: psb_ipk_, psb_d_vect_type
+      implicit none
+      type(psb_d_vect_type), intent(inout)    :: b
+      character(len=*), intent(in)            :: header
+      integer(psb_ipk_), intent(out)          :: info
+      integer(psb_ipk_), optional, intent(in) :: iunit
+      character(len=*), optional, intent(in)  :: filename
+    end subroutine mm_dvect_write
+  end interface
+
+#if ! defined(PSB_HAVE_BUGGY_GENERICS)
+  interface mm_vet_write
+    procedure mm_dvet1_write, mm_dvet2_write
+  end interface
+#endif
+
+  interface mm_mat_read
+    subroutine dmm_mat_read(a, info, iunit, filename)   
+      import :: psb_dspmat_type, psb_ipk_
+      implicit none
+      type(psb_dspmat_type), intent(out)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine dmm_mat_read
+    subroutine ldmm_mat_read(a, info, iunit, filename)   
+      import :: psb_ldspmat_type, psb_ipk_
+      implicit none
+      type(psb_ldspmat_type), intent(out)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine ldmm_mat_read
+  end interface
+
+  interface mm_mat_write
+    subroutine dmm_mat_write(a,mtitle,info,iunit,filename)
+      import :: psb_dspmat_type, psb_ipk_ 
+      implicit none
+      type(psb_dspmat_type), intent(in)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      character(len=*), intent(in) :: mtitle
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine dmm_mat_write
+    subroutine ldmm_mat_write(a,mtitle,info,iunit,filename)
+      import :: psb_ldspmat_type, psb_ipk_ 
+      implicit none
+      type(psb_ldspmat_type), intent(in)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      character(len=*), intent(in) :: mtitle
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine ldmm_mat_write
+  end interface
+#if 0
+#if ! defined(PSB_HAVE_BUGGY_GENERICS)
+  public mm_vet_read, mm_vet_write
+#endif
+#endif
+end module psb_d_mmio_mod
diff --git a/util/psb_i_mmio_mod.F90 b/util/psb_i_mmio_mod.F90
new file mode 100644
index 000000000..57fb503f3
--- /dev/null
+++ b/util/psb_i_mmio_mod.F90
@@ -0,0 +1,160 @@
+!   
+!                Parallel Sparse BLAS  version 3.5
+!      (C) Copyright 2006-2018
+!        Salvatore Filippone    
+!        Alfredo Buttari      
+!   
+!    Redistribution and use in source and binary forms, with or without
+!    modification, are permitted provided that the following conditions
+!    are met:
+!      1. Redistributions of source code must retain the above copyright
+!         notice, this list of conditions and the following disclaimer.
+!      2. Redistributions in binary form must reproduce the above copyright
+!         notice, this list of conditions, and the following disclaimer in the
+!         documentation and/or other materials provided with the distribution.
+!      3. The name of the PSBLAS group or the names of its contributors may
+!         not be used to endorse or promote products derived from this
+!         software without specific prior written permission.
+!   
+!    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+!    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
+!    TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+!    PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
+!    BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+!    CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+!    SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+!    INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+!    CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+!    ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+!    POSSIBILITY OF SUCH DAMAGE.
+!   
+!    
+module psb_i_mmio_mod
+
+
+  use psb_base_mod, only :  psb_ipk_, psb_lpk_,&
+       & psb_i_vect_type, psb_l_vect_type
+
+  !public mm_mat_read, mm_mat_write, mm_array_read, mm_array_write
+
+  interface mm_array_read
+    subroutine mm_ivet_read(b, info, iunit, filename)   
+      import :: psb_ipk_
+      implicit none
+      integer(psb_ipk_), allocatable, intent(out)  :: b(:)
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_ivet_read
+    subroutine mm_ivet2_read(b, info, iunit, filename)   
+      import :: psb_ipk_
+      implicit none
+      integer(psb_ipk_), allocatable, intent(out)  :: b(:,:)
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_ivet2_read
+    subroutine mm_ivect_read(b, info, iunit, filename)   
+      import :: psb_ipk_, psb_i_vect_type
+      implicit none
+      type(psb_i_vect_type), intent(inout)  :: b
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in) :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_ivect_read
+    subroutine mm_lvect_read(b, info, iunit, filename)   
+      import :: psb_ipk_, psb_l_vect_type
+      implicit none
+      type(psb_l_vect_type), intent(inout)  :: b
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in) :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_lvect_read
+#if defined(PSB_IPK4) && defined(PSB_LPK8) 
+    subroutine mm_lvet_read(b, info, iunit, filename)   
+      import :: psb_ipk_, psb_lpk_
+      implicit none
+      integer(psb_lpk_), allocatable, intent(out)  :: b(:)
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_lvet_read
+    subroutine mm_lvet2_read(b, info, iunit, filename)   
+      import :: psb_ipk_, psb_lpk_
+      implicit none
+      integer(psb_lpk_), allocatable, intent(out)  :: b(:,:)
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_lvet2_read
+  end interface mm_array_read
+#endif
+
+  interface mm_array_write
+    subroutine mm_ivet2_write(b, header, info, iunit, filename)   
+      import :: psb_ipk_
+      implicit none
+      integer(psb_ipk_), intent(in)  :: b(:,:)
+      character(len=*), intent(in) :: header
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_ivet2_write
+    subroutine mm_ivet1_write(b, header, info, iunit, filename)   
+      import :: psb_ipk_
+      implicit none
+      integer(psb_ipk_), intent(in)  :: b(:)
+      character(len=*), intent(in) :: header
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_ivet1_write
+#if defined(PSB_IPK4) && defined(PSB_LPK8) 
+    subroutine mm_lvet2_write(b, header, info, iunit, filename)   
+      import :: psb_ipk_, psb_lpk_
+      implicit none
+      integer(psb_lpk_), intent(in)  :: b(:,:)
+      character(len=*), intent(in) :: header
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_lvet2_write
+    subroutine mm_lvet1_write(b, header, info, iunit, filename)   
+      import :: psb_ipk_, psb_lpk_
+      implicit none
+      integer(psb_lpk_), intent(in)  :: b(:)
+      character(len=*), intent(in) :: header
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_lvet1_write
+    subroutine mm_ivect_write(b, header, info, iunit, filename)   
+      import :: psb_ipk_,psb_i_vect_type
+      implicit none
+      type(psb_i_vect_type), intent(inout)  :: b
+      character(len=*), intent(in) :: header
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_ivect_write
+    subroutine mm_lvect_write(b, header, info, iunit, filename)   
+      import :: psb_ipk_,psb_l_vect_type
+      implicit none
+      type(psb_l_vect_type), intent(inout)  :: b
+      character(len=*), intent(in) :: header
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_lvect_write
+  end interface mm_array_write
+#endif
+  interface mm_vet_write
+    procedure mm_ivet1_write, mm_ivet2_write
+  end interface
+
+#if 0
+#if ! defined(PSB_HAVE_BUGGY_GENERICS)
+  public mm_vet_read, mm_vet_write
+#endif
+#endif
+end module psb_i_mmio_mod
diff --git a/util/psb_mmio_mod.F90 b/util/psb_mmio_mod.F90
index bc0630c28..0eaf97db7 100644
--- a/util/psb_mmio_mod.F90
+++ b/util/psb_mmio_mod.F90
@@ -30,497 +30,9 @@
 !   
 !    
 module psb_mmio_mod
-
-
-  use psb_base_mod, only :  psb_ipk_, psb_lpk_, psb_spk_, psb_dpk_,&
-       & psb_s_vect_type, psb_d_vect_type, psb_i_vect_type, psb_l_vect_type,&
-       & psb_c_vect_type, psb_z_vect_type, &
-       & psb_sspmat_type, psb_cspmat_type, &
-       & psb_dspmat_type, psb_zspmat_type, &
-       & psb_lsspmat_type, psb_lcspmat_type, &
-       & psb_ldspmat_type, psb_lzspmat_type
-
-  public mm_mat_read, mm_mat_write, mm_array_read, mm_array_write
-
-#if ! defined(PSB_HAVE_BUGGY_GENERICS)
-  public mm_vet_read, mm_vet_write
-#endif
-
-  interface mm_array_read
-    subroutine mm_svet_read(b, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_
-      implicit none
-      real(psb_spk_), allocatable, intent(out)  :: b(:)
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_svet_read
-    subroutine mm_dvet_read(b, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_
-      implicit none
-      real(psb_dpk_), allocatable, intent(out)  :: b(:)
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_dvet_read
-    subroutine mm_cvet_read(b, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_
-      implicit none
-      complex(psb_spk_), allocatable, intent(out)  :: b(:)
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_cvet_read
-    subroutine mm_zvet_read(b, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_
-      implicit none
-      complex(psb_dpk_), allocatable, intent(out)  :: b(:)
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_zvet_read
-    subroutine mm_svet2_read(b, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_
-      implicit none
-      real(psb_spk_), allocatable, intent(out)  :: b(:,:)
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_svet2_read
-    subroutine mm_dvet2_read(b, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_
-      implicit none
-      real(psb_dpk_), allocatable, intent(out)  :: b(:,:)
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_dvet2_read
-    subroutine mm_cvet2_read(b, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_
-      implicit none
-      complex(psb_spk_), allocatable, intent(out)  :: b(:,:)
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_cvet2_read
-    subroutine mm_zvet2_read(b, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_
-      implicit none
-      complex(psb_dpk_), allocatable, intent(out)  :: b(:,:)
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_zvet2_read
-    subroutine mm_ivet_read(b, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_
-      implicit none
-      integer(psb_ipk_), allocatable, intent(out)  :: b(:)
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_ivet_read
-    subroutine mm_ivet2_read(b, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_
-      implicit none
-      integer(psb_ipk_), allocatable, intent(out)  :: b(:,:)
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_ivet2_read
-#if defined(PSB_IPK4) && defined(PSB_LPK8) 
-    subroutine mm_lvet_read(b, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_, psb_lpk_
-      implicit none
-      integer(psb_lpk_), allocatable, intent(out)  :: b(:)
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_lvet_read
-    subroutine mm_lvet2_read(b, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_, psb_lpk_
-      implicit none
-      integer(psb_lpk_), allocatable, intent(out)  :: b(:,:)
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_lvet2_read
-#endif
-    subroutine mm_svect_read(b, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_,psb_s_vect_type
-      implicit none
-      type(psb_s_vect_type), intent(inout)  :: b
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in) :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_svect_read
-    subroutine mm_dvect_read(b, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_,psb_d_vect_type
-      implicit none
-      type(psb_d_vect_type), intent(inout)    :: b
-      integer(psb_ipk_), intent(out)          :: info
-      integer(psb_ipk_), optional, intent(in) :: iunit
-      character(len=*), optional, intent(in)  :: filename
-    end subroutine mm_dvect_read
-    subroutine mm_cvect_read(b, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_,psb_c_vect_type
-      implicit none
-      type(psb_c_vect_type), intent(inout)    :: b
-      integer(psb_ipk_), intent(out)          :: info
-      integer(psb_ipk_), optional, intent(in) :: iunit
-      character(len=*), optional, intent(in)  :: filename
-    end subroutine mm_cvect_read
-    subroutine mm_zvect_read(b, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_,psb_z_vect_type
-      implicit none
-      type(psb_z_vect_type), intent(inout)    :: b
-      integer(psb_ipk_), intent(out)          :: info
-      integer(psb_ipk_), optional, intent(in) :: iunit
-      character(len=*), optional, intent(in)  :: filename
-    end subroutine mm_zvect_read
-    subroutine mm_ivect_read(b, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_,psb_i_vect_type
-      implicit none
-      type(psb_i_vect_type), intent(inout)  :: b
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in) :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_ivect_read
-    subroutine mm_lvect_read(b, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_,psb_l_vect_type
-      implicit none
-      type(psb_l_vect_type), intent(inout)  :: b
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in) :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_lvect_read
-  end interface
-
-
-#if ! defined(PSB_HAVE_BUGGY_GENERICS)
-  interface mm_vet_read
-    procedure mm_svet_read, mm_dvet_read, mm_cvet_read,&
-         & mm_zvet_read, mm_svet2_read, mm_dvet2_read, &
-         & mm_cvet2_read, mm_zvet2_read
-  end interface
-#endif
-
-
-  interface mm_array_write
-    subroutine mm_svet2_write(b, header, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_
-      implicit none
-      real(psb_spk_), intent(in)  :: b(:,:)
-      character(len=*), intent(in) :: header
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_svet2_write
-    subroutine mm_svet1_write(b, header, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_
-      implicit none
-      real(psb_spk_), intent(in)  :: b(:)
-      character(len=*), intent(in) :: header
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_svet1_write
-    subroutine mm_dvet2_write(b, header, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_
-      implicit none
-      real(psb_dpk_), intent(in)  :: b(:,:)
-      character(len=*), intent(in) :: header
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_dvet2_write
-    subroutine mm_dvet1_write(b, header, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_
-      implicit none
-      real(psb_dpk_), intent(in)  :: b(:)
-      character(len=*), intent(in) :: header
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_dvet1_write
-    subroutine mm_cvet2_write(b, header, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_
-      implicit none
-      complex(psb_spk_), intent(in)  :: b(:,:)
-      character(len=*), intent(in) :: header
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_cvet2_write
-    subroutine mm_cvet1_write(b, header, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_
-      implicit none
-      complex(psb_spk_), intent(in)  :: b(:)
-      character(len=*), intent(in) :: header
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_cvet1_write
-    subroutine mm_zvet2_write(b, header, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_
-      implicit none
-      complex(psb_dpk_), intent(in)  :: b(:,:)
-      character(len=*), intent(in) :: header
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_zvet2_write
-    subroutine mm_zvet1_write(b, header, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_
-      implicit none
-      complex(psb_dpk_), intent(in)  :: b(:)
-      character(len=*), intent(in) :: header
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_zvet1_write
-    subroutine mm_ivet2_write(b, header, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_
-      implicit none
-      integer(psb_ipk_), intent(in)  :: b(:,:)
-      character(len=*), intent(in) :: header
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_ivet2_write
-    subroutine mm_ivet1_write(b, header, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_
-      implicit none
-      integer(psb_ipk_), intent(in)  :: b(:)
-      character(len=*), intent(in) :: header
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_ivet1_write
-#if defined(PSB_IPK4) && defined(PSB_LPK8) 
-    subroutine mm_lvet2_write(b, header, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_, psb_lpk_
-      implicit none
-      integer(psb_lpk_), intent(in)  :: b(:,:)
-      character(len=*), intent(in) :: header
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_lvet2_write
-    subroutine mm_lvet1_write(b, header, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_, psb_lpk_
-      implicit none
-      integer(psb_lpk_), intent(in)  :: b(:)
-      character(len=*), intent(in) :: header
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_lvet1_write
-#endif
-    subroutine mm_svect_write(b, header, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_,psb_s_vect_type
-      implicit none
-      type(psb_s_vect_type), intent(inout)    :: b
-      character(len=*), intent(in)            :: header
-      integer(psb_ipk_), intent(out)          :: info
-      integer(psb_ipk_), optional, intent(in) :: iunit
-      character(len=*), optional, intent(in)  :: filename
-    end subroutine mm_svect_write  
-    subroutine mm_dvect_write(b, header, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_,psb_d_vect_type
-      implicit none
-      type(psb_d_vect_type), intent(inout)    :: b
-      character(len=*), intent(in)            :: header
-      integer(psb_ipk_), intent(out)          :: info
-      integer(psb_ipk_), optional, intent(in) :: iunit
-      character(len=*), optional, intent(in)  :: filename
-    end subroutine mm_dvect_write
-    subroutine mm_cvect_write(b, header, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_,psb_c_vect_type
-      implicit none
-      type(psb_c_vect_type), intent(inout)    :: b
-      character(len=*), intent(in)            :: header
-      integer(psb_ipk_), intent(out)          :: info
-      integer(psb_ipk_), optional, intent(in) :: iunit
-      character(len=*), optional, intent(in)  :: filename
-    end subroutine mm_cvect_write
-    subroutine mm_zvect_write(b, header, info, iunit, filename)   
-      import :: psb_dpk_, psb_ipk_,psb_z_vect_type
-      implicit none
-      type(psb_z_vect_type), intent(inout)    :: b
-      character(len=*), intent(in)            :: header
-      integer(psb_ipk_), intent(out)          :: info
-      integer(psb_ipk_), optional, intent(in) :: iunit
-      character(len=*), optional, intent(in)  :: filename
-    end subroutine mm_zvect_write
-    subroutine mm_ivect_write(b, header, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_,psb_i_vect_type
-      implicit none
-      type(psb_i_vect_type), intent(inout)  :: b
-      character(len=*), intent(in) :: header
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_ivect_write
-    subroutine mm_lvect_write(b, header, info, iunit, filename)   
-      import :: psb_spk_, psb_ipk_,psb_l_vect_type
-      implicit none
-      type(psb_l_vect_type), intent(inout)  :: b
-      character(len=*), intent(in) :: header
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine mm_lvect_write
-  end interface
-
-#if ! defined(PSB_HAVE_BUGGY_GENERICS)
-  interface mm_vet_write
-    procedure mm_svet1_write, mm_dvet1_write, mm_cvet1_write,&
-         & mm_zvet1_write, mm_svet2_write, mm_dvet2_write, &
-         & mm_cvet2_write, mm_zvet2_write, &
-         & mm_ivet1_write, mm_ivet2_write
-  end interface
-#endif
-
-  interface mm_mat_read
-    subroutine smm_mat_read(a, info, iunit, filename)   
-      import :: psb_sspmat_type, psb_ipk_
-      implicit none
-      type(psb_sspmat_type), intent(out)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine smm_mat_read
-    subroutine dmm_mat_read(a, info, iunit, filename)   
-      import :: psb_dspmat_type, psb_ipk_
-      implicit none
-      type(psb_dspmat_type), intent(out)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine dmm_mat_read
-    subroutine cmm_mat_read(a, info, iunit, filename)   
-      import :: psb_cspmat_type, psb_ipk_
-      implicit none
-      type(psb_cspmat_type), intent(out)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine cmm_mat_read
-    subroutine zmm_mat_read(a, info, iunit, filename)   
-      import :: psb_zspmat_type, psb_ipk_
-      implicit none
-      type(psb_zspmat_type), intent(out)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine zmm_mat_read
-    subroutine lsmm_mat_read(a, info, iunit, filename)   
-      import :: psb_lsspmat_type, psb_ipk_
-      implicit none
-      type(psb_lsspmat_type), intent(out)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine lsmm_mat_read
-    subroutine ldmm_mat_read(a, info, iunit, filename)   
-      import :: psb_ldspmat_type, psb_ipk_
-      implicit none
-      type(psb_ldspmat_type), intent(out)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine ldmm_mat_read
-    subroutine lcmm_mat_read(a, info, iunit, filename)   
-      import :: psb_lcspmat_type, psb_ipk_
-      implicit none
-      type(psb_lcspmat_type), intent(out)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine lcmm_mat_read
-    subroutine lzmm_mat_read(a, info, iunit, filename)   
-      import :: psb_lzspmat_type, psb_ipk_
-      implicit none
-      type(psb_lzspmat_type), intent(out)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine lzmm_mat_read
-  end interface
-
-  interface mm_mat_write
-    subroutine smm_mat_write(a,mtitle,info,iunit,filename)
-      import :: psb_sspmat_type, psb_ipk_ 
-      implicit none
-      type(psb_sspmat_type), intent(in)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      character(len=*), intent(in) :: mtitle
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine smm_mat_write
-    subroutine dmm_mat_write(a,mtitle,info,iunit,filename)
-      import :: psb_dspmat_type, psb_ipk_
-      implicit none
-      type(psb_dspmat_type), intent(in)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      character(len=*), intent(in) :: mtitle
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine dmm_mat_write
-    subroutine cmm_mat_write(a,mtitle,info,iunit,filename)
-      import :: psb_cspmat_type, psb_ipk_
-      implicit none
-      type(psb_cspmat_type), intent(in)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      character(len=*), intent(in) :: mtitle
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine cmm_mat_write
-    subroutine zmm_mat_write(a,mtitle,info,iunit,filename)
-      import :: psb_zspmat_type, psb_ipk_
-      implicit none
-      type(psb_zspmat_type), intent(in)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      character(len=*), intent(in) :: mtitle
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine zmm_mat_write
-    subroutine lsmm_mat_write(a,mtitle,info,iunit,filename)
-      import :: psb_lsspmat_type, psb_ipk_ 
-      implicit none
-      type(psb_lsspmat_type), intent(in)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      character(len=*), intent(in) :: mtitle
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine lsmm_mat_write
-    subroutine ldmm_mat_write(a,mtitle,info,iunit,filename)
-      import :: psb_ldspmat_type, psb_ipk_
-      implicit none
-      type(psb_ldspmat_type), intent(in)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      character(len=*), intent(in) :: mtitle
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine ldmm_mat_write
-    subroutine lcmm_mat_write(a,mtitle,info,iunit,filename)
-      import :: psb_lcspmat_type, psb_ipk_
-      implicit none
-      type(psb_lcspmat_type), intent(in)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      character(len=*), intent(in) :: mtitle
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine lcmm_mat_write
-    subroutine lzmm_mat_write(a,mtitle,info,iunit,filename)
-      import :: psb_lzspmat_type, psb_ipk_
-      implicit none
-      type(psb_lzspmat_type), intent(in)  :: a
-      integer(psb_ipk_), intent(out)        :: info
-      character(len=*), intent(in) :: mtitle
-      integer(psb_ipk_), optional, intent(in)          :: iunit
-      character(len=*), optional, intent(in) :: filename
-    end subroutine lzmm_mat_write
-  end interface
-
+  use psb_i_mmio_mod
+  use psb_s_mmio_mod
+  use psb_d_mmio_mod
+  use psb_c_mmio_mod
+  use psb_z_mmio_mod
 end module psb_mmio_mod
diff --git a/util/psb_s_mmio_mod.F90 b/util/psb_s_mmio_mod.F90
new file mode 100644
index 000000000..3eb4cf521
--- /dev/null
+++ b/util/psb_s_mmio_mod.F90
@@ -0,0 +1,155 @@
+!   
+!                Parallel Sparse BLAS  version 3.5
+!      (C) Copyright 2006-2018
+!        Salvatore Filippone    
+!        Alfredo Buttari      
+!   
+!    Redistribution and use in source and binary forms, with or without
+!    modification, are permitted provided that the following conditions
+!    are met:
+!      1. Redistributions of source code must retain the above copyright
+!         notice, this list of conditions and the following disclaimer.
+!      2. Redistributions in binary form must reproduce the above copyright
+!         notice, this list of conditions, and the following disclaimer in the
+!         documentation and/or other materials provided with the distribution.
+!      3. The name of the PSBLAS group or the names of its contributors may
+!         not be used to endorse or promote products derived from this
+!         software without specific prior written permission.
+!   
+!    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+!    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
+!    TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+!    PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
+!    BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+!    CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+!    SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+!    INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+!    CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+!    ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+!    POSSIBILITY OF SUCH DAMAGE.
+!   
+!    
+module psb_s_mmio_mod
+
+
+  use psb_base_mod, only :  psb_ipk_, psb_lpk_, psb_spk_,&
+       & psb_s_vect_type, psb_sspmat_type, psb_lsspmat_type
+
+  !public mm_mat_read, mm_mat_write, mm_array_read, mm_array_write
+  interface mm_array_read
+    subroutine mm_svet_read(b, info, iunit, filename)   
+      import :: psb_spk_, psb_ipk_
+      implicit none
+      real(psb_spk_), allocatable, intent(out)  :: b(:)
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_svet_read
+    subroutine mm_svet2_read(b, info, iunit, filename)   
+      import :: psb_spk_, psb_ipk_
+      implicit none
+      real(psb_spk_), allocatable, intent(out)  :: b(:,:)
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_svet2_read
+    subroutine mm_svect_read(b, info, iunit, filename)   
+      import :: psb_ipk_,psb_s_vect_type
+      implicit none
+      type(psb_s_vect_type), intent(inout)  :: b
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in) :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_svect_read
+  end interface mm_array_read
+
+
+#if 0
+#if ! defined(PSB_HAVE_BUGGY_GENERICS)
+  interface mm_vet_read
+    procedure mm_svet_read, mm_svet2_read
+  end interface
+#endif
+#endif
+  
+  interface mm_array_write
+    subroutine mm_svet2_write(b, header, info, iunit, filename)   
+      import :: psb_spk_, psb_ipk_
+      implicit none
+      real(psb_spk_), intent(in)  :: b(:,:)
+      character(len=*), intent(in) :: header
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_svet2_write
+    subroutine mm_svet1_write(b, header, info, iunit, filename)   
+      import :: psb_spk_, psb_ipk_
+      implicit none
+      real(psb_spk_), intent(in)  :: b(:)
+      character(len=*), intent(in) :: header
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_svet1_write
+    subroutine mm_svect_write(b, header, info, iunit, filename)   
+      import :: psb_ipk_, psb_s_vect_type
+      implicit none
+      type(psb_s_vect_type), intent(inout)    :: b
+      character(len=*), intent(in)            :: header
+      integer(psb_ipk_), intent(out)          :: info
+      integer(psb_ipk_), optional, intent(in) :: iunit
+      character(len=*), optional, intent(in)  :: filename
+    end subroutine mm_svect_write
+  end interface
+
+#if ! defined(PSB_HAVE_BUGGY_GENERICS)
+  interface mm_vet_write
+    procedure mm_svet1_write, mm_svet2_write
+  end interface
+#endif
+
+  interface mm_mat_read
+    subroutine smm_mat_read(a, info, iunit, filename)   
+      import :: psb_sspmat_type, psb_ipk_
+      implicit none
+      type(psb_sspmat_type), intent(out)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine smm_mat_read
+    subroutine lsmm_mat_read(a, info, iunit, filename)   
+      import :: psb_lsspmat_type, psb_ipk_
+      implicit none
+      type(psb_lsspmat_type), intent(out)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine lsmm_mat_read
+  end interface
+
+  interface mm_mat_write
+    subroutine smm_mat_write(a,mtitle,info,iunit,filename)
+      import :: psb_sspmat_type, psb_ipk_ 
+      implicit none
+      type(psb_sspmat_type), intent(in)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      character(len=*), intent(in) :: mtitle
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine smm_mat_write
+    subroutine lsmm_mat_write(a,mtitle,info,iunit,filename)
+      import :: psb_lsspmat_type, psb_ipk_ 
+      implicit none
+      type(psb_lsspmat_type), intent(in)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      character(len=*), intent(in) :: mtitle
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine lsmm_mat_write
+  end interface
+#if 0
+#if ! defined(PSB_HAVE_BUGGY_GENERICS)
+  public mm_vet_read, mm_vet_write
+#endif
+#endif
+end module psb_s_mmio_mod
diff --git a/util/psb_z_mmio_mod.F90 b/util/psb_z_mmio_mod.F90
new file mode 100644
index 000000000..e03d029da
--- /dev/null
+++ b/util/psb_z_mmio_mod.F90
@@ -0,0 +1,155 @@
+!   
+!                Parallel Sparse BLAS  version 3.5
+!      (C) Copyright 2006-2018
+!        Salvatore Filippone    
+!        Alfredo Buttari      
+!   
+!    Redistribution and use in source and binary forms, with or without
+!    modification, are permitted provided that the following conditions
+!    are met:
+!      1. Redistributions of source code must retain the above copyright
+!         notice, this list of conditions and the following disclaimer.
+!      2. Redistributions in binary form must reproduce the above copyright
+!         notice, this list of conditions, and the following disclaimer in the
+!         documentation and/or other materials provided with the distribution.
+!      3. The name of the PSBLAS group or the names of its contributors may
+!         not be used to endorse or promote products derived from this
+!         software without specific prior written permission.
+!   
+!    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+!    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
+!    TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+!    PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
+!    BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+!    CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+!    SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+!    INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+!    CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+!    ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+!    POSSIBILITY OF SUCH DAMAGE.
+!   
+!    
+module psb_z_mmio_mod
+
+
+  use psb_base_mod, only :  psb_ipk_, psb_lpk_, psb_dpk_,&
+       & psb_z_vect_type, psb_zspmat_type, psb_lzspmat_type
+
+  !public mm_mat_read, mm_mat_write, mm_array_read, mm_array_write
+  interface mm_array_read
+    subroutine mm_zvet_read(b, info, iunit, filename)   
+      import :: psb_dpk_, psb_ipk_
+      implicit none
+      complex(psb_dpk_), allocatable, intent(out)  :: b(:)
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_zvet_read
+    subroutine mm_zvet2_read(b, info, iunit, filename)   
+      import :: psb_dpk_, psb_ipk_
+      implicit none
+      complex(psb_dpk_), allocatable, intent(out)  :: b(:,:)
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_zvet2_read
+    subroutine mm_zvect_read(b, info, iunit, filename)   
+      import :: psb_ipk_,psb_z_vect_type
+      implicit none
+      type(psb_z_vect_type), intent(inout)  :: b
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in) :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_zvect_read
+  end interface mm_array_read
+
+
+#if 0
+#if ! defined(PSB_HAVE_BUGGY_GENERICS)
+  interface mm_vet_read
+    procedure mm_zvet_read, mm_zvet2_read
+  end interface
+#endif
+#endif
+  
+  interface mm_array_write
+    subroutine mm_zvet2_write(b, header, info, iunit, filename)   
+      import :: psb_dpk_, psb_ipk_
+      implicit none
+      complex(psb_dpk_), intent(in)  :: b(:,:)
+      character(len=*), intent(in) :: header
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_zvet2_write
+    subroutine mm_zvet1_write(b, header, info, iunit, filename)   
+      import :: psb_dpk_, psb_ipk_
+      implicit none
+      complex(psb_dpk_), intent(in)  :: b(:)
+      character(len=*), intent(in) :: header
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine mm_zvet1_write
+    subroutine mm_zvect_write(b, header, info, iunit, filename)   
+      import :: psb_ipk_, psb_z_vect_type
+      implicit none
+      type(psb_z_vect_type), intent(inout)    :: b
+      character(len=*), intent(in)            :: header
+      integer(psb_ipk_), intent(out)          :: info
+      integer(psb_ipk_), optional, intent(in) :: iunit
+      character(len=*), optional, intent(in)  :: filename
+    end subroutine mm_zvect_write
+  end interface
+
+#if ! defined(PSB_HAVE_BUGGY_GENERICS)
+  interface mm_vet_write
+    procedure mm_zvet1_write, mm_zvet2_write
+  end interface
+#endif
+
+  interface mm_mat_read
+    subroutine zmm_mat_read(a, info, iunit, filename)   
+      import :: psb_zspmat_type, psb_ipk_
+      implicit none
+      type(psb_zspmat_type), intent(out)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine zmm_mat_read
+    subroutine lzmm_mat_read(a, info, iunit, filename)   
+      import :: psb_lzspmat_type, psb_ipk_
+      implicit none
+      type(psb_lzspmat_type), intent(out)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine lzmm_mat_read
+  end interface
+
+  interface mm_mat_write
+    subroutine zmm_mat_write(a,mtitle,info,iunit,filename)
+      import :: psb_zspmat_type, psb_ipk_ 
+      implicit none
+      type(psb_zspmat_type), intent(in)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      character(len=*), intent(in) :: mtitle
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine zmm_mat_write
+    subroutine lzmm_mat_write(a,mtitle,info,iunit,filename)
+      import :: psb_lzspmat_type, psb_ipk_ 
+      implicit none
+      type(psb_lzspmat_type), intent(in)  :: a
+      integer(psb_ipk_), intent(out)        :: info
+      character(len=*), intent(in) :: mtitle
+      integer(psb_ipk_), optional, intent(in)          :: iunit
+      character(len=*), optional, intent(in) :: filename
+    end subroutine lzmm_mat_write
+  end interface
+#if 0
+#if ! defined(PSB_HAVE_BUGGY_GENERICS)
+  public mm_vet_read, mm_vet_write
+#endif
+#endif
+end module psb_z_mmio_mod