From bb97f1977bbf475ec8f1e9f8a52854ccecff1dac Mon Sep 17 00:00:00 2001
From: Theophile Terraz <theophile.terraz@inrae.fr>
Date: Fri, 1 Dec 2023 15:23:59 +0100
Subject: [PATCH] debug
---
src/Model/Geometry/ProfileXYZ.py | 20 ++++++++++----------
1 file changed, 10 insertions(+), 10 deletions(-)
diff --git a/src/Model/Geometry/ProfileXYZ.py b/src/Model/Geometry/ProfileXYZ.py
index 81fa89dc..ad3ebb5e 100644
--- a/src/Model/Geometry/ProfileXYZ.py
+++ b/src/Model/Geometry/ProfileXYZ.py
@@ -461,23 +461,23 @@ class ProfileXYZ(Profile, SQLSubModel):
if self.point(i).z <= z:
ird = i
# interpolation des points ptX et ptY
- if (irg > 0):
+ if (irg < self.number_points-1):
x=np.interp(z,
- [self.point(irg-1).z,self.point(irg).z],
- [self.point(irg-1).x,self.point(irg).x])
+ [self.point(irg).z,self.point(irg+1).z],
+ [self.point(irg).x,self.point(irg+1).x])
y=np.interp(z,
- [self.point(irg-1).z,self.point(irg).z],
- [self.point(irg-1).y,self.point(irg).y])
+ [self.point(irg).z,self.point(irg+1).z],
+ [self.point(irg).y,self.point(irg+1).y])
ptX=PointXYZ(x,y,z)
else:
ptX = self.point(0)
- if (ird < self.number_points-1):
+ if (ird > 0):
x=np.interp(z,
- [self.point(ird).z,self.point(ird+1).z],
- [self.point(ird).x,self.point(ird+1).x])
+ [self.point(ird-1).z,self.point(ird).z],
+ [self.point(ird-1).x,self.point(ird).x])
y=np.interp(z,
- [self.point(ird).z,self.point(ird+1).z],
- [self.point(ird).y,self.point(ird+1).y])
+ [self.point(ird).z,self.point(ird-1).z],
+ [self.point(ird).y,self.point(ird-1).y])
ptY=PointXYZ(x,y,z)
else:
ptY = self.point(self.number_points-1)
--
GitLab