Codice che va in loop infinito

[jonny78]

Salve a tutti, sono nuovo del forum è non sono un programmatore, ma un autodidatta; ho il seguente problema ho trovato un codice per l'algoritmo del simplesso (minimizzare/massimizzare una funzione in più variabili soggetta a dei vincoli), ho provato a compilare con il FreePascal con l'opzione "compatibilità con Turbo Pascal" ed inserendo sia l'esempio dell'autore che i valori numerici che m'interessano il programma va in loop ininito, ho provato a contattare l'autore ma non ho avuto risposta, visto che sul suo sito vi era anche il file com compilato con TP3, ho provato gli esempi sotto dosbox è funzionava alla grande, ho installato quindi TP7_1 ho compilato il codice modificando solo alcune istruzioni WRITELN "che portavano un istruzione CON, penso console" ed il programma funzionava, ho cercato di impostare gli stessi tipi (reali e interi) del programma originale in FreePascal ma il programma andava sempre in loop infinito, allora ho provato ad inserire i tipi Double e longint in TP7_1 ed a compilare ed il programma mi diceva di inserire l'opzione 8087 mentre standard era scelta numeric processing 'emulation', ebbene selezionando 8087/80287 il programma si comportava come nel FreePascal, cioè andava in loop infinito. Grazie in anticipo per qualsiasi suggerimento/aiuto
riporto sia la pagina web dell'autore che il listato con le modifiche da me effettuate.
http://www.calcna.ab.ca/~sfgale/simplex.html

[Avogadro]

nel codice vedo :

procedure optimize( trow : integer; var error : integer);
var xrow, xcol, iterate : integer;
begin
   repeat
      price( xcol, trow, error );
      if ( error = 0 ) then
         leave( xrow, xcol, error );
      if ( error = 0 ) then
         pivot( xrow, xcol );
   until ( error <> 0 )
end;

Non è che è qui che va in loop  perchè magari non trova mai la condizione di "stop" ? Banalmente, magari dopo 100 iterazioni lo si puo' far fermare => i codici delle procedure numeriche "autoconsitenti" fanno così, per evitare di girare  all' infinito quando non convergono verso una soluzione.

In ogni caso, hai seguito step by step il codice con il debugger ?

Infine, se puo servire, a pag 526 di "Numerical recipes in Pascal: the art of scientific computing" parlano dell'algoritmo del simplesso:

https://pdfs.semanticscholar.org/e795/26ba140a32c3bdcc59142818efc90327b87f.pdf

Io ho la copia cartacea , in inglese, comprata qualche anno fa, ma si dovrebbe trovare facilmente sul web  in pdf .


Sempre se puo' servire

https://wiki.freepascal.org/DMath

"The program includes routines for:

    Elementary functions
    Special functions
    Probability distributions
    Linear algebra
    Optimisation
...
"

con l' ultima versione di lazarus si installa in un attimo

Ciao

[jonny78]

Grazie mille per la risposta, il loop è proprio quello, cercherò di adattare il codice di numerical recipies la libreria DMATH già la conoscevo quando si chiamava TPMATH, essa comprende la procedura del simplesso solo che non supporta tutte le condizioni "<=>" di vincolo, la cosa che non capisco è perche sotto dos senza attivare quella opzione del compilatore 8087/80287 funziona e se l'attivo non funziona va in loop. Cmq ancora grazie per l'altro testo.

[Avogadro]

Non lo so perchè accade, non sono ferrato in linguaggio macchina/assempler,  al massimo mi sono cimentato con il lingiaggio "C" quando c'era da fare qualcosa fuori dall'ordinario del classico Fortran in tema di manipolazione di entità numeriche (sul VAX della Digital con potevi fare routine con diversi linguaggi e poi unirle con il linker in un unico eseguibile, altri tempi).

L'algoritmo del simplesso è interessante anche se non non ho mai avuto il tempo di approfondire; in tema di ottimizzazione quando è stato necessario mi sono limitato al DOE (es. https://www.qualityi.it/doe.html)  che per le cose che devo fare io si puo' tranquillamente implementare su un foglio elettronico .

Per vedere perchè il codice non gira e va in loop Lazarus ha un un buon debugger, puoi seguire il codice e vedere che valore assumono man mano le variabili .

[Avogadro]

Se puo' servire

"LINEAR PROGRAMMING IN PASCAL "

http://jean-pierre.moreau.pagesperso-orange.fr/p_linear.html

[SB]

Errore subdolo...
Allego il programma corretto e funzionante in Lazarus senza modificare impostazioni

Il problema si trovava nella funzione price nella penultima riga
   if ( quant = -eqzero ) then
Mai testare l'uguaglianza tra numeri perchè ci sono problemi di rappresentazione interna
Meglio
   if ( abs(quant - -eqzero) < 1e-10 ) then

Puoi cambiare 1e-10 con un numero più piccolo consentito dalla precisione

Non ho guardato se ci siano altre parti da modificare, ma è probabile

Ovviamente è valido anche il suggerimento di Avogadro sul controllo del numero di cicli

[Avogadro]

Giusta osservazione.

Se puo' servire in lazarus c'è già una funzione ad hoc, "iszero"

https://www.freepascal.org/docs-html/rtl/math/iszero.html

"IsZero checks whether the float value A is zero, up to a precision of Epsilon. It returns True if Abs(A) is less than Epsilon.

The default value for Epsilon depends on the type of the argument: it is 1E-4 for Single, 1E-12 for Double and 1E-16 for extended."

[Avogadro]

Se puo' servire, una pagina web sull'algoritmo del simplesso e tante altre cosette interessanti *

https://www.zweigmedia.com/RealWorld/tutorialsf4/frames4_3.html

*
https://www.zweigmedia.com/RealWorld/dietProblem/diet.html

[SB]

Sito interessante
Grazie della segnalazione

[doc]

Esistono anche le librerie numeriche per diversi problemi di calcolo numerico. Dateci una occhiata: https://wiki.freepascal.org/NumLib
Ciaooooo

[Avogadro]

In merito alla interessante libreria NumLib *  leggo la sua storia:

https://wiki.freepascal.org/NumLib

"Durante i primi anni novanta, ho fatto le mie lezioni di matematica numerica usando questa libreria. Circa un decennio dopo (2003?), a causa di un contatto esistente con il signor Tom Verhoeff, gli ho chiesto se la lib esisteva ancora. Ha contattato gli autori originali ed è stato in grado di fornire una copia del manuale e un backup della directory di lavoro dell'autore (l'autore era in pensione iirc), oltre al permesso verbale di rilasciare il pacchetto sotto la licenza LGPL di FPC.

Ho ripulito la maggior parte dei sorgenti nel 2003-2004 e nel 2005 ho ricevuto feedback, patch e aiuto con la documentazione da Darius Blaszijk (Jedi-Math). Comunque sia finito, l'aiuto di Darius è stato inestimabile, dal momento che ho ricevuto fonti critiche da lui, dopo che la mia casa è stata scassinata e il mio computer è stato rubato.

All'inizio del 2010 ho menzionato il pacchetto sul lazarusforum tedesco e ho deciso di dargli una spinta.

Sfortunatamente non avevo più la copia scansionata del manuale, quindi sto rifacendo quel lavoro poco a poco (e questa volta lo metto sul server FPC) "




Morale: c'è tanto codice interessante (e valido) in giro  ma con il passare del tempo sta andando perduto .

Ai primi tempi del web si trovava in giro codice d'ogni tipo, adesso si trovano solo scatole vuote, il web è diventato solo un manifesto pubblicitario infinito, vuoto d'ogni contenuto.
 
Che tristezza.




*
qui la documentazione

https://wiki.freepascal.org/NumLib_Documentation

[Avogadro]

In merito alla interessante libreria NumLib *  leggo la sua storia:

https://wiki.freepascal.org/NumLib

"Durante i primi anni novanta, ho fatto le mie lezioni di matematica numerica usando questa libreria. Circa un decennio dopo (2003?), a causa di un contatto esistente con il signor Tom Verhoeff, gli ho chiesto se la lib esisteva ancora. Ha contattato gli autori originali ed è stato in grado di fornire una copia del manuale e un backup della directory di lavoro dell'autore (l'autore era in pensione iirc), oltre al permesso verbale di rilasciare il pacchetto sotto la licenza LGPL di FPC.

Ho ripulito la maggior parte dei sorgenti nel 2003-2004 e nel 2005 ho ricevuto feedback, patch e aiuto con la documentazione da Darius Blaszijk (Jedi-Math). Comunque sia finito, l'aiuto di Darius è stato inestimabile, dal momento che ho ricevuto fonti critiche da lui, dopo che la mia casa è stata scassinata e il mio computer è stato rubato.

All'inizio del 2010 ho menzionato il pacchetto sul lazarusforum tedesco e ho deciso di dargli una spinta.

Sfortunatamente non avevo più la copia scansionata del manuale, quindi sto rifacendo quel lavoro poco a poco (e questa volta lo metto sul server FPC) "



Ancora: "dsl: Unità sconosciuta. Non esiste alcuna documentazione per esso, non è commentato e non può riconoscere direttamente l'algoritmo. I suffissi delle procedure sembrano indicare alcune caratteristiche del tipo a matrice (dall'unità SLE). Probabilmente qualche matrice pivot?"


Morale: c'è tanto codice interessante (e valido) in giro  ma con il passare del tempo sta andando perduto .

Ai primi tempi del web si trovava in giro codice d'ogni tipo, adesso si trovano solo scatole vuote, il web è diventato solo un manifesto pubblicitario infinito, vuoto d'ogni contenuto.
 
Che tristezza.




*
qui la documentazione

https://wiki.freepascal.org/NumLib_Documentation

[doc]

Purtroppo è così, molto materiale "storico" ma utile è andato perso.

[DragoRosso]

Non concordo completamente con quanto dite, a livello di conoscenza e di know how, e di relativo codice e documentazione, sono convinto che ci sia ancora molto che "cresce".

Certo "antiche" conoscenze, o meglio il codice, sono andate perse: ma altro arriva.

C'è tantissimo codice disponibile, su tanti fronti delle scienze umane. Arricchirlo spetta a noi tutti, intendo noi  programmatori, ricercatori, insegnanti, professionisti, etc ...

Ricordo che circa 35 anni fà mi venne chiesto di effettuare un progetto di ricerca per un sistema di scansione dei documenti e di riconoscimento ottico. Allora c'erano solo prodotti a pagamento (praticamente uno o due), ora un riconoscimento OCR lo si fà con un qualsiasi apparato informatico anche hobbistico e gratuitamente con sorgenti e quant'altro.

Tutte ciò che non c'è più a livello di codice si può ricostruire, anche tenendo conto delle nuove tecnologie: una volta per fare una operazione matematica ci si impiegava molto tempo, il codice era complesso .... ora esistono tecnologie come SSE, AVX e altro e il codice andrebbe comunque riscritto.

Ci sono librerie di tutti i generi, molto spesso a pagamento e super ottimizzate ed efficienti (magari senza codice sorgente), non direi che la conoscenza è andata persa.

Io sono positivo, perchè ricordo i tempi in cui occorreva entrare nelle BBS (via modem e chiamata telefonica) per trovare qualche cosa oppure avere la fortuna di conoscere qualche gruppo che condivideva le conoscenze.

No, direi che l'informatica si è decisamente evoluta e come nell'era moderna esistono le "fake news" e bisogna destreggiarsi per riconoscere le cose "buone", così ritengo che esistano come dite voi le scatole vuote  o anche di peggio e bisogna destreggiarsi per riconoscere quanto di buono esiste.

Ciao