An algorithm for the solution of the Fermi-Pasta-Ulam problem on a vector processor

E' presentato un algoritmo per it calcolo della soluzione del problema Fermi-Pasta-Ulam (una stringa monodimensionale di particelle accoppiate con forze non lineari) su elaboratori vettoriali. L'algoritmo è stato provato su un Array Processor FPS-100 utilizzando il linguaggio CPFORTRAN con chiamate alle routine vettoriali. Un'analisi dettagliata dell'algoritmo mostra che quando l'ordine del sistema di equazioni differenziali ordinarie aumenta l'efficienza diventa unitaria. Risulta che la dimensione per la quale si ottiene la metà della velocità massima é n_1/2= 120 e che si supera il 90% della velocità massima del processore (8 Mega FLOPS) se l'ordine del sistema è maggiore di 1800. L'algoritmo può essere facilmente portato su altri elaboratori vettoriali ed il metodo può essere usato per risolvere problemi simili.