Il 23/01/2014 19:45, Paolo Russo ha scritto:
> [marcofuics:]
>> suppongo che tu voglia soluzioni CAS
>
> Uhm. Column address strobe, no; compare and swap, no;
> complicazione affari semplici, preferirei evitarla. :-)
>
Ciao, di programmi in C o C++ non ne conosco, in fortran si'. In ogni
caso, come fatto notare da marcofuics, la complessita' del programma
dipende fortemente da quali funzioni supponi siano disponibili `a
gratis', in altre parole se ti va bene usare librerie (le boost in C)
per cose tipo trovare autovalori/autovettori o se vuoi un programma del
tutto self-contained. In ogni caso il metodo dell'articolo citato (che
ho letto molto velocemente) e' semplice concettualmente ma NON
efficiente per potenziali generici! calcolare f(X) dove X e' una matrice
costa un sacco di operazioni e ci sono metodo molto piu' efficienti. In
particolare, il metodo citato mi pare ad occhio una riscoperta del
metodo DVR (discrete variable representation) cosi' come era stato
scoperto negli anni '60. Se consideri disponibile una routine per
autovalori/vettori di una matrice uno dei metodi migliori IMHO e' il
cosiddetto DVR Colbert-Miller (formula A7 di
http://xbeams.chem.yale.edu/~batista/v572/ColbertMiller.pdf , e aggiungi
il potenziale sulla diagonale della matrice), altrimenti se vuoi fare
tutto tu da solo considera raffinamenti del classico metodo di Numerov,
vedi ad es.
http://scitation.aip.org/content/aip/journal/jcp/88/6/10.1063/1.453880
In ogni caso, a cosa ti serve? Qual'e' la forma del potenziale? Ha
singolarita'? Ti interessano solo gli stati legati? Risonanze?
L.
Received on Tue Jan 28 2014 - 23:14:38 CET