We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 603ab45 commit 16184bcCopy full SHA for 16184bc
stagpy/plates.py
@@ -37,12 +37,12 @@ def detect_plates_vzcheck(step, vz_thres_ratio=0):
37
dvphi = np.diff(vphi_surf) / (r_c[-1] * dphi)
38
dvx_thres = 16 * step.timeinfo.loc['vrms']
39
40
- limits = []
41
- for phi in range(0, nphi):
42
- mark = all(abs(dvphi[i % nphi]) <= abs(dvphi[phi])
43
- for i in range(phi - nphi // 33, phi + nphi // 33))
44
- if mark and abs(dvphi[phi]) >= dvx_thres:
45
- limits.append(phi)
+ limits = [
+ phi for phi in range(nphi)
+ if (abs(dvphi[phi]) >= dvx_thres and
+ all(abs(dvphi[i % nphi]) <= abs(dvphi[phi])
+ for i in range(phi - nphi // 33, phi + nphi // 33)))
+ ]
46
47
# verifying vertical velocity
48
vz_thres = 0
0 commit comments