Programmi

Vettore residuo

Subroutine per il calcolo del vettore residuo

C====================================================================
C Subroutine per il calcolo del vettore residuo
C Utilizza la subroutine matxvet (prodotto matrice vettore)
C====================================================================
      SUBROUTINE residuo (A, b, x, n, maxdim, r)
      IMPLICIT NONE
C
C Variabili INGRESSO:
C     A       matrice di dimensione n
C     b	      vettore termine noto di dimensione n
C     x       vettore iterata di dimensione n
C     n       dimensione di A, di b, di x e di r
C     maxdim  numero massimo di righe di A (come da dimensionamento)
C
C Variabili USCITA:
C     r       vettore residuo di dimensione n
C
C====================
C Parte dichiarativa
C====================
      INTEGER n, maxdim
      INTEGER i
      REAL*8 A(1:maxdim,*), b(1:*), x(1:*), r(1:*)
C====================
C Parte esecutiva
C====================
C
C Calcolo di A*x
C
      CALL matxvet (A, x, n, n, maxdim, r)
C
C Calcolo del vettore residuo
C
      do i = 1, n
        r(i) = b(i) - r(i)
      end do
      return
      end

• » Invia tramite EMAIL
• » Versione per la STAMPA
• » Le vostre opinioni