Skip to content

Commit a2db490

Browse files
JWock82JWock82
authored
Merge pull request #205 from connorferster/fixes/numpy-itemset-deprecation
fix: remove .itemset from codebase
2 parents 33d3fad + 290c217 commit a2db490

File tree

2 files changed

+57
-57
lines changed

2 files changed

+57
-57
lines changed

‎PyNite/FixedEndReactions.py

Lines changed: 30 additions & 30 deletions
Original file line numberDiff line numberDiff line change
@@ -33,15 +33,15 @@ def FER_PtLoad(P, x, L, Direction):
3333

3434
# Populate the fixed end reaction vector
3535
if Direction == "Fy":
36-
FER.itemset((1, 0), -P*b**2*(L+2*x)/L**3)
37-
FER.itemset((5, 0), -P*x*b**2/L**2)
38-
FER.itemset((7, 0), -P*x**2*(L+2*b)/L**3)
39-
FER.itemset((11, 0), P*x**2*b/L**2)
36+
FER[1, 0] = -P*b**2*(L+2*x)/L**3
37+
FER[5, 0] = -P*x*b**2/L**2
38+
FER[7, 0] = -P*x**2*(L+2*b)/L**3
39+
FER[11, 0] = P*x**2*b/L**2
4040
elif Direction == "Fz":
41-
FER.itemset((2, 0), -P*b**2*(L+2*x)/L**3)
42-
FER.itemset((4, 0), P*x*b**2/L**2)
43-
FER.itemset((8, 0), -P*x**2*(L+2*b)/L**3)
44-
FER.itemset((10, 0), -P*x**2*b/L**2)
41+
FER[2, 0] = -P*b**2*(L+2*x)/L**3
42+
FER[4, 0] = P*x*b**2/L**2
43+
FER[8, 0] = -P*x**2*(L+2*b)/L**3
44+
FER[10, 0] = -P*x**2*b/L**2
4545

4646
return FER
4747

@@ -70,15 +70,15 @@ def FER_Moment(M, x, L, Direction):
7070

7171
# Populate the fixed end reaction vector
7272
if Direction == "Mz":
73-
FER.itemset((1, 0), 6*M*x*b/L**3)
74-
FER.itemset((5, 0), M*b*(2*x-b)/L**2)
75-
FER.itemset((7, 0), -6*M*x*b/L**3)
76-
FER.itemset((11, 0), M*x*(2*b-x)/L**2)
73+
FER[1, 0] = 6*M*x*b/L**3
74+
FER[5, 0] = M*b*(2*x-b)/L**2
75+
FER[7, 0] = -6*M*x*b/L**3
76+
FER[11, 0] = M*x*(2*b-x)/L**2
7777
elif Direction == "My":
78-
FER.itemset((2, 0), -6*M*x*b/L**3)
79-
FER.itemset((4, 0), M*b*(2*x-b)/L**2)
80-
FER.itemset((8, 0), 6*M*x*b/L**3)
81-
FER.itemset((10, 0), M*x*(2*b-x)/L**2)
78+
FER[2, 0] = -6*M*x*b/L**3
79+
FER[4, 0] = M*b*(2*x-b)/L**2
80+
FER[8, 0] = 6*M*x*b/L**3
81+
FER[10, 0] = M*x*(2*b-x)/L**2
8282
return FER
8383

8484
# %%
@@ -90,15 +90,15 @@ def FER_LinLoad(w1, w2, x1, x2, L, Direction):
9090

9191
# Populate the fixed end reaction vector
9292
if Direction == 'Fy':
93-
FER.itemset((1, 0), (x1 - x2)*(10*L**3*w1 + 10*L**3*w2 - 15*L*w1*x1**2 - 10*L*w1*x1*x2 - 5*L*w1*x2**2 - 5*L*w2*x1**2 - 10*L*w2*x1*x2 - 15*L*w2*x2**2 + 8*w1*x1**3 + 6*w1*x1**2*x2 + 4*w1*x1*x2**2 + 2*w1*x2**3 + 2*w2*x1**3 + 4*w2*x1**2*x2 + 6*w2*x1*x2**2 + 8*w2*x2**3)/(20*L**3))
94-
FER.itemset((5, 0), (x1 - x2)*(20*L**2*w1*x1 + 10*L**2*w1*x2 + 10*L**2*w2*x1 + 20*L**2*w2*x2 - 30*L*w1*x1**2 - 20*L*w1*x1*x2 - 10*L*w1*x2**2 - 10*L*w2*x1**2 - 20*L*w2*x1*x2 - 30*L*w2*x2**2 + 12*w1*x1**3 + 9*w1*x1**2*x2 + 6*w1*x1*x2**2 + 3*w1*x2**3 + 3*w2*x1**3 + 6*w2*x1**2*x2 + 9*w2*x1*x2**2 + 12*w2*x2**3)/(60*L**2))
95-
FER.itemset((7, 0), -(x1 - x2)*(-15*L*w1*x1**2 - 10*L*w1*x1*x2 - 5*L*w1*x2**2 - 5*L*w2*x1**2 - 10*L*w2*x1*x2 - 15*L*w2*x2**2 + 8*w1*x1**3 + 6*w1*x1**2*x2 + 4*w1*x1*x2**2 + 2*w1*x2**3 + 2*w2*x1**3 + 4*w2*x1**2*x2 + 6*w2*x1*x2**2 + 8*w2*x2**3)/(20*L**3))
96-
FER.itemset((11, 0), (x1 - x2)*(-15*L*w1*x1**2 - 10*L*w1*x1*x2 - 5*L*w1*x2**2 - 5*L*w2*x1**2 - 10*L*w2*x1*x2 - 15*L*w2*x2**2 + 12*w1*x1**3 + 9*w1*x1**2*x2 + 6*w1*x1*x2**2 + 3*w1*x2**3 + 3*w2*x1**3 + 6*w2*x1**2*x2 + 9*w2*x1*x2**2 + 12*w2*x2**3)/(60*L**2))
93+
FER[1, 0] = (x1 - x2)*(10*L**3*w1 + 10*L**3*w2 - 15*L*w1*x1**2 - 10*L*w1*x1*x2 - 5*L*w1*x2**2 - 5*L*w2*x1**2 - 10*L*w2*x1*x2 - 15*L*w2*x2**2 + 8*w1*x1**3 + 6*w1*x1**2*x2 + 4*w1*x1*x2**2 + 2*w1*x2**3 + 2*w2*x1**3 + 4*w2*x1**2*x2 + 6*w2*x1*x2**2 + 8*w2*x2**3)/(20*L**3)
94+
FER[5, 0] = (x1 - x2)*(20*L**2*w1*x1 + 10*L**2*w1*x2 + 10*L**2*w2*x1 + 20*L**2*w2*x2 - 30*L*w1*x1**2 - 20*L*w1*x1*x2 - 10*L*w1*x2**2 - 10*L*w2*x1**2 - 20*L*w2*x1*x2 - 30*L*w2*x2**2 + 12*w1*x1**3 + 9*w1*x1**2*x2 + 6*w1*x1*x2**2 + 3*w1*x2**3 + 3*w2*x1**3 + 6*w2*x1**2*x2 + 9*w2*x1*x2**2 + 12*w2*x2**3)/(60*L**2)
95+
FER[7, 0] = -(x1 - x2)*(-15*L*w1*x1**2 - 10*L*w1*x1*x2 - 5*L*w1*x2**2 - 5*L*w2*x1**2 - 10*L*w2*x1*x2 - 15*L*w2*x2**2 + 8*w1*x1**3 + 6*w1*x1**2*x2 + 4*w1*x1*x2**2 + 2*w1*x2**3 + 2*w2*x1**3 + 4*w2*x1**2*x2 + 6*w2*x1*x2**2 + 8*w2*x2**3)/(20*L**3)
96+
FER[11, 0] = (x1 - x2)*(-15*L*w1*x1**2 - 10*L*w1*x1*x2 - 5*L*w1*x2**2 - 5*L*w2*x1**2 - 10*L*w2*x1*x2 - 15*L*w2*x2**2 + 12*w1*x1**3 + 9*w1*x1**2*x2 + 6*w1*x1*x2**2 + 3*w1*x2**3 + 3*w2*x1**3 + 6*w2*x1**2*x2 + 9*w2*x1*x2**2 + 12*w2*x2**3)/(60*L**2)
9797
elif Direction == 'Fz':
98-
FER.itemset((2, 0), (x1 - x2)*(10*L**3*w1 + 10*L**3*w2 - 15*L*w1*x1**2 - 10*L*w1*x1*x2 - 5*L*w1*x2**2 - 5*L*w2*x1**2 - 10*L*w2*x1*x2 - 15*L*w2*x2**2 + 8*w1*x1**3 + 6*w1*x1**2*x2 + 4*w1*x1*x2**2 + 2*w1*x2**3 + 2*w2*x1**3 + 4*w2*x1**2*x2 + 6*w2*x1*x2**2 + 8*w2*x2**3)/(20*L**3))
99-
FER.itemset((4, 0), -(x1 - x2)*(20*L**2*w1*x1 + 10*L**2*w1*x2 + 10*L**2*w2*x1 + 20*L**2*w2*x2 - 30*L*w1*x1**2 - 20*L*w1*x1*x2 - 10*L*w1*x2**2 - 10*L*w2*x1**2 - 20*L*w2*x1*x2 - 30*L*w2*x2**2 + 12*w1*x1**3 + 9*w1*x1**2*x2 + 6*w1*x1*x2**2 + 3*w1*x2**3 + 3*w2*x1**3 + 6*w2*x1**2*x2 + 9*w2*x1*x2**2 + 12*w2*x2**3)/(60*L**2))
100-
FER.itemset((8, 0), -(x1 - x2)*(-15*L*w1*x1**2 - 10*L*w1*x1*x2 - 5*L*w1*x2**2 - 5*L*w2*x1**2 - 10*L*w2*x1*x2 - 15*L*w2*x2**2 + 8*w1*x1**3 + 6*w1*x1**2*x2 + 4*w1*x1*x2**2 + 2*w1*x2**3 + 2*w2*x1**3 + 4*w2*x1**2*x2 + 6*w2*x1*x2**2 + 8*w2*x2**3)/(20*L**3))
101-
FER.itemset((10, 0), -(x1 - x2)*(-15*L*w1*x1**2 - 10*L*w1*x1*x2 - 5*L*w1*x2**2 - 5*L*w2*x1**2 - 10*L*w2*x1*x2 - 15*L*w2*x2**2 + 12*w1*x1**3 + 9*w1*x1**2*x2 + 6*w1*x1*x2**2 + 3*w1*x2**3 + 3*w2*x1**3 + 6*w2*x1**2*x2 + 9*w2*x1*x2**2 + 12*w2*x2**3)/(60*L**2))
98+
FER[2, 0] = (x1 - x2)*(10*L**3*w1 + 10*L**3*w2 - 15*L*w1*x1**2 - 10*L*w1*x1*x2 - 5*L*w1*x2**2 - 5*L*w2*x1**2 - 10*L*w2*x1*x2 - 15*L*w2*x2**2 + 8*w1*x1**3 + 6*w1*x1**2*x2 + 4*w1*x1*x2**2 + 2*w1*x2**3 + 2*w2*x1**3 + 4*w2*x1**2*x2 + 6*w2*x1*x2**2 + 8*w2*x2**3)/(20*L**3)
99+
FER[4, 0] = -(x1 - x2)*(20*L**2*w1*x1 + 10*L**2*w1*x2 + 10*L**2*w2*x1 + 20*L**2*w2*x2 - 30*L*w1*x1**2 - 20*L*w1*x1*x2 - 10*L*w1*x2**2 - 10*L*w2*x1**2 - 20*L*w2*x1*x2 - 30*L*w2*x2**2 + 12*w1*x1**3 + 9*w1*x1**2*x2 + 6*w1*x1*x2**2 + 3*w1*x2**3 + 3*w2*x1**3 + 6*w2*x1**2*x2 + 9*w2*x1*x2**2 + 12*w2*x2**3)/(60*L**2)
100+
FER[8, 0] = -(x1 - x2)*(-15*L*w1*x1**2 - 10*L*w1*x1*x2 - 5*L*w1*x2**2 - 5*L*w2*x1**2 - 10*L*w2*x1*x2 - 15*L*w2*x2**2 + 8*w1*x1**3 + 6*w1*x1**2*x2 + 4*w1*x1*x2**2 + 2*w1*x2**3 + 2*w2*x1**3 + 4*w2*x1**2*x2 + 6*w2*x1*x2**2 + 8*w2*x2**3)/(20*L**3)
101+
FER[10, 0] = -(x1 - x2)*(-15*L*w1*x1**2 - 10*L*w1*x1*x2 - 5*L*w1*x2**2 - 5*L*w2*x1**2 - 10*L*w2*x1*x2 - 15*L*w2*x2**2 + 12*w1*x1**3 + 9*w1*x1**2*x2 + 6*w1*x1*x2**2 + 3*w1*x2**3 + 3*w2*x1**3 + 6*w2*x1**2*x2 + 9*w2*x1*x2**2 + 12*w2*x2**3)/(60*L**2)
102102

103103
return FER
104104

@@ -110,8 +110,8 @@ def FER_AxialPtLoad(P, x, L):
110110
FER = zeros((12, 1))
111111

112112
# Populate the fixed end reaction vector
113-
FER.itemset((0, 0), -P*(L-x)/L)
114-
FER.itemset((6, 0), -P*x/L)
113+
FER[0, 0] = -P*(L-x)/L
114+
FER[6, 0] = -P*x/L
115115

116116
return FER
117117

@@ -123,8 +123,8 @@ def FER_AxialLinLoad(p1, p2, x1, x2, L):
123123
FER = zeros((12, 1))
124124

125125
# Populate the fixed end reaction vector
126-
FER.itemset((0, 0), 1/(6*L)*(x1-x2)*(3*L*p1+3*L*p2-2*p1*x1-p1*x2-p2*x1-2*p2*x2))
127-
FER.itemset((6, 0), 1/(6*L)*(x1-x2)*(2*p1*x1+p1*x2+p2*x1+2*p2*x2))
126+
FER[0, 0] = 1/(6*L)*(x1-x2)*(3*L*p1+3*L*p2-2*p1*x1-p1*x2-p2*x1-2*p2*x2)
127+
FER[6, 0] = 1/(6*L)*(x1-x2)*(2*p1*x1+p1*x2+p2*x1+2*p2*x2)
128128

129129
return FER
130130

@@ -144,7 +144,7 @@ def FER_Torque(T, x, L):
144144
FER = zeros((12, 1))
145145

146146
# Populate the fixed end reaction vector
147-
FER.itemset((3, 0), -T*(L - x)/L)
148-
FER.itemset((9, 0), -T*x/L)
147+
FER[3, 0] = -T*(L - x)/L
148+
FER[9, 0] = -T*x/L
149149

150150
return FER

‎PyNite/Plate3D.py

Lines changed: 27 additions & 27 deletions
Original file line numberDiff line numberDiff line change
@@ -404,33 +404,33 @@ def D(self, combo_name='Combo 1'):
404404
D = zeros((24, 1))
405405

406406
# Read in the global displacements from the nodes
407-
D.itemset((0, 0), self.i_node.DX[combo_name])
408-
D.itemset((1, 0), self.i_node.DY[combo_name])
409-
D.itemset((2, 0), self.i_node.DZ[combo_name])
410-
D.itemset((3, 0), self.i_node.RX[combo_name])
411-
D.itemset((4, 0), self.i_node.RY[combo_name])
412-
D.itemset((5, 0), self.i_node.RZ[combo_name])
413-
414-
D.itemset((6, 0), self.j_node.DX[combo_name])
415-
D.itemset((7, 0), self.j_node.DY[combo_name])
416-
D.itemset((8, 0), self.j_node.DZ[combo_name])
417-
D.itemset((9, 0), self.j_node.RX[combo_name])
418-
D.itemset((10, 0), self.j_node.RY[combo_name])
419-
D.itemset((11, 0), self.j_node.RZ[combo_name])
420-
421-
D.itemset((12, 0), self.m_node.DX[combo_name])
422-
D.itemset((13, 0), self.m_node.DY[combo_name])
423-
D.itemset((14, 0), self.m_node.DZ[combo_name])
424-
D.itemset((15, 0), self.m_node.RX[combo_name])
425-
D.itemset((16, 0), self.m_node.RY[combo_name])
426-
D.itemset((17, 0), self.m_node.RZ[combo_name])
427-
428-
D.itemset((18, 0), self.n_node.DX[combo_name])
429-
D.itemset((19, 0), self.n_node.DY[combo_name])
430-
D.itemset((20, 0), self.n_node.DZ[combo_name])
431-
D.itemset((21, 0), self.n_node.RX[combo_name])
432-
D.itemset((22, 0), self.n_node.RY[combo_name])
433-
D.itemset((23, 0), self.n_node.RZ[combo_name])
407+
D[0, 0] = self.i_node.DX[combo_name]
408+
D[1, 0] = self.i_node.DY[combo_name]
409+
D[2, 0] = self.i_node.DZ[combo_name]
410+
D[3, 0] = self.i_node.RX[combo_name]
411+
D[4, 0] = self.i_node.RY[combo_name]
412+
D[5, 0] = self.i_node.RZ[combo_name]
413+
414+
D[6, 0] = self.j_node.DX[combo_name]
415+
D[7, 0] = self.j_node.DY[combo_name]
416+
D[8, 0] = self.j_node.DZ[combo_name]
417+
D[9, 0] = self.j_node.RX[combo_name]
418+
D[10, 0] = self.j_node.RY[combo_name]
419+
D[11, 0] = self.j_node.RZ[combo_name]
420+
421+
D[12, 0] = self.m_node.DX[combo_name]
422+
D[13, 0] = self.m_node.DY[combo_name]
423+
D[14, 0] = self.m_node.DZ[combo_name]
424+
D[15, 0] = self.m_node.RX[combo_name]
425+
D[16, 0] = self.m_node.RY[combo_name]
426+
D[17, 0] = self.m_node.RZ[combo_name]
427+
428+
D[18, 0] = self.n_node.DX[combo_name]
429+
D[19, 0] = self.n_node.DY[combo_name]
430+
D[20, 0] = self.n_node.DZ[combo_name]
431+
D[21, 0] = self.n_node.RX[combo_name]
432+
D[22, 0] = self.n_node.RY[combo_name]
433+
D[23, 0] = self.n_node.RZ[combo_name]
434434

435435
# Return the global displacement vector
436436
return D

0 commit comments

Comments
 (0)