Commande Fortran - Maximum d'occurence d'un nombre dans un vecteur

Version imprimable

11/12/2013, 11h46
nicolas823

Commande Fortran - Maximum d'occurence d'un nombre dans un vecteur

Bonjour à tous,

J'aurai besoin de votre aide sur l'existence d'une commande Fortran qui pourrait m'aider à simplifier mon problème.

Voici le problème:

Je dispose d'un vecteur contenant par exemple =[1,2,2,6,8,2]

Je voudrais repérer le nombre dont le nombre d'occurrence est maximal, ici, 2.

Existe-il une commande qui pourrait m'aider ? (même moyennant l'écriture d'une petite boucle)

Comment gérer le cas d'égalité ?

Cordialement,
Nicolas823
12/12/2013, 17h19
wafiflicha
Bonsoir,

essai de faire ça :
Code:

1 2 A=[1,2,2,6,8,2] MAXVAL(A,dim=1,mask=A=2)
:D

Est-ce que la solution suivante rencontre le critère de « commande Fortran » ?

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
   integer Valeur,NbrVal
   integer i,Nbr
   ...
   NbrVal = 0
   do i = 1, size(N)
      Nbr = count(N(i:) == N(i))
      if (Nbr > NbrVal) then
         NbrVal = Nbr
         Valeur = N(i)
      endif
   enddo
 
   write(*,*) Valeur, NbrVal

Salut,
je propose ma version :

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
30
31
32
33
34
35
36
37
program main
    implicit none
    integer, allocatable :: arr(:)
    integer :: maxi
 
    arr = (/ 1, 2, 2, 6, 8, 2 /)
    print '(99i4)', arr
    maxi = val_occ_max (arr)
    print '(i4)', maxi
 
contains
 
function val_occ_max (v1)
    integer, intent(in) :: v1(:)
    integer, allocatable :: v2(:)
    integer :: value, occur, val_occ_max, nb_occ_max
 
    allocate(v2(size(v1)))
    v2(:) = v1(:)
    val_occ_max = 0
    nb_occ_max = 0
    do
        if (size(v2) <= 0) exit
        value = maxval(v2)
        occur = count(v2 == value)
        print '(3i4)', value, occur, size(v2)
        if (occur > nb_occ_max) then
            val_occ_max = value
            nb_occ_max = occur
        endif
        v2 = pack(v1, v1 < value)
    enddo
    deallocate(v2)
 
end function val_occ_max
 
end program main

13/12/2013, 09h13
nicolas823

Merci à tous.

Effectivement, je me suis débrouillé avec maxval et suis arrivé à mes fins.