psblas3:

base/modules/psb_base_mat_mod.f03 base/serial/f03/psb_c_coo_impl.f03 base/serial/f03/psb_d_coo_impl.f03 base/serial/f03/psb_s_coo_impl.f03 base/serial/f03/psb_z_coo_impl.f03 Fixed idiotic bug in mv_coo_from_coo.
2026-10-20 21:02:26 +00:00 · 2009-10-15 13:25:01 +00:00
parent b74413cd69
commit 1c49f3ed5e
5 changed files with 3 additions and 28 deletions
@@ -483,10 +483,8 @@ contains
    a%unitd     = b%unitd
    a%upper     = b%upper
    a%sorted    = b%sorted
-    if (allocated(b%aux)) then 
-      allocate(a%aux(size(b%aux)))
-      a%aux(:)       = b%aux(:)
-    end if
+    a%aux       = b%aux
+
    return

  end subroutine base_cp_from
@@ -510,10 +508,7 @@ contains
    a%unitd     = b%unitd
    a%upper     = .not.b%upper
    a%sorted    = .false.
-    if (allocated(b%aux)) then 
-      allocate(a%aux(size(b%aux)))
-      a%aux(:)       = b%aux(:)
-    end if
+    a%aux       = b%aux
    
    return
    
@@ -2580,11 +2580,6 @@ subroutine c_mv_coo_from_coo_impl(a,b,info)
  call move_alloc(b%ja , a%ja   )
  call move_alloc(b%val, a%val )
  call b%free()
-
-  a%ia(:)  = b%ia(:)
-  a%ja(:)  = b%ja(:)
-  a%val(:) = b%val(:)
-
  call a%fix(info)

  if (info /= 0) goto 9999
@@ -2380,11 +2380,6 @@ subroutine d_mv_coo_from_coo_impl(a,b,info)
  call move_alloc(b%ja , a%ja   )
  call move_alloc(b%val, a%val )
  call b%free()
-
-  a%ia(:)  = b%ia(:)
-  a%ja(:)  = b%ja(:)
-  a%val(:) = b%val(:)
-
  call a%fix(info)

  if (info /= 0) goto 9999
@@ -2378,11 +2378,6 @@ subroutine s_mv_coo_from_coo_impl(a,b,info)
  call move_alloc(b%ja , a%ja   )
  call move_alloc(b%val, a%val )
  call b%free()
-
-  a%ia(:)  = b%ia(:)
-  a%ja(:)  = b%ja(:)
-  a%val(:) = b%val(:)
-
  call a%fix(info)

  if (info /= 0) goto 9999
@@ -2580,11 +2580,6 @@ subroutine z_mv_coo_from_coo_impl(a,b,info)
  call move_alloc(b%ja , a%ja   )
  call move_alloc(b%val, a%val )
  call b%free()
-
-  a%ia(:)  = b%ia(:)
-  a%ja(:)  = b%ja(:)
-  a%val(:) = b%val(:)
-
  call a%fix(info)

  if (info /= 0) goto 9999