Added average for the (unlikely) case of an overlap in the base data distribution.

mlprec/mld_dprec_aply.f90

@@ -99,7 +99,7 @@ subroutine mld_dprec_aply(prec,x,y,desc_data,info,trans,work)
   call psb_info(ictxt, me, np)
 
   if (present(trans)) then 
-    trans_=trans
+    trans_=toupper(trans)
   else
     trans_='N'
   end if
@@ -139,6 +139,10 @@ subroutine mld_dprec_aply(prec,x,y,desc_data,info,trans,work)
     goto 9999
   endif
 
+  ! If the original distribution has an overlap we should fix that. 
+  call psb_ovrl(y,desc_data,info,update=psb_avg_)
+
+
   if (present(work)) then 
   else
     deallocate(work_)

mlprec/mld_zprec_aply.f90

@@ -99,7 +99,7 @@ subroutine mld_zprec_aply(prec,x,y,desc_data,info,trans,work)
   call psb_info(ictxt, me, np)
 
   if (present(trans)) then 
-    trans_=trans
+    trans_=toupper(trans)
   else
     trans_='N'
   end if
@@ -139,6 +139,9 @@ subroutine mld_zprec_aply(prec,x,y,desc_data,info,trans,work)
     goto 9999
   endif
 
+  ! If the original distribution has an overlap we should fix that. 
+  call psb_ovrl(y,desc_data,info,update=psb_avg_)
+
   if (present(work)) then 
   else
     deallocate(work_)