-
Notifications
You must be signed in to change notification settings - Fork 0
/
particle.py
46 lines (36 loc) · 934 Bytes
/
particle.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
#!/usr/bin/env python
# -*- coding: utf-8 -*-
from math import log
import random
class Particle:
def __init__(self, element , r):
self.elem = element
self.x = r[0]
self.y = r[1]
self.z = r[2]
self.setForce([0,0,0])
def r(self):
return [self.x, self.y, self.z]
def p(self):
return [self.px, self.py, self.pz]
def F(self):
return [self.Fx, self.Fy, self.Fz]
def setPosition(self, r):
self.x = r[0]
self.y = r[1]
self.z = r[2]
def setMomentum(self, momentum = [0,0,0], maxwell = False):
if maxwell:
self.px= log(random.uniform(0,1))
self.py= log(random.uniform(0,1))
self.pz= log(random.uniform(0,1))
else :
self.px= momentum[0]
self.py= momentum[1]
self.pz= momentum[2]
def setForce(self, force):
self.Fx = force[0]
self.Fy = force[1]
self.Fz = force[2]
def __str__ (self):
return self.elem + " " + str(self.x) + " " + str(self.y) + " " + str(self.z) + "\n"