* * * *
19 Visitatori, 1 Utente
 

Autore Topic: Codice che va in loop infinito  (Letto 245 volte)

jonny78

  • Newbie
  • *
  • Post: 3
  • Karma: +0/-0
Codice che va in loop infinito
« il: 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

Avogadro

  • Full Member
  • ***
  • Post: 105
  • Karma: +0/-0
Re:Codice che va in loop infinito
« Risposta #1 il: 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
« Ultima modifica: Settembre 08, 2019, 10:26:35 pm da Avogadro »

jonny78

  • Newbie
  • *
  • Post: 3
  • Karma: +0/-0
Re:Codice che va in loop infinito
« Risposta #2 il: 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.

Avogadro

  • Full Member
  • ***
  • Post: 105
  • Karma: +0/-0
Re:Codice che va in loop infinito
« Risposta #3 il: 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 .
« Ultima modifica: Settembre 08, 2019, 10:27:21 pm da Avogadro »

Avogadro

  • Full Member
  • ***
  • Post: 105
  • Karma: +0/-0
Re:Codice che va in loop infinito
« Risposta #4 il: Settembre 09, 2019, 04:09:37 pm »
Se puo' servire

"LINEAR PROGRAMMING IN PASCAL "

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


SB

  • Scrittore
  • Full Member
  • *****
  • Post: 201
  • Karma: +0/-0
Re:Codice che va in loop infinito
« Risposta #5 il: 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





Avogadro

  • Full Member
  • ***
  • Post: 105
  • Karma: +0/-0
Re:Codice che va in loop infinito
« Risposta #6 il: 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."



 

Recenti

How To

Trascinare un file nel programma da xinyiman
DB concetti fondamentali e ZeosLib da xinyiman
Recuperare codice HTML da pagina web da xinyiman
Mandare mail con Lazarus da xinyiman
Stabilire il sistema operativo da xinyiman
Esempio lista in pascal da xinyiman
File INI da xinyiman
Codice di attivazione da xinyiman
Realizzare programmi multilingua da xinyiman
Lavorare con le directory da xinyiman
Utilizzare Unità esterne da Loryea
TTreeView da xinyiman
TTreeview e Menu da xinyiman
Generare controlli RUN-TIME da xinyiman
LazReport, PDF ed immagini da xinyiman
Intercettare tasti premuti da xinyiman
Ampliare Lazarus da xinyiman
Lazarus e la crittografia da xinyiman
System Tray con Lazarus da xinyiman
UIB: Unified Interbase da Microges2000
Il file: questo sconosciuto da Microges2000
Conferma di chiusura di un applicazione da xinyiman
Liste e puntatori da Microges2000
Overload di funzioni da Microges2000
Funzioni a parametri variabili da Microges2000
Proprietà da Microges2000
Conversione numerica da Microges2000
TImage su Form e Panel da Maverich
Indy gestiore server FTP lato Client da Maverich
PopUpMenu sotto Pulsante (TSpeedButton) da Maverich
Direttiva $macro da Microges2000
Toolbar da xinyiman
Evidenziare voci TreeView da Maverich
Visualizzare un file Html esterno da Maverich
StatusBar - aggirare l'errore variabile duplicata da Maverich
Da DataSource a Excel da xinyiman
Le permutazioni da xinyiman
Brute force da xinyiman
Indy 10 - Invio email con allegati da Maverich
La gestione degli errori in Lazarus da xinyiman
Pascal Script da xinyiman
Linux + Zeos + Firebird da xinyiman
Dataset virtuale da xinyiman
Overload di operatori da Microges2000
Lavorare con file in formato JSON con Lazarus da nomorelogic
Zeos ... dietro le quinte (prima parte) da Stilgar
Disporre le finestre in un blocco unico (come Delphi) da Maverich
Aspetto retrò (Cmd Line) da xinyiman
Come interfacciare periferica twain da Narciso
Ubuntu - aggiornare free pascal e lazarus da xinyiman
fpcup: installazioni parallele di lazarus e fpc da nomorelogic
Free Pascal e Lazarus sul Raspberry Pi da nomorelogic
Cifratura: breve guida all'uso dell'algoritmo BlowFish con lazarus e free pascal. da nomorelogic
Creare un server multithread da xinyiman
guida all'installazione di fpc trunk da subversion in linux gentoo da nomorelogic
Indice da nomorelogic
DB concetti fondamentali e connessioni standard da xinyiman
Advanced Record Syntax da nomorelogic
DB concetti fondamentali e DBGrid da xinyiman
DB concetti fondamentali e TDBEdit, TDBMemo e TDBText da xinyiman
Advanced Record Syntax: un esempio pratico da nomorelogic
Superclasse form base per programmi gestionali (e non) da nomorelogic
Superclasse form base per programmi gestionali (e non) #2 - log, exception call stack, application toolbox da nomorelogic
Superclasse form base per programmi gestionali (e non) #3 - traduzione delle form da nomorelogic
Superclasse form base per programmi gestionali (e non) #4 - wait animation da nomorelogic
Un dialog per la connessione al database:TfmSimpleDbConnectionDialog da nomorelogic
Installare lazarus su mac osx sierra da xinyiman
immagine docker per lavorare con lazarus e free pascal da nomorelogic
Utenti
  • Utenti in totale: 636
  • Latest: @diego
Stats
  • Post in totale: 12668
  • Topic in totale: 1578
  • Online Today: 22
  • Online Ever: 520
  • (Novembre 04, 2018, 12:00:31 pm)
Utenti Online
Users: 1
Guests: 19
Total: 20

Disclaimer:

Questo blog non rappresenta una testata giornalistica poiché viene aggiornato senza alcuna periodicità. Non può pertanto considerarsi un prodotto editoriale ai sensi della legge n. 62/2001.