From 4d7d70260697179bb90dfa9a637914bcb3e11908 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Franck=20P=C3=A9rignon?= <franck.perignon@imag.fr>
Date: Tue, 21 Oct 2014 16:50:58 +0200
Subject: [PATCH] rien
---
HySoP/hysop/operator/discrete/drag_and_lift.py | 16 ++++++++--------
1 file changed, 8 insertions(+), 8 deletions(-)
diff --git a/HySoP/hysop/operator/discrete/drag_and_lift.py b/HySoP/hysop/operator/discrete/drag_and_lift.py
index efe88ecb1..caec11784 100644
--- a/HySoP/hysop/operator/discrete/drag_and_lift.py
+++ b/HySoP/hysop/operator/discrete/drag_and_lift.py
@@ -339,9 +339,9 @@ class NocaForces(Forces):
for j in i_t:
buff[...] = vd[j1[j]][sl] * wd[j2[j]][sl]\
- vd[j2[j]][sl] * wd[j1[j]][sl]
- res[i_n] -= self._coeff * normal[i_n] * npw.real_sum(coords[j]
- * buff)
- res[j] -= self._coeff * normal[i_n] * coords[i_n] *\
+ res[i_n] -= self._coeff * normal * npw.real_sum(coords[j]
+ * buff)
+ res[j] -= self._coeff * normal * coords[i_n] *\
npw.real_sum(buff)
# Last part
@@ -350,18 +350,18 @@ class NocaForces(Forces):
self._laplacian.fd_scheme.computeIndices(sl)
for j in i_t:
self._rwork[...] = self._laplacian(vd[j], self._rwork)
- res[i_n] += self._coeff * self.nu * normal[i_n]\
+ res[i_n] += self._coeff * self.nu * normal\
* npw.real_sum(coords[j] * self._rwork[sl])
- res[j] -= self._coeff * self.nu * normal[i_n] * coords[i_n] * \
+ res[j] -= self._coeff * self.nu * normal * coords[i_n] * \
npw.real_sum(self._rwork[sl])
self._fd_scheme.computeIndices(sl)
self._fd_scheme.compute(vd[i_n], i_n, self._rwork)
- res[i_n] += 2.0 * normal[i_n] * self.nu * npw.real_sum(self._rwork[sl])
+ res[i_n] += 2.0 * normal * self.nu * npw.real_sum(self._rwork[sl])
for j in i_t:
self._fd_scheme.compute(vd[i_n], j, self._rwork)
- res[j] += normal[i_n] * self.nu * npw.real_sum(self._rwork[sl])
+ res[j] += normal * self.nu * npw.real_sum(self._rwork[sl])
self._fd_scheme.compute(vd[j], i_n, self._rwork)
- res[j] += normal[i_n] * self.nu * npw.real_sum(self._rwork[sl])
+ res[j] += normal * self.nu * npw.real_sum(self._rwork[sl])
res *= dsurf
return res
--
GitLab