mld2p4-2:

mlprec/impl/level/mld_c_base_onelev_cseti.F90
 mlprec/impl/level/mld_c_base_onelev_seti.F90
 mlprec/impl/level/mld_d_base_onelev_cseti.F90
 mlprec/impl/level/mld_d_base_onelev_seti.F90
 mlprec/impl/level/mld_s_base_onelev_cseti.F90
 mlprec/impl/level/mld_s_base_onelev_seti.F90
 mlprec/impl/level/mld_z_base_onelev_cseti.F90
 mlprec/impl/level/mld_z_base_onelev_seti.F90
 mlprec/mld_c_mumps_solver.F90
 mlprec/mld_d_mumps_solver.F90
 mlprec/mld_s_mumps_solver.F90
 mlprec/mld_z_mumps_solver.F90

Fixed default for second smoother.
get_fmt in mumps_solver.

mlprec/impl/level/mld_c_base_onelev_cseti.F90

@@ -132,7 +132,12 @@ subroutine mld_c_base_onelev_cseti(lv,what,val,info,pos)
       ! Do nothing and hope for the best :) 
       !
     end select
-    if (allocated(lv%sm)) call lv%sm%default()
+    if (ipos_==mld_pre_smooth_) then 
+      if (allocated(lv%sm)) call lv%sm%default()
+    else if (ipos_==mld_post_smooth_) then
+      if (allocated(lv%sm2a)) call lv%sm2a%default()
+    end if
+    
 
   case('SUB_SOLVE')
     select case (val) 

mlprec/impl/level/mld_c_base_onelev_seti.F90

@@ -133,7 +133,12 @@ subroutine mld_c_base_onelev_seti(lv,what,val,info,pos)
       ! Do nothing and hope for the best :) 
       !
     end select
-    if (allocated(lv%sm)) call lv%sm%default()
+    if (ipos_==mld_pre_smooth_) then 
+      if (allocated(lv%sm)) call lv%sm%default()
+    else if (ipos_==mld_post_smooth_) then
+      if (allocated(lv%sm2a)) call lv%sm2a%default()
+    end if
+    
 
   case(mld_sub_solve_)
     select case (val) 

mlprec/impl/level/mld_d_base_onelev_cseti.F90

@@ -138,7 +138,12 @@ subroutine mld_d_base_onelev_cseti(lv,what,val,info,pos)
       ! Do nothing and hope for the best :) 
       !
     end select
-    if (allocated(lv%sm)) call lv%sm%default()
+    if (ipos_==mld_pre_smooth_) then 
+      if (allocated(lv%sm)) call lv%sm%default()
+    else if (ipos_==mld_post_smooth_) then
+      if (allocated(lv%sm2a)) call lv%sm2a%default()
+    end if
+    
 
   case('SUB_SOLVE')
     select case (val) 

mlprec/impl/level/mld_d_base_onelev_seti.F90

@@ -139,7 +139,12 @@ subroutine mld_d_base_onelev_seti(lv,what,val,info,pos)
       ! Do nothing and hope for the best :) 
       !
     end select
-    if (allocated(lv%sm)) call lv%sm%default()
+    if (ipos_==mld_pre_smooth_) then 
+      if (allocated(lv%sm)) call lv%sm%default()
+    else if (ipos_==mld_post_smooth_) then
+      if (allocated(lv%sm2a)) call lv%sm2a%default()
+    end if
+    
 
   case(mld_sub_solve_)
     select case (val) 

mlprec/impl/level/mld_s_base_onelev_cseti.F90

@@ -132,7 +132,12 @@ subroutine mld_s_base_onelev_cseti(lv,what,val,info,pos)
       ! Do nothing and hope for the best :) 
       !
     end select
-    if (allocated(lv%sm)) call lv%sm%default()
+    if (ipos_==mld_pre_smooth_) then 
+      if (allocated(lv%sm)) call lv%sm%default()
+    else if (ipos_==mld_post_smooth_) then
+      if (allocated(lv%sm2a)) call lv%sm2a%default()
+    end if
+    
 
   case('SUB_SOLVE')
     select case (val) 

mlprec/impl/level/mld_s_base_onelev_seti.F90

@@ -133,7 +133,12 @@ subroutine mld_s_base_onelev_seti(lv,what,val,info,pos)
       ! Do nothing and hope for the best :) 
       !
     end select
-    if (allocated(lv%sm)) call lv%sm%default()
+    if (ipos_==mld_pre_smooth_) then 
+      if (allocated(lv%sm)) call lv%sm%default()
+    else if (ipos_==mld_post_smooth_) then
+      if (allocated(lv%sm2a)) call lv%sm2a%default()
+    end if
+    
 
   case(mld_sub_solve_)
     select case (val) 

mlprec/impl/level/mld_z_base_onelev_cseti.F90

@@ -138,7 +138,12 @@ subroutine mld_z_base_onelev_cseti(lv,what,val,info,pos)
       ! Do nothing and hope for the best :) 
       !
     end select
-    if (allocated(lv%sm)) call lv%sm%default()
+    if (ipos_==mld_pre_smooth_) then 
+      if (allocated(lv%sm)) call lv%sm%default()
+    else if (ipos_==mld_post_smooth_) then
+      if (allocated(lv%sm2a)) call lv%sm2a%default()
+    end if
+    
 
   case('SUB_SOLVE')
     select case (val) 

mlprec/impl/level/mld_z_base_onelev_seti.F90

@@ -139,7 +139,12 @@ subroutine mld_z_base_onelev_seti(lv,what,val,info,pos)
       ! Do nothing and hope for the best :) 
       !
     end select
-    if (allocated(lv%sm)) call lv%sm%default()
+    if (ipos_==mld_pre_smooth_) then 
+      if (allocated(lv%sm)) call lv%sm%default()
+    else if (ipos_==mld_post_smooth_) then
+      if (allocated(lv%sm2a)) call lv%sm2a%default()
+    end if
+    
 
   case(mld_sub_solve_)
     select case (val) 

mlprec/mld_c_mumps_solver.F90

@@ -80,6 +80,7 @@ module mld_c_mumps_solver
     procedure, pass(sv) :: cseti    =>c_mumps_solver_cseti
     procedure, pass(sv) :: csetr    => c_mumps_solver_csetr
     procedure, pass(sv) :: default  => c_mumps_solver_default
+    procedure, nopass   :: get_fmt  => c_mumps_get_fmt
 #if defined(HAVE_FINAL) 
 
     final               :: c_mumps_solver_finalize
@@ -477,5 +478,13 @@ contains
     return
   end function c_mumps_solver_sizeof
 #endif
+
+  function c_mumps_get_fmt() result(val)
+    implicit none 
+    character(len=32)  :: val
+
+    val = "MUMPS solver"
+  end function c_mumps_get_fmt
+
 end module mld_c_mumps_solver
 

mlprec/mld_d_mumps_solver.F90

@@ -80,6 +80,7 @@ module mld_d_mumps_solver
     procedure, pass(sv) :: cseti    =>d_mumps_solver_cseti
     procedure, pass(sv) :: csetr    => d_mumps_solver_csetr
     procedure, pass(sv) :: default  => d_mumps_solver_default
+    procedure, nopass   :: get_fmt  => d_mumps_get_fmt
 #if defined(HAVE_FINAL) 
 
     final               :: d_mumps_solver_finalize
@@ -477,5 +478,13 @@ contains
     return
   end function d_mumps_solver_sizeof
 #endif
+
+  function d_mumps_get_fmt() result(val)
+    implicit none 
+    character(len=32)  :: val
+
+    val = "MUMPS solver"
+  end function d_mumps_get_fmt
+
 end module mld_d_mumps_solver
 

mlprec/mld_s_mumps_solver.F90

@@ -80,6 +80,7 @@ module mld_s_mumps_solver
     procedure, pass(sv) :: cseti    =>s_mumps_solver_cseti
     procedure, pass(sv) :: csetr    => s_mumps_solver_csetr
     procedure, pass(sv) :: default  => s_mumps_solver_default
+    procedure, nopass   :: get_fmt  => s_mumps_get_fmt
 #if defined(HAVE_FINAL) 
 
     final               :: s_mumps_solver_finalize
@@ -477,5 +478,13 @@ contains
     return
   end function s_mumps_solver_sizeof
 #endif
+
+  function s_mumps_get_fmt() result(val)
+    implicit none 
+    character(len=32)  :: val
+
+    val = "MUMPS solver"
+  end function s_mumps_get_fmt
+
 end module mld_s_mumps_solver
 

mlprec/mld_z_mumps_solver.F90

@@ -80,6 +80,7 @@ module mld_z_mumps_solver
     procedure, pass(sv) :: cseti    =>z_mumps_solver_cseti
     procedure, pass(sv) :: csetr    => z_mumps_solver_csetr
     procedure, pass(sv) :: default  => z_mumps_solver_default
+    procedure, nopass   :: get_fmt  => z_mumps_get_fmt
 #if defined(HAVE_FINAL) 
 
     final               :: z_mumps_solver_finalize
@@ -477,5 +478,13 @@ contains
     return
   end function z_mumps_solver_sizeof
 #endif
+
+  function z_mumps_get_fmt() result(val)
+    implicit none 
+    character(len=32)  :: val
+
+    val = "MUMPS solver"
+  end function z_mumps_get_fmt
+
 end module mld_z_mumps_solver