Risoluzione numerica dell'equazione scalare delle onde mediante ...
Risoluzione numerica dell'equazione scalare delle onde mediante ...
Risoluzione numerica dell'equazione scalare delle onde mediante ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
PDE Laboratory @ Crs4 8!---------------------------------------------------------------------! Calcolo uu(t) n+1!---------------------------------------------------------------------do i = 1,nnoduu(i) = k2*(rhs(i) - KUn(i))/mm(i) + hh(i)enddoSui punti della frontiera di Dirichlet viene imposta la soluzione analitica del problema.! Correction for Dirichlet nodesdo ie=1,nbardo i=1,2irig=con_edge(ie,i)uu(irig) = sol_an(xx(irig),yy(irig),t)uu0(irig) = sol_an(xx(irig),yy(irig),t-dt)enddoenddoLa valutazione dell’errore in norma infinito è svolta dal codice seguente:! Error evaluation; <strong>numerica</strong>l and exact solutions are written in ’sol.d’!---------------------------------------------------------------------allocate (uu_an(nnod))do i=1,nnoduu_an(i)=sol_an(xx(i),yy(i),t_fin)enddonum=0.d0; den=0.d0do i=1,nnoddiff=dabs(uu(i)-uu_an(i))if (diff.gt.num) thennum=diff4 luglio 2005 Propagazione <strong>delle</strong> <strong>onde</strong>
Change language