Italian community of Lazarus and Free Pascal

Programmazione => Generale => Topic aperto da: DragoRosso - Maggio 23, 2024, 11:55:44 am

Titolo: Divisione per interi ... una giungla
Inserito da: DragoRosso - Maggio 23, 2024, 11:55:44 am
Mi sono imbattuto in una diatriba matematica, irrisolta e irrisolvibile probabilmente.
Per adesso ho imposto la mia idea, ossia il calcolo in PASCAL (potevo farlo essendo il capo progetto) e la questione operativa è chiusa ... ma non è per niente risolta la problematica pratica.

Eseguendo una operazione matematica in pascal, tramite interi e resti (div e mod), i risultati delle operazioni sono a mio avviso chiare e non ho mai avuto dubbi, sino a quando non mi sono interfacciato con dei moduli in python (e mi dicono che anche con moduli google si ottengo risultati identici al python) in cui i risultati differiscono dal Pascal (ma anche dal C)

Esempio:
Codice: [Seleziona]
//Questa è equivalente per entrambe
var i, j: INT64;
begin
  i := 9;
  j := 5;
  ShowMessage(IntToString(i div j)) // ==> 1
 //Python: (9 // 5) ==> 1
end;

//Qui iniziano i dolori !!!
var i, j: INT64;
begin
  i := -9;
  j := 5;
  ShowMessage(IntToString(i div j)) // ==> -1
 //Python: (-9 // 5) ==> -2                                                    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!
end;

Per il modulo vale ovviamente la stessa differenziazione (quindi i risultati sono differenti in base al fatto che i risultati siano negativi o positivi).
Ciò perchè le operazioni vengono effettuate sempre in float per alcuni linguaggi e i risultati vengono "arrotondati" per difetto o per eccesso (floor o ceil) a seconda delle convenzioni del linguaggio.

Avete informazioni su ciò? Ovvero avete notizia di come funzioni in altri linguaggi ?
Ulteriormente, secondo voi è più corretto una / altra o entrambe (ogni linguaggio si gestisce a modo proprio) ?

La mia opinione è che non ci dovrebbero essere dubbi su queste cose (e io non ce ne ho, sono convinto della mia idea magari antiquata), invece vedo che siamo ancora in alto mare ...

N.B.: lasciando perdere l'arrotondamento, il problema principale e che i risultati e le successive operazioni possono essere notevolmente dissimili tra due procedure effettuate con linguaggi differenti, ad esempio nel caso si debbano usare per vari motivi moduli diversi.
Titolo: Re:Divisione per interi ... una giungla
Inserito da: bonmario - Maggio 23, 2024, 01:11:13 pm
Ciao,
premesso che non avevo mai analizzato questa problematica prima, secondo me dipende da cosa rappresenta quel numero ed il motivo per cui lo stai dividendo.

P.S. Ho provato a fare -9/5 con la calcolatrice di Windows 11, e mi da come risultato -1.8. Quindi, in alcune condizioni, (-9 div 5) dovrebbe dare -1 come risultato !


Ciao, Mario
Titolo: Re:Divisione per interi ... una giungla
Inserito da: nomorelogic - Maggio 23, 2024, 01:27:54 pm
In effetti il risultato è emblematico.

L'operatore // in python (come si può vedere anche nel link sotto)
https://www.w3schools.com/python/trypython.asp?filename=demo_oper_floordiv (https://www.w3schools.com/python/trypython.asp?filename=demo_oper_floordiv)
in realtà arrotonda il risultato all'intero più vicino.

Bene, però qualcosa non torna.
9/5, come dice bonmario fa 1.8
però mentre è corretto a questo punto che -9/5 dia come risultato -2, è errato che 9/5 dia come risultato 1!  Dovrebbe dare 2.

Riguardo l'operatore div, secondo me il pascal risponde correttamente.
Ho voluto simulare quello che fa il python con Round (la Round del pascal arrotonda all'intero più vicino).
Infatti eseguendo il codice sotto, vedrete che i risultati sono esatti:


Codice: [Seleziona]
program Hello;
var i, j: INT64;
begin
  writeln ('a)', 9 div 5);
  writeln ('b)', -9 div 5);
  writeln ('c)', Round(9 / 5));
  writeln ('d)', Round(-9 / 5));
end.

Direi che è chiaro che bisogna sapere che operazione si vuole svolgere, ma il python secondo me in questo caso non da risultati corretti.


Edit:
in effetti credo di aver preso una stecca  ::)
l'operatore // del python "rounds the result down to the nearest whole number"
che in effetti significa che arrotonda all'interno minore

quindi l'intero minore di 1.8 è 1
mentre l'intero minore di -1.8 è -2

l'inghippo è che non si tratta solo di una divisione ma di una divisione con un arrotondamento all'intero minore
che comunque sono due operazioni che in effetti vengono svolte correttamente

NB:
la Round del pascal arrotonda all'interno più vicino
Titolo: Re:Divisione per interi ... una giungla
Inserito da: DragoRosso - Maggio 23, 2024, 03:41:51 pm
L'operazione che fà Python è corretta dal suo punto di vista, lo afferma la definizione dell'operatore "//" ma secondo me illogica anche se dal punto di vista prettamente matematico sembra che ci siano delle prevalenze verso il Python.

Dal punto di vista informatico invece sembra che ci sia preponderanza verso il calcolo Pascal, per cui di fatto una divisione intera ha come risultato una divisione troncata con segno (o come meglio si definisce arrotondata allo ZERO).

Python invece arrotonda con "floor" o "ceil"  in funzione al "segno" del risultato (appunto all'intero minore, dove il limite non è lo ZERO ma sono bensì il + / - INFINITO). Quindi il minore è il valore verso il - INFINITO e il maggiore è il valore verso il + INFINITO.
Tendenzialmente l'informatico il valore "minore" in termini di arrotondamento lo intende come valore verso lo ZERO.

Per intenderci, un risultato in Python di questa divisione (e i matematici sembrano prevalere per questo ?) sarebbe è:

(-9 // 5) = -2                 (9 // -5) = -2         (-9 // -5) = 1           (9 // 5) = 1
(-9 % 5) = 1                 (9 % -5) = -1        (-9 % -5) = -4         (9 % 5) = 4


in Pascal (e gli informatici sono più concordi su ciò):

(-9 div 5) = -1              (9 div -5) = -1         (-9 div -5) = 1         (9 div 5) = 1
(-9 mod 5) = -4            (9 mod -5) = 4       (-9 mod -5) = -4      (9 mod 5) = 4   

Entrambi dal punto di vista pratico danno un risultato equivalente,  ma potete capire che la rappresentazione non è proprio così equivalente. A me sembrano più coerenti i risultati Pascal, però sicuramente i matematici hanno le loro ragioni.

Comunque per rispondere a @BonMario, il campo di applicazione da dove ciò è comparso e l'applicazione di analisi spettrografica, dove la rappresentazione era legata ad alcuni algoritmi che a me (in Pascal) davano i risultati aspettati mentre in Python no.
E sinceramente mi ci è voluto un pò per riuscire a fare quadrato, perchè mentre in Pascal è tutto come ti aspetti (quando cambiano i segni degli operandi cambiano i segni dei risultati ma i valori "assoluti" sono sempre gli stessi, in Python non è così).

E' come un pò dire che (3 * 5) = (10+5) oppure  (3 * 5) = (20 - 5), entrambi sono corretti però dal punto di vista logico il primo mi appare un pò più ... come si può dire ... logico ...  ;D
Titolo: Re:Divisione per interi ... una giungla
Inserito da: nomorelogic - Maggio 23, 2024, 06:21:57 pm
Per intenderci, un risultato in Python di questa divisione (e i matematici sembrano prevalere per questo ?) sarebbe:

(-9 // 5) = -2                 (9 // -5) = -2         (-9 // -5) = 1           (9 // 5) = 1
(-9 % 5) = 1                 (9 % -5) = -1        (-9 % -5) = -4         (9 % 5) = 4


onestamente questa cosa mi destabilizza un po'...   :o
non mi pare tanto normale, sarà che sono un informatico e non un matematico ma, questi calcoli, farei fatica a farmeli calare

Titolo: Re:Divisione per interi ... una giungla
Inserito da: DragoRosso - Maggio 23, 2024, 06:50:25 pm
Per intenderci, un risultato in Python di questa divisione (e i matematici sembrano prevalere per questo ?) sarebbe:

(-9 // 5) = -2                 (9 // -5) = -2         (-9 // -5) = 1           (9 // 5) = 1
(-9 % 5) = 1                 (9 % -5) = -1        (-9 % -5) = -4         (9 % 5) = 4


onestamente questa cosa mi destabilizza un po'...   :o
non mi pare tanto normale, sarà che sono un informatico e non un matematico ma, questi calcoli, farei fatica a farmeli calare

A guarda, ti quoto e sottoscrivo in pieno, ma è quello che Python fà. E non solo, ma Google a ruota ... e prova la calcolatrice di Windows (ha solo il MOD) e vedrai gli stessi risultati.
Fortunatamente i fogli di calcolo lavorano ancora come in Pascal e come logica comanda.

Giusto per confermare la mia logica (anche se magari un pò demenziale  :o ) :
tralasciando lo ZERO se ho un qualcosa tipo -nx , sarà che il resto (ciò che resta  ::) ) sia comunque un -ax, dove a è minore dir n ovviamente, ma minore in senso numerico assoluto.
Se ho -9 mi avanzerà qualcosa da -1 a -9 (tralasciando che può rimanere anche ZERO). E infatti in Pascal il resto mantiene il segno del dividendo, mentre il risultato "troncato" ha la combinazione dei segni del dividendo e del divisore.

Mi sembra sicuramente più ovvio e gestibile così. Non avrei mai pensato a qualcosa di diverso. Oltretutto, che senso ha che (9 div 5) o (9 mod 5) abbiano valori assoluti differenti con segni differenti ?

Mah, rimango basito soprattutto perchè ho sempre più spesso a che fare con terze parti e nascono continue sorprese.

Mandi

P.S.: non è che (-6 // 5) in Python abbia un risultato diverso, fà sempre -2 ... però con resto di 4 ....
Titolo: Re:Divisione per interi ... una giungla
Inserito da: quack - Maggio 23, 2024, 07:16:19 pm
Credo che il "mod" di phyton faccia cosi:

-6//5 = -1.2 e lo arrotonda a -2

poi calcola l'intero più vicino come -2*5 = -10

trova il resto con -6 -(-10) = 4

Altro esempio

-9 // 5 = -1.8 arrotondato a -2

poi calcola l'intero più vicino come  -2*5 = -10

trova il resto con -9 -(-10) = 1

Per tornare al numero di partenza -2*5 + 1 = -9

Anche la calcolatrice di linux e openoffice danno gli stessi risultati.
Titolo: Re:Divisione per interi ... una giungla
Inserito da: nomorelogic - Maggio 23, 2024, 07:47:20 pm
ok, giusto per curiosità
come calcolatrice di linux cosa usi?

io uso dal un bel po' la mia fida speedcrunch
e il risultato è quello Pascal - style
Titolo: Re:Divisione per interi ... una giungla
Inserito da: quack - Maggio 23, 2024, 07:51:40 pm
Ho usato quella stock di gnome 46.
Titolo: Re:Divisione per interi ... una giungla
Inserito da: nomorelogic - Maggio 23, 2024, 08:10:16 pm
è ufficiale, ci sono 2 scuole di pensiero...

forse c'è il modulo dei boomer e quello della gen Z  :D
Titolo: Re:Divisione per interi ... una giungla
Inserito da: quack - Maggio 23, 2024, 08:40:07 pm
si mette male... anche su wolframalpha -9 mod 5 = 1

 :'(   
Titolo: Re:Divisione per interi ... una giungla
Inserito da: DragoRosso - Maggio 23, 2024, 10:57:51 pm
Bhè, non siamo soli.

Divisione intera con segno processore Intel in assembler (idiv): -9 div 5 = -1 quoziente, -4 resto  :D :D :D

E fino qui tutto bene: o il Pascal è l'unico che esegue calcoli in interi, insieme ad Intel oppure sfugge qualcosa.

P.S.: il foglio di calcolo in LibreOffice esegue il "resto" come Python.
Titolo: Re:Divisione per interi ... una giungla
Inserito da: nomorelogic - Maggio 24, 2024, 12:17:42 am
https://www.omnicalculator.com/it/matematica/modulo (https://www.omnicalculator.com/it/matematica/modulo)

Provate questo calcolatore di modulo.
Prima con 9 mod 5, poi con -9 mod 5.

Col nr. negativo verranno restituiti 2 risultati.
E una breve spiegazione...
Titolo: Re:Divisione per interi ... una giungla
Inserito da: quack - Maggio 24, 2024, 08:29:29 am
Se cercate "Divisione Euclidea" su Wikipedia c'è la spiegazione. Pascal fa la divisione con troncamento phyton quella con parte intera.

Ciao
Titolo: Re:Divisione per interi ... una giungla
Inserito da: DragoRosso - Maggio 24, 2024, 08:44:01 am
https://www.omnicalculator.com/it/matematica/modulo (https://www.omnicalculator.com/it/matematica/modulo)

Provate questo calcolatore di modulo.
Prima con 9 mod 5, poi con -9 mod 5.

Col nr. negativo verranno restituiti 2 risultati.
E una breve spiegazione...

Spiegazione corta (ma anche quella più lunga) efficiente e coerente tutto sommato, ma che comunque non si addice ai risultati di Python e company, dove il resto è anche negativo.
L'unica regola che posso comprendere è che il resto matematicamente sia compreso tra ZERO e il valore del divisore (minore o maggiore di quest'ultimo in base al segno  :o .....).

Comunque diatriba senza fine probabilmente, perchè una buona parte dei programmatori (se non tutti) non accetteranno mai che una divisione intera tipo (-9 / 5) o (-6 / 5) dia -2 come risultato, come i matematici non accetteranno mai l'inverso.

Rimango della mia idea, e se mi accorgessi del cambiamento repentino in favore di quest'ultima soluzione (matematica), bhè farei le mie routine in assembler o userei altri metodi per comunque avere i risultati alla "vecchia" maniera.

Ho verificato ora sui mie sorgenti: divisioni e mod su interi ne uso una sbarcata, e una simulazione banale su quello che succederebbe ad applicare il calcolo come fà Python ... bhè disastro totale ... almeno nel mio mondo dell'industria.

Nel panorama industriale la stragrande maggioranza dei calcoli (a parte il machine learing) coinvolge dati "interi", quindi usare la logica del float e degli arrotondamenti "fantasiosi" ( ;) ) non ha molto senso.

Sempre in un mio progetto ho conteggiato le dichiarazioni di variabile DOUBLE e INTEGER (senza contare i nativeint, i cardinal, le word, gli INTxx e gli UINTxx) quante volte comparivano: 87% integer e 13% double.

Ciao
Titolo: Re:Divisione per interi ... una giungla
Inserito da: DragoRosso - Maggio 24, 2024, 09:04:55 am
Se cercate "Divisione Euclidea" su Wikipedia c'è la spiegazione. Pascal fa la divisione con troncamento phyton quella con parte intera.
Ciao

Si, ma non torna comunque con i calcoli di Python, dove i valori di resto sono anche negativi (nel calcolo Euclideo il resto può essere solo positivo). Quindi Python non applica quel tipo di calcolo. Applica semplicemente un calcolo float dove l'arrotondamento e una scelta progettuale.

Parto dal presupposto che è ovvio e chiaro:
Citazione
"La divisione euclidea o divisione con resto è intuitivamente quell'operazione che si fa quando si suddivide un numero a di oggetti in gruppi di b oggetti ciascuno e quindi si conta quanti gruppi sono stati formati e quanti oggetti sono rimasti".

Quindi non ci possono essere più oggetti formati come risultato (quoziente) rispetto ai gruppo suddivisi: (-6 div 5) fa -1 cioè gli informatici ritengo che i gruppi siano presi in valore assoluto e il segno si applica poi. Inoltre il minore e maggiore è sempre rilevante verso le 0 per capire cosa si stà facendo. Dire che (-6 div 5) fà -2 perchè -2 è il numero di gruppi"minore" che si può ottenere (-1 è comunque maggiore di -2 e quindi la divisione da come risultato -2 ovviamente) mi pare ...  :-X (meglio che non mi esprima).

Cioè in poche parole per gli informatici (magari non tutti), che sono uomini di azione e vedono il concreto al di là delle teorie filosofiche, il risultati vanno presi in assoluto e poi si applicano i segni (questo detto in parole misere).

P.S.: comunque essere supportati dalle istruzioni di un processore INTEL ... bhè questa la dice tutta  :D :D :D :D :D :D.
Titolo: Re:Divisione per interi ... una giungla
Inserito da: nomorelogic - Maggio 24, 2024, 09:24:38 am
ho chiesto ad un mio collega matematico e mi ha liquidato dicendo che in matematica il modulo per come lo intendiamo noi informatici... non esiste

nel senso che non c'è un'operazione chiamata modulo ma ci sono dei procedimenti che possono calcolarlo, se abbiamo bisogno di quel calcolo

una risposta sicuramente discutibile ma che tutto sommato rende l'idea
il tutto sta nel tipo di calcolo che si necessita

probabilmente gli informatici hanno avuto quella necessità e lo hanno implementato così
nella versione python si è voluto essere aderenti ad una certa dimostrazione matematica

come si dice: non si finisce mai di imparare :)

Edit:
è comunque mia opinione che quelli del python avrebbero dovuto lasciare il calcolo del modulo
come già era in altri linguaggi e magari introdurne un secondo tipo
infatti quello che ci rimane in mano è 1 concetto con 2 metodi di calcolo...
Titolo: Re:Divisione per interi ... una giungla
Inserito da: DragoRosso - Maggio 24, 2024, 09:26:41 am
Voglio riportarvi un esempio pratico e banale:

un encoder, contaimpulsi che "conta" in positivo o negativo a seconda se ruota in senso orario o antiorario.
Non è che se questo fà -1000 impulsi e ogni giro viene effettuato con 750 impulsi fà -2 giri   :-X :-X :-X :

Pascal: (-1000 div 750) = -1 con resto xxx ==> deduco che ha fatto un (UNO) giro in senso antiorario ad esempio.
Python: (-1000 // 750) = -2 con resto xxx ==> deduco che ha fatto due (DUE) giri in senso antiorario..... (è barrato perchè non voglio che venga preso come esempio applicabile !!!)

In Pascal senza analizzare il resto ho già i miei dati, con Python deve farmi mille elucubrazioni sul quoziente e sul resto .... Questo è un esempio pratico e banale. Poi uno potrebbe aggiustarlo come vuole con mille elucubrazioni, algoritmi e fantasie varie: ma quello indicate è il modo più veloce, naturale e efficacie che non lascia nulla alla fantasia interpretativa.
Titolo: Re:Divisione per interi ... una giungla
Inserito da: DragoRosso - Maggio 24, 2024, 09:30:14 am
ho chiesto ad un mio collega matematico e mi ha liquidato dicendo che in matematica il modulo per come lo intendiamo noi informatici... non esiste

nel senso che non c'è un'operazione chiamata modulo ma ci sono dei procedimenti che possono calcolarlo, se abbiamo bisogno di quel calcolo

una risposta sicuramente discutibile ma che tutto sommato rende l'idea
il tutto sta nel tipo di calcolo che si necessita

probabilmente gli informatici hanno avuto quella necessità e lo hanno implementato così
nella versione python si è voluto essere aderenti ad una certa dimostrazione matematica

come si dice: non si finisce mai di imparare :)

Edit:
è comunque mia opinione che quelli del python avrebbero dovuto lasciare il calcolo del modulo
come già era in altri linguaggi e magari introdurne un secondo tipo
infatti quello che ci rimane in mano è 1 concetto con 2 metodi di calcolo...

Hi, hi, l'hanno detto anche a me, un fisico, che si è messo a ridere e ha concluso che noi informatici ci inventiamo molto spesso la matematica "per scansare la fatica del ragionamento"

P.S.: Comunque, il fatto che Intel (e penso anche AMD ovviamente) abbia implementato proprio in HARDWARE una istruzione aderente al ragionamento che si addice a noi informatici la dice lunga.
Titolo: Re:Divisione per interi ... una giungla
Inserito da: nomorelogic - Maggio 24, 2024, 11:49:21 am
P.S.: Comunque, il fatto che Intel (e penso anche AMD ovviamente) abbia implementato proprio in HARDWARE una istruzione aderente al ragionamento che si addice a noi informatici la dice lunga.

si, diciamo che è una scuola di pensiero radicata  :)
Titolo: Re:Divisione per interi ... una giungla
Inserito da: DragoRosso - Maggio 24, 2024, 03:34:25 pm
Per chi vuole, questo è il codice per fare con processori INTEL X64 in tutti i sistemi operativi teoricamente il calcolo tra numeri interi in modo "informatico" (ciò che fà il Pascal con il DIV e il MOD.

Poi farò anche i vari overload per i tipi e per x86.

Codice: [Seleziona]
{$IFDEF CPUX64}
/// <summary> Esegue una divisione per interi fornendo come risultato quoziente e resto </summary>
/// <remarks> Esegue il calcolo da "informatico" stile tutto chiaro !!! </remarks>
function DivandMod(Dividendo: int64; Divisore: int64; var Resto: int64): int64; register;
begin
  {$ASMMODE intel}
  asm
  //Dividendo, in RCX
  //Divisore, in RDX
  //Resto, memoria puntata da R8
  //Risultato (quoziente) in RAX
  mov RAX, RCX
  mov RCX, RDX
  CQO
  //idiv esegue la divisione intera e riporta il quoziente in RAX e in RDX il resto
  idiv RCX
  //SONO EQUIVALENTI LE DUE SUCCESSIVE IN DELPHI, IN LAZARUS NO !!!
  mov QWORD PTR [R8], RDX
  //mov QWORD PTR [Resto], RDX     //SOLO X DELPHI

  end;
end;
{$ENDIF}

procedure TForm1.Button1Click(Sender: TObject);
var quoziente: integer;
    resto: int64;
    Dividendo: int64;
    Divisore: int64;
begin
  Dividendo := 9;
  Divisore := -5;
  resto := 15;
  quoziente := DivandMod(Dividendo, Divisore, resto);
  ShowMessage('Quoziente= '+quoziente.ToString+' Resto= '+resto.ToString); 
end;[/code