*** empty log message ***

2026-10-20 12:52:14 +00:00 · 2005-12-05 12:23:27 +00:00
parent 1b7d51a008
commit b76b6bb0e0
2 changed files with 17 additions and 7 deletions
@@ -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