84:
192:
182:
158:
74:
53:
22:
998:. In fact, they seem to be identical if F is not a function of v - in that case they only differ in the labeling of the nuisance variable v. Yet Leapfrog is a second-order method, while Euler-Cromer is first order. So something is strange here. I'm sure other readers than I may be confused, so a discussion of this might be nice to have in the article.
1017:
by Hairer et al. on the history and multitude of names of the Newton-Verlet-Stoermer-... method. Also on how to interpret the methods in terms of operator/vector field splittings. It seems that you are right, the time-shift of the velocity is the only difference, and it leads to a global O(h) error
343:
The article states that the leapfrog method "conserves a (slightly modified) energy of dynamical systems". This is a widely spread misconception and is simply not true, except in very special cases - most notably linear systems (harm. osc.). The symplecticity implies phase-space area preservation,
498:
The existence of a conserved "shadow
Hamiltonian" which converges to the exact Hamiltonian in the limit of vanishing time step is quite general. A rather straightforward proof is presented in section 13 of chapter 3 of Tuckerman's Statistical Mechanics textbook, which I have now referenced in the
990:
767:
1103:
Using RK methods for n-body gravitational simulators is not very useful as they do not conserve energy, so after a few orbits all bodies begin to spiral either inwards or outwards when N is larger than 3. Leapfrog conserves energy better, so provides much better long term simulations.
503:
mean that there is non-conservation. Most physical dynamical systems beyond of few particles exhibit chaotic dynamics, and this is precisely why symplectic integrators are used in dynamical simulations; they help to minimize energy drift by sticking to the shadow
284:
This term can also be found in the field of electronic filter realisation. A standard passive filter can be actively simulated (called functional simulation) by interchanging current and voltage in a "leapfrog" way in the equations describing the circuit.
776:
553:
781:
558:
1166:
1204:
1081:
262:
418:
1199:
480:
441:
140:
1214:
252:
985:{\displaystyle {\begin{aligned}a_{i}&=F(x_{i})\\v_{i+1/2}&=v_{i-1/2}+a_{i}\,\Delta t,\\x_{i+1}&=x_{i}+v_{i+1/2}\,\Delta t,\\\end{aligned}}}
762:{\displaystyle {\begin{aligned}x_{i}&=x_{i-1}+v_{i-1/2}\,\Delta t,\\a_{i}&=F(x_{i})\\v_{i+1/2}&=v_{i-1/2}+a_{i}\,\Delta t,\end{aligned}}}
1194:
130:
344:
yes - but not necessarily
Hamiltonian conservation. Instead, typical systems even in one dim are non-integrable, and display more or less chaos.
214:
365:
Please propose a better formulation. If I read Hairer right, it should at least be true that a modified
Hamiltonian is preserved to the order
83:
1209:
356:
238:
1189:
1137:
1111:
324:
1170:
106:
529:
205:
163:
301:
97:
58:
33:
482:. KAM style chaos is perfectly compatible with Hamiltionian preservation, no energy is being destroyed or created.--
995:
1021:
1141:
1115:
352:
328:
1018:
in the velocity. The position in Euler-Cromer should still be O(h²) for situations that can be formulated as
368:
21:
348:
39:
191:
1162:
1133:
1107:
320:
289:
446:
1130:
This method gives me different results to the midstep method even with constant timestep. Flaws?
213:
on
Knowledge. If you would like to participate, please visit the project page, where you can join
105:
on
Knowledge. If you would like to participate, please visit the project page, where you can join
1014:
197:
181:
157:
1003:
526:
297:
423:
1088:
487:
244:
89:
1183:
1157:
999:
293:
210:
73:
52:
1174:
1145:
1119:
1092:
1007:
491:
332:
305:
1084:
483:
187:
79:
443:), and that the difference between original and modified Hamiltonian is
102:
247:
in the banner shell. Please resolve this conflict if possible.
243:
This article has been given a rating which conflicts with the
15:
523:
Statistical
Mechanics: Theory and Molecular Simulation
1024:
779:
556:
449:
426:
371:
209:, a collaborative effort to improve the coverage of
101:, a collaborative effort to improve the coverage of
499:article. A system exhibiting chaotic dynamics does
1158:http://www.artcompsci.org/vol_1/v1_web/node34.html
1075:
984:
761:
474:
435:
412:
1205:C-Class articles with conflicting quality ratings
317:Should it be stated explicitly that v = dx/dt?
525:(1 ed.). Oxford University Press. pp. 121–124.
1200:Stub-Class physics articles of Low-importance
1076:{\displaystyle {\ddot {x}}(t)=\nabla P(x(t))}
8:
1160:
1131:
1105:
152:
47:
19:
1026:
1025:
1023:
957:
947:
934:
911:
886:
869:
859:
838:
828:
811:
788:
780:
778:
738:
721:
711:
690:
680:
663:
640:
611:
601:
582:
565:
557:
555:
463:
448:
425:
401:
382:
370:
548:The algorithm described in the article,
1167:2A02:120B:2C67:1960:1C5D:82C3:B66A:106F
967:
892:
744:
621:
513:
154:
49:
413:{\displaystyle O(e^{Lt}*\Delta t^{4})}
7:
203:This article is within the scope of
95:This article is within the scope of
420:(time symmetry kills odd powers of
38:It is of interest to the following
1049:
969:
894:
746:
623:
456:
427:
394:
245:project-independent quality rating
14:
1215:Low-priority mathematics articles
223:Knowledge:WikiProject Mathematics
226:Template:WikiProject Mathematics
190:
180:
156:
82:
72:
51:
20:
1195:Low-importance physics articles
475:{\displaystyle O(\Delta t^{2})}
257:This article has been rated as
135:This article has been rated as
1070:
1067:
1061:
1055:
1043:
1037:
817:
804:
669:
656:
469:
453:
407:
375:
347:— Preceding comment added by
1:
1146:13:19, 19 February 2019 (UTC)
1120:12:24, 19 February 2019 (UTC)
280:Use in electronic engineering
217:and see a list of open tasks.
115:Knowledge:WikiProject Physics
109:and see a list of open tasks.
1210:C-Class mathematics articles
1175:08:32, 13 January 2022 (UTC)
306:18:05, 15 January 2009 (UTC)
118:Template:WikiProject Physics
1190:Stub-Class physics articles
521:Tuckerman, Mark E. (2010).
1231:
1156:This reference is broken:
994:looks very similar to the
544:Contrast with Euler-Cromer
492:10:53, 10 March 2014 (UTC)
141:project's importance scale
1093:14:57, 23 June 2014 (UTC)
1008:12:50, 23 June 2014 (UTC)
333:19:40, 30 July 2009 (UTC)
256:
242:
175:
134:
67:
46:
436:{\displaystyle \Delta t}
359:) 20:14, 7 February 2013
263:project's priority scale
206:WikiProject Mathematics
1126:Kick-Drift-Kick method
1077:
986:
763:
476:
437:
414:
28:This article is rated
1078:
987:
764:
477:
438:
415:
1022:
777:
554:
447:
424:
369:
229:mathematics articles
1013:Read the very good
996:Euler-Cromer method
98:WikiProject Physics
1073:
982:
980:
968:
893:
759:
757:
745:
622:
472:
433:
410:
198:Mathematics portal
34:content assessment
1177:
1165:comment added by
1148:
1136:comment added by
1122:
1110:comment added by
1034:
361:
339:conservation of H
323:comment added by
309:
292:comment added by
277:
276:
273:
272:
269:
268:
151:
150:
147:
146:
1222:
1152:Broken reference
1082:
1080:
1079:
1074:
1036:
1035:
1027:
991:
989:
988:
983:
981:
966:
965:
961:
939:
938:
922:
921:
891:
890:
878:
877:
873:
847:
846:
842:
816:
815:
793:
792:
771:or equivalently
768:
766:
765:
760:
758:
743:
742:
730:
729:
725:
699:
698:
694:
668:
667:
645:
644:
620:
619:
615:
593:
592:
570:
569:
536:
535:
518:
481:
479:
478:
473:
468:
467:
442:
440:
439:
434:
419:
417:
416:
411:
406:
405:
390:
389:
360:
345:
335:
308:
286:
231:
230:
227:
224:
221:
200:
195:
194:
184:
177:
176:
171:
168:
160:
153:
123:
122:
121:physics articles
119:
116:
113:
92:
87:
86:
76:
69:
68:
63:
55:
48:
31:
25:
24:
16:
1230:
1229:
1225:
1224:
1223:
1221:
1220:
1219:
1180:
1179:
1154:
1128:
1101:
1020:
1019:
979:
978:
943:
930:
923:
907:
904:
903:
882:
855:
848:
824:
821:
820:
807:
794:
784:
775:
774:
756:
755:
734:
707:
700:
676:
673:
672:
659:
646:
636:
633:
632:
597:
578:
571:
561:
552:
551:
546:
541:
540:
539:
532:
520:
519:
515:
459:
445:
444:
422:
421:
397:
378:
367:
366:
349:130.235.189.247
346:
341:
318:
315:
313:definition of v
287:
282:
228:
225:
222:
219:
218:
196:
189:
169:
166:
120:
117:
114:
111:
110:
88:
81:
61:
32:on Knowledge's
29:
12:
11:
5:
1228:
1226:
1218:
1217:
1212:
1207:
1202:
1197:
1192:
1182:
1181:
1153:
1150:
1138:131.227.74.101
1127:
1124:
1112:131.227.74.101
1100:
1099:Superior to RK
1097:
1096:
1095:
1072:
1069:
1066:
1063:
1060:
1057:
1054:
1051:
1048:
1045:
1042:
1039:
1033:
1030:
977:
974:
971:
964:
960:
956:
953:
950:
946:
942:
937:
933:
929:
926:
924:
920:
917:
914:
910:
906:
905:
902:
899:
896:
889:
885:
881:
876:
872:
868:
865:
862:
858:
854:
851:
849:
845:
841:
837:
834:
831:
827:
823:
822:
819:
814:
810:
806:
803:
800:
797:
795:
791:
787:
783:
782:
754:
751:
748:
741:
737:
733:
728:
724:
720:
717:
714:
710:
706:
703:
701:
697:
693:
689:
686:
683:
679:
675:
674:
671:
666:
662:
658:
655:
652:
649:
647:
643:
639:
635:
634:
631:
628:
625:
618:
614:
610:
607:
604:
600:
596:
591:
588:
585:
581:
577:
574:
572:
568:
564:
560:
559:
545:
542:
538:
537:
530:
512:
511:
507:
506:
505:
495:
494:
471:
466:
462:
458:
455:
452:
432:
429:
409:
404:
400:
396:
393:
388:
385:
381:
377:
374:
340:
337:
325:128.100.76.206
314:
311:
281:
278:
275:
274:
271:
270:
267:
266:
255:
249:
248:
241:
235:
234:
232:
215:the discussion
202:
201:
185:
173:
172:
161:
149:
148:
145:
144:
137:Low-importance
133:
127:
126:
124:
107:the discussion
94:
93:
90:Physics portal
77:
65:
64:
62:Low‑importance
56:
44:
43:
37:
26:
13:
10:
9:
6:
4:
3:
2:
1227:
1216:
1213:
1211:
1208:
1206:
1203:
1201:
1198:
1196:
1193:
1191:
1188:
1187:
1185:
1178:
1176:
1172:
1168:
1164:
1159:
1151:
1149:
1147:
1143:
1139:
1135:
1125:
1123:
1121:
1117:
1113:
1109:
1098:
1094:
1090:
1086:
1064:
1058:
1052:
1046:
1040:
1031:
1028:
1016:
1012:
1011:
1010:
1009:
1005:
1001:
997:
992:
975:
972:
962:
958:
954:
951:
948:
944:
940:
935:
931:
927:
925:
918:
915:
912:
908:
900:
897:
887:
883:
879:
874:
870:
866:
863:
860:
856:
852:
850:
843:
839:
835:
832:
829:
825:
812:
808:
801:
798:
796:
789:
785:
772:
769:
752:
749:
739:
735:
731:
726:
722:
718:
715:
712:
708:
704:
702:
695:
691:
687:
684:
681:
677:
664:
660:
653:
650:
648:
641:
637:
629:
626:
616:
612:
608:
605:
602:
598:
594:
589:
586:
583:
579:
575:
573:
566:
562:
549:
543:
533:
531:9780198525264
528:
524:
517:
514:
510:
502:
497:
496:
493:
489:
485:
464:
460:
450:
430:
402:
398:
391:
386:
383:
379:
372:
364:
363:
362:
358:
354:
350:
338:
336:
334:
330:
326:
322:
312:
310:
307:
303:
299:
295:
291:
279:
264:
260:
254:
251:
250:
246:
240:
237:
236:
233:
216:
212:
208:
207:
199:
193:
188:
186:
183:
179:
178:
174:
165:
162:
159:
155:
142:
138:
132:
129:
128:
125:
108:
104:
100:
99:
91:
85:
80:
78:
75:
71:
70:
66:
60:
57:
54:
50:
45:
41:
35:
27:
23:
18:
17:
1161:— Preceding
1155:
1132:— Preceding
1129:
1106:— Preceding
1102:
993:
773:
770:
550:
547:
522:
516:
508:
500:
342:
316:
283:
259:Low-priority
258:
204:
170:Low‑priority
136:
96:
40:WikiProjects
319:—Preceding
288:—Preceding
220:Mathematics
211:mathematics
164:Mathematics
1184:Categories
509:References
30:Stub-class
504:manifold.
1163:unsigned
1134:unsigned
1108:unsigned
357:contribs
321:unsigned
302:contribs
290:unsigned
1000:Amaurea
294:Aphexer
261:on the
167:C‑class
139:on the
112:Physics
103:Physics
59:Physics
36:scale.
1085:LutzL
1015:paper
484:LutzL
1171:talk
1142:talk
1116:talk
1089:talk
1004:talk
527:ISBN
488:talk
353:talk
329:talk
298:talk
1083:.--
501:not
253:Low
131:Low
1186::
1173:)
1144:)
1118:)
1091:)
1050:∇
1032:¨
1006:)
970:Δ
895:Δ
864:−
747:Δ
716:−
624:Δ
606:−
587:−
490:)
457:Δ
428:Δ
395:Δ
392:∗
355:•
331:)
304:)
300:•
1169:(
1140:(
1114:(
1087:(
1071:)
1068:)
1065:t
1062:(
1059:x
1056:(
1053:P
1047:=
1044:)
1041:t
1038:(
1029:x
1002:(
976:,
973:t
963:2
959:/
955:1
952:+
949:i
945:v
941:+
936:i
932:x
928:=
919:1
916:+
913:i
909:x
901:,
898:t
888:i
884:a
880:+
875:2
871:/
867:1
861:i
857:v
853:=
844:2
840:/
836:1
833:+
830:i
826:v
818:)
813:i
809:x
805:(
802:F
799:=
790:i
786:a
753:,
750:t
740:i
736:a
732:+
727:2
723:/
719:1
713:i
709:v
705:=
696:2
692:/
688:1
685:+
682:i
678:v
670:)
665:i
661:x
657:(
654:F
651:=
642:i
638:a
630:,
627:t
617:2
613:/
609:1
603:i
599:v
595:+
590:1
584:i
580:x
576:=
567:i
563:x
534:.
486:(
470:)
465:2
461:t
454:(
451:O
431:t
408:)
403:4
399:t
387:t
384:L
380:e
376:(
373:O
351:(
327:(
296:(
265:.
239:C
143:.
42::
Text is available under the Creative Commons Attribution-ShareAlike License. Additional terms may apply.