Моделирование квантовой пены

[MallSerg]

Монте карло ваще ни к селу ни к городу. )

[MakarovDs]

Решил зайти с тыла нашел 2d лавовую лампу у одно чела, и решил её переписать в 3d но получилась лампа с мушками:

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

class Ball:
    def __init__(self, x, y, z, vx, vy, vz, radius):
        self.x = x
        self.y = y
        self.z = z
        self.vx = vx
        self.vy = vy
        self.vz = vz
        self.radius = radius

    def update(self):
        self.x += self.vx
        self.y += self.vy
        self.z += self.vz

        if self.x < 0 or self.x > 100:
            self.vx *= -1
        if self.y < 0 or self.y > 100:
            self.vy *= -1
        if self.z < 0 or self.z > 100:
            self.vz *= -1

    def get_position(self):
        return self.x, self.y, self.z

class LavaLamp:
    def __init__(self, num_balls):
        self.num_balls = num_balls
        self.balls = []

        for i in range(num_balls):
            x = np.random.uniform(0, 100)
            y = np.random.uniform(0, 100)
            z = np.random.uniform(0, 100)
            vx = np.random.uniform(-1, 1)
            vy = np.random.uniform(-1, 1)
            vz = np.random.uniform(-1, 1)
            radius = np.random.uniform(1, 5)

            self.balls.append(Ball(x, y, z, vx, vy, vz, radius))

    def update(self):
        for ball in self.balls:
            ball.update()

    def get_positions(self):
        positions = []

        for ball in self.balls:
            positions.append(ball.get_position())

        return positions

def main():
    num_balls = 100
    lamp = LavaLamp(num_balls)

    fig = plt.figure()
    ax = fig.add_subplot(111, projection='3d')

    for i in range(100):
        lamp.update()
        positions = lamp.get_positions()

        ax.clear()
        ax.set_xlim(0, 100)
        ax.set_ylim(0, 100)
        ax.set_zlim(0, 100)

        for position in positions:
            ax.scatter(position[0], position[1], position[2], c='r', marker='o', s=10)

        plt.pause(0.01)

    plt.show()

if __name__ == "__main__":
    main()

[MakarovDs]

Давайте исправим этот недостаток:

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

class Ball:
    def __init__(self, x, y, z, vx, vy, vz, radius):
        self.x = x
        self.y = y
        self.z = z
        self.vx = vx
        self.vy = vy
        self.vz = vz
        self.radius = radius

    def update(self):
        self.x += self.vx
        self.y += self.vy
        self.z += self.vz

        if self.x < 0 or self.x > 100:
            self.vx *= -1
        if self.y < 0 or self.y > 100:
            self.vy *= -1
        if self.z < 0 or self.z > 100:
            self.vz *= -1

    def get_position(self):
        return self.x, self.y, self.z

class LavaLamp:
    def __init__(self, num_balls):
        self.num_balls = num_balls
        self.balls = []
        self.size_factors = [1, 2, 3, 4, 5]
        self.probabilities = [0.4, 0.3, 0.2, 0.08, 0.02]  # Корректируйте эти вероятности по мере необходимости

        for i in range(num_balls):
            x = np.random.uniform(0, 100)
            y = np.random.uniform(0, 100)
            z = np.random.uniform(0, 100)
            vx = np.random.uniform(-1, 1)
            vy = np.random.uniform(-1, 1)
            vz = np.random.uniform(-1, 1)

            size_factor = np.random.choice(self.size_factors, p=self.probabilities)
            radius = size_factor

            self.balls.append(Ball(x, y, z, vx, vy, vz, radius))

    def update(self):
        for ball in self.balls:
            ball.update()

    def get_positions(self):
        positions = []

        for ball in self.balls:
            positions.append(ball.get_position())

        return positions

def main():
    num_balls = 100
    lamp = LavaLamp(num_balls)

    fig = plt.figure()
    ax = fig.add_subplot(111, projection='3d')

    for i in range(100):
        lamp.update()
        positions = lamp.get_positions()

        ax.clear()
        ax.set_xlim(0, 100)
        ax.set_ylim(0, 100)
        ax.set_zlim(0, 100)

        for position in positions:
            ax.scatter(position[0], position[1], position[2], c='r', marker='o', s=lamp.balls[positions.index(position)].radius ** 2 * 10)

        plt.pause(0.01)

    plt.show()

if __name__ == "__main__":
    main()

Теперь это больше похоже на лавовую лампу

[MakarovDs]

Короче решил сделать бесконечную процедурно генерируемую лавовую лампу на Babylon.js вот что у меня получилось:

ССЫЛКА

Как вам моя лавовая лампа? Перемещение с помощью стрелок.

Вот вторая версия

[MakarovDs]

Я В БЕСКОНЕЧНОЙ ЛАВОВОЙ ЛАМПЕ UwU



Это уже третья версия лавовой лампы!



Как вам!?