Comment initialiser les tableaux à deux dimensions en Fortran

En C, vous pouvez facilement initialiser un tableau en utilisant les accolades de la syntaxe, si je me souviens bien:

int* a = new int[] { 1, 2, 3, 4 };

Comment pouvez-vous faire de même en Fortran pour les tableaux à deux dimensions quand vous voulez initialiser une matrice spécifique avec des valeurs de test de mathématiques fins? (Sans avoir à doublement de l'indice de chaque élément sur des relevés distincts)

Le tableau est défini par

real, dimension(3, 3) :: a

real, dimension(:), allocatable :: a

InformationsquelleAutor Fludlu McBorry | 2010-09-14

51

Vous pouvez le faire en utilisant remodeler et forme intrinsèques. Quelque chose comme:
```
INTEGER, DIMENSION(3, 3) :: array
array = reshape((/1, 2, 3, 4, 5, 6, 7, 8, 9 /), shape(array))
```
Mais n'oubliez pas la colonne ordre majeur. Le tableau sera
```
1   4   7
2   5   8
3   6   9
```
après remodelage.

De sorte à obtenir:
```
1   2   3
4   5   6
7   8   9
```
vous avez également besoin d' transposer intrinsèque:
```
array = transpose(reshape((/1, 2, 3, 4, 5, 6, 7, 8, 9 /), shape(array)))
```
Pour plus d'exemple général (affectables tableau 2D avec des dimensions différentes), on a besoin taille intrinsèque:
```
PROGRAM main

  IMPLICIT NONE

  INTEGER, DIMENSION(:, :), ALLOCATABLE :: array

  ALLOCATE (array(2, 3))

  array = transpose(reshape((/1, 2, 3, 4, 5, 6 /),                            &
    (/size(array, 2), size(array, 1) /)))

  DEALLOCATE (array)

END PROGRAM main
```
- 1) la Plupart des compilateurs acceptent désormais le Fortran 2003 notation [] pour initialiser les tableaux, au lieu de la un peu maladroit (/ /). 2) Pour les cas simples, vous pouvez omettre les transposer en fournissant les valeurs dans la colonne ordre majeur: array = remodeler ( [1, 4, 7, 2, 5, 8, 3, 6, 9 ], forme (array) )
- J'ai oublié de mentionner que nous sommes tenus de travailler en Fortran 90.
InformationsquelleAutor Wildcat
19

Pour multidimensionnelle (rang>1) des tableaux, le Fortran façon pour l'initialisation diffère de la C de la solution parce que dans C les tableaux multidimensionnels sont juste des tableaux de tableaux de, etc.

En Fortran, à chaque rang correspond à un attribut différent de la modification du type de données. Mais il y a un seul constructeur array, pour le rang-1 tableaux. À partir de ces deux raisons, l'initialisation par le biais de la matrice de constructeur requiert la RESHAPE propre fonction.

En plus de ce qui a déjà été répondu, il y a une façon plus directe de la saisie de la valeur d'une matrice par ligne plutôt que par colonne: remodeler a un argument optionnel ORDER qui peut être utilisé pour modifier l'ordre de remplissage de l'élément du tableau multidimensionnel avec les entrées de la matrice de constructeur.

Par exemple, dans le cas de l'exemple dans la première réponse, on pourrait écrire:
```
INTEGER, DIMENSION(3, 3) :: array=reshape( (/1, 2, 3, &
                                              4, 5, 6, &
                                              7, 8, 9 /), &
                                           shape(array), order=(/2,1/) )
```
obtenir le remplissage de la matrice exactement dans l'ordre indiqué par les lignes de code.

Le tableau (/2, 1/) les forces de l'index de colonne (2) avoir préséance sur l'index de ligne (1), donnant l'effet désiré.

InformationsquelleAutor GiorgioP
7

Initialisation de tableau peut être fait dans le tableau de déclaration elle-même, comme indiqué ci-dessous:
```
program test
 real:: x(3) = (/1,2,3/)
 real:: y(3,3) = reshape((/1,2,3,4,5,6,7,8,9/), (/3,3/))
 integer:: i(3,2,2) = reshape((/1,2,3,4,5,6,7,8,9,10,11,12/), (/3,2,2/))

end program test
```
Ça m'étonne que
```
 real:: y(3,3) = (/(/1,2,3/),(/4,5,6/),(/7,8,9/)/)
```
est pas accepté par le compilateur (essayé g95, gfortran). Il s'avère que le shape de
(/(/1,2,3/),(/4,5,6/),(/7,8,9/)/) est 9 et pas 3 3!

InformationsquelleAutor Sourav Mukherjee

Vous devez vous connecter pour publier un commentaire.