From b76b6bb0e0dc1231cbfee7024c07b7bfae8f0d14 Mon Sep 17 00:00:00 2001
From: Alfredo Buttari <alfredo.buttari@enseeiht.fr>
Date: Mon, 5 Dec 2005 12:23:27 +0000
Subject: [PATCH] *** empty log message ***

---
 src/internals/psi_crea_index.f90 |  6 ++++--
 src/prec/psb_dbldaggrmat.f90     | 18 +++++++++++++-----
 2 files changed, 17 insertions(+), 7 deletions(-)

diff --git a/src/internals/psi_crea_index.f90 b/src/internals/psi_crea_index.f90
index e21ad1ce..72e61c23 100644
--- a/src/internals/psi_crea_index.f90
+++ b/src/internals/psi_crea_index.f90
@@ -72,9 +72,11 @@ subroutine psi_crea_index(desc_a,index_in,index_out,glob_idx,info)
   endif
 
   ! allocate dependency list
-!  call psi_compute_size(desc_a%matrix_data, index_in, dl_lda, info)
-!  dl_lda=dl_lda+3
+  !  call psi_compute_size(desc_a%matrix_data, index_in, dl_lda, info)
+  ! This should be computed more efficiently to save space when
+  ! the number of processors becomes very high
   dl_lda=np+1
+
   allocate(dep_list(max(1,dl_lda),0:np),length_dl(0:np))
   ! ...extract dependence list (ordered list of identifer process
   !    which every process must communcate with...
diff --git a/src/prec/psb_dbldaggrmat.f90 b/src/prec/psb_dbldaggrmat.f90
index e0ebf493..b98131bd 100644
--- a/src/prec/psb_dbldaggrmat.f90
+++ b/src/prec/psb_dbldaggrmat.f90
@@ -382,11 +382,19 @@ contains
       goto 9999
     end if
 
-    where (p%dorig /= zero) 
-      p%dorig = one / p%dorig
-    elsewhere
-      p%dorig = one
-    end where
+    do i=1,size(p%dorig)
+       if (p%dorig(i) /= zero) then
+          p%dorig(i) = one / p%dorig(i)
+       else
+          p%dorig(i) = one
+       end if
+    end do
+
+!     where (p%dorig /= zero) 
+!       p%dorig = one / p%dorig
+!     elsewhere
+!       p%dorig = one
+!     end where
 
 
     ! 1. Allocate Ptilde in sparse matrix form