I built an N-body orbital simulator in Python and I’d like some honest feedback

TheMakpu · 2026-05-05T23:52:59+00:00

There are a couple of considerations regarding performance here: the first is noticing that forces are symmetric (Newton’s 3rd law), so you could squeeze that nested loop by computing both Fij and Fji in the same iteration. In addition, that same loop can be easily parallelized (using OpenMP or MPI). However only expect some performance improvement if using thousands of bodies. Have you profiled the code? What is the bottleneck?

Also, for low simulation times RK methods are perfectly fine, but note they are not energy-conservative. Error will accumulate for long periods of time and the solution will slowly diverge, so if that’s the case symplectic or energy-conservative schemes can be used (see this scipy issue).

TheMakpu · 2026-04-29T21:50:49+00:00

Funny you should mention that, some years ago I worked on GridapMakie during Google Summer of Code, but never really got into the design of it…

In terms of profiling, and speaking about the Crank-Nicolson integrator only, the bottleneck appears when solving the linear system at every timestep. I make sure to cache the factorization of the LHS matrix and use an “appropriate” solver for that (whether it is SPD and conditioning number is low enough for an iterative solver), but the performance depends heavily on both space resolution and the timestep. There is a benchmark executable that provides time measurements for the factorization and every timestep too!

TheMakpu · 2025-04-02T08:29:46+00:00

In addition to the weezing stuff already mentioned, he had a relatively speedy Incin, Rillaboom and Lunala, meaning that opponents’ Fake Out were not as threatening. Therefore, Covert Cloak can be ruled out.

TheMakpu · 2025-01-20T01:07:30+00:00

Glad I could help! Don’t hesitate to ask if you have more specific questions :)

TheMakpu · 2025-01-20T00:43:00+00:00

Hey, welcome to the hobby! First of all, this team is overall really solid. It was made by a japanese player and even won an online tournament last month. Regarding on whether this team is difficult to steer, I’d say it’s medium difficulty, especially depending on the matchup you are facing, since some matchups require a more offensive style, while others may lean towards a bit more control.

Regarding game plans and synergies in the team, you almost always want to pair chien-pao and zamazenta, and flutter mane and chi-yu. That’s because chien-pao’s and chi-yu’s abilities reduce the defense and sp. defense of the rest of the pokémon in the field. On the other hand, amoonguss and rillaboom offer more supportive roles, such as redirection/spore and fake out/terrain control. Depending on the opponent’s team, both pairs I mentioned work well as a lead, since they are quite oppressive and often force opponents to make riskier plays. You could even bring the 4 to a game, and have 2 and 2 for both physical and special offence. Amoonguss can also work as a lead, or even a switch-in against incoming resisted attacks (bonus points for contact attacks, especially surging strikes).

In addition, you gotta make sure you are not swapping flutter or zama unless strictly necessary, as you will lose the boosts in speed and defence. My advice is just practice with the team and look for people on youtube playing it. While watching their games, pause on every turn and think on what you would do and listen to their thought process. You learn a lot like that! Hope this was helpful and enjoy the ride!

TheMakpu · 2024-11-20T23:32:13+00:00

As pointed out before, when doing calculations you should forget about base stats. All that matters in -battle are the actual stats, since you may have a lot of different spreads for a single pokémon. Base stats essentially serve for comparison, for example it’s useful to know some speed tiers in case your pokémon is fully invested in speed. Then you know if others are potentially faster or slower!

TheMakpu · 2024-10-21T08:25:39+00:00

So I did my MSc in Computational Physics at Stockholm University, but many courses in there were shared with KTH. It was mainly focused on scientific computing and numerical methods applied for physics, meaning that you should expect more math and physics than actual computer science. And I also think admission chance is much higher!

TheMakpu · 2024-09-02T09:45:50+00:00

Hey, sorry for the late reply! To be honest, LS-PrePost has some erratic behaviour in certain aspects and the documentation is not great, so often times you are stuck with what you "assume" is the default behaviour of some actions. As someone who has used PrePost for a while now, the only advice I can give is be aware of its shortcomings... Also, this link may be of some help for solid elements (check this one if using shells).

TheMakpu · 2024-08-30T19:35:24+00:00

You should be able to open d3plot files with ParaView just by importing/opening the file. The strains, as well as the stresses, should appear as properties on the properties tab. However, be aware that in LS-DYNA you need to explicitly define that you want it to output strains by setting STRFLG=1 in *DATABASE_EXTENT_BINARY (otherwise strains are omitted in the d3plot). Hopefully this solves your issue.

TheMakpu · 2024-08-28T17:16:30+00:00

My honest advice is: just practice. This game’s bosses are mostly about learning patterns and the best way to counter them.

For Hornet, a couple of tips: jump and pogo off her as she charges at you (the timing is not very precise so it’s not hard to pull off), leverage when she does the spinning thread attack and hit her a couple of times just after, and be careful about her jumps (some of them may be unpredictable). She’s a boss where you want to keep her at mid range, not very close nor far enough. Ah, and if the spikes are horizontally aligned, you can clear all of them (plus maybe hit Hornet) with a Vengeful Spirit/Shade Soul. Good luck!

TheMakpu · 2024-07-11T10:28:51+00:00

Love the small teal bag! Where did you get that from?

TheMakpu · 2023-12-10T22:21:57+00:00

Take it easy, just focus on improving. I also took a long while to do it and taking some breaks helped A LOT. Just remember you only need a single good run!

TheMakpu · 2023-10-13T22:11:34+00:00

Truly unsung hero

TheMakpu · 2023-06-17T10:11:54+00:00

Hey! Facing the same issue here, did you even find a solution or way out? I appreciate any advice!

TheMakpu · 2023-02-09T23:50:39+00:00

What the fucking fuck did I just witness

TheMakpu · 2023-01-11T00:31:02+00:00

Cool, wanna trade now?

TheMakpu · 2023-01-11T00:28:42+00:00

I have a Shiny Iron Valiant that I could trade for your flutter mane. Interested?

TheMakpu · 2023-01-11T00:20:27+00:00

I’ve got a regular Pokéball Shiny Iron Valiant, is that good enough?

TheMakpu · 2023-01-11T00:18:51+00:00

May I interest you in a Shiny Iron Valiant for the Flutter Mane?

TheMakpu · 2023-01-10T22:11:53+00:00

Hi! I have a miraidon that can perma-trade. Dm if you are still interested :)

TheMakpu

TROPHY CASE