From 04d794524106885cb6b7fc01d2c31a8989c32622 Mon Sep 17 00:00:00 2001
From: Guy Scharf <guy.oscar@moxis.com>
Date: Sat, 17 Jul 2021 16:22:40 -0700
Subject: [PATCH] Use rolling average when computing TV from flow graph Also,
 start with a guessed-at reasonable value to try to minimized spike artifacts
 resulting from poor calculation of initial breath volume.

---
 oscar/SleepLib/calcs.cpp | 10 ++++++++--
 1 file changed, 8 insertions(+), 2 deletions(-)

diff --git a/oscar/SleepLib/calcs.cpp b/oscar/SleepLib/calcs.cpp
index a7f5b196..dec607f1 100644
--- a/oscar/SleepLib/calcs.cpp
+++ b/oscar/SleepLib/calcs.cpp
@@ -480,7 +480,7 @@ void FlowParser::calc(bool calcResp, bool calcTv, bool calcTi, bool calcTe, bool
     quint32 *tv_tptr = nullptr;
     EventStoreType *tv_dptr = nullptr;
     int tv_count = 0;
-    double tvlast, tvlast2, tvlast3;
+    double tvlast = 0, tvlast2 = 0, tvlast3 = 0;
 
     if (calcTv) {
         TV = m_session->AddEventList(CPAP_TidalVolume, EVL_Event);
@@ -600,8 +600,14 @@ void FlowParser::calc(bool calcResp, bool calcTv, bool calcTi, bool calcTe, bool
             //val2=x;
 
             // Average TV over last three data points
-            if (tv_count == 0)
+            if (tv_count == 0) {
+                if (tv > 800)           // Very much a Q&D patch to handle the first TV value not being calculated well
+                    tv = 300;           // If unreasonable, just set to a "reasonable" number.
                 tvlast = tvlast2 = tvlast3 = tv;
+            }
+            if (tv_count < 4) {
+                qDebug() << "tv" << tv << tvlast << tvlast2 << tvlast3 << "avg" << (tvlast + tvlast2 + tvlast3 + tv*2)/5;
+            }
             tv = (tvlast + tvlast2 + tvlast3 + tv*2)/5;
             tvlast3 = tvlast2;
             tvlast2 = tvlast;