Skip to content

Commit 025e286

Browse files
committed
fix: refactor fault ellipsoid plotter to not use fault builder
1 parent 6bd359f commit 025e286

File tree

1 file changed

+11
-6
lines changed

1 file changed

+11
-6
lines changed

LoopStructural/modelling/features/fault/_fault_segment.py

Lines changed: 11 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -124,12 +124,17 @@ def displacementfeature(self):
124124
def fault_ellipsoid(self, **kwargs):
125125
try:
126126
import pyvista as pv
127-
128-
fault_ellipsoid = pv.PolyData(
129-
self.model.rescale(self.fault_centre[None, :], inplace=False)
130-
)
131-
fault_ellipsoid["norm"] = self.builder.fault_normal_vector[None, :]
132-
127+
if self.model is None:
128+
pts = self.fault_centre[None, :]
129+
else:
130+
pts = self.model.rescale(self.fault_centre[None, :], inplace=False)
131+
# pts = self.fault_centre[None, :]
132+
fault_ellipsoid = pv.PolyData(pts)
133+
# fault_ellipsoid = pv.PolyData(
134+
# self.model.rescale(self.fault_centre[None, :], inplace=False)
135+
# )
136+
fault_ellipsoid["norm"] = self.fault_normal_vector[None, :]
137+
fault_ellipsoid['norm'] /= np.linalg.norm(fault_ellipsoid['norm'], axis=1)[:, None]
133138
geom = pv.ParametricEllipsoid(
134139
self.fault_minor_axis,
135140
self.fault_major_axis,

0 commit comments

Comments
 (0)