Italian community of Lazarus and Free Pascal

Programmazione => Generale => Topic aperto da: jonny78 - Settembre 05, 2019, 07:16:17 pm

Titolo: Codice che va in loop infinito
Inserito da: jonny78 - Settembre 05, 2019, 07:16:17 pm
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 (http://www.calcna.ab.ca/~sfgale/simplex.html)
Titolo: Re:Codice che va in loop infinito
Inserito da: Avogadro - Settembre 05, 2019, 11:23:17 pm
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
Titolo: Re:Codice che va in loop infinito
Inserito da: jonny78 - Settembre 06, 2019, 12:18:45 am
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.
Titolo: Re:Codice che va in loop infinito
Inserito da: Avogadro - Settembre 08, 2019, 10:20:00 pm
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 .
Titolo: Re:Codice che va in loop infinito
Inserito da: Avogadro - Settembre 09, 2019, 04:09:37 pm
Se puo' servire

"LINEAR PROGRAMMING IN PASCAL "

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

Titolo: Re:Codice che va in loop infinito
Inserito da: SB - Settembre 15, 2019, 09:34:08 am
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




Titolo: Re:Codice che va in loop infinito
Inserito da: Avogadro - Settembre 22, 2019, 01:16:23 am
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."