←back to thread

60 points melector | 1 comments | | HN request time: 0s | source
Show context
Evidlo ◴[] No.45512953[source]
Anyone else feel like the pendulum motion seems off? Maybe the default mass settings are weird, but the movement just does not look physical to me.
replies(5): >>45513075 #>>45513084 #>>45513483 #>>45514149 #>>45538570 #
magicalhippo ◴[] No.45514149[source]
It's way off. My first guess was that there was something wrong with the physics code, but after carefully checking against this[1] derivation of the Hamiltonian it seemed fine, and once I wrapped my head around the JavaScript the RK4 integration[2] checked out as well.

So, what else might be wrong I wondered. Well, it seems to move in the wrong direction... so I checked how the pendulum is displayed. And sure enough, I think there's a sign error:

  getUpperBob() {
     const { x0, y0, ang0, l0 } = this;
     const { x, y } = this.calculateBobPosition(x0, y0, ang0, l0);
     return { x, y };
   }

   getLowerBob() {
     const upperBobPos = this.getUpperBob();
     const { ang1, l1 } = this;
     const { x, y } = this.calculateBobPosition(
    upperBobPos.x,
    upperBobPos.y,
    -ang1,
    l1
     );
     return { x, y };
   }
Note how the upper bob uses ang0 while the lower one has -ang1. Meanwhile the physics derivation assumes both angles are against the vertical, so have same sign.

Changing -ang1 to ang1 does indeed make the pendulum move in a natural way, except now dragging it is flipped. Ie you drag it left and it moves right. Another sign error in setLowerBobPos. Fixing that as well it now works as I'd expect.

[1]: https://dassencio.org/46

[2]: https://lpsa.swarthmore.edu/NumInt/NumIntFourth.html##sectio...

replies(1): >>45518965 #
Evidlo ◴[] No.45518965[source]
Just made a PR
replies(2): >>45524146 #>>45537196 #
cabirum ◴[] No.45537196[source]
Hn shows this story is 3 hours old, yet your GitHub PR was made 2 days ago, linking here?
replies(1): >>45538147 #
1. pansa2 ◴[] No.45538147[source]
The story has been reposted to HN via the “second-chance pool”, which resets the dates.

https://news.ycombinator.com/item?id=26998308