Help with LU factorization

General Usage
1

Is that possible to Implement ๐ฟ๐‘ˆ factorization without row pivoting ? can someone give me an example?

2

lu(A,Val(false))

1 Like
3

I donโ€™t want to use LinearAlgebra package , I want to implement the algorithm from scratch

4

Well, you can look at the source code for the Pivot == false case. But if this is a learning exercise (or a homework problem), you might want to try to implement it yourself from a textbook description.

2 Likes
5

Than you .Can you tell me the name of textbook that can describe how to implement this code ?

6

If this really is a homework assignment, you should do some of the work yourself. As for LA resources, check out e.g. http://ulaff.net/.

2 Likes
7

I tried it like that :

using LinearAlgebra
using SparseArrays
function LU2(A::Union{Array{real},Array{Float64}})
n = size(A,1); # Assuming A is nxn.
U = copy(A);
L = eye(eltype(A),n);
for k=2:n
for i=k:n
L[i,k-1] = U[i,k-1]/U[k-1,k-1];
for j=k-1:n
U[i,j] = U[i,j] - L[i,k-1]*U[k-1,j];
end
end
end
return L,U
end

but the problem now is when I test it i got this error ( UndefVarError: eye not defined)

8

You can get the n \times n identity matrix using Matrix{eltype(A)}(I,(n,n)).

1 Like
9

I got UndefVarError: eltpe not defined

10

Sorry, typo. Fixed now.

11

thanks it works

12

Deleted code

13

You probably want abs.(array).

14

I did it like this

pivot = indmax(abs.(U[(k-1):n,k-1])) + k - 2;

but i got UndefVarError: indmax not defined

15

The function that used to be called.indmax is now called findmax.

3 Likes
16

Thanks.it works but i got different error
MethodError: no method matching +(::Tuple{Float64,Int64}, ::Int64)

17

Did you try running the findmax function to see what kind of object it returns? And you can read the documentation of the function using ?findmax from the REPL.

18

Ah thereโ€™s also argmax.

1 Like
19

Yes i did. it works but when I test it on a matrix I got this strange error
( MethodError: no method matching +(::Tuple{Float64,Int64}, ::Int64) )
I tried to figure out the problem in the code again but i didnโ€™t note whatโ€™s the problem

20

findmax returns a tuple.
It is telling you that you are trying to add that tuple to an integer.

2 Likes