Home My Page Projects Code Snippets Project Openings diderot
Summary Activity Tracker Tasks SCM

SCM Repository

[diderot] Diff of /branches/fem/test/lic.diderot
ViewVC logotype

Diff of /branches/fem/test/lic.diderot

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1115, Thu May 5 04:42:18 2011 UTC revision 1330, Mon Jun 13 00:13:54 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);

Legend:
Removed from v.1115  
changed lines
  Added in v.1330

root@smlnj-gforge.cs.uchicago.edu
ViewVC Help
Powered by ViewVC 1.0.0