Fix jac%pa handling

mlprec/impl/smoother/mld_c_jac_smoother_bld.f90

@@ -76,7 +76,7 @@ subroutine mld_c_jac_smoother_bld(a,desc_a,sm,info,amold,vmold,imold)
   select type (smsv => sm%sv)
   class is (mld_c_diag_solver_type)
     call sm%nd%free()
-      sm%nd_nnz_tot = nztota
+    sm%nd_nnz_tot = nztota
     call psb_sum(ictxt,sm%nd_nnz_tot)
     call sm%sv%build(a,desc_a,info,amold=amold,vmold=vmold)
 

mlprec/impl/smoother/mld_c_l1_jac_smoother_bld.f90

@@ -72,12 +72,11 @@ subroutine mld_c_l1_jac_smoother_bld(a,desc_a,sm,info,amold,vmold,imold)
   nrow_a = a%get_nrows()
   nztota = a%get_nzeros()
 
-  if( sm%checkres ) sm%pa => a
+  sm%pa => a
 
   select type (smsv => sm%sv)
   class is (mld_c_diag_solver_type)
     call sm%nd%free()
-    sm%pa => a
     sm%nd_nnz_tot = nztota
     
     call psb_sum(ictxt,sm%nd_nnz_tot)

mlprec/impl/smoother/mld_d_jac_smoother_bld.f90

@@ -76,7 +76,7 @@ subroutine mld_d_jac_smoother_bld(a,desc_a,sm,info,amold,vmold,imold)
   select type (smsv => sm%sv)
   class is (mld_d_diag_solver_type)
     call sm%nd%free()
-      sm%nd_nnz_tot = nztota
+    sm%nd_nnz_tot = nztota
     call psb_sum(ictxt,sm%nd_nnz_tot)
     call sm%sv%build(a,desc_a,info,amold=amold,vmold=vmold)
 

mlprec/impl/smoother/mld_d_l1_jac_smoother_bld.f90

@@ -72,12 +72,11 @@ subroutine mld_d_l1_jac_smoother_bld(a,desc_a,sm,info,amold,vmold,imold)
   nrow_a = a%get_nrows()
   nztota = a%get_nzeros()
 
-  if( sm%checkres ) sm%pa => a
+  sm%pa => a
 
   select type (smsv => sm%sv)
   class is (mld_d_diag_solver_type)
     call sm%nd%free()
-    sm%pa => a
     sm%nd_nnz_tot = nztota
     
     call psb_sum(ictxt,sm%nd_nnz_tot)

mlprec/impl/smoother/mld_s_jac_smoother_bld.f90

@@ -76,7 +76,7 @@ subroutine mld_s_jac_smoother_bld(a,desc_a,sm,info,amold,vmold,imold)
   select type (smsv => sm%sv)
   class is (mld_s_diag_solver_type)
     call sm%nd%free()
-      sm%nd_nnz_tot = nztota
+    sm%nd_nnz_tot = nztota
     call psb_sum(ictxt,sm%nd_nnz_tot)
     call sm%sv%build(a,desc_a,info,amold=amold,vmold=vmold)
 

mlprec/impl/smoother/mld_s_l1_jac_smoother_bld.f90

@@ -72,12 +72,11 @@ subroutine mld_s_l1_jac_smoother_bld(a,desc_a,sm,info,amold,vmold,imold)
   nrow_a = a%get_nrows()
   nztota = a%get_nzeros()
 
-  if( sm%checkres ) sm%pa => a
+  sm%pa => a
 
   select type (smsv => sm%sv)
   class is (mld_s_diag_solver_type)
     call sm%nd%free()
-    sm%pa => a
     sm%nd_nnz_tot = nztota
     
     call psb_sum(ictxt,sm%nd_nnz_tot)

mlprec/impl/smoother/mld_z_jac_smoother_bld.f90

@@ -76,7 +76,7 @@ subroutine mld_z_jac_smoother_bld(a,desc_a,sm,info,amold,vmold,imold)
   select type (smsv => sm%sv)
   class is (mld_z_diag_solver_type)
     call sm%nd%free()
-      sm%nd_nnz_tot = nztota
+    sm%nd_nnz_tot = nztota
     call psb_sum(ictxt,sm%nd_nnz_tot)
     call sm%sv%build(a,desc_a,info,amold=amold,vmold=vmold)
 

mlprec/impl/smoother/mld_z_l1_jac_smoother_bld.f90

@@ -72,12 +72,11 @@ subroutine mld_z_l1_jac_smoother_bld(a,desc_a,sm,info,amold,vmold,imold)
   nrow_a = a%get_nrows()
   nztota = a%get_nzeros()
 
-  if( sm%checkres ) sm%pa => a
+  sm%pa => a
 
   select type (smsv => sm%sv)
   class is (mld_z_diag_solver_type)
     call sm%nd%free()
-    sm%pa => a
     sm%nd_nnz_tot = nztota
     
     call psb_sum(ictxt,sm%nd_nnz_tot)