axpym.f

Go to the documentation of this file.

C                  CVM Class Library
C                  http://cvmlib.com
C
C          Copyright Sergei Nikolaev 1992-2013
C Distributed under the Boost Software License, Version 1.0.
C    (See accompanying file LICENSE_1_0.txt or copy at
C          http://www.boost.org/LICENSE_1_0.txt)
C
C
C     Matrix axpy routines
C
C     The ?axpy routines perform a vector-vector operation defined as
C     y := a*x + y
C
C
C     Input/Output parameters:
C
C     M   - rows (int)(input)
C     N   - columns (int)(input)
C     A   - multiplier (real)(input)
C     X   - source matrix (real)(input)
C     LDX - leading dimesion of X (int)(input)
C     Y   - destination matrix (real)(output)
C     LDY - leading dimesion of Y (int)(input)

      SUBROUTINE saxpym (M, N, A, X, LDX, Y, LDY)
CDEC$ IF DEFINED (FTN_EXPORTS)
CDEC$     ATTRIBUTES DLLEXPORT::SAXPYM
CDEC$ ENDIF
      INTEGER m, n, ldx, ldy
      REAL a, x(ldx*n), y(ldy*n)
      INTEGER i

      IF (m .EQ. ldx .AND. m .EQ. ldy) THEN
          CALL saxpy(m * n, a, x, 1, y, 1)
      ELSE
          DO 10 i = 0, n-1
              CALL saxpy(m, a, x(i*ldx+1), 1, y(i*ldy+1), 1)
10        CONTINUE
      END IF
      RETURN
      END !SUBROUTINE SAXPYM


      SUBROUTINE daxpym (M, N, A, X, LDX, Y, LDY)
CDEC$ IF DEFINED (FTN_EXPORTS)
CDEC$     ATTRIBUTES DLLEXPORT::DAXPYM
CDEC$ ENDIF
      INTEGER m, n, ldx, ldy
      DOUBLE PRECISION a, x(ldx*n), y(ldy*n)
      INTEGER i

      IF (m .EQ. ldx .AND. m .EQ. ldy) THEN
          CALL daxpy(m * n, a, x, 1, y, 1)
      ELSE
          DO 20 i = 0, n-1
              CALL daxpy(m, a, x(i*ldx+1), 1, y(i*ldy+1), 1)
20        CONTINUE
      END IF
      RETURN
      END !SUBROUTINE DAXPYM


      SUBROUTINE caxpym (M, N, A, X, LDX, Y, LDY)
CDEC$ IF DEFINED (FTN_EXPORTS)
CDEC$     ATTRIBUTES DLLEXPORT::CAXPYM
CDEC$ ENDIF
      INTEGER m, n, ldx, ldy
      COMPLEX a, x(ldx*n), y(ldy*n)
      INTEGER i

      IF (m .EQ. ldx .AND. m .EQ. ldy) THEN
          CALL caxpy(m * n, a, x, 1, y, 1)
      ELSE
          DO 30 i = 0, n-1
              CALL caxpy(m, a, x(i*ldx+1), 1, y(i*ldy+1), 1)
30        CONTINUE
      END IF
      RETURN
      END !SUBROUTINE CAXPYM


      SUBROUTINE zaxpym (M, N, A, X, LDX, Y, LDY)
CDEC$ IF DEFINED (FTN_EXPORTS)
CDEC$     ATTRIBUTES DLLEXPORT::ZAXPYM
CDEC$ ENDIF
      INTEGER m, n, ldx, ldy
      DOUBLE COMPLEX a, x(ldx*n), y(ldy*n)
      INTEGER i

      IF (m .EQ. ldx .AND. m .EQ. ldy) THEN
          CALL zaxpy(m * n, a, x, 1, y, 1)
      ELSE
          DO 40 i = 0, n-1
              CALL zaxpy(m, a, x(i*ldx+1), 1, y(i*ldy+1), 1)
40        CONTINUE
      END IF
      RETURN
      END !SUBROUTINE ZAXPYM