Avoid integer overflow when computing nnz.

6 years ago · 331e57ed35
parent a0343aa36a
commit 331e57ed35
4 changed files with 6 additions and 6 deletions
--- a/tests/pdegen/mld_d_pde2d.f90
+++ b/tests/pdegen/mld_d_pde2d.f90
@ -247,8 +247,8 @@ contains

    m   = idim*idim
    n   = m
-    nnz = ((n*7)/(np))
-    if(iam == psb_root_) write(psb_out_unit,'("Generating Matrix (size=",i0,")...")')n
+    nnz = 7*((n+np-1)/np)
+    if (iam == psb_root_) write(psb_out_unit,'("Generating Matrix (size=",i0,")...")')n
    t0 = psb_wtime()
    select case(partition_)
    case(1)
--- a/tests/pdegen/mld_d_pde3d.f90
+++ b/tests/pdegen/mld_d_pde3d.f90
@ -263,7 +263,7 @@ contains
    
    m   = idim*idim*idim
    n   = m
-    nnz = ((n*7)/(np))
+    nnz = 7*((n+np-1)/np)
    if(iam == psb_root_) write(psb_out_unit,'("Generating Matrix (size=",i0,")...")')n
    t0 = psb_wtime()
    select case(partition_)
--- a/tests/pdegen/mld_s_pde2d.f90
+++ b/tests/pdegen/mld_s_pde2d.f90
@ -247,8 +247,8 @@ contains

    m   = idim*idim
    n   = m
-    nnz = ((n*7)/(np))
-    if(iam == psb_root_) write(psb_out_unit,'("Generating Matrix (size=",i0,")...")')n
+    nnz = 7*((n+np-1)/np)
+    if (iam == psb_root_) write(psb_out_unit,'("Generating Matrix (size=",i0,")...")')n
    t0 = psb_wtime()
    select case(partition_)
    case(1)
--- a/tests/pdegen/mld_s_pde3d.f90
+++ b/tests/pdegen/mld_s_pde3d.f90
@ -263,7 +263,7 @@ contains
    
    m   = idim*idim*idim
    n   = m
-    nnz = ((n*7)/(np))
+    nnz = 7*((n+np-1)/np)
    if(iam == psb_root_) write(psb_out_unit,'("Generating Matrix (size=",i0,")...")')n
    t0 = psb_wtime()
    select case(partition_)