Fix bounds computation in set_vect. Bug report by Alexandre Silva Lopes

6 years ago · 1a5ee8c46a
parent a87bf0ad7f
commit 1a5ee8c46a
5 changed files with 15 additions and 10 deletions
--- a/base/modules/serial/psb_c_base_vect_mod.f90
+++ b/base/modules/serial/psb_c_base_vect_mod.f90
@ -749,9 +749,10 @@ contains
        
    integer(psb_ipk_) :: info, first_, last_, nr

-    first_=1
-    last_=min(psb_size(x%v),size(val))
+
+    first_                     = 1
    if (present(first)) first_ = max(1,first)
+    last_                      = min(psb_size(x%v),first_+size(val)-1)
    if (present(last))  last_  = min(last,last_)

    if (allocated(x%v)) then 
--- a/base/modules/serial/psb_d_base_vect_mod.f90
+++ b/base/modules/serial/psb_d_base_vect_mod.f90
@ -749,9 +749,10 @@ contains
        
    integer(psb_ipk_) :: info, first_, last_, nr

-    first_=1
-    last_=min(psb_size(x%v),size(val))
+
+    first_                     = 1
    if (present(first)) first_ = max(1,first)
+    last_                      = min(psb_size(x%v),first_+size(val)-1)
    if (present(last))  last_  = min(last,last_)

    if (allocated(x%v)) then 
--- a/base/modules/serial/psb_i_base_vect_mod.f90
+++ b/base/modules/serial/psb_i_base_vect_mod.f90
@ -718,9 +718,10 @@ contains
        
    integer(psb_ipk_) :: info, first_, last_, nr

-    first_=1
-    last_=min(psb_size(x%v),size(val))
+
+    first_                     = 1
    if (present(first)) first_ = max(1,first)
+    last_                      = min(psb_size(x%v),first_+size(val)-1)
    if (present(last))  last_  = min(last,last_)

    if (allocated(x%v)) then 
--- a/base/modules/serial/psb_s_base_vect_mod.f90
+++ b/base/modules/serial/psb_s_base_vect_mod.f90
@ -749,9 +749,10 @@ contains
        
    integer(psb_ipk_) :: info, first_, last_, nr

-    first_=1
-    last_=min(psb_size(x%v),size(val))
+
+    first_                     = 1
    if (present(first)) first_ = max(1,first)
+    last_                      = min(psb_size(x%v),first_+size(val)-1)
    if (present(last))  last_  = min(last,last_)

    if (allocated(x%v)) then 
--- a/base/modules/serial/psb_z_base_vect_mod.f90
+++ b/base/modules/serial/psb_z_base_vect_mod.f90
@ -749,9 +749,10 @@ contains
        
    integer(psb_ipk_) :: info, first_, last_, nr

-    first_=1
-    last_=min(psb_size(x%v),size(val))
+
+    first_                     = 1
    if (present(first)) first_ = max(1,first)
+    last_                      = min(psb_size(x%v),first_+size(val)-1)
    if (present(last))  last_  = min(last,last_)

    if (allocated(x%v)) then