In general_distance.cpp tutti i file sono caricati, analizzati ed e' delegato al resto. Il file piu' interessante e' partitions.cpp, in cui la classe general_partition e' contenuta - tutte le operazioni tra partizioni, altamente ottimizzate, sono li' definite. In ising_simulation.cpp sono definite le classi per le simulazioni con generiche adiacenze. In adj_handler.cpp i file per caricare e gestire le possibili topologie.

Un paio di note:

I file in input sono in formato binario, intendo un file di interi a 32 bit, in cui ogni intero rappresenta un dato della configurazione, senza nessun spazio ne "a capo" tra le diverse configurazioni. In Matlab si ottiene molto semplicemente:

fid = fopen('configurazioni.bin','wb')\n fwrite(fid,matricione,'int32')

L'output, se sono richieste le matrici delle distanze, sono matrici binarie di double: -# output-shan.bin: distanza di Rohlin, usando entropia di Shannon, non ridotta -# output-shan_r.bin: come sopra, ma con la riduzione specificata (partizione comune o eliminazione diretta, epsilon 0 o no) -# output-top.bin: distanza di Rohlin, usando entropia topologica, non ridotta -# output-top_r.bin: come sopra, con la riduzione specificata, anche se non ha senso per epsilon > 0

fid = fopen('output-shan_r.bin','rb')\n dist=fread(fid,[dimensione,dimensione],'double');

Esempi:

distanze_generiche -torus 100 -file input.bin -num 1000 Legge dal file "input.bin" 1000 configurazioni, crea le partizioni sapendo che rappresentano delle configurazioni su un toro 100x100, calcola le complete matrici delle distanze Se il numero di configurazioni da leggere non e' specificato, controlla che il file "input.bin" abbia il numero giusto di dati, e ne carica tante quanto il file contiene.

distanze_generiche -sequence 566 -file input.bin -num 1400 -epsilon 3 Legge stavolta 1400 configurazioni, sempre da un file binario, in cui ogni stato e' rappresentato da un intero, considerate come catene lunghe 566. Le distanze sono calcolate con riduzione "a meno di epsilon uguale o minore a 3".

distanze_generiche -random -num 50 -symbols 2 -sierpinski 5 -common Crea 50 configurazioni random, usando 2 simboli (equivalente ad un Ising), su un sierpinski di generazione 5, utilizzando come metodo di riduzione, quello del confronto con una partizione comune

distanze_generiche -adj colonne.bin righe.bin -file configurazioni.bin -num 100 Carica le informazioni sull'adiacenza dai due file "colonne.bin" e "righe.bin", contenenti gli elementi nonnulli della matrice di adiacenza, legge 100 configurazioni e le trasforma in partizioni utilizzando le informazioni dell'adiacenza, e ne calcola tutte le distanze, usando come riduzione il metodo di default, la riduzione diretta tra atomi identici nelle due partizioni

distanze_generiche -sierpinski 8 -microcanonical -num 5000 -sweeps 10 -skip 1000 -beta 0.44 Su un Sierpinski gasket di generazione 8, esegue 5000 iterazioni di Ising con la dinamica di tipo microcanonico, con 10 full sweep del reticolo per istante di tempo, saltando 1000 istanti di tempo all'inizio, calcolando le distanze tra una configurazione e la successiva, ad un beta (approssimativamente) di 0.44

distanze_generiche -sierpinski 8 -microcanonical -num 5000 -sweeps 10 -skip 1000 -beta 0.44 -v -v Come prima, ma generando molto output, nei file "states.bin" e "partitions.bin" ci saranno 5000 configurazioni e partizioni rispettivamente, codificate come int32. Il file "output.txt" conterra' lo stato delle variabili iterazione per iterazione, "medie.txt" e "varianze.txt" le medie e varianze delle osservabili calcolate durante tutto il run del programma.

distanze_generiche -sierpinski 8 -microcanonical -num 5000 -sweeps 10 -skip 1000 -beta 0.3,1.0 -v -v Come prima, solo che un lato del gasket ha beta pari a 0.3, l'altro 1.0

distanze_generiche -sierpinski 8 -file states.bin -num 200 Calcola la matrici complete delle distanze tra le prime 200 configurazioni in output dal run precedente (ci mette davvero un attimo)

distanze_generiche -demo -square 1000 Mostra statistiche sulla velocita' per operazioni eseguite su un reticolo con periodicita' cilindrica, 1000x1000