Fix check on descriptor pointer

5 years ago · 5c2eef76ab
parent c9e1719098
commit 5c2eef76ab
4 changed files with 4 additions and 4 deletions
--- a/mlprec/impl/mld_c_hierarchy_bld.f90
+++ b/mlprec/impl/mld_c_hierarchy_bld.f90
@ -443,7 +443,7 @@ subroutine mld_c_hierarchy_bld(a,desc_a,prec,info)
    iszv = newsz
    ! Fix the pointers, but the level 1 should
    ! be treated differently
-    if (.not.associated(prec%precv(i)%base_desc,desc_a)) then
+    if (.not.associated(prec%precv(1)%base_desc,desc_a)) then
      prec%precv(1)%base_desc => prec%precv(1)%desc_ac
    end if
    do i=2, iszv 
--- a/mlprec/impl/mld_d_hierarchy_bld.f90
+++ b/mlprec/impl/mld_d_hierarchy_bld.f90
@ -443,7 +443,7 @@ subroutine mld_d_hierarchy_bld(a,desc_a,prec,info)
    iszv = newsz
    ! Fix the pointers, but the level 1 should
    ! be treated differently
-    if (.not.associated(prec%precv(i)%base_desc,desc_a)) then
+    if (.not.associated(prec%precv(1)%base_desc,desc_a)) then
      prec%precv(1)%base_desc => prec%precv(1)%desc_ac
    end if
    do i=2, iszv 
--- a/mlprec/impl/mld_s_hierarchy_bld.f90
+++ b/mlprec/impl/mld_s_hierarchy_bld.f90
@ -443,7 +443,7 @@ subroutine mld_s_hierarchy_bld(a,desc_a,prec,info)
    iszv = newsz
    ! Fix the pointers, but the level 1 should
    ! be treated differently
-    if (.not.associated(prec%precv(i)%base_desc,desc_a)) then
+    if (.not.associated(prec%precv(1)%base_desc,desc_a)) then
      prec%precv(1)%base_desc => prec%precv(1)%desc_ac
    end if
    do i=2, iszv 
--- a/mlprec/impl/mld_z_hierarchy_bld.f90
+++ b/mlprec/impl/mld_z_hierarchy_bld.f90
@ -443,7 +443,7 @@ subroutine mld_z_hierarchy_bld(a,desc_a,prec,info)
    iszv = newsz
    ! Fix the pointers, but the level 1 should
    ! be treated differently
-    if (.not.associated(prec%precv(i)%base_desc,desc_a)) then
+    if (.not.associated(prec%precv(1)%base_desc,desc_a)) then
      prec%precv(1)%base_desc => prec%precv(1)%desc_ac
    end if
    do i=2, iszv