Hello,

I have nowhere seen a proposal for Systems of Linear Ecuations solving.

I know they can be solved with LU descomposition for quadratic systems.

Eg: 

Ax = b;

A = LU with LU transform;

LUx = b; note Ux = y, then 

Ly = b => y = LTRIS(L, b), then

Ux = y => x = UTRIS(U, y);

Note: LTRIS and UTRIS are methods for solving triangular lower/upper systems of ecuations.

LTRIS

x = b;

for i = 1 : n

for j = 1 : i - 1

x(i) = x(i) - L(i,j)x(j)

end

x(i) = x(i) / L(i,i);

end

UTRIS:

for i = n : -1: 1

for j = i + 1 : n

x(i) = x(i) - U(i,j)x(j)

end

x(i) = x(i) / U(i,i);

end

Can this thing be proposed in Boost uBLAS?

Thank you,

Iulian Calciu