Avoid integer overflow when computing nnz.

6 years ago · 532ad262b0
parent 5848b592bd
commit 532ad262b0
4 changed files with 4 additions and 4 deletions
--- a/tests/pdegen/mld_d_pde2d.f90
+++ b/tests/pdegen/mld_d_pde2d.f90
@ -247,7 +247,7 @@ contains

    m   = 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_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,7 +247,7 @@ contains

    m   = 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_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_)