laboratorio di calcolo II AA 2003/04

1. UNIVERSITA’ DEGLI STUDI ROMA TRE DIPARTIMENTO DI FISICA “E. AMALDI” laboratorio di calcolo IIAA 2003/04 a cura di Domizia Orestano Dipartimento di Fisica Stanza 159 - tel. (06 5517) 7281 www.fis.uniroma3.it/~orestano orestano@fis.uniroma3.it

2. Introduzione: obiettivi del corso Introduzione alla Programmazione ad Oggetti Accento posto sulle “Idee” e non sugli aspetti tecnici “Architetti” e non “Esperti del Linguaggio C/C++” Imparare a progettare un approccio al Calcolo Restano fuori Elementi di Analisi Numerica Aspetti di Grafica o di Tecnologie Multimediali

3. Contenuti Introduzione al linguaggio C++ e alla Programmazione ad Oggetti Elementi di linguaggio UML Elementi di strutture riutilizzabili Metodo Comprensione del problema, ricerca della soluzione in ambito OO, implementazione in C++ Lezioni Frontali (circa 40%) e Laboratorio (circa 60%) Utilizzo di un problema guida (Sistema Solare) Sintassi in Appendice a ciascuna lezione Introduzione: organizzazione del corso

4. Introduzione: orario del corso Indicativamente... Abbiamo 10 settimane • 10 o 11 lezioni frontali da 2 ore • 0.5 - Introduzione • 8.5 - Linguaggio C++ • 1.5 - UML e strutture riutilizzabili • 0.5 - Conclusione Si svolgeranno il lunedì dalle 14 alle 16 • 9 o 10 esercitazioni di laboratorio da 4 ore Si svolgeranno presso i laboratorio didattici il martedì (studenti del II anno) e il giovedì (studenti del I anno) dalle 9 alle 13

5. Introduzione: strumenti di lavoro • Dispense • Esercitazioni di Laboratorio (in gruppi di 2 persone) • Manuale C++ • Come riferimento per la sintassi • Attenzione all’impostazione (partiamo direttamente dal C++, e non dal C!) • Alcuni suggerimenti nella prossima trasparenza • Presenza attiva in Laboratorio e in Aula

6. Introduzione: manuali (1)

7. Introduzione: manuali (2)

8. Introduzione: svolgimento degli esami • Prova Pratica • Prova in Laboratorio, in linguaggio C++ • Prova Orale • Discussione su un problema scelto dallo studente e preparato personalmente durante il corso • Ogni studente sceglie un problema di Fisica (semplice!) concordato o suggerito dal docente • Lo formalizza nella logica della programmazione ad oggetti • Lo scrive in C++ e lo discute nel corso della Prova Orale

9. Problematiche del calcolo moderno Applicazione del calcolo a sistemi sempre più complessi per quantità di dati da manipolare, numero di CPU utilizzate, distribuzione geografica dei soggetti coinvolti Esempi: esperimenti di fisica delle alte energie e di astroparticelle, biologia molecolare e progetto genoma, ma anche telecomunicazioni, sistemi bancari, reti di trasporto.

10. Requisiti per il software moderno • Robustezzaprotezioni nell’accesso ai dati • Possibilità di ri-utilizzo del codice • economia di risorse umane ed economiche • maggiore affidabilità • Portabilità • verso sistemi operativi diversi • verso diverse versioni di uno stesso sistema • Flessibilità e Organizzazione del Codice semplicità di gestione e sviluppo successivo del codice

11. Procedurale: si sceglie l’algoritmo che risolva il problema dato e il linguaggio permette di implementare l'algoritmo scelto Modulare: si decide quali moduli utilizzare e si separano i dati in modo danasconderli entro ciascun modulo Definita dall’utente: si decide quali tipi utilizzare e si genera un insieme di operazioni per i vari tipi Orientata agli Oggetti: si decide quali classi utilizzare, si genera un insieme di operazioni per le varie classi e si raggruppano le proprietà comuni tramite il meccanismodell'ereditarietà Generica: si identificano una serie di algoritmi e si parametrizzano in modo che possanoessereutilizzatiperpiùclassiotipidefiniti dall'utente. Paradigmi di programmazione

12. Paradigmi e linguaggi di programmazione I diversi linguaggi di programmazione possono essere più o meno indicati ad esprimere un tipo di programmazione (paradigma). Il C++ permette di implementare unaprogrammazione definita dall'utente, consentendo la definizione di tipi di oggetti e delle operazioni effettuabili con essi, con possibilità di mettere in evidenza i collegamenti e le similitudini tra di essi mediante il meccanismo dell’ereditarietà (ovvero una programmazione Orientata agli Oggetti) ma consente anche di realizzare una programmazione Generica, creando codice utilizzabile in situazioni formalmente simili anche in assenza di un collegamento concettuale.

13. Caratteristiche della programmazione OO • Incapsulamento:i dati dell’oggetto sono nascosti ad altri oggetti ed è possibile accedervi solo attraverso modalita’ ben definite Da questa caratteristica derivano robustezza e flessibilità del codice • Ereditarietà:gli oggetti complessi possono essere costruiti a partire da oggetti più semplici. Gli oggetti complessi derivano tutto o parte del loro comportamento dagli oggetti a partire daiquali sono stati generati Ne derivano ri-utilizzabilità e possibilità di organizzazione del codice • Polimorfismo:oggetti simili possono essere trattati, in alcuni casi, come se fossero dello stesso tipo, senza la necessità di implementare trattamenti specifici per distinguere tra le varie tipologie Rende il codice più flessibile e aumenta l’efficienza dell’ereditarietà La Portabilità invece è una caratteristica del linguaggio ed in particolare del C++ standard (ANSi-C++)

14. Progettazione di un programma • Scelta degli Oggetti • Identificazione delle relazioni tra gli Oggetti (proprietà statiche o “geometriche”) • Esame dei casi d’uso (use case) • Identificazione delle interazioni tra gli Oggetti (proprietà dinamiche) • Eventuale ridefinizione delle responsabilita’ delle classi NB gli Oggetti vengono descritti mediante classi, ovvero tipi definiti dall’utente che realizzano l’incapsulamento dei dati

15. La simulazione del Sistema Solare Il problema fisico: Determinare l’evoluzione temporale di un certo numero di corpi sottoposti a reciproca attrazione gravitazionale F12= m a F12= - G m1 m2r / r3 G= 6.673e-11Unità mks

16. Scelta degli Oggetti Quali “Oggetti” entrano nel problema? corpi celesti, sonde, pianeti, stelle, satelliti … Quali proprieta’ hanno tali oggetti? massa, posizione, velocità…

17. Relazioni tra gli Oggetti Quali sono le relazioni tra gli Oggetti? Possono essere raggruppati? Quali proprieta’ sono condivise da Oggetti simili? Quali proprieta’ permettono di distinguere Oggetti simili tra loro? un pianeta è un corpo celeste, una sonda è un satellite… … ed è anche un corpo celeste tutti gli oggetti possono essere ricondotti a Corpi Celesti o a Sonde, dove la Sonda è un Corpo Celeste con qualche proprietà in più Quali oggetti sono superflui? Quali oggetti mancano? manca il Sistema Solare … e un Sistema Solare è composto da dei Corpi Celesti!

18. UML Universal Modeling Language

19. Casi d’uso A cosa serve esattamente il programma? simulazione del moto dei vari Corpi Celesti inseriti nel Sistema Solare In quali situazioni concrete tali oggetti dovranno essere utilizzati? • Si vuole conoscere lo stato del sistema (posizione di ciascun Corpo Celeste) dopo un tempo T dall’istante iniziale • Si suddivide l’intervallo T in n intervallini dT nei quali si possano trascurare variazioni di forza e di velocità • Per ogni intervallo dT • Nota la posizione di ciascun Corpo Celeste all’inizio dell’intervallo • La si usa per determinare la forza agente su ciascun Corpo Celeste • Nota la forza si determina l’accelerazione di ogni Corpo Celeste • Si calcola la velocità del Corpo Celeste, nota la velocità all’inizio dell’intervallo, assumendo che l’accelerazione resti costante • Si calcola la nuova posizione, assumendo che la velocità resti costante

20. Scenari Oltre ai casi d’uso si possono introdurre gli scenari. Questi descrivono in maniera piu’ dettagliata situazioni particolari che possono presentarsi e devono essere previste dal programma. Nel nostro caso possiamo ad esempio prevedere di voler lanciare una sonda dalla superficie di un altro Corpo Celeste X ad un istante t compreso tra 0 e T. Il programma dovrà allora anche controllare quando si raggiunga il tempo t e aggiungere al Sistema Solare un Oggetto Sonda con una certa velocità iniziale rispetto al Corpo Celeste X, assegnandole la posizione di X ed una velocità ottenuta dalla combinazione della velocità iniziale relativa della Sonda rispetto a X e della velocità del Corpo Celeste X.

21. Interazioni tra gli Oggetti Non è sufficiente capire la relazione “geometrica” tra gli oggetti, è necessario anche valutare come gli oggetti interagiscono “dinamicamente” tra loro Quali compiti ha ciascuna classe? Chi determina le nuove coordinate del sistema? Chi invoca i metodi per far evolvere il sistema? Una possibile scelta: • CorpoCeleste ha il metodo CalcolaPosizione • SistemaSolare ha il metodo Evolvi che determina il nuovo stato del sistema dopo un tempo T • Il metodo Evolvi (di SistemaSolare) si occupa di invocare i metodi CalcolaPosizione di ciascun pianeta del sistema solare.

22. Responsabilità delle classi è chiaro di chi sia la responsabilita’ di un punto chiave del programma? Una classe ha troppe responsabilita’? Una classe è irrilevante e può essere “assorbita” da un’altra? • Distribuire le responsabilità tra le varie classi • Sfruttare un procedimento a cascata • Limitare le responsabilità di una classe ad una o due

23. Interaction Diagram (semplificato)

24. Sessione CRC CRC: Classi, Responsabilità, Collaborazione capire come funziona dinamicamente il programma che si vuole realizzare • Giuoco di Ruolo: • ciascuno dei partecipanti è un Oggetto del programma • uno dei partecipanti è il Programma "Main • Ogni partecipante deve conoscere le caratteristiche di ciascuna Classe: • "come è fatta" • "quali azionì può compiere” • Durante la sessione tutto si svolge come nel corso dell'esecuzione del Programma, • in questo modo emergono casi non previsti ed eventuali problemi.

25. Nome Classe Attributi Metodi Costruzione di una classe • Attributi • public • private • protected • Metodi • public • private • protected • Obbligatorio il tipo • (“nessun tipo” void ) • Accesso agli attributi di una classe • (tipo Set e tipo Get) Set Get

26. La dichiarazione di una classe: l’header file (.h) CorpoCeleste Nome (stringa) m (num. reale) x (num. reale) y (num. reale) vx (num. reale) vy (num. reale) CalcolaPosizione(forza, dt) StampaVelocita’() StampaPosizione() M() X() Y() Vx() Vy() Contiene la dichiarazione degli Attributi e dei Metodi Nomeclasse.h

27. L’implementazione di una Classe: l’implementation file ( .cc) Contiene l’implementazione degli Attributi e dei Metodi Nomeclasse.cc

28. #ifndef CORPOCELESTE_H #define CORPOCELESTE_H class CorpoCeleste { private: char *Nome; double m; double x; double y; double vx; double vy; public: CorpoCeleste() ; CorpoCeleste (const char *nomeCorpo, float mass,float xpos, float ypos, float vxi,float vyi); ~CorpoCeleste() ; void calcolaPosizione(float fx, float fy, float t); void stampaPosizione(); void stampaVelocita(); const char *nome() ; double M() ; double X() ; double Y() ; double Vx() ; double Vy() ; }; #endif CorpoCeleste.h

29. #ifndef CORPOCELESTE_H #define CORPOCELESTE_H class CorpoCeleste { protected: char *Nome; double m; double x; double y; double vx; double vy; public: CorpoCeleste() ; CorpoCeleste (const char *nomeCorpo, float mass,float xpos, float ypos, float vxi,float vyi); ~CorpoCeleste() ; virtual void calcolaPosizione(float fx, float fy, float t); void stampaPosizione(); void stampaVelocita(); const char *nome() ; double M() ; double X() ; double Y() ; double Vx() ; double Vy() ; }; #endif CorpoCeleste.h con modifiche per ereditarietà

30. #include "CorpoCeleste.h“ CorpoCeleste::CorpoCeleste( ) { } CorpoCeleste::CorpoCeleste (const char *nomeCorpo, float mass, float xpos, float ypos, float vxi, float vyi) { } CorpoCeleste::~CorpoCeleste() { } void CorpoCeleste::calcolaPosizione(float fx, float fy, float t) { } void CorpoCeleste::stampaPosizione() { } void CorpoCeleste::stampaVelocita() { } const char* CorpoCeleste::nome() { } double CorpoCeleste::M() { } double CorpoCeleste::X() { } double CorpoCeleste::Y(){ } double CorpoCeleste::Vx(){ } double CorpoCeleste::Vy(){ } CorpoCeleste.cc Schema vuoto

31. CorpoCeleste.cc Prima parte #include "CorpoCeleste.h" #include <string.h> #include <iostream.h> #include <iomanip.h> CorpoCeleste::CorpoCeleste() { } CorpoCeleste::CorpoCeleste (const char *nomeCorpo, float mass, float xpos, float ypos, float vxi, float vyi) { Nome = new char[strlen(nomeCorpo)]; strcpy(Nome, nomeCorpo); m = mass; x = xpos; y = ypos; vx = vxi; vy = vyi; } void CorpoCeleste::calcolaPosizione(float fx, float fy, float t) { double ax = fx/m; double ay = fy/m; vx += ax*t; vy += ay*t; x += vx*t; y += vy*t; }

32. void CorpoCeleste::stampaPosizione() { cout.setf(ios::fixed); cout.setf(ios::showpos); cout << " " << setprecision(4) << setw(9) << x*1.e-11 << " " << setprecision(4) << setw(9) << y*1e-11 ; } void CorpoCeleste::stampaVelocita() { cout.setf(ios::fixed); cout.setf(ios::showpos); cout << " " << vx << " " << vy ; } CorpoCeleste::~CorpoCeleste() { } const char* CorpoCeleste::nome() {return Nome; } double CorpoCeleste::M() { return m; } double CorpoCeleste::X() { return x; } double CorpoCeleste::Y(){return y; } double CorpoCeleste::Vx(){return vx; } double CorpoCeleste::Vy(){return vy; } CorpoCeleste.cc Seconda parte

33. #ifndef CORPOCELESTE_H #define CORPOCELESTE_H class CorpoCeleste { protected: char *Nome; double m; double x; double y; double vx; double vy; public: CorpoCeleste() ; CorpoCeleste (const char *nomeCorpo, float mass,float xpos, float ypos, float vxi,float vyi); ~CorpoCeleste() { } ; virtual void calcolaPosizione(float fx, float fy, float t); void stampaPosizione(); void stampaVelocita(); const char *nome() {return Nome;}; double M() {return m;}; double X() {return x;}; double Y() {return y;}; double Vx() {return vx;}; double Vy() {return vy;}; }; #endif CorpoCeleste.h con impletazione di alcuni metodi “inline”