MultiPy/tests.py at master · songbird175/MultiPy · GitHub

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
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
"""
A script solely for testing
"""

import multipy as mp
import numpy as np
from sympy import *

r_ = mp.Vector(5,6,7)
rnaught = mp.Vector(1,2,3)
v_ = mp.Vector(2,0,4)
w_ = mp.Vector(10,12,14)
a_ = mp.Vector(8,9)
A = mp.Point(1, 2, 3)
B = mp.Point(4, 0, 2)
AB = mp.Vector(3, -2, -1)
paramAB = mp.Parametric(A, AB)
pnaught = mp.Point(1, 3, -2)
n = mp.Vector(2, 1, -1)
planeP = mp.Plane(pnaught, n)
r_of_t = mp.VVF(cos(Symbol('t')), sin(Symbol('t')), Symbol('t'))


# print(r_of_t.secant_vector().x)

# rt2 = mp.VVF((Symbol('t')**2), Symbol('t'), Symbol('t')**3)
# print(rt2.derivative())
# print(r_of_t.derivative())

# print(r_of_t.continuous(2*np.pi))

# example1 = mp.VVF((3*cos(2*Symbol('t'))), (4*sin(2*Symbol('t'))), Symbol('t'))
# example1.vis()

# r_of_t.vis(np.arange(-5,5,0.2))
# print(r_of_t.x)

# print(paramAB.r().x)

# print(planeP.gen_eq)
# print(planeP.x_intercept())
# print(planeP.y_intercept())
# print(planeP.z_intercept())

# print(planeP.equation)

# print(paramAB.eq_x)

# print(paramAB.xy_intercept())
# print(paramAB.yz_intercept())
# print(paramAB.xz_intercept())

# print(paramAB.equations())
# print(paramAB.symmetric_x, paramAB.symmetric_y, paramAB.symmetric_z)

# print(r_.ihat, r_.jhat, r_.khat)

# print(r_.magnitude())

# print(rnaught.direction())

# print(mp.vector_add(r_, rnaught))
# print(mp.vector_sub(rnaught, v_))

# print(mp.scalar_mult(rnaught, 0))
# print(r_.magnitude() + v_.magnitude())

# print(mp.dot_prod(rnaught, v_))

# print(mp.is_mult(r_, w_))

# a_ = mp.find_v(5, 30)
# print(a_.khat)

# print(mp.theta_between(r_, rnaught))

# print(mp.projection(r_, rnaught))
# print(mp.projection_theta(r_, 30))

# print(mp.cross_prod(r_, a_))
# print(a_.z)