mld2p4-dev:

Fixes for possible empty local row set.
16 years ago · f51126b977
parent d8284f038a
commit f51126b977
6 changed files with 6 additions and 6 deletions
--- a/examples/pdegen/mld_dexample_1lev.f90
+++ b/examples/pdegen/mld_dexample_1lev.f90
@ -322,7 +322,7 @@ contains
    ! Using a simple BLOCK distribution.
    !
    nt = (m+np-1)/np
-    nr = min(nt,m-(iam*nt))
+    nr = max(0,min(nt,m-(iam*nt)))

    nt = nr
    call psb_sum(ictxt,nt) 
--- a/examples/pdegen/mld_dexample_ml.f90
+++ b/examples/pdegen/mld_dexample_ml.f90
@ -360,7 +360,7 @@ contains
    ! Using a simple BLOCK distribution.
    !
    nt = (m+np-1)/np
-    nr = min(nt,m-(iam*nt))
+    nr = max(0,min(nt,m-(iam*nt)))

    nt = nr
    call psb_sum(ictxt,nt) 
--- a/examples/pdegen/mld_sexample_1lev.f90
+++ b/examples/pdegen/mld_sexample_1lev.f90
@ -323,7 +323,7 @@ contains
    ! Using a simple BLOCK distribution.
    !
    nt = (m+np-1)/np
-    nr = min(nt,m-(iam*nt))
+    nr = max(0,min(nt,m-(iam*nt)))

    nt = nr
    call psb_sum(ictxt,nt) 
--- a/examples/pdegen/mld_sexample_ml.f90
+++ b/examples/pdegen/mld_sexample_ml.f90
@ -361,7 +361,7 @@ contains
    ! Using a simple BLOCK distribution.
    !
    nt = (m+np-1)/np
-    nr = min(nt,m-(iam*nt))
+    nr = max(0,min(nt,m-(iam*nt)))

    nt = nr
    call psb_sum(ictxt,nt) 
--- a/tests/pdegen/ppde.f90
+++ b/tests/pdegen/ppde.f90
@ -456,7 +456,7 @@ contains
    ! Using a simple BLOCK distribution.
    !
    nt = (m+np-1)/np
-    nr = min(nt,m-(iam*nt))
+    nr = max(0,min(nt,m-(iam*nt)))

    nt = nr
    call psb_sum(ictxt,nt) 
--- a/tests/pdegen/spde.f90
+++ b/tests/pdegen/spde.f90
@ -456,7 +456,7 @@ contains
    ! Using a simple BLOCK distribution.
    !
    nt = (m+np-1)/np
-    nr = min(nt,m-(iam*nt))
+    nr = max(0,min(nt,m-(iam*nt)))

    nt = nr
    call psb_sum(ictxt,nt)