@@ -32,21 +32,18 @@ function searchmax_pivot!(piv_idx::AbstractArray{Int}, piv_val::AbstractArray{T}
3232 max_idx = argmax (abs .(A[:]))
3333 piv_idx[1 ] = offset+ max_idx
3434 piv_val[1 ] = A[max_idx]
35- println (" searchmax_pivot: " , piv_idx[1 ], " \n " , abs (piv_val[1 ]))
3635end
3736
3837function update_ipiv! (ipivl, piv_idx:: AbstractArray{Int} , piv_val:: AbstractArray{T} , k:: Int , nb:: Int ) where T
3938 max_piv_idx = argmax (abs .(piv_val))
4039 ipivl[1 ] = (max_piv_idx+ k- 2 )* nb + piv_idx[max_piv_idx]
41- println (" update_ipiv: " , ipivl[1 ])
4240end
4341
4442function swaprows_panel! (A:: AbstractArray{T} , M:: AbstractArray{T} , ipivl:: AbstractVector{Int} , m:: Int , p:: Int , nb:: Int ) where T
4543 q = div (ipivl[1 ]- 1 ,nb) + 1
4644 r = (ipivl[1 ]- 1 )% nb+ 1
4745 if m == q
4846 A[p,:], M[r,:] = M[r,:], A[p,:]
49- println (" swaprows_panel: " , imag .(A[p,:]), " \n " , imag .(M[r,:]))
5047 end
5148end
5249
@@ -62,7 +59,6 @@ function swaprows_trail!(A::AbstractArray{T}, M::AbstractArray{T}, ipiv::Abstrac
6259 r = (ipiv[p]- 1 )% nb+ 1
6360 if m == q
6461 A[p,:], M[r,:] = M[r,:], A[p,:]
65- println (" swaprows_trail: " , imag .(A[p,:]), " \n " , imag .(M[r,:]))
6662 end
6763 end
6864end
0 commit comments