Dovrei (sorry ) aver corretto tutti gli errori, ho impostato anche la pressione in bar ... in doppia precisione va in overflow con indice nitertemp>23 . E' possibile far uscire il log su file, e sono pronto per generare un foglio excel alla fine dell'esecuzione e con la math precision...ti gira l'eseguibile con la dll sotto?

ma prova a mettere deltat =0,000001

non ho la possibilità...c'è ancora qualche errore?
la velocità rispetto il foglio?
cmq la precisione del double vb fa veramente pena...oltre l'esponente volevo vedere di quanto miglioravano i risultati...a discapito sicuramente della velocità...ti gira MathPrecision.exe?

C'è ancora il problema della memoria ( bombicchia...prova ad aumentare gli indici...)

l'eseguibile gira sembra giusto ora lo controllo numer per numero


comunque rispetto a calc (openoffice.org) la velocità è x200 o più, rispetto ad excel non so un ipoteticamente X 50

ok ho controllato è come il mio quindi o sono entrambi giusti o entrambi sbagliati


comunque se p =0 i soldatini marciano indisturbati OK

ora si tratta di vedere se per un tempo abbastanza lungo la colonna accelera all'indietro secondo le usuali leggi ma bisogna simulare almeno per 0,5sec...

implemeto la sub evoluzione N=


quindi il tempo di interazione totale sarà niter*N*deltat


è ridicolo come dei computer ritenuti velocissimi "mostrano subito le corde" se usati per una cosa un minimo seria come le simulazioni numeriche

bene bene,
mi sarei aspettato un risultato migliore di 200x cmq son curioso di capire quanto migliora con VC6...
oggi sono rimasto incredulo sulla precisione dei double in vb ...si impiccia terribilmente...sto andando a fondo per capire se questo errore è presente su tutta la piatttaforma visual...cmq è pronta:
1: la versione con precisione quadrupla ( o anche più con precisione a scelta)...non ho capito se ti gira questo eseguibile:LINK per me è importante capirlo visto che usa una dll esterna
2: genera (ops dovrebbe generare) un file excel tale e quale a quello che genera il tuo codice (l'utilità è che dovrebbe permettere di fare calcoli grafici etc)
3: lo 'spezzettamento' della memoria allocato dinamicamente


CITAZIONE

Ho allunghigghiato il vettore ...questo è il valore massimo consecutivo di memoria che riesco ad allocare cioè 10000x10000...poi passo a spezzettarlo e riidicizzarlo. Funzica l'allocazione di memoria anche su altri pc?
Ho spostato nitertemp ...l'ho messo sotto il vettore psi...ho eliminato la parte testuale ( non ce la fa la textbox devo passare alla richbox...mi occorre sapere se gira la dll importata) ...il log esce su file ...(può essere importato su excel fino a che non lo genera l'exe stesso).
Sul mio pc l'elaborazione di default impiega 25 secondi...considerate che un pc molto veloce...

for I=1 to 9999

psi(0,I)= L0+psi(0,I-1) 'prima colonna nel foglio "A"
AAA(I)=psi(0,I)
next I
for I=0 to 9999

BBB(I)= AAA(I)-V0*deltat

next I
'girotondo 0

CCC(0) = 2*BBB(0)-AAA(0)+1/Rho*(deltat/L0)^2*(P*L0+B*(BBB(1)-BBB(0)-L0))
for I=1 to 9998
CCC(I) = 2*BBB(I)-AAA(I)+B/Rho*(deltat/L0)^2*(BBB(I+1)-2*BBB(I)+BBB(I-1))
next I
CCC(9999) = 2*BBB(9999)-AAA(9999)+B/Rho*(deltat/L0)^2*(BBB(9999)-BBB(9998)-L0)

for j=1 to nitertemp
'girotondo 1

For K = 0 to N-2

AAA(0) = 2*CCC(0)-BBB(0)+1/Rho*(deltat/L0)^2*(P*L0+B*(CCC(1)-CCC(0)-L0))
for I=1 to 9998
AAA(I) = 2*CCC(I)-BBB(I)+B/Rho*(deltat/L0)^2*(CCC(I+1)-2*CCC(I)+CCC(I-1))
next I
AAA(9999) = 2*CCC(9999)-BBB(9999)+B/Rho*(deltat/L0)^2*(CCC(9999)-CCC(9998)-L0)

'girotondo 2

BBB(0) = 2*AAA(0)-CCC(0)+1/Rho*(deltat/L0)^2*(P*L0+B*(AAA(1)-AAA(0)-L0))
for I=1 to 9998
BBB(I) = 2*AAA(I)-CCC(I)+B/Rho*(deltat/L0)^2*(AAA(I+1)-2*AAA(I)+AAA(I-1))
next I
BBB(9999) = 2*AAA(9999)-CCC(9999)+B/Rho*(deltat/L0)^2*(AAA(9999)-AAA(9998)-L0)
'girotondo extra

CCC(0) = 2*BBB(0)-AAA(0)+1/Rho*(deltat/L0)^2*(P*L0+B*(BBB(1)-BBB(0)-L0))
for I=1 to 9998
CCC(I) = 2*BBB(I)-AAA(I)+B/Rho*(deltat/L0)^2*(BBB(I+1)-2*BBB(I)+BBB(I-1))
next I
CCC(9999) = 2*BBB(9999)-AAA(9999)+B/Rho*(deltat/L0)^2*(BBB(9999)-BBB(9998)-L0)


next k

for I=0 to 9999
psi(j,I)=CCC(I)
next I

If j=1 Then N=N+1


next j

azz...hai già fatto?

eccomi...versione implementata del girotondo (Go girotondo!), ho modificato alcuni parametri di default per evitare l'overflow, funge sempre il GO (double)!.
Inoltre, a fine generazione , se il flag LogFile Excel è settato, chiede se si vuole la generazione del file LogFile.xls ... va lanciato da windows e apre il foglio...mi serve un riscontro da parte di qualcuno se gira su un pc 'normale'...

frena frena non ho ancora fatto in tempo a testarlo


il girotondo serve associato al numero N a fare molte interazioni ma salvare sulla matrice solo quelle necessarie


ad es se si mette N=5000

nintertemp=500


e deltat =0,000 000 010 secondi

si ha una memorizzazione ogni 15 000 deltat


quindi ogni 0,15ms

per 75ms in totale


CITAZIONE

azz incredibile azz azz azz FUNZIONA FUNZIONA


nel senso che i risultati con il girotondo sono esatti ma va in overflow ... bisogna indagare

sarebbe bene stampare sul txt anche l'ultima particella comunque faccio anche altri controlli

miiiiiizzeca...

senti, ma riesci a fare la prova se ti gira MathPrecision.exe? (LINK)se non lo vuoi fare...non ci sono problemi, sul prossimo eseguibile faccio un sotto processo che lancia le operazioni, salva il risultato su una stringa e la prima volta che ti connetti lo invia al server del forum...dai...mi fai quella prova?

A parte gli scherzi, la precisione dei double sui prodotti microsoft fa pietà e misericordia...se nel foglio ottieni gli stessi valori dell'eseguibile...non può andare...troppa tolleranza, più si aumentano le operazioni e più il risultato va a cavolfiore...

A proposito...il foglio excel lo genera? funzica?

ok per l'ultima particella (di sodio)!

CITAZIONE (qqcreafis, 04/02/2015 23:39:18 )

Versione 1.6 ...viene intercettato l'errore di overflow ( se necessario si può aumentare l'estensione dell'esponente facendo fare calcoli in precisione maggiore) ed aggiunta l'ultima particella sia sul file txt che sul file excel.

Screenshot:

CITAZIONE (qqcreafis, 05/02/2015 12:21:26 )