sistlspf.f90

subroutine sistlspf(a, b, u, ip)

!Resolución del sistema triangular inferior, con unos en la diagonal, salvo la permutación de filas ip
!Entrada:	a (triangular inferior: L y triangular superior (sin la diagonal de 1): U)
!		b (términos independientes)
!		ip (puntero de las permutaciones)
!Salida:	u (solución de L*u=Pb)

implicit none
real(8), intent (in) :: a(:,:), b(:)
integer, intent (in) :: ip(:)
real(8), intent (out) :: u(:)
real(8) :: sum
integer :: i, j, n

n = size(b)

do i = 1, n
	sum = 0d0
	do j = 1, i-1
		sum = sum+a(ip(i), j)*u(j)
	end do
	u(i) = (b(ip(i))-sum)
end do

return

end subroutine sistlspf