/algorithms_for_tapenade/fortran_code/dense_inverse.f

http://github.com/b45ch1/hpsc_hanoi_2009_walter · FORTRAN Legacy · 25 lines · 21 code · 3 blank · 1 comment · 0 complexity · 80de0a4491d8840de49338f8f08f9dec MD5 · raw file 

    

  1. subroutine inv(A,QT,R,NA)
    2. 

  2. 	implicit none
    3. 

  3. 	integer*4 NA
    4. 

  4. 	real*8 A(NA*NA), QT(NA*NA), R(NA*NA)
    5. 

  5. 	integer*4 n,m,k
    6. 

  6. 	real*8 tmp
    7. 

  7. 	call qr(A,QT,R,NA)
    8. 

  8. 	
    9. 

  9. 	do n= NA-1,0,-1
    10. 

  10. 		tmp = R(n*NA + n + 1)
    11. 

  11. 		do m = 0, NA-1
    12. 

  12. 			R (n*NA + m + 1) = R (n*NA + m + 1)/ tmp
    13. 

  13. 			QT(n*NA + m + 1) = QT(n*NA + m + 1)/ tmp
    14. 

  14. 		enddo
    15. 

  15. 		do m=n+1,NA-1
    16. 

  16. 			tmp = R(n*NA + m + 1)
    17. 

  17. 			R( n*NA + m + 1) = 0
    18. 

  18. 			do k=0,NA-1
    19. 

  19. 				QT(n*NA + k + 1) = QT(n*NA + k + 1) - QT( m*NA + k + 1) * tmp
    20. 

  20. 			enddo
    21. 

  21. 		enddo
    22. 

  22. 	enddo
    23. 

  23. 	
    24. 

  24. 	end
    25.