resueltos - producto vectorial fortran

indexar múltiples elementos de matriz no adyacentes fortran

Question

¿Hay alguna manera en fortran de acceder a muchos elementos de una matriz fuera de un bucle? Por ejemplo, una matriz dada de 100 elementos

real(100) :: a

¿Puedo hacer algo como esto para acceder a los elementos 1,4,7,54,81 que no siguen un paso regular?

a(1,4,7,54,81)= 3.21423

Answer 1

RESHAPE y DONDE valen la pena mirar. Si está determinando qué elementos ''extraer'', entonces ASIGNAR una nueva variable y rellenar los elementos de A en B. Tal vez algo como esto:

REAL,    DIMENSION(100)            :: A
LOGICAL, DIMENSION(100)            :: A_Mask
INTEGER                            :: SizeB
REAL,    DIMENSION(:), ALLOCATABLE :: B
!...
A_Mask = .FALSE.
WHERE(A > 1.0)
  A_Mask = .TRUE.
ENDWHERE
SizeB = SUM(A_Mask)
!... then allocate B and fill it.

Answer 2

podrías usar un subíndice vectorial: a( (/1,4,7,54,81/) )= 3.21423