←back to thread

60 points melector | 1 comments | | HN request time: 0.2s | 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 #
1. magicalhippo ◴[] No.45524146[source]
I considered it, but had to go to work and I figured since the last commit was over 4 years ago it might just sit there idly, so I skipped it.

But no, PR was merged in short order. Lesson learned.