code: def CAXPY(N, CA, CX, INCX, CY, INCY):
      #  Purpose
      #  =======
      #
      #     CAXPY constant times a vector plus a vector.
      #
      #  Further Details
      #  ===============
      #
      #     jack dongarra, linpack, 3/11/78.
      #     modified 12/3/93, array(1) declarations changed to array(#)
      #
      #  =====================================================================
      #
      if N <= 0:
          return
      if abs(CA) == 0.0e0:
          return
      if INCX == 1 and INCY == 1:
          "omp parallel for"
          for I in range(N):
              CY[I] = CY[I] + CA * CX[I]
      #
      #        code for both increments equal to 1
      #
      else:
          #
          #        code for unequal increments or equal increments
          #          not equal to 1
          #
          IX = 0
          IY = 0
          if INCX < 0:
              IX = (-N + 1) * INCX
          if INCY < 0:
              IY = (-N + 1) * INCY
          for I in range(N):
              CY[IY] = CY[IY] + CA * CX[IX]
              IX = IX + INCX
              IY = IY + INCY
      return CY
  
label_desc: Use image node_id 1 for calling a global function with example usage: CAXPY(N, CA, CX, INCX, CY, INCY) and returns: CY