Questo sito contribuisce alla audience di

Operazione saxpy nel Fortran

Come realizzare l'operazione saxpy in un programma fortran

fortran 90,intel fortran,visual fortran,compilatore fortran

Un’operazione vettoriale di base è l’operazione saxpy, si tratta in realtà di una combinazione tra una moltiplicazione vettoriale e un’addizione tra vettori. L’espressione che ci permette di eseguire una tale operazione è la seguente:

y: = alfa * x + y

dove alfa è uno scalare, ma x e y sono vettori.

Ecco una subroutine semplice per questa operazione:

subroutine saxpy (n, alpha, x, y)
      integer n
      real alpha, x(*), y(*)
c
c Saxpy: Compute y := alpha*x + y,
c where x and y are vectors of length n (at least).
c
c Local variables
      integer i
c
      do 10 i = 1, n
         y(i) = alpha*x(i) + y(i)
   10 continue
c
      return
      end