From 1413321d13ddedcc050e86fbb4ea3d2811296959 Mon Sep 17 00:00:00 2001
From: dtgaebe <dtgaebe@sandia.gov>
Date: Tue, 12 Dec 2023 13:54:49 -0800
Subject: [PATCH] corrected sign convention in tut4

---
 examples/tutorial_4_Pioneer.ipynb | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/examples/tutorial_4_Pioneer.ipynb b/examples/tutorial_4_Pioneer.ipynb
index 455fbedb..e1ef6055 100644
--- a/examples/tutorial_4_Pioneer.ipynb
+++ b/examples/tutorial_4_Pioneer.ipynb
@@ -772,11 +772,11 @@
     "P_radiated = ((1/2)*(Rad_res * np.abs(Vel)**2 ).squeeze().sum('omega').item()) # after Falnes Eq. 6.4\n",
     "P_excited= (1/4)*(Fex*np.conjugate(Vel) + np.conjugate(Fex)*Vel ).squeeze().sum('omega').item() # after Falnes Eq. 6.3\n",
     "\n",
-    "power_flows = {'Optimal Excitation' : 2* (np.real(P_max_absorbable)),    #positive because the only inflow\n",
-    "            'Radiated': -1*(np.real(P_radiated)), #negative because \"out\"flow\n",
-    "            'Actual Excitation': -1*(np.real(P_excited)), #negative because \"out\"flow\n",
+    "power_flows = {'Optimal Excitation' : -2* (np.real(P_max_absorbable)),    \n",
+    "            'Radiated': -1*(np.real(P_radiated)), \n",
+    "            'Actual Excitation': -1*(np.real(P_excited)), \n",
     "}\n",
-    "power_flows['Unused Potential'] =  -1*power_flows['Optimal Excitation'] - power_flows['Actual Excitation']\n",
+    "power_flows['Unused Potential'] =  power_flows['Optimal Excitation'] - power_flows['Actual Excitation']\n",
     "power_flows['Absorbed'] =  power_flows['Actual Excitation'] - power_flows['Radiated']  # after Falnes Eq. 6.2 \n",
     "#absorbed by WEC-PTO (difference between actual excitation power and actual radiated power needs to be absorbed by PTO)\n",
     "\n"