Crivello di Eratostene

[Top Fuel]

A 3 miliardi ti ha dato out of bound perchè il limite di elementi è 2^32-1.

[Gianluigi]

A 3 miliardi ti ha dato out of bound perchè il limite di elementi è 2^32-1.

Il limite di elementi di cosa?

[vuott]

Se non erro il codice di Vuott non fa uso di vettori e questo potrebbe spiegare il fatto che è veloce il doppio.

...c'è un Puntatore ad un Intero:

Codice: [Seleziona]

int *list

Non vi sono Moltiplicazioni.

[Gianluigi]

Nel mio codice quello che rallenta tanto è il secondo loop, la prima parte individua i numeri non primi mettendoli a True in 7/8 secondi, tutto il resto se lo posta via questo:

Codice: [Seleziona]

For i = 2 To n
    If v[i] = False Then
      Print i
    Endif
  Next

Purtroppo io non conosco il C, tu non potresti tradurre in gambese il codice scritto in C per vedere se funziona meglio 

 

PS Chissà Top Fuel cosa voleva dire...

[vuott]

...potresti tradurre in gambese il codice scritto in C ....

...gambese ?   
http://www.jokisaari.net/gambia/pics/gambia-map.gif

Ad ogni modo ecco la traduzione:

Codice: [Seleziona]

Public Sub Main()

  Dim limite, dime, lim, i, j, z As Integer
  Dim lista As Integer[]
  
    limite = 1000000
    
    dime = limite / 2 + 1
    
    lista = New Integer[dime]
    
    Print 2
    
    lim = Sqr(CFloat(limite))

    j = 1
    i = 3
      
    While i <= lim
      If lista[j] = 0 Then
        z = j + i
        While z < dime
          lista[z] = -1
          z += i
        Wend
        Print i
      Endif
      i += 2
      Inc j
    Wend

    While i < limite - 1
      If Not lista[j] Then Print i
      i += 2
      Inc j
    Wend

End

[Gianluigi]

...gambese ?   
http://www.jokisaari.net/gambia/pics/gambia-map.gif

Ho notato che anche la console rallenta, infatti nel codice C si ottiene il result a seguire, pertanto se noi inseriamo dietro a Print nei codici in gambase  la virgola (il doppio punto e virgola è più lento) riduciamo il gap col C, questi i risultati su 100 milioni nel mio computer:

Codice: [Seleziona]

Crivello in C vuott   22371,8590736389 msec / 22567,2857761383 msec
Setaccio Vuott        24636,3999843597 msec / 24362,3239994049 msec
Setaccio Top Fuel     26364,4378185272 msec / 25966,2098884583 msec
Setaccio Gianluigi    30971,825838089  msec / 29990,3969764709 msec

Il nuovo codice su 2 miliardi mi da memoria terminata.
 

[vuott]

Nei programmi Gambas puro hai usato Fast ?

[Gianluigi]

Nei programmi Gambas puro hai usato Fast ?

Si in tutti escluso il tuo con la libreria C, senza Fast il tuo codice impiega 30520,5111503601 msec invece di 24636,3999843597 msec

[Top Fuel]

A 3 miliardi ti ha dato out of bound perchè il limite di elementi è 2^32-1.

Il limite di elementi di cosa?

Il limite di elementi che può avere un array. Dall'errore si deduce che Gambas usa un intero di 4 byte per conteggiare il numero di elementi.

[Gianluigi]

A 3 miliardi ti ha dato out of bound perchè il limite di elementi è 2^32-1.

Il limite di elementi di cosa?

Il limite di elementi che può avere un array. Dall'errore si deduce che Gambas usa un intero di 4 byte per conteggiare il numero di elementi.

Si è vero, hai ragione.
Grazie 

[vuott]

http://www.lescienze.it/news/2016/09/30/news/nuovo_metodo_calcolare_numeri_computer-3251946/

[Top Fuel]

Interessante, ma come fai in Gambas a leggere la cache della CPU? E sopratutto come fai a sapere in quale cache sono scritti i numeri? L1, L2, L3, core 0, core 1...

[vuott]

...come fai in Gambas a leggere la cache della CPU?

A mio modesto parere il problema non è se con Gambas sia possibile leggere la cache della CPU, ma più in generale "come", e soprattutto "con quali risorse di quale libreria" sia possibile.



eloaders aveva mostrato una libreria esterna (ce n'è comunque anche qualche altra) per ottenere semplicemente alcune informazioni generali sulla CPU:
http://www.gambas-it.org/smf/index.php?topic=5197.0
Non so se questa o altre consentano di leggere dentro la cache.