729:
nature of the technique draws upon theoretical concepts in the domain of physics known as analytical dynamics. In "modern times" advanced analytical dynamics has become merged with differential geometry, topology, functional analysis, and related topics and at its higher realization seems to be in the process of being absorbed by what is called geometric mechanics. Einstein's general relativity is an early example of the merger of dynamics with differential geometry, but it took a few decades before this took off as a general trend. At about the time it did take off, the higher reaches of mathematics had started to transition from proofs dominated by set theoretic concepts to to being ever more reliant upon the reconceptualization of everything in terms of differential geometry, topology and other related disciplines which have been ascendant for some time. As a result, there is nowadays an extensive reliance upon multiple advanced vocabularies and concept sets which can, if used, render the simplest ideas quite incomprehensible to most readers. It is my personal observation (I avoid calling it merely a belief) that projects like
Knowledge have been negatively impacted by over use of esoteric and excessively abstract representations of things. You should not need multiple PhDs to extract useful content from every article on physics in particular. For a very on-point, and but harsher, statement of this position addressed to a professional audience, see the nicely worded Introduction in this paper: David Hestenes, Hamiltonian Mechanics with Geometric Calculus In: Z. Oziewicz et al (eds.), Spinors, Twistors, Clifford Algebras and Quantum Deformations, Kluwer: Dordercht/Boston (1993), 203–214. The article in isolation was found at
719:
a discussion on the plague as such at this time. As it is, what I will say about it may well inspire well meaning devotees to characterize my remarks as demonstrating too much Point of View. I'll take that risk for three reasons: (1) I am a PhD physicist and so I have a POV that is more professional than personal,(2) the problem I am going to address is a growing quality and policy problem among a large set of math, physics, and similar articles and (3) the question leading off this section is not unlike others I have seen and reflects the same comprehensibility issues. Despite my education, I too have found many articles to be incomprehensible. I am one guy with a modest amount of time on my hands and I can only insert myself into a few articles. Perhaps by improving this one and explaining why and what policy ideas guided it, I can inspire others similarly situated to do the same.
84:
744:). Anyway, while I doubt it will be particularly necessary to improve this article, I welcome the collaboration of anyone with a stronger background in differential geometry and/or geometric mechanics. Having spent serious time in front of a blackboard teaching physics to undergraduates, however, I am confident that I can distinguish incomprehensible explanations from comprehensible ones and that things do not always have to be dumbed down to be made comprehensible. Pedagogy is about being clear and complete and knowing when an advanced concept is really helpful versus justified only to refer the reader to sources of deeper insight.
189:
179:
158:
308:. I must admit that the form used to describe the algorithm here is elegant and short. It is quite likely the prefered form for mathematicians. But in order to broaden the audience I would ask anyone reading this comment and being able to translate this article into more traditional and ancient algorithm language to perform this translation. Assuming such a translation is possible. If not I beg for a short explanation of the equation
74:
53:
22:
743:
program. Like pretty much any PhD graduate, I nevertheless have enough training to proceed further as needed and unravel things as necessary. (P.S. my advisor (condensed matter theory), brilliant though he was, did not support my taking modern math courses though I did sneak off an take at least one<g: -->
945:
Note: I now at least inserted a remark in the 'Examples' section stating that the index is always traversed in decreasing order. It would be, however, much better, if the operator statements were rephrased so they have the ordering with the lowest index value to the far right (applied first), because
718:
I have only been examining this particular
Knowledge article since somewhat over a month ago. In my humble -- but I think widely shared -- perception it is an example, though a relatively mild one, of a plague that has descended upon articles at this level throughout Knowledge. I don't want to start
708:
Its quite a few years since the above question was asked. I think this article is one of many I have seen which are less comprehensible than they should be. In the spirit of
Knowledge have decided to be bold and clean up this one. I am first going to insert some personal insights that motivate what I
834:
The application "In Plasma
Physics" seems oddly specific. Especially with involvement of Yoshida's trick used to construct higher order splitting. Could Yoshida's trick be discussed in a separate section? Could there perhaps be some example of a practical implementation (pseudo-code)? Could there be
742:
I am a non-specialized theoretical physicist. Despite that in graduate school my department included a mathematical physicist of some renown, I was only exposed to differential geometry in the context of general relativity and had no prior background in it prior to my late-in-life entry into a PhD
896:
As stated in the page, the proposed second-order example should match the Verlet method. However, by doing the math, it seems that this is only true with c1 = 1 and c2 = 0, and not the other way around (d1 and d2 are OK). I tried to use c1 = 0 and c2 = 1 (as it is currently stated in the article),
728:
Symplectic
Integrators are really an applied physics topic or, more specifically, computational physics, rather than one considered within theoretical physics. Mathematicians also have plenty to say about symplectic spaces, but its physicists and research engineers who use those integrators. The
870:
Need to add: partitioned Runge-Kutta, Gauss-Legendre, history, connection with generating functions, advantages of symplectic integrators like KAM / Nekhoroshev results, approximate energy preservation, linear error growth. Problems with variable step size. Mention that there are more splitting
784:
I am glad that someone said what needed to be said. Yes, this is a plague. I too have a PhD in physics. Why do I have to go on a scavenger hunt of obscure terms like "wedge product" and "cohomology" to figure out what an article is saying? Knowledge of all places should default to a lower level
812:
I am a mere pusher of 1s and 0s, and don't really know what the advantage is of mostly conserving the original
Hamiltonian. In addition to the statement above, I would love to see a more "practical" statement of the virtues. Maybe something like: "the error is therefore of the order (dt^n)" or
901:
I was also confused by this
However, a look at previous equations reveals that the higher order indices are applied first! Since they are operators, they act from right to left. Therefore, Verlet's method is recovered by first setting c=1, d=1/2, then again with c=0, d=1/2. Therefore, the
303:
This article is supposed to be of "mid-importance within physics" however it is written with "advanced" mathematic notations. As a physician, I can not understand it and I find this quite unfortunate. Conversely it is perfeclty easy to read and implement algorithm as describe in page like
276:
What happens if you try to use a symplectic integrator whose time evolution cannot be derived from a hamiltonian (eg. if the system also includes a dissipative force)? I think this should be mentioned in the article.
287:
A symplectic integrator is defined in terms of a
Hamiltonian, H. If there is no H (as you say, the system cannot be derived from a Hamiltonian) then there is simply no integrator to iterate. Does that make sense?
245:
665:
902:
coefficients are ok, but the "simplified equations" are, in my opinion, wrong. They should regressive equations, giving e.g. v_{i-1} in terms of x_i and v_i, to be run from i=k to i=1. --
582:
987:
485:
523:
344:
918:
This seems an important **open issue** in this page. I ran across it when implementing the fourth-order version, and running into trouble. I then saw this discussion:
426:
685:
446:
399:
381:
is defined to name the set of position coordinate q and momentum coordinate p. For example, if you are thinking about the motion of a particle along a coordinate
379:
140:
813:"numeric stability is improved significantly over non-SI methods", or whatever will help me compare this usefully to say Runge-Kutta or other non-SI methods.
919:
997:
809:"Symplectic integrators possess as a conserved quantity a Hamiltonian which is slightly perturbed from the original one. By virtue of these advantages".
235:
753:
I start soon and will go slowly. Comments on the above will be cheerfully received, but lets spend most of our time working on article content. Thanks
982:
130:
946:
that would then make everything below very easy to follow. It would also require a change in the order of the listed coefficients further below.
211:
83:
992:
897:
but this results in integration steps that I can clearly see that do not hold for constant accelerations. Am I interpreting something wrong?
820:
786:
347:
922:...which agrees with your statement that the order should be reversed. How to best do this without changing too much at the same time?
687:
was introduced only to simplify the expression so that the essense of the derivation become clear. I added an explanation in the text.
953:
929:
106:
977:
202:
163:
864:, these subpages are now deprecated. The comments may be irrelevant or outdated; if so, please feel free to remove this section.
920:
https://scicomp.stackexchange.com/questions/20533/test-of-3rd-order-vs-4th-order-symplectic-integrator-with-strange-result
97:
58:
33:
587:
824:
790:
351:
957:
933:
876:
528:
39:
188:
840:
487:. If you are thinking about a particle in a 3D space and representing its position by the coordinates
949:
925:
816:
693:
21:
907:
305:
210:
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
836:
785:
language. Why is it that every physics article is written like a graduate-level math textbook???
688:
194:
178:
157:
278:
872:
766:
758:
451:
490:
311:
404:
273:"Symplectic integrators are designed for the numerical solution of Hamilton's equations"
903:
89:
670:
431:
384:
364:
971:
888:
Last edited at 13:54, 2 August 2007 (UTC). Substituted at 07:29, 30 April 2016 (UTC)
861:
709:
want to accomplish and which address the concerns expressed in the original question.
762:
754:
961:
937:
911:
880:
844:
828:
794:
770:
697:
355:
292:
289:
281:
207:
73:
52:
730:
184:
79:
102:
15:
835:
some application/examples involving non-autonomous systems?
856:
673:
590:
531:
493:
454:
434:
407:
387:
367:
346:
it hinders my understanding of those lines. regards.
314:
206:, a collaborative effort to improve the coverage of
101:, a collaborative effort to improve the coverage of
679:
659:
576:
517:
479:
440:
420:
393:
373:
338:
860:, and are posted here for posterity. Following
854:The comment(s) below were originally left at
8:
892:problem in "A second-order example" section?
731:http://geocalc.clas.asu.edu/pdf/Hamilton.pdf
660:{\displaystyle z=(x,y,z,p_{x},p_{y},p_{z})}
988:C-Class physics articles of Mid-importance
947:
923:
152:
47:
672:
648:
635:
622:
589:
565:
552:
539:
530:
492:
468:
453:
433:
412:
406:
386:
366:
313:
154:
49:
19:
761:) 20:31, 24 March 2013 (UTC) (revised
7:
200:This article is within the scope of
95:This article is within the scope of
857:Talk:Symplectic integrator/Comments
577:{\displaystyle (p_{x},p_{y},p_{z})}
38:It is of interest to the following
14:
998:Mid-priority mathematics articles
862:several discussions in past years
805:A Request for more practical info
220:Knowledge:WikiProject Mathematics
448:just means a pair of variables,
223:Template:WikiProject Mathematics
187:
177:
156:
82:
72:
51:
20:
983:Mid-importance physics articles
240:This article has been rated as
135:This article has been rated as
912:11:29, 17 September 2021 (UTC)
654:
597:
571:
532:
512:
494:
474:
455:
356:07:16, 24 September 2008 (UTC)
333:
321:
1:
962:11:55, 15 November 2023 (UTC)
938:11:41, 15 November 2023 (UTC)
845:19:22, 2 September 2019 (UTC)
214:and see a list of open tasks.
115:Knowledge:WikiProject Physics
109:and see a list of open tasks.
993:C-Class mathematics articles
293:22:16, 27 October 2007 (UTC)
282:12:48, 25 October 2007 (UTC)
118:Template:WikiProject Physics
829:18:43, 31 August 2011 (UTC)
1014:
881:13:54, 2 August 2007 (UTC)
771:21:01, 26 March 2013 (UTC)
698:05:54, 1 August 2010 (UTC)
141:project's importance scale
869:
480:{\displaystyle (x,p_{x})}
239:
172:
134:
67:
46:
978:C-Class physics articles
795:03:10, 20 May 2021 (UTC)
270:It says in the article:
246:project's priority scale
518:{\displaystyle (x,y,z)}
401:and if its momentum is
339:{\displaystyle z=(p,q)}
203:WikiProject Mathematics
681:
661:
578:
519:
481:
442:
422:
395:
375:
340:
28:This article is rated
733:, accessed 3/25/2013.
682:
662:
579:
520:
482:
443:
423:
421:{\displaystyle p_{x}}
396:
376:
341:
671:
588:
529:
525:and its momentum by
491:
452:
432:
405:
385:
365:
312:
226:mathematics articles
306:Runge-Kutta_methods
98:WikiProject Physics
850:Assessment comment
677:
657:
574:
515:
477:
438:
418:
391:
371:
336:
195:Mathematics portal
34:content assessment
964:
952:comment added by
940:
928:comment added by
886:
885:
819:comment added by
680:{\displaystyle z}
441:{\displaystyle z}
394:{\displaystyle x}
374:{\displaystyle z}
299:An other question
260:
259:
256:
255:
252:
251:
151:
150:
147:
146:
1005:
867:
866:
859:
831:
701:
686:
684:
683:
678:
666:
664:
663:
658:
653:
652:
640:
639:
627:
626:
583:
581:
580:
575:
570:
569:
557:
556:
544:
543:
524:
522:
521:
516:
486:
484:
483:
478:
473:
472:
447:
445:
444:
439:
427:
425:
424:
419:
417:
416:
400:
398:
397:
392:
380:
378:
377:
372:
345:
343:
342:
337:
228:
227:
224:
221:
218:
197:
192:
191:
181:
174:
173:
168:
160:
153:
123:
122:
121:physics articles
119:
116:
113:
92:
87:
86:
76:
69:
68:
63:
55:
48:
31:
25:
24:
16:
1013:
1012:
1008:
1007:
1006:
1004:
1003:
1002:
968:
967:
894:
855:
852:
814:
807:
691:
669:
668:
644:
631:
618:
586:
585:
561:
548:
535:
527:
526:
489:
488:
464:
450:
449:
430:
429:
408:
403:
402:
383:
382:
363:
362:
310:
309:
301:
268:
225:
222:
219:
216:
215:
193:
186:
166:
120:
117:
114:
111:
110:
88:
81:
61:
32:on Knowledge's
29:
12:
11:
5:
1011:
1009:
1001:
1000:
995:
990:
985:
980:
970:
969:
966:
965:
942:
941:
915:
914:
893:
890:
884:
883:
851:
848:
821:76.104.205.238
806:
803:
802:
801:
800:
799:
798:
797:
787:74.140.199.156
777:
776:
775:
774:
748:
747:
746:
745:
737:
736:
735:
734:
723:
722:
721:
720:
713:
712:
711:
710:
703:
702:
696:comment added
676:
656:
651:
647:
643:
638:
634:
630:
625:
621:
617:
614:
611:
608:
605:
602:
599:
596:
593:
573:
568:
564:
560:
555:
551:
547:
542:
538:
534:
514:
511:
508:
505:
502:
499:
496:
476:
471:
467:
463:
460:
457:
437:
415:
411:
390:
370:
348:192.93.101.133
335:
332:
329:
326:
323:
320:
317:
300:
297:
296:
295:
267:
264:
262:
258:
257:
254:
253:
250:
249:
238:
232:
231:
229:
212:the discussion
199:
198:
182:
170:
169:
161:
149:
148:
145:
144:
137:Mid-importance
133:
127:
126:
124:
107:the discussion
94:
93:
90:Physics portal
77:
65:
64:
62:Mid‑importance
56:
44:
43:
37:
26:
13:
10:
9:
6:
4:
3:
2:
1010:
999:
996:
994:
991:
989:
986:
984:
981:
979:
976:
975:
973:
963:
959:
955:
951:
944:
943:
939:
935:
931:
927:
921:
917:
916:
913:
909:
905:
900:
899:
898:
891:
889:
882:
878:
874:
868:
865:
863:
858:
849:
847:
846:
842:
838:
832:
830:
826:
822:
818:
810:
804:
796:
792:
788:
783:
782:
781:
780:
779:
778:
772:
768:
764:
760:
756:
752:
751:
750:
749:
741:
740:
739:
738:
732:
727:
726:
725:
724:
717:
716:
715:
714:
707:
706:
705:
704:
699:
695:
690:
674:
667:. The symbol
649:
645:
641:
636:
632:
628:
623:
619:
615:
612:
609:
606:
603:
600:
594:
591:
566:
562:
558:
553:
549:
545:
540:
536:
509:
506:
503:
500:
497:
469:
465:
461:
458:
435:
413:
409:
388:
368:
361:In the text,
360:
359:
358:
357:
353:
349:
330:
327:
324:
318:
315:
307:
298:
294:
291:
286:
285:
284:
283:
280:
274:
271:
265:
263:
247:
243:
237:
234:
233:
230:
213:
209:
205:
204:
196:
190:
185:
183:
180:
176:
175:
171:
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:
954:92.225.65.73
948:— Preceding
930:92.225.65.73
924:— Preceding
895:
887:
873:Jitse Niesen
871:methods. --
853:
833:
815:— Preceding
811:
808:
302:
275:
272:
269:
261:
242:Mid-priority
241:
201:
167:Mid‑priority
136:
96:
40:WikiProjects
692:—Preceding
217:Mathematics
208:mathematics
164:Mathematics
972:Categories
266:A question
950:unsigned
926:unsigned
817:unsigned
763:scanyon
755:scanyon
694:undated
584:, then
428:, then
244:on the
139:on the
112:Physics
103:Physics
59:Physics
30:C-class
904:Daniel
837:Egeris
290:Jcandy
36:scale.
689:Norio
958:talk
934:talk
908:talk
877:talk
841:talk
825:talk
791:talk
767:talk
759:talk
352:talk
279:Ossi
236:Mid
131:Mid
974::
960:)
936:)
910:)
879:)
843:)
827:)
793:)
769:)
354:)
956:(
932:(
906:(
875:(
839:(
823:(
789:(
773:)
765:(
757:(
700:.
675:z
655:)
650:z
646:p
642:,
637:y
633:p
629:,
624:x
620:p
616:,
613:z
610:,
607:y
604:,
601:x
598:(
595:=
592:z
572:)
567:z
563:p
559:,
554:y
550:p
546:,
541:x
537:p
533:(
513:)
510:z
507:,
504:y
501:,
498:x
495:(
475:)
470:x
466:p
462:,
459:x
456:(
436:z
414:x
410:p
389:x
369:z
350:(
334:)
331:q
328:,
325:p
322:(
319:=
316:z
248:.
143:.
42::
Text is available under the Creative Commons Attribution-ShareAlike License. Additional terms may apply.
