rien

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