home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / fortran / library / ssp / matops / cadd.for next >
Encoding:
Text File  |  1985-11-29  |  1.6 KB  |  60 lines
  1. C
  2. C     ..................................................................
  3. C
  4. C        SUBROUTINE CADD
  5. C
  6. C        PURPOSE
  7. C           ADD COLUMN OF ONE MATRIX TO COLUMN OF ANOTHER MATRIX
  8. C
  9. C        USAGE
  10. C           CALL CADD(A,ICA,R,ICR,N,M,MS,L)
  11. C
  12. C        DESCRIPTION OF PARAMETERS
  13. C           A   - NAME OF INPUT MATRIX
  14. C           ICA - COLUMN IN MATRIX A TO BE ADDED TO COLUMN ICR OF R
  15. C           R   - NAME OF OUTPUT MATRIX
  16. C           ICR - COLUMN IN MATRIX R WHERE SUMMATION IS DEVELOPED
  17. C           N   - NUMBER OF ROWS IN A AND R
  18. C           M   - NUMBER OF COLUMNS IN A
  19. C           MS  - ONE DIGIT NUMBER FOR STORAGE MODE OF MATRIX A
  20. C                  0 - GENERAL
  21. C                  1 - SYMMETRIC
  22. C                  2 - DIAGONAL
  23. C           L   - NUMBER OF COLUMNS IN R
  24. C
  25. C        REMARKS
  26. C           MATRIX R MUST BE A GENERAL MATRIX
  27. C           MATRIX R CANNOT BE IN THE SAME LOCATION AS MATRIX A UNLESS
  28. C           A IS GENERAL
  29. C
  30. C        SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
  31. C           LOC
  32. C
  33. C        METHOD
  34. C           EACH ELEMENT OF COLUMN ICA OF MATRIX A IS ADDED TO
  35. C           CORRESPONDING ELEMENT OF COLUMN ICR OF MATRIX R
  36. C
  37. C     ..................................................................
  38. C
  39.       SUBROUTINE CADD(A,ICA,R,ICR,N,M,MS,L)
  40.       DIMENSION A(1),R(1)
  41. C
  42.       IR=N*(ICR-1)
  43.       DO 2 I=1,N
  44.       IR=IR+1
  45. C
  46. C        LOCATE INPUT ELEMENT FOR ANY MATRIX STORAGE MODE
  47. C
  48.       CALL LOC(I,ICA,IA,N,M,MS)
  49. C
  50. C        TEST FOR ZERO ELEMENT IN DIAGONAL MATRIX
  51. C
  52.       IF(IA) 1,2,1
  53. C
  54. C        ADD ELEMENTS
  55. C
  56.     1 R(IR)=R(IR)+A(IA)
  57.     2 CONTINUE
  58.       RETURN
  59.       END
  60.