Retrait de certaines cases dans un tableau Fortran90

Version imprimable

07/07/2016, 12h07
delphine2202

Retrait de certaines cases dans un tableau Fortran90
Bonjour à tous,

Je débute en Fortran et j'ai quelques difficultés à réaliser ce que je souhaite...
Alors, je dispose d'un tableau de réels du type :
0 1 2 3
1 0 4 5
2 4 0 6
3 5 6 0
Tableau dist

Je souhaiterais faire "remonter" mes valeurs pour remplir les cases où j'ai des 0, ou plus simplement, retirer les 0 de mon tableau pour obtenir un second tableau de réels du type :
1 1 2 3
2 4 4 5
3 5 6 6
Tableau temp

J'ai essayé plusieurs approches qui compilent mais qui ne me donnent pas ce que je souhaite (j'adapte ici mon code à l'exemple ci-dessus, en réalité les dimensions des deux tableaux sont variables mais toujours de la forme (n,n) pour dist et (n-1,n) pour temp).
Code:

1 2 3 4 real(kind = dp), parameter :: eps = 1.0E-08 integer(kind = sp) :: i,j real(kind = dp), dimension(4,4) :: dist real(kind = dp), dimension(3,4) :: temp
Essai 1 :
Code:

1 2 3 4 5 6 7 Do i = 1,4 Do j = 1,4 if(abs(dist(j,i)) > eps) then temp(j,i) = dist(j,i) Endif Enddo Enddo
Rq : le "abs" est ici superflu car mon tableau ne contient que des valeurs >= 0 mais je l'ai ajouté en prévision de possibles changements. Cela ne modifie en rien mon code normalement..

Essai 2 :
Code:

1 2 3 4 5 6 Do i = 1,4 Do j = 1,4 if(abs(dist(j,i)) < eps) cycle temp(j,i) = dist(j,i) Enddo Enddo
Essai 3 :
Code:

1 2 3 4 5 6 7 Do i = 1,4 Do j = 1,4 if(i /= j) then temp(j,i) = dist(j,i) Endif Enddo Enddo
Aucune de ces solutions ne fonctionne ... Il y a manifestement un problème au niveau de ma condition qui n'est pas prise en compte mais je ne sais pas pour quelle raison..
Je précise également que je souhaite effectuer cette manipulation pour ensuite déterminer la valeur minimale pour chaque colonne à l'aide de la fonction minval. Mon tableau dist est un tableau de distance entre points.
Les 0 correspondent en fait à la distance entre un point et lui-même. Afin de trouver le point le plus proche, il faudrait donc que j'enlève ces 0 pour obtenir la véritable valeur minimale..

Je vous remercie par avance pour vos réponses et si vous avez besoin de compléments d'information, n'hésitez pas ..!

Bonne journée à tous

Bonjour.
Le problème est mal posé. Faut-il supprimer les zéros ou la diagonale principale ?
S'il s'agit du premier cas, Fortran 90 possède des fonctions intrinsèques permettant des traiter les tableaux. Une solution avec pack, transpose et reshape :
Code:

1 2 3 4 5 6 implicit none integer, dimension(:,:), allocatable :: m, t m = reshape( (/ 0, 1, 2, 3, 1, 0, 4, 5, 2, 4, 0, 6, 3, 5, 6, 0 /), (/4, 4/) ) t = transpose(m) t = reshape( pack(t, t /= 0), (/3, 4/) ) print '(4(i3,x))', transpose(t)

Hello,

Pourquoi ne pas plus simplement changer les valeurs "0" en des valeurs arbitrairement grandes ?

Alors le minval te reverras les valeurs les plus petites autres que "0".

Tu peux sinon toujours faire le minval à la main en rajoutant une condition d'exclusion sur zero:
Code:

1 2 3 4 5 6 7 8 9 10 REAL(KIND=8) :: MIN_VAL MIN_VAL=1d9 !valeur arbitraire d'initialisation de la valeur minimale DO J=1,N !ou J=1,M si N/=M DO I=1,N IF(TAB(I,J)<MIN_VAL .and. TAB(I,J)/=0d0) MIN_VAL=TAB(I,J) END DO END DO
La même version gérant le cas où tu veux le min pour chacune des colonnes :
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 REAL(KIND=8),ALLOCATABLE,DIMENSION(:)::MIN_VAL ALLOCATE(MIN_VAL(N)) MIN_VAL=1d9 !valeur arbitraire d'initialisation de la valeur minimale DO J=1,N !ou J=1,M si N/=M DO I=1,N IF(TAB(I,J)<MIN_VAL(J) .and. TAB(I,J)/=0d0) MIN_VAL(J)=TAB(I,J) END DO END DO
En espérant avoir répondu à ta question,

Marlan

08/07/2016, 09h56
delphine2202

Bonjour à vous et merci pour vos réponses.

_dardanos_ : dans mon cas, enlever les zéros est équivalent à supprimer la diagonale principale car les 0 ne peuvent se trouver que sur cette diagonale. J'avais effectivement trouvé la fonction pack dans des docs mais je ne savais pas bien comment l'utiliser. Si je comprends bien d'après votre exemple, on reverse dans t les valeurs de t différentes de 0 et ce sur 3 lignes et 4 colonnes ..?
J'ai essayé de l'inclure dans mon code, cela compile mais cela m'indique "Segmentation fault" en sortie lorsque je l'exécute, auriez-vous une idée sur la raison de cet affichage en sortie ?

Marlan : J'avais pensé à la solution des valeurs arbitrairement grandes aussi et je pense que c'est celle-ci que je vais mettre en œuvre car plus facile ..! Cela fonctionne bien à première vue en tout cas ! Je garde en tête votre exemple du minval à la main qui peut m'être utile également :)

Merci en tout cas pour vos réponses, plusieurs solutions pour un même problème c'est toujours mieux :)

Bonne journée à tous

Citation:

Envoyé par delphine2202

J'ai essayé de l'inclure dans mon code, cela compile mais cela m'indique "Segmentation fault" en sortie lorsque je l'exécute, auriez-vous une idée sur la raison de cet affichage en sortie ?

Difficile à dire sans voir le code. Voici un exemple :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
program delete_diag
    implicit none
    integer, dimension(:,:), allocatable :: m, t
 
    m = reshape( (/ 0, 1, 2, 3, 1, 0, 4, 5, &
        2, 4, 0, 6, 3, 5, 6, 0 /), (/4, 4/) )
    call display(m)
    t = clean(m)
    call display(t)
 
contains
 
function clean(a)
    integer, dimension(:,:), allocatable :: a
    integer, dimension(:,:), allocatable :: clean
    clean = reshape( pack(a, a /= 0), (/ubound(a,1) - 1, ubound(a,2)/) )
    clean = transpose(clean)
end
 
subroutine display(a)
    integer, dimension(:,:), allocatable :: a
    integer :: i, j
    do j = lbound(a, 2), ubound(a, 2)
        print '(100(i2,x))', (a(i, j), i = lbound(a, 1), ubound(a, 1))
    enddo
    print *
end
 
end program

Code:

1
2
3
4
5
6
7
8
9
10
11
12
$ gfortran -Wall delete_diag.f90 
$  ./a.out 
 0  1  2  3
 1  0  4  5
 2  4  0  6
 3  5  6  0
 
 1  1  2  3
 2  4  4  5
 3  5  6  6
 
$

08/07/2016, 19h23
Sylvain Bergeron

@Dardanos: Ton code requiert un compilateur F2003 (clean = ...), ce qui n'est peut-être pas son cas.