SCM Repository
| [diderot] / branches / pure-cfg / test / lic.diderot |
Diff of /branches/pure-cfg/test/lic.diderot
Parent Directory
| 
Revision Log
| 
Patch
| revision 1381, Thu Jun 23 19:20:18 2011 UTC | revision 1382, Thu Jun 23 20:03:05 2011 UTC | |
|---|---|---|
| # | Line 28 | Line 28 |
| 28 | // forw = forw + h*V(forw); | // forw = forw + h*V(forw); |
| 29 | // back = back - h*V(back); | // back = back - h*V(back); |
| 30 | // Midpoint method step | // Midpoint method step |
| 31 | forw = forw + h*V(forw + 0.5*h*V(forw)); | forw += h*V(forw + 0.5*h*V(forw)); |
| 32 | back = back - h*V(back - 0.5*h*V(back)); | back -= h*V(back - 0.5*h*V(back)); |
| 33 | sum = sum + R(forw) + R(back); | sum += R(forw) + R(back); |
| 34 | step = step + 1; | step += 1; |
| 35 | if (step == stepNum) { | if (step == stepNum) { |
| 36 | // modulate output by velocity at initial position | // modulate output by velocity at initial position |
| 37 | sum = |V(pos0)|*sum/real(1 + 2*stepNum); | sum = |V(pos0)|*sum/real(1 + 2*stepNum); |
|
||||||||
| root@smlnj-gforge.cs.uchicago.edu | ViewVC Help |
| Powered by ViewVC 1.0.0 |
