Very nice work. Just spotted this though: the length() function will use a sqrt() which is time consuming in itself. You have accidentally if LastPos.... inside the for loop. Move this and you gain a lot of extra cpu time. Even a much longer tick list.

really interesting stuff! I've been toying with similar stuff for a long time just for fun. it's really interesting to see how you're approaching it!

Looking forward to the delicious links. :P

im sure python will handle 100,000 particles just fine 😅

I just want to say before even starting, I literally have spent an embarrassing amount of time playing Dwarf Fortress, Rogue, Nethack, CDDA, NeoScav, Rimworld, Warsim:KoA So, so many of these games. I genuinely can't believe that me trying to learn how to do this myself completely raw no experience leads my first video to being this, and you open with that. I can't explain the coincidence here, I just rambled with my fingers.

I'm not gonna lie. I would love to see the code of this.

Bro the fact my reason to do this too was because of dwarf fortress as well 💀

Very cool!

My guy, you are a neat dude. You know how to do everything I want to learn. This is looking like some Wargames type of effect.

This is great! I really enjoyed this video and your procedural generation video. You have a very direct way of speaking and I think that is crucial in coding tutorials. I really hope you do more videos. I'm still learning python and you help out alot. I would like to see a bit more on noise and how to manipulate the maps to achieve a desired look (mostly lakes similar to the game Intergalactic Fishing). I also want to find a way to make that star map from Starbound. Its basically a procedural map of stars that you could hover over a star and see the random planets and moons inside. Anyways, keep it up! you have something here.

You can just make your own smooth noise function and use that as a base for an fbm loop. You could have also explored various cellular automata solutions, or the newfangled wave function collapse. Still, it works.

I like it, keep it up!

Your channel is dope

How do you expand the radius of orbit in the window ?

Beautifully done! I have been working out a similar Python problem. I will share mine - once I have completed a solution.

Excellent vid

This is very nice but the result is not realistic at all. Particles behave very different than they should. I think n-body is usefull only when one particle is is much larger than all the others.

wow, helped me a lot, I knew how to get the velocity, but I couldn't figure out how to get the direction :) thanks

If you import sqrt by writing from math import sqrt, it is significantly faster. Don’t ask me why but it is.

You should record the video in 720p/1080p

Thank you very much for the detailed explanation. I have compared your code with other options. yours is much shorter. the fewer lines the better. I tried to write down the same code, exactly. can you send me your code via mail or share github? Because my code doesn't run. I'll attach my code if anything from tkinter import * import numpy as np import random import math fig = Tk() fig.wm_title('N body problem') canvas = Canvas(fig, width =700, height=700, bg='black') canvas.grid(row=0, column=0) G= 1 bodies= [] class body: def __init__(self, Mass, Px, Py, V): self.Mass = Mass self.Px = Px self.Py =Py self.P = np.array([self.Px, self.Py], dtype = float) self.V = np.array(V, dtype = float) self.Logx = [] self.Logy = [] self.R = (self.Mass / 3.14)**(1/3) def Fg(self, other): if self == other: self.dV = np.array([0,0], dtype= float) else: self.dist =math.sqrt((self.P[0] - other.P[0])**2 + (self.P[1] - other.P[1])**2) self.dV = np.array([0,0], dtype= float) if self.dist == 0: self.dV = 0 else: self.Fg = [-G * self.Mass * other.Mass / self.dist**2] * (self.P - other.P) / self.dist self.dv = self.Fg / self.Mass if self.dist < self.R + other.R +5: self.dV *= -0.5 self.v +=self.dv if self.P[0] < 10 + self.R: self.V[0] *= -0.9 self.P[0] +=1 if self.P[0] > 698 - self.R: self.V[0] *= -0.9 self.P[0] -=1 if self.P[1] < 10 + self.R: self.V[1] *= -0.9 self.P[1] +=1 if self.P[1] > 698 - self.R: self.V[1] *= -0.9 self.P[1] -=1 BODIES = 25 for i in range (0, BODIES): BODIES.append(Body(random.randrange(10, 500), random.randrange(50, 650), random.randrange(50, 650), [(random.randrange(-10,10)/20), random.randrange(-10,10)/20])) T = 0 Tx = True while Tx ==True: T+=1 canvas.delete ('all') for Body1 in Bodies: for Body2 in Bodies: Body1.Fg[Body2] Body1.P +=Body1.V if Body1.P[1] < 5: Body1.V[1] *= -0.5 if Body1.P[1] > 695: Body1.V[1] *=-0.5 if Body1.P[0] < 5: Body1.V[0] *= -0.9 if Body1.P[0] > 695: Body1.V[0] *=-0.9 canvas.create_oval(Body1.P[0] -Body1.R, Body1.P[1] - Body1.R, Body1.P[0] + Body1.R, Body1.P[1] + Body1.R, fill = 'yellow') canvas.update() if T ==1000000: Tx = False mainloop()

Body() takes no arguments

This awesome and very powerful, I appreciate it all.

hello! how do you create another circle around the earth orbital? I try scrolling and nothing is working

This is pretty excellent

Did you ever end up doing the 4D sim? Sounds interesting

Awesome!

I don't feel like this is properly addressing the N-body problem, surely you would find the center of gravity/mass of the field of particles then for each "particle" you would find it's acceleration towards that singular point. Your process is iterating through each particle individually which is kind of wrong

good video! but why not 720p or above?

good vid mate, but can i ask something. how can we make the result reproducible?

w...what is that abomination? i dont get the reference

how to visualize using this simulation code?

hey man can you show this with 3d plotting maybe opencl or something thanks :D

this is technically my whole fyp lol

Very nice mate, liked how you didnt type the code as you explained but rather just explained it. Very informative. Look forward 2 more stuff