FuckAuthority
19-03-2010, 09:35
Ho un problemino che magari qualche ingegnere un pò esperto in calcolo numerico sa meglio di me :D
Sto sviluppando un codice in python che calcola pressioni, portate e wss in una rete vascolare (generalmente del braccio) dato come ingresso un segnale di un flusso pulsatile modellato con una serie di Fourier. Come condizioni al contorno impongo pressioni transmurali dei vasi nulli, e ogni vaso è modellato come un circuito RCL (womersley model).
Il mio problema è che facendo girare il codice su pc a 32 bit e su pc a 64 bit ottengo risultati leggermente diversi. Le medie sono praticamente molto simili, ma i valori cambiano (soprattutto nei picchi).
E' un problema di "rumore numerico" dovuto alle diverse architetture, ma sto impazzendo e non so come gestirlo.
Il calcolo numerico risolve semplicemente un sistema differenziale matriciale con 3 matrici (ordine zero, uno, due) tramite il metodo del trapezio del primo ordine ( o equivalentemente con Newmark, non cambia nulla...)
help, per anubi!
Sto sviluppando un codice in python che calcola pressioni, portate e wss in una rete vascolare (generalmente del braccio) dato come ingresso un segnale di un flusso pulsatile modellato con una serie di Fourier. Come condizioni al contorno impongo pressioni transmurali dei vasi nulli, e ogni vaso è modellato come un circuito RCL (womersley model).
Il mio problema è che facendo girare il codice su pc a 32 bit e su pc a 64 bit ottengo risultati leggermente diversi. Le medie sono praticamente molto simili, ma i valori cambiano (soprattutto nei picchi).
E' un problema di "rumore numerico" dovuto alle diverse architetture, ma sto impazzendo e non so come gestirlo.
Il calcolo numerico risolve semplicemente un sistema differenziale matriciale con 3 matrici (ordine zero, uno, due) tramite il metodo del trapezio del primo ordine ( o equivalentemente con Newmark, non cambia nulla...)
help, per anubi!