Skip to content

Instantly share code, notes, and snippets.

@stavlor

stavlor/current_patch.diff

Created August 28, 2017 17:05
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save stavlor/83b693277c44d22eb94ad2563ccab20a to your computer and use it in GitHub Desktop.
Save stavlor/83b693277c44d22eb94ad2563ccab20a to your computer and use it in GitHub Desktop.
Download ZIP
Raw
current_patch.diff
diff --git a/app/src/main/assets/OpenAPSAMA/basal-set-temp.js b/app/src/main/assets/OpenAPSAMA/basal-set-temp.js
index 97458691..75ea8835 100644
--- a/app/src/main/assets/OpenAPSAMA/basal-set-temp.js
+++ b/app/src/main/assets/OpenAPSAMA/basal-set-temp.js
@@ -19,7 +19,7 @@ tempBasalFunctions.setTempBasal = function setTempBasal(rate, duration, profile,
//var maxSafeBasal = Math.min(profile.max_basal, 3 * profile.max_daily_basal, 4 * profile.current_basal);
var maxSafeBasal = tempBasalFunctions.getMaxSafeBasal(profile);
-var round_basal = require('./round-basal');
+var round_basal = require('../round-basal');
if (rate < 0) {
rate = 0;
@@ -29,8 +29,8 @@ var round_basal = require('./round-basal');
}
var suggestedRate = round_basal(rate, profile);
- if (typeof(currenttemp) !== 'undefined' && typeof(currenttemp.duration) !== 'undefined' && typeof(currenttemp.rate) !== 'undefined' && currenttemp.duration > 20 && suggestedRate <= currenttemp.rate * 1.2 && suggestedRate >= currenttemp.rate * 0.8) {
- rT.reason += ", but "+currenttemp.duration+"m left and " + currenttemp.rate + " ~ req " + suggestedRate + "U/hr: no action required";
+ if (typeof(currenttemp) !== 'undefined' && typeof(currenttemp.duration) !== 'undefined' && typeof(currenttemp.rate) !== 'undefined' && currenttemp.duration > (duration-10) && currenttemp.duration <= 120 && suggestedRate <= currenttemp.rate * 1.2 && suggestedRate >= currenttemp.rate * 0.8) {
+ rT.reason += " "+currenttemp.duration+"m left and " + currenttemp.rate + " ~ req " + suggestedRate + "U/hr: no temp required";
return rT;
}
@@ -58,4 +58,4 @@ var round_basal = require('./round-basal');
}
};
-module.exports = tempBasalFunctions;
+module.exports = tempBasalFunctions;
\ No newline at end of file
diff --git a/app/src/main/assets/OpenAPSAMA/determine-basal.js b/app/src/main/assets/OpenAPSAMA/determine-basal.js
index 08991568..4e05b618 100644
--- a/app/src/main/assets/OpenAPSAMA/determine-basal.js
+++ b/app/src/main/assets/OpenAPSAMA/determine-basal.js
@@ -1,9 +1,7 @@
/*
Determine Basal
-
Released under MIT license. See the accompanying LICENSE.txt file for
full terms and conditions
-
THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
@@ -19,6 +17,7 @@ var round_basal = require('../round-basal')
// Rounds value to 'digits' decimal places
function round(value, digits)
{
+ if (! digits) { digits = 0; }
var scale = Math.pow(10, digits);
return Math.round(value * scale) / scale;
}
@@ -43,43 +42,48 @@ function convert_bg(value, profile)
}
else
{
- return value.toFixed(0);
+ return Math.round(value);
}
}
-var determine_basal = function determine_basal(glucose_status, currenttemp, iob_data, profile, autosens_data, meal_data, tempBasalFunctions) {
- var rT = { //short for requestedTemp
- };
+var determine_basal = function determine_basal(glucose_status, currenttemp, iob_data, profile, autosens_data, meal_data, tempBasalFunctions, microBolusAllowed, reservoir_data) {
+ var rT = {}; //short for requestedTemp
+
+ var deliverAt = new Date();
if (typeof profile === 'undefined' || typeof profile.current_basal === 'undefined') {
rT.error ='Error: could not get current basal rate';
return rT;
}
- var basal = profile.current_basal;
+ var profile_current_basal = round_basal(profile.current_basal, profile);
+ var basal = profile_current_basal;
if (typeof autosens_data !== 'undefined' ) {
basal = profile.current_basal * autosens_data.ratio;
basal = round_basal(basal, profile);
- if (basal != profile.current_basal) {
- console.error("Adjusting basal from "+profile.current_basal+" to "+basal);
+ if (basal != profile_current_basal) {
+ console.error("Autosens adjusting basal from "+profile_current_basal+" to "+basal+"; ");
+ } else {
+ console.error("Basal unchanged: "+basal+"; ");
}
}
var bg = glucose_status.glucose;
- // TODO: figure out how to use raw isig data to estimate BG
if (bg < 39) { //Dexcom is in ??? mode or calibrating
rT.reason = "CGM is calibrating or in ??? state";
if (basal <= currenttemp.rate * 1.2) { // high temp is running
- rT.reason += "; setting current basal of " + basal + " as temp";
+ rT.reason += "; setting current basal of " + basal + " as temp. ";
+ rT.deliverAt = deliverAt;
+ rT.temp = 'absolute';
return tempBasalFunctions.setTempBasal(basal, 30, profile, rT, currenttemp);
} else { //do nothing.
- rT.reason += ", temp " + currenttemp.rate + " <~ current basal " + basal + "U/hr";
+ rT.reason += ", temp " + currenttemp.rate + " <~ current basal " + basal + "U/hr. ";
return rT;
}
}
var max_iob = profile.max_iob; // maximum amount of non-bolus IOB OpenAPS will ever deliver
- // if target_bg is set, great. otherwise, if min and max are set, then set target to their average
+ // if min and max are set, then set target to their average
var target_bg;
var min_bg;
var max_bg;
@@ -89,36 +93,33 @@ var determine_basal = function determine_basal(glucose_status, currenttemp, iob_
if (typeof profile.max_bg !== 'undefined') {
max_bg = profile.max_bg;
}
- if (typeof profile.target_bg !== 'undefined') {
- target_bg = profile.target_bg;
+ if (typeof profile.min_bg !== 'undefined' && typeof profile.max_bg !== 'undefined') {
+ target_bg = (profile.min_bg + profile.max_bg) / 2;
} else {
- if (typeof profile.min_bg !== 'undefined' && typeof profile.max_bg !== 'undefined') {
- target_bg = (profile.min_bg + profile.max_bg) / 2;
- } else {
- rT.error ='Error: could not determine target_bg';
- return rT;
- }
+ rT.error ='Error: could not determine target_bg. ';
+ return rT;
}
// adjust min, max, and target BG for sensitivity, such that 50% increase in ISF raises target from 100 to 120
if (typeof autosens_data !== 'undefined' && profile.autosens_adjust_targets) {
if (profile.temptargetSet) {
- console.error("Temp Target set, not adjusting with autosens");
+ console.error("Temp Target set, not adjusting with autosens; ");
} else {
- min_bg = Math.round((min_bg - 60) / autosens_data.ratio) + 60;
- max_bg = Math.round((max_bg - 60) / autosens_data.ratio) + 60;
- new_target_bg = Math.round((target_bg - 60) / autosens_data.ratio) + 60;
+ // with a target of 100, default 0.7-1.2 autosens min/max range would allow a 93-117 target range
+ min_bg = round((min_bg - 60) / autosens_data.ratio) + 60;
+ max_bg = round((max_bg - 60) / autosens_data.ratio) + 60;
+ new_target_bg = round((target_bg - 60) / autosens_data.ratio) + 60;
if (target_bg == new_target_bg) {
- console.error("target_bg unchanged:", new_target_bg);
+ console.error("target_bg unchanged: "+new_target_bg+"; ");
} else {
- console.error("Adjusting target_bg from", target_bg, "to", new_target_bg);
+ console.error("target_bg from "+target_bg+" to "+new_target_bg+"; ");
}
target_bg = new_target_bg;
}
}
if (typeof iob_data === 'undefined' ) {
- rT.error ='Error: iob_data undefined';
+ rT.error ='Error: iob_data undefined. ';
return rT;
}
@@ -129,7 +130,7 @@ var determine_basal = function determine_basal(glucose_status, currenttemp, iob_
}
if (typeof iob_data.activity === 'undefined' || typeof iob_data.iob === 'undefined' ) {
- rT.error ='Error: iob_data missing some property';
+ rT.error ='Error: iob_data missing some property. ';
return rT;
}
@@ -140,50 +141,89 @@ var determine_basal = function determine_basal(glucose_status, currenttemp, iob_
} else {
tick = round(glucose_status.delta,0);
}
- var minDelta = Math.min(glucose_status.delta, glucose_status.short_avgdelta, glucose_status.long_avgdelta);
+ //var minDelta = Math.min(glucose_status.delta, glucose_status.short_avgdelta, glucose_status.long_avgdelta);
+ var minDelta = Math.min(glucose_status.delta, glucose_status.short_avgdelta);
var minAvgDelta = Math.min(glucose_status.short_avgdelta, glucose_status.long_avgdelta);
+ var profile_sens = round(profile.sens,1)
var sens = profile.sens;
if (typeof autosens_data !== 'undefined' ) {
sens = profile.sens / autosens_data.ratio;
sens = round(sens, 1);
- if (sens != profile.sens) {
- console.error("Adjusting sens from "+profile.sens+" to "+sens);
+ if (sens != profile_sens) {
+ console.error("sens from "+profile_sens+" to "+sens);
+ } else {
+ console.error("sens unchanged: "+sens);
}
+ console.error(" (autosens ratio "+autosens_data.ratio+")");
}
+ console.error("");
//calculate BG impact: the amount BG "should" be rising or falling based on insulin activity alone
var bgi = round(( -iob_data.activity * sens * 5 ), 2);
// project deviations for 30 minutes
- var deviation = Math.round( 30 / 5 * ( minDelta - bgi ) );
+ var deviation = round( 30 / 5 * ( minDelta - bgi ) );
// don't overreact to a big negative delta: use minAvgDelta if deviation is negative
if (deviation < 0) {
- deviation = Math.round( (30 / 5) * ( minAvgDelta - bgi ) );
+ deviation = round( (30 / 5) * ( minAvgDelta - bgi ) );
}
// calculate the naive (bolus calculator math) eventual BG based on net IOB and sensitivity
if (iob_data.iob > 0) {
- var naive_eventualBG = Math.round( bg - (iob_data.iob * sens) );
+ var naive_eventualBG = round( bg - (iob_data.iob * sens) );
} else { // if IOB is negative, be more conservative and use the lower of sens, profile.sens
- var naive_eventualBG = Math.round( bg - (iob_data.iob * Math.min(sens, profile.sens) ) );
+ var naive_eventualBG = round( bg - (iob_data.iob * Math.min(sens, profile.sens) ) );
}
// and adjust it for the deviation above
var eventualBG = naive_eventualBG + deviation;
// calculate what portion of that is due to bolussnooze
var bolusContrib = iob_data.bolussnooze * sens;
// and add it back in to get snoozeBG, plus another 50% to avoid low-temping at mealtime
- var naive_snoozeBG = Math.round( naive_eventualBG + 1.5 * bolusContrib );
+ var naive_snoozeBG = round( naive_eventualBG + 1.5 * bolusContrib );
// adjust that for deviation like we did eventualBG
var snoozeBG = naive_snoozeBG + deviation;
+ // adjust target BG range if needed to safely bring down high BG faster without causing lows
+ if ( bg > max_bg && profile.adv_target_adjustments ) {
+ // with target=100, as BG rises from 100 to 160, adjustedTarget drops from 100 to 80
+ var adjustedMinBG = round(Math.max(80, min_bg - (bg - min_bg)/3 ),0);
+ var adjustedTargetBG =round( Math.max(80, target_bg - (bg - target_bg)/3 ),0);
+ var adjustedMaxBG = round(Math.max(80, max_bg - (bg - max_bg)/3 ),0);
+ // if eventualBG, naive_eventualBG, and target_bg aren't all above adjustedMinBG, don’t use it
+ //console.error("naive_eventualBG:",naive_eventualBG+", eventualBG:",eventualBG);
+ if (eventualBG > adjustedMinBG && naive_eventualBG > adjustedMinBG && min_bg > adjustedMinBG) {
+ console.error("Adjusting targets for high BG: min_bg from "+min_bg+" to "+adjustedMinBG+"; ");
+ min_bg = adjustedMinBG;
+ } else {
+ console.error("min_bg unchanged: "+min_bg+"; ");
+ }
+ // if eventualBG, naive_eventualBG, and target_bg aren't all above adjustedTargetBG, don’t use it
+ if (eventualBG > adjustedTargetBG && naive_eventualBG > adjustedTargetBG && target_bg > adjustedTargetBG) {
+ console.error("target_bg from "+target_bg+" to "+adjustedTargetBG+"; ");
+ target_bg = adjustedTargetBG;
+ } else {
+ console.error("target_bg unchanged: "+target_bg+"; ");
+ }
+ // if eventualBG, naive_eventualBG, and max_bg aren't all above adjustedMaxBG, don’t use it
+ if (eventualBG > adjustedMaxBG && naive_eventualBG > adjustedMaxBG && max_bg > adjustedMaxBG) {
+ console.error("max_bg from "+max_bg+" to "+adjustedMaxBG);
+ max_bg = adjustedMaxBG;
+ } else {
+ console.error("max_bg unchanged: "+max_bg);
+ }
+ }
+
var expectedDelta = calculate_expected_delta(profile.dia, target_bg, eventualBG, bgi);
if (typeof eventualBG === 'undefined' || isNaN(eventualBG)) {
- rT.error ='Error: could not calculate eventualBG';
+ rT.error ='Error: could not calculate eventualBG. ';
return rT;
}
- // min_bg of 90 -> threshold of 70, 110 -> 80, and 130 -> 90
- var threshold = min_bg - 0.5*(min_bg-50);
+ // min_bg of 90 -> threshold of 65, 100 -> 70 110 -> 75, and 130 -> 85
+ var threshold = min_bg - 0.5*(min_bg-40);
+
+ //console.error(reservoir_data);
+ var deliverAt = new Date();
rT = {
'temp': 'absolute'
@@ -191,80 +231,178 @@ var determine_basal = function determine_basal(glucose_status, currenttemp, iob_
, 'tick': tick
, 'eventualBG': eventualBG
, 'snoozeBG': snoozeBG
+ , 'insulinReq': 0
+ , 'reservoir' : reservoir_data // The expected reservoir volume at which to deliver the microbolus (the reservoir volume from immediately before the last pumphistory run)
+ , 'deliverAt' : deliverAt // The time at which the microbolus should be delivered
+ , 'minPredBG' : 999
};
- var basaliob;
- if (iob_data.basaliob) { basaliob = iob_data.basaliob; }
- else { basaliob = iob_data.iob - iob_data.bolussnooze; }
+ var basaliob = iob_data.basaliob;
+ //if (iob_data.basaliob) { basaliob = iob_data.basaliob; }
+ //else { basaliob = iob_data.iob - iob_data.bolussnooze; }
+ var bolusiob = iob_data.iob - basaliob;
// generate predicted future BGs based on IOB, COB, and current absorption rate
var COBpredBGs = [];
var aCOBpredBGs = [];
var IOBpredBGs = [];
+ var UAMpredBGs = [];
COBpredBGs.push(bg);
aCOBpredBGs.push(bg);
IOBpredBGs.push(bg);
+ UAMpredBGs.push(bg);
+
+ // enable SMB whenever we have COB or UAM is enabled
+ // SMB is diabled by default, unless explicitly enabled in preferences.json
+ var enableSMB=false;
+ // disable SMB when a high temptarget is set
+ if (profile.temptargetSet && target_bg > 100) {
+ enableSMB=false;
+ // enable SMB/UAM (if enabled in preferences) for DIA hours after bolus
+ } else if (profile.enableSMB_with_bolus && bolusiob > 0.1) {
+ enableSMB=true;
+ // enable SMB/UAM (if enabled in preferences) while we have COB
+ } else if (profile.enableSMB_with_COB && meal_data.mealCOB) {
+ enableSMB=true;
+ // enable SMB/UAM (if enabled in preferences) if a low temptarget is set
+ } else if (profile.enableSMB_with_temptarget && (profile.temptargetSet && target_bg < 100)) {
+ enableSMB=true;
+ // enable SMB/UAM (if enabled in preferences) for a full 6 hours after any carb entry
+ // (6 hours is defined in carbWindow in lib/meal/total.js)
+ } else if (profile.enableSMB_after_carbs && meal_data.carbs) {
+ enableSMB=true;
+ }
+ // enable UAM (if enabled in preferences) for DIA hours after bolus, or if SMB is enabled
+ var enableUAM=(profile.enableUAM && (bolusiob > 0.1 || enableSMB));
+
+
//console.error(meal_data);
// carb impact and duration are 0 unless changed below
var ci = 0;
var cid = 0;
// calculate current carb absorption rate, and how long to absorb all carbs
// CI = current carb impact on BG in mg/dL/5m
- ci = Math.round((minDelta - bgi)*10)/10;
- if (meal_data.mealCOB * 2 > meal_data.carbs) {
- // set ci to a minimum of 3mg/dL/5m (default) if less than half of carbs have absorbed
- ci = Math.max(profile.min_5m_carbimpact, ci);
- }
+ ci = round((minDelta - bgi),1);
+ uci = round((minAvgDelta - bgi),1);
+ // ISF (mg/dL/U) / CR (g/U) = CSF (mg/dL/g)
+ var csf = sens / profile.carb_ratio
+ // set meal_carbimpact high enough to absorb all meal carbs over 6 hours
+ // total_impact (mg/dL) = CSF (mg/dL/g) * carbs (g)
+ //console.error(csf * meal_data.carbs);
+ // meal_carbimpact (mg/dL/5m) = CSF (mg/dL/g) * carbs (g) / 6 (h) * (1h/60m) * 5 (m/5m) * 2 (for linear decay)
+ //var meal_carbimpact = round((csf * meal_data.carbs / 6 / 60 * 5 * 2),1)
+ // calculate the number of carbs absorbed over 4h at current CI
+ // CI (mg/dL/5m) * (5m)/5 (m) * 60 (min/hr) * 4 (h) / 2 (linear decay factor) = total carb impact (mg/dL)
+ var totalCI = Math.max(0, ci / 5 * 60 * 4 / 2);
+ // totalCI (mg/dL) / CSF (mg/dL/g) = total carbs absorbed (g)
+ var totalCA = totalCI / csf;
+ // exclude the last 1/3 of carbs from remainingCarbs, and then cap it at 90
+ var remainingCarbsCap = 90; // default to 90
+ var remainingCarbsFraction = 1;
+ if (profile.remainingCarbsCap) { remainingCarbsCap = Math.min(90,profile.remainingCarbsCap); }
+ if (profile.remainingCarbsFraction) { remainingCarbsFraction = Math.min(1,profile.remainingCarbsFraction); }
+ var remainingCarbsIgnore = 1 - remainingCarbsFraction;
+ var remainingCarbs = Math.max(0, meal_data.mealCOB - totalCA - meal_data.carbs*remainingCarbsIgnore);
+ remainingCarbs = Math.min(remainingCarbsCap,remainingCarbs);
+ // assume remainingCarbs will absorb over 4h
+ // remainingCI (mg/dL/5m) = remainingCarbs (g) * CSF (mg/dL/g) * 5 (m/5m) * 1h/60m / 4 (h)
+ var remainingCI = remainingCarbs * csf * 5 / 60 / 4;
+ //console.error(profile.min_5m_carbimpact,ci,totalCI,totalCA,remainingCarbs,remainingCI);
+ //if (meal_data.mealCOB * 3 > meal_data.carbs) { }
+
+ // calculate peak deviation in last hour, and slope from that to current deviation
+ var minDeviationSlope = round(meal_data.minDeviationSlope,2);
+ //console.error(minDeviationSlope);
+
aci = 10;
//5m data points = g * (1U/10g) * (40mg/dL/1U) / (mg/dL/5m)
- cid = meal_data.mealCOB * ( sens / profile.carb_ratio ) / ci;
- acid = meal_data.mealCOB * ( sens / profile.carb_ratio ) / aci;
- console.error("Carb Impact:",ci,"mg/dL per 5m; CI Duration:",Math.round(10*cid/6)/10,"hours");
- console.error("Accel. Carb Impact:",aci,"mg/dL per 5m; ACI Duration:",Math.round(10*acid/6)/10,"hours");
- var minPredBG = 999;
- var maxPredBG = bg;
+ // duration (in 5m data points) = COB (g) * CSF (mg/dL/g) / ci (mg/dL/5m)
+ cid = Math.max(0, meal_data.mealCOB * csf / ci );
+ acid = Math.max(0, meal_data.mealCOB * csf / aci );
+ // duration (hours) = duration (5m) * 5 / 60 * 2 (to account for linear decay)
+ console.error("Carb Impact:",ci,"mg/dL per 5m; CI Duration:",round(cid*5/60*2,1),"hours; remaining 4h+ CI:",round(remainingCI,1),"mg/dL per 5m");
+ console.error("Accel. Carb Impact:",aci,"mg/dL per 5m; ACI Duration:",round(acid*5/60*2,1),"hours");
+ var minIOBPredBG = 999;
+ var minCOBPredBG = 999;
+ var minUAMPredBG = 999;
+ var minPredBG;
+ var avgPredBG;
+ var IOBpredBG = eventualBG;
+ var maxIOBPredBG = bg;
+ var maxCOBPredBG = bg;
+ var maxUAMPredBG = bg;
+ //var maxPredBG = bg;
var eventualPredBG = bg;
+ var lastIOBpredBG;
+ var lastCOBpredBG;
+ var lastUAMpredBG;
+ var UAMduration = 0;
try {
iobArray.forEach(function(iobTick) {
//console.error(iobTick);
predBGI = round(( -iobTick.activity * sens * 5 ), 2);
- // predicted deviation impact drops linearly from current deviation down to zero
+ // for IOBpredBGs, predicted deviation impact drops linearly from current deviation down to zero
// over 60 minutes (data points every 5m)
predDev = ci * ( 1 - Math.min(1,IOBpredBGs.length/(60/5)) );
IOBpredBG = IOBpredBGs[IOBpredBGs.length-1] + predBGI + predDev;
//IOBpredBG = IOBpredBGs[IOBpredBGs.length-1] + predBGI;
- // predicted carb impact drops linearly from current carb impact down to zero
+ // for COBpredBGs, predicted carb impact drops linearly from current carb impact down to zero
// eventually accounting for all carbs (if they can be absorbed over DIA)
- predCI = Math.max(0, ci * ( 1 - COBpredBGs.length/Math.max(cid*2,1) ) );
- predACI = Math.max(0, aci * ( 1 - COBpredBGs.length/Math.max(acid*2,1) ) );
- COBpredBG = COBpredBGs[COBpredBGs.length-1] + predBGI + Math.min(0,predDev) + predCI;
+ predCI = Math.max(0, Math.max(0,ci) * ( 1 - COBpredBGs.length/Math.max(cid*2,1) ) );
+ predACI = Math.max(0, Math.max(0,aci) * ( 1 - COBpredBGs.length/Math.max(acid*2,1) ) );
+ // if any carbs aren't absorbed after 4 hours, assume they'll absorb at a constant rate for next 4h
+ COBpredBG = COBpredBGs[COBpredBGs.length-1] + predBGI + Math.min(0,predDev) + predCI + remainingCI;
+ // stop adding remainingCI after 4h
+ if (COBpredBGs.length > 4 * 60 / 5) { remainingCI = 0; }
aCOBpredBG = aCOBpredBGs[aCOBpredBGs.length-1] + predBGI + Math.min(0,predDev) + predACI;
- //console.error(predBGI, predCI, predBG);
- IOBpredBGs.push(IOBpredBG);
- COBpredBGs.push(COBpredBG);
- aCOBpredBGs.push(aCOBpredBG);
- // wait 45m before setting minPredBG
- if ( COBpredBGs.length > 9 && (COBpredBG < minPredBG) ) { minPredBG = COBpredBG; }
- if ( COBpredBG > maxPredBG ) { maxPredBG = COBpredBG; }
+ // for UAMpredBGs, predicted carb impact drops at minDeviationSlope
+ // calculate predicted CI from UAM based on minDeviationSlope
+ predUCIslope = Math.max(0, uci + ( UAMpredBGs.length*minDeviationSlope ) );
+ // if minDeviationSlope is too flat, predicted deviation impact drops linearly from
+ // current deviation down to zero over DIA (data points every 5m)
+ predUCIdia = Math.max(0, uci * ( 1 - UAMpredBGs.length/Math.max(profile.dia*60/5,1) ) );
+ //console.error(predUCIslope, predUCIdia);
+ // predicted CI from UAM is the lesser of CI based on deviationSlope or DIA
+ predUCI = Math.min(predUCIslope, predUCIdia);
+ if(predUCI>0) {
+ //console.error(UAMpredBGs.length,minDeviationSlope, predUCI);
+ UAMduration=round((UAMpredBGs.length+1)*5/60,1);
+ }
+ UAMpredBG = UAMpredBGs[UAMpredBGs.length-1] + predBGI + Math.min(0, predDev) + predUCI;
+ //console.error(predBGI, predCI, predUCI);
+ // truncate all BG predictions at 3.5 hours
+ if ( IOBpredBGs.length < 42) { IOBpredBGs.push(IOBpredBG); }
+ if ( COBpredBGs.length < 42) { COBpredBGs.push(COBpredBG); }
+ if ( aCOBpredBGs.length < 42) { aCOBpredBGs.push(aCOBpredBG); }
+ if ( UAMpredBGs.length < 42) { UAMpredBGs.push(UAMpredBG); }
+ // wait 90m before setting minIOBPredBG
+ if ( IOBpredBGs.length > 18 && (IOBpredBG < minIOBPredBG) ) { minIOBPredBG = round(IOBpredBG); }
+ if ( IOBpredBG > maxIOBPredBG ) { maxIOBPredBG = IOBpredBG; }
+ // wait 60m before setting COB and UAM minPredBGs
+ if ( (cid || remainingCI > 0) && COBpredBGs.length > 12 && (COBpredBG < minCOBPredBG) ) { minCOBPredBG = round(COBpredBG); }
+ if ( (cid || remainingCI > 0) && COBpredBG > maxIOBPredBG ) { maxCOBPredBG = COBpredBG; }
+ if ( enableUAM && UAMpredBGs.length > 12 && (UAMpredBG < minUAMPredBG) ) { minUAMPredBG = round(UAMpredBG); }
+ if ( enableUAM && UAMpredBG > maxIOBPredBG ) { maxUAMPredBG = UAMpredBG; }
});
// set eventualBG to include effect of carbs
//console.error("PredBGs:",JSON.stringify(predBGs));
} catch (e) {
- console.error("Problem with iobArray. Optional feature Advanced Meal Assist disabled.");
+ console.error("Problem with iobArray. Optional feature Advanced Meal Assist disabled:",e);
}
rT.predBGs = {};
IOBpredBGs.forEach(function(p, i, theArray) {
- theArray[i] = Math.round(Math.min(401,Math.max(39,p)));
+ theArray[i] = round(Math.min(401,Math.max(39,p)));
});
for (var i=IOBpredBGs.length-1; i > 12; i--) {
if (IOBpredBGs[i-1] != IOBpredBGs[i]) { break; }
else { IOBpredBGs.pop(); }
}
rT.predBGs.IOB = IOBpredBGs;
+ lastIOBpredBG=round(IOBpredBGs[IOBpredBGs.length-1]);
if (meal_data.mealCOB > 0) {
aCOBpredBGs.forEach(function(p, i, theArray) {
- theArray[i] = Math.round(Math.min(401,Math.max(39,p)));
+ theArray[i] = round(Math.min(401,Math.max(39,p)));
});
for (var i=aCOBpredBGs.length-1; i > 12; i--) {
if (aCOBpredBGs[i-1] != aCOBpredBGs[i]) { break; }
@@ -272,32 +410,168 @@ var determine_basal = function determine_basal(glucose_status, currenttemp, iob_
}
rT.predBGs.aCOB = aCOBpredBGs;
}
- if (meal_data.mealCOB > 0 && ci > 0 ) {
+ if (meal_data.mealCOB > 0 && ( ci > 0 || remainingCI > 0 )) {
COBpredBGs.forEach(function(p, i, theArray) {
- theArray[i] = Math.round(Math.min(401,Math.max(39,p)));
+ theArray[i] = round(Math.min(401,Math.max(39,p)));
});
for (var i=COBpredBGs.length-1; i > 12; i--) {
if (COBpredBGs[i-1] != COBpredBGs[i]) { break; }
else { COBpredBGs.pop(); }
}
rT.predBGs.COB = COBpredBGs;
- eventualBG = Math.max(eventualBG, Math.round(COBpredBGs[COBpredBGs.length-1]) );
+ lastCOBpredBG=round(COBpredBGs[COBpredBGs.length-1]);
+ eventualBG = Math.max(eventualBG, round(COBpredBGs[COBpredBGs.length-1]) );
+ }
+ if (ci > 0 || remainingCI > 0) {
+ if (enableUAM) {
+ UAMpredBGs.forEach(function(p, i, theArray) {
+ theArray[i] = round(Math.min(401,Math.max(39,p)));
+ });
+ for (var i=UAMpredBGs.length-1; i > 12; i--) {
+ if (UAMpredBGs[i-1] != UAMpredBGs[i]) { break; }
+ else { UAMpredBGs.pop(); }
+ }
+ rT.predBGs.UAM = UAMpredBGs;
+ lastUAMpredBG=round(UAMpredBGs[UAMpredBGs.length-1]);
+ eventualBG = Math.max(eventualBG, round(UAMpredBGs[UAMpredBGs.length-1]) );
+ }
+
+ // set eventualBG and snoozeBG based on COB or UAM predBGs
rT.eventualBG = eventualBG;
- minPredBG = Math.min(minPredBG, eventualBG);
- // set snoozeBG to minPredBG
- snoozeBG = Math.round(Math.max(snoozeBG,minPredBG));
- rT.snoozeBG = snoozeBG;
}
+ console.error("UAM Impact:",uci,"mg/dL per 5m; UAM Duration:",UAMduration,"hours");
+
+ minIOBPredBG = Math.max(39,minIOBPredBG);
+ minCOBPredBG = Math.max(39,minCOBPredBG);
+ minUAMPredBG = Math.max(39,minUAMPredBG);
+ minPredBG = round(minIOBPredBG);
+
+ // if we have COB and UAM is enabled, average all three
+ if ( minUAMPredBG < 400 && minCOBPredBG < 400 ) {
+ avgPredBG = round( (IOBpredBG + UAMpredBG + COBpredBG)/3 );
+ // if UAM is disabled, average IOB and COB
+ } else if ( minCOBPredBG < 400 ) {
+ avgPredBG = round( (IOBpredBG + COBpredBG)/2 );
+ // if carbs are expired, use IOB instead of COB
+ } else if ( meal_data.carbs && minUAMPredBG < 400 ) {
+ avgPredBG = round( (2*IOBpredBG + UAMpredBG)/3 );
+ // in pure UAM mode, just average IOB and UAM
+ } else if ( minUAMPredBG < 400 ) {
+ avgPredBG = round( (IOBpredBG + UAMpredBG)/2 );
+ } else {
+ avgPredBG = round( IOBpredBG );
+ }
+
+ // if any carbs have been entered recently
+ if (meal_data.carbs) {
+ // average the minIOBPredBG and minUAMPredBG if available
+ if ( minUAMPredBG < 400 ) {
+ avgMinPredBG = round( (minIOBPredBG+minUAMPredBG)/2 );
+ } else {
+ avgMinPredBG = minIOBPredBG;
+ }
+
+ // if UAM is disabled, use max of minIOBPredBG, minCOBPredBG
+ if ( ! enableUAM && minCOBPredBG < 400 ) {
+ minPredBG = round(Math.max(minIOBPredBG, minCOBPredBG));
+ // if we have COB, use minCOBPredBG, or blendedMinPredBG if it's higher
+ } else if ( minCOBPredBG < 400 ) {
+ // calculate blendedMinPredBG based on how many carbs remain as COB
+ fractionCarbsLeft = meal_data.mealCOB/meal_data.carbs;
+ blendedMinPredBG = fractionCarbsLeft*minCOBPredBG + (1-fractionCarbsLeft)*avgMinPredBG;
+ // if blendedMinPredBG > minCOBPredBG, use that instead
+ minPredBG = round(Math.max(minIOBPredBG, minCOBPredBG, blendedMinPredBG));
+ // if carbs have been entered, but have expired, use avg of minIOBPredBG and minUAMPredBG
+ } else {
+ minPredBG = avgMinPredBG;
+ }
+ // in pure UAM mode, use the higher of minIOBPredBG,minUAMPredBG
+ } else if ( enableUAM ) {
+ minPredBG = round(Math.max(minIOBPredBG,minUAMPredBG));
+ }
+
+ // make sure minPredBG isn't higher than avgPredBG
+ minPredBG = Math.min( minPredBG, avgPredBG );
+
+ console.error("minPredBG: "+minPredBG+" minIOBPredBG: "+minIOBPredBG);
+ if (minCOBPredBG < 400) {
+ console.error(" minCOBPredBG: "+minCOBPredBG);
+ }
+ if (minUAMPredBG < 400) {
+ console.error(" minUAMPredBG: "+minUAMPredBG);
+ }
+ console.error(" avgPredBG:",avgPredBG,"COB:",meal_data.mealCOB,"carbs:",meal_data.carbs);
+ // But if the COB line falls off a cliff, don't trust UAM too much:
+ // use maxCOBPredBG if it's been set and lower than minPredBG
+ if ( maxCOBPredBG > bg ) {
+ minPredBG = Math.min(minPredBG, maxCOBPredBG);
+ }
+ // set snoozeBG to minPredBG if it's higher
+ if (minPredBG < 400) {
+ snoozeBG = round(Math.max(snoozeBG,minPredBG));
+ }
+ rT.snoozeBG = snoozeBG;
+ //console.error(minPredBG, minIOBPredBG, minUAMPredBG, minCOBPredBG, maxCOBPredBG, snoozeBG);
+
rT.COB=meal_data.mealCOB;
rT.IOB=iob_data.iob;
- rT.reason="COB: " + meal_data.mealCOB + ", Dev: " + deviation + ", BGI: " + bgi + ", ISF: " + convert_bg(sens, profile) + ", Target: " + convert_bg(target_bg, profile) + "; ";
- if (typeof autosens_data !== 'undefined' && profile.autosens_adjust_targets && autosens_data.ratio != 1)
- rT.reason += "Autosens: " + autosens_data.ratio + "; ";
- if (bg < threshold) { // low glucose suspend mode: BG is < ~80
+ rT.reason="COB: " + meal_data.mealCOB + ", Dev: " + deviation + ", BGI: " + bgi + ", ISF: " + convert_bg(sens, profile) + ", Target: " + convert_bg(target_bg, profile) + ", minPredBG " + convert_bg(minPredBG, profile) + ", IOBpredBG " + convert_bg(lastIOBpredBG, profile);
+ if (lastCOBpredBG > 0) {
+ rT.reason += ", COBpredBG " + convert_bg(lastCOBpredBG, profile);
+ }
+ if (lastUAMpredBG > 0) {
+ rT.reason += ", UAMpredBG " + convert_bg(lastUAMpredBG, profile)
+ }
+ rT.reason += "; ";
+ var bgUndershoot = target_bg - Math.max( naive_eventualBG, eventualBG, lastIOBpredBG );
+ // calculate how long until COB (or IOB) predBGs drop below min_bg
+ var minutesAboveMinBG = 240;
+ if (meal_data.mealCOB > 0 && ( ci > 0 || remainingCI > 0 )) {
+ for (var i=0; i<COBpredBGs.length; i++) {
+ //console.error(COBpredBGs[i], min_bg);
+ if ( COBpredBGs[i] < min_bg ) {
+ minutesAboveMinBG = 5*i;
+ break;
+ }
+ }
+ } else {
+ for (var i=0; i<IOBpredBGs.length; i++) {
+ //console.error(IOBpredBGs[i], min_bg);
+ if ( IOBpredBGs[i] < min_bg ) {
+ minutesAboveMinBG = 5*i;
+ break;
+ }
+ }
+ }
+
+ console.error("BG projected to remain above",min_bg,"for",minutesAboveMinBG,"minutes");
+ // include at least minutesAboveMinBG worth of zero temps in calculating carbsReq
+ // always include at least 30m worth of zero temp (carbs to 80, low temp up to target)
+ var zeroTempDuration = Math.max(30,minutesAboveMinBG);
+ // BG undershoot, minus effect of zero temps until hitting min_bg, converted to grams, minus COB
+ var zeroTempEffect = profile.current_basal*sens*zeroTempDuration/60;
+ var carbsReq = (bgUndershoot - zeroTempEffect) / csf - meal_data.mealCOB;
+ zeroTempEffect = round(zeroTempEffect);
+ carbsReq = round(carbsReq);
+ console.error("bgUndershoot:",bgUndershoot,"zeroTempDuration:",zeroTempDuration,"zeroTempEffect:",zeroTempEffect,"carbsReq:",carbsReq);
+ if ( carbsReq > 0 ) {
+ rT.carbsReq = carbsReq;
+ rT.reason += carbsReq + " add'l carbs req + " + minutesAboveMinBG + "m zero temp; ";
+ }
+ // don't low glucose suspend if IOB is already super negative and BG is rising faster than predicted
+ if (bg < threshold && iob_data.iob < -profile.current_basal*20/60 && minDelta > 0 && minDelta > expectedDelta) {
+ rT.reason += "IOB "+iob_data.iob+" < " + round(-profile.current_basal*20/60,2);
+ rT.reason += " and minDelta " + minDelta + " > " + "expectedDelta " + expectedDelta + "; ";
+ }
+ // low glucose suspend mode: BG is < ~80
+ else if (bg < threshold) {
rT.reason += "BG " + convert_bg(bg, profile) + "<" + convert_bg(threshold, profile);
if ((glucose_status.delta <= 0 && minDelta <= 0) || (glucose_status.delta < expectedDelta && minDelta < expectedDelta) || bg < 60 ) {
// BG is still falling / rising slower than predicted
+ if ( minDelta < expectedDelta ) {
+ rT.reason += ", minDelta " + minDelta + " < " + "expectedDelta " + expectedDelta + "; ";
+ }
return tempBasalFunctions.setTempBasal(0, 30, profile, rT, currenttemp);
}
if (glucose_status.delta > minDelta) {
@@ -306,10 +580,10 @@ var determine_basal = function determine_basal(glucose_status, currenttemp, iob_
rT.reason += ", min delta " + minDelta.toFixed(2) + ">0";
}
if (currenttemp.duration > 15 && (round_basal(basal, profile) === round_basal(currenttemp.rate, profile))) {
- rT.reason += ", temp " + currenttemp.rate + " ~ req " + basal + "U/hr";
+ rT.reason += ", temp " + currenttemp.rate + " ~ req " + basal + "U/hr. ";
return rT;
} else {
- rT.reason += "; setting current basal of " + basal + " as temp";
+ rT.reason += "; setting current basal of " + basal + " as temp. ";
return tempBasalFunctions.setTempBasal(basal, 30, profile, rT, currenttemp);
}
}
@@ -318,16 +592,21 @@ var determine_basal = function determine_basal(glucose_status, currenttemp, iob_
rT.reason += "Eventual BG " + convert_bg(eventualBG, profile) + " < " + convert_bg(min_bg, profile);
// if 5m or 30m avg BG is rising faster than expected delta
if (minDelta > expectedDelta && minDelta > 0) {
+ // if naive_eventualBG < 40, set a 30m zero temp (oref0-pump-loop will let any longer SMB zero temp run)
+ if (naive_eventualBG < 40) {
+ rT.reason += ", naive_eventualBG < 40. ";
+ return tempBasalFunctions.setTempBasal(0, 30, profile, rT, currenttemp);
+ }
if (glucose_status.delta > minDelta) {
- rT.reason += ", but Delta " + tick + " > Exp. Delta " + expectedDelta;
+ rT.reason += ", but Delta " + tick + " > expectedDelta " + expectedDelta;
} else {
rT.reason += ", but Min. Delta " + minDelta.toFixed(2) + " > Exp. Delta " + expectedDelta;
}
if (currenttemp.duration > 15 && (round_basal(basal, profile) === round_basal(currenttemp.rate, profile))) {
- rT.reason += ", temp " + currenttemp.rate + " ~ req " + basal + "U/hr";
+ rT.reason += ", temp " + currenttemp.rate + " ~ req " + basal + "U/hr. ";
return rT;
} else {
- rT.reason += "; setting current basal of " + basal + " as temp";
+ rT.reason += "; setting current basal of " + basal + " as temp. ";
return tempBasalFunctions.setTempBasal(basal, 30, profile, rT, currenttemp);
}
}
@@ -335,14 +614,17 @@ var determine_basal = function determine_basal(glucose_status, currenttemp, iob_
if (eventualBG < min_bg) {
// if we've bolused recently, we can snooze until the bolus IOB decays (at double speed)
if (snoozeBG > min_bg) { // if adding back in the bolus contribution BG would be above min
- rT.reason += ", bolus snooze: eventual BG range " + convert_bg(eventualBG, profile) + "-" + convert_bg(snoozeBG, profile);
- //console.error(currenttemp, basal );
- if (currenttemp.duration > 15 && (round_basal(basal, profile) === round_basal(currenttemp.rate, profile))) {
- rT.reason += ", temp " + currenttemp.rate + " ~ req " + basal + "U/hr";
- return rT;
- } else {
- rT.reason += "; setting current basal of " + basal + " as temp";
- return tempBasalFunctions.setTempBasal(basal, 30, profile, rT, currenttemp);
+ // If we're not in SMB mode with COB, or lastCOBpredBG > target_bg, bolus snooze
+ if (! (microBolusAllowed && rT.COB) || lastCOBpredBG > target_bg) {
+ rT.reason += ", bolus snooze: eventual BG range " + convert_bg(eventualBG, profile) + "-" + convert_bg(snoozeBG, profile);
+ //console.error(currenttemp, basal );
+ if (currenttemp.duration > 15 && (round_basal(basal, profile) === round_basal(currenttemp.rate, profile))) {
+ rT.reason += ", temp " + currenttemp.rate + " ~ req " + basal + "U/hr. ";
+ return rT;
+ } else {
+ rT.reason += "; setting current basal of " + basal + " as temp. ";
+ return tempBasalFunctions.setTempBasal(basal, 30, profile, rT, currenttemp);
+ }
}
} else {
// calculate 30m low-temp required to get projected BG up to target
@@ -350,6 +632,9 @@ var determine_basal = function determine_basal(glucose_status, currenttemp, iob_
// multiply by 2 to low-temp faster for increased hypo safety
var insulinReq = 2 * Math.min(0, (snoozeBG - target_bg) / sens);
insulinReq = round( insulinReq , 2);
+ // calculate naiveInsulinReq based on naive_eventualBG
+ var naiveInsulinReq = Math.min(0, (naive_eventualBG - target_bg) / sens);
+ naiveInsulinReq = round( naiveInsulinReq , 2);
if (minDelta < 0 && minDelta > expectedDelta) {
// if we're barely falling, newinsulinReq should be barely negative
rT.reason += ", Snooze BG " + convert_bg(snoozeBG, profile);
@@ -362,21 +647,45 @@ var determine_basal = function determine_basal(glucose_status, currenttemp, iob_
rate = round_basal(rate, profile);
// if required temp < existing temp basal
var insulinScheduled = currenttemp.duration * (currenttemp.rate - basal) / 60;
- if (insulinScheduled < insulinReq - basal*0.3) { // if current temp would deliver a lot (30% of basal) less than the required insulin, raise the rate
- rT.reason += ", "+currenttemp.duration + "m@" + (currenttemp.rate - basal).toFixed(3) + " = " + insulinScheduled.toFixed(3) + " < req " + insulinReq + "-" + basal*0.3;
+ // if current temp would deliver a lot (30% of basal) less than the required insulin,
+ // by both normal and naive calculations, then raise the rate
+ var minInsulinReq = Math.min(insulinReq,naiveInsulinReq);
+ if (insulinScheduled < minInsulinReq - basal*0.3) {
+ rT.reason += ", "+currenttemp.duration + "m@" + (currenttemp.rate).toFixed(2) + " is a lot less than needed. ";
return tempBasalFunctions.setTempBasal(rate, 30, profile, rT, currenttemp);
}
if (typeof currenttemp.rate !== 'undefined' && (currenttemp.duration > 5 && rate >= currenttemp.rate * 0.8)) {
- rT.reason += ", temp " + currenttemp.rate + " ~< req " + rate + "U/hr";
+ rT.reason += ", temp " + currenttemp.rate + " ~< req " + rate + "U/hr. ";
return rT;
} else {
- rT.reason += ", setting " + rate + "U/hr";
+ // calculate a long enough zero temp to eventually correct back up to target
+ if ( rate < 0 ) {
+ var bgUndershoot = target_bg - naive_eventualBG;
+ var worstCaseInsulinReq = bgUndershoot / sens;
+ var durationReq = round(60*worstCaseInsulinReq / profile.current_basal);
+ if (durationReq < 0) {
+ durationReq = 0;
+ // don't set a temp longer than 120 minutes
+ } else {
+ durationReq = round(durationReq/30)*30;
+ durationReq = Math.min(120,Math.max(0,durationReq));
+ }
+ //console.error(durationReq);
+ //rT.reason += "insulinReq " + insulinReq + "; "
+ if (durationReq > 0) {
+ rT.reason += ", setting " + durationReq + "m zero temp. ";
+ return tempBasalFunctions.setTempBasal(rate, durationReq, profile, rT, currenttemp);
+ }
+ } else {
+ rT.reason += ", setting " + rate + "U/hr. ";
+ }
return tempBasalFunctions.setTempBasal(rate, 30, profile, rT, currenttemp);
}
}
}
}
+ /*
var minutes_running;
if (typeof currenttemp.duration == 'undefined' || currenttemp.duration == 0) {
minutes_running = 30;
@@ -386,70 +695,82 @@ var determine_basal = function determine_basal(glucose_status, currenttemp, iob_
} else {
minutes_running = 30 - currenttemp.duration;
}
-
// if there is a low-temp running, and eventualBG would be below min_bg without it, let it run
if (round_basal(currenttemp.rate, profile) < round_basal(basal, profile) ) {
var lowtempimpact = (currenttemp.rate - basal) * ((30-minutes_running)/60) * sens;
var adjEventualBG = eventualBG + lowtempimpact;
- if ( adjEventualBG < min_bg ) {
+ // don't return early if microBolusAllowed etc.
+ if ( adjEventualBG < min_bg && ! (microBolusAllowed && enableSMB)) {
rT.reason += "letting low temp of " + currenttemp.rate + " run.";
return rT;
}
}
+ */
// if eventual BG is above min but BG is falling faster than expected Delta
if (minDelta < expectedDelta) {
- if (glucose_status.delta < minDelta) {
- rT.reason += "Eventual BG " + convert_bg(eventualBG, profile) + " > " + convert_bg(min_bg, profile) + " but Delta " + tick + " < Exp. Delta " + expectedDelta;
- } else {
- rT.reason += "Eventual BG " + convert_bg(eventualBG, profile) + " > " + convert_bg(min_bg, profile) + " but Min. Delta " + minDelta.toFixed(2) + " < Exp. Delta " + expectedDelta;
- }
- if (currenttemp.duration > 15 && (round_basal(basal, profile) === round_basal(currenttemp.rate, profile))) {
- rT.reason += ", temp " + currenttemp.rate + " ~ req " + basal + "U/hr";
- return rT;
- } else {
- rT.reason += "; setting current basal of " + basal + " as temp";
- return tempBasalFunctions.setTempBasal(basal, 30, profile, rT, currenttemp);
+ // if in SMB mode, don't cancel SMB zero temp
+ if (! (microBolusAllowed && enableSMB)) {
+ if (glucose_status.delta < minDelta) {
+ rT.reason += "Eventual BG " + convert_bg(eventualBG, profile) + " > " + convert_bg(min_bg, profile) + " but Delta " + tick + " < Exp. Delta " + expectedDelta;
+ } else {
+ rT.reason += "Eventual BG " + convert_bg(eventualBG, profile) + " > " + convert_bg(min_bg, profile) + " but Min. Delta " + minDelta.toFixed(2) + " < Exp. Delta " + expectedDelta;
+ }
+ if (currenttemp.duration > 15 && (round_basal(basal, profile) === round_basal(currenttemp.rate, profile))) {
+ rT.reason += ", temp " + currenttemp.rate + " ~ req " + basal + "U/hr. ";
+ return rT;
+ } else {
+ rT.reason += "; setting current basal of " + basal + " as temp. ";
+ return tempBasalFunctions.setTempBasal(basal, 30, profile, rT, currenttemp);
+ }
}
}
- // eventualBG or snoozeBG (from minPredBG) is below max_bg
- if (eventualBG < max_bg || snoozeBG < max_bg) {
+ // eventualBG, snoozeBG, or minPredBG is below max_bg
+ if (Math.min(eventualBG,snoozeBG,minPredBG) < max_bg) {
// if there is a high-temp running and eventualBG > max_bg, let it run
- if (eventualBG > max_bg && round_basal(currenttemp.rate, profile) > round_basal(basal, profile) ) {
- rT.reason += ", " + eventualBG + " > " + max_bg + ": no action required (letting high temp of " + currenttemp.rate + " run)."
+ if (eventualBG > max_bg && round_basal(currenttemp.rate, profile) > round_basal(basal, profile) && currenttemp.duration > 5 ) {
+ rT.reason += eventualBG + " > " + max_bg + ": no temp required (letting high temp of " + currenttemp.rate + " run). "
return rT;
}
- rT.reason += convert_bg(eventualBG, profile)+"-"+convert_bg(snoozeBG, profile)+" in range: no temp required";
- if (currenttemp.duration > 15 && (round_basal(basal, profile) === round_basal(currenttemp.rate, profile))) {
- rT.reason += ", temp " + currenttemp.rate + " ~ req " + basal + "U/hr";
- return rT;
- } else {
- rT.reason += "; setting current basal of " + basal + " as temp";
- return tempBasalFunctions.setTempBasal(basal, 30, profile, rT, currenttemp);
+ // if in SMB mode, don't cancel SMB zero temp
+ if (! (microBolusAllowed && enableSMB )) {
+ rT.reason += convert_bg(eventualBG, profile)+"-"+convert_bg(Math.min(minPredBG,snoozeBG), profile)+" in range: no temp required";
+ if (currenttemp.duration > 15 && (round_basal(basal, profile) === round_basal(currenttemp.rate, profile))) {
+ rT.reason += ", temp " + currenttemp.rate + " ~ req " + basal + "U/hr. ";
+ return rT;
+ } else {
+ rT.reason += "; setting current basal of " + basal + " as temp. ";
+ return tempBasalFunctions.setTempBasal(basal, 30, profile, rT, currenttemp);
+ }
}
}
- // eventual BG is at/above target:
+ // eventual BG is at/above target (or bolus snooze disabled for SMB)
// if iob is over max, just cancel any temps
var basaliob;
if (iob_data.basaliob) { basaliob = iob_data.basaliob; }
else { basaliob = iob_data.iob - iob_data.bolussnooze; }
- rT.reason += "Eventual BG " + convert_bg(eventualBG, profile) + " >= " + convert_bg(max_bg, profile) + ", ";
+ // if we're not here because of SMB, eventual BG is at/above target
+ if (! (microBolusAllowed && rT.COB)) {
+ rT.reason += "Eventual BG " + convert_bg(eventualBG, profile) + " >= " + convert_bg(max_bg, profile) + ", ";
+ }
if (basaliob > max_iob) {
rT.reason += "basaliob " + round(basaliob,2) + " > max_iob " + max_iob;
if (currenttemp.duration > 15 && (round_basal(basal, profile) === round_basal(currenttemp.rate, profile))) {
- rT.reason += ", temp " + currenttemp.rate + " ~ req " + basal + "U/hr";
+ rT.reason += ", temp " + currenttemp.rate + " ~ req " + basal + "U/hr. ";
return rT;
} else {
- rT.reason += "; setting current basal of " + basal + " as temp";
+ rT.reason += "; setting current basal of " + basal + " as temp. ";
return tempBasalFunctions.setTempBasal(basal, 30, profile, rT, currenttemp);
}
} else { // otherwise, calculate 30m high-temp required to get projected BG down to target
- // insulinReq is the additional insulin required to get down to max bg:
- // if in meal assist mode, check if snoozeBG is lower, as eventualBG is not dependent on IOB
- var insulinReq = round( (Math.min(snoozeBG,eventualBG) - target_bg) / sens, 2);
+ // insulinReq is the additional insulin required to get minPredBG down to target_bg
+ //console.error(minPredBG,snoozeBG,eventualBG);
+ var insulinReq = round( (Math.min(minPredBG,snoozeBG,eventualBG) - target_bg) / sens, 2);
+ // when dropping, but not as fast as expected, reduce insulinReq proportionally
+ // to the what fraction of expectedDelta we're dropping at
if (minDelta < 0 && minDelta > expectedDelta) {
var newinsulinReq = round(( insulinReq * (1 - (minDelta / expectedDelta)) ), 2);
//console.error("Reducing insulinReq from " + insulinReq + " to " + newinsulinReq);
@@ -464,8 +785,75 @@ var determine_basal = function determine_basal(glucose_status, currenttemp, iob_
// rate required to deliver insulinReq more insulin over 30m:
var rate = basal + (2 * insulinReq);
rate = round_basal(rate, profile);
+ insulinReq = round(insulinReq,3);
+ rT.insulinReq = insulinReq;
+ rT.minPredBG = minPredBG;
+ //console.error(iob_data.lastBolusTime);
+ // minutes since last bolus
+ var lastBolusAge = round(( new Date().getTime() - iob_data.lastBolusTime ) / 60000,1);
+ //console.error(lastBolusAge);
+ //console.error(profile.temptargetSet, target_bg, rT.COB);
+ // only allow microboluses with COB or low temp targets, or within DIA hours of a bolus
+ // only microbolus if 0.1U SMB represents 20m or less of basal (0.3U/hr or higher)
+ if (microBolusAllowed && enableSMB && profile.current_basal >= 0.3 && bg > threshold) {
+ // never bolus more than 30m worth of basal
+ maxBolus = round(profile.current_basal/2,1);
+ // bolus 1/3 the insulinReq, up to maxBolus
+ microBolus = round(Math.min(insulinReq/3,maxBolus),1);
+
+ // calculate a long enough zero temp to eventually correct back up to target
+ var smbTarget = target_bg;
+ var worstCaseInsulinReq = (smbTarget - (naive_eventualBG + minIOBPredBG)/2 ) / sens;
+ var durationReq = round(60*worstCaseInsulinReq / profile.current_basal);
+
+ // if no microBolus required, snoozeBG > target_bg, and lastCOBpredBG > target_bg, don't set a zero temp
+ if (microBolus < 0.1 && snoozeBG > target_bg && lastCOBpredBG > target_bg) {
+ durationReq = 0;
+ }
-// var maxSafeBasal = Math.min(profile.max_basal, 3 * profile.max_daily_basal, 4 * basal);
+ if (durationReq < 0) {
+ durationReq = 0;
+ // don't set a temp longer than 120 minutes
+ } else {
+ durationReq = round(durationReq/30)*30;
+ durationReq = Math.min(120,Math.max(0,durationReq));
+ }
+ rT.reason += " insulinReq " + insulinReq;
+ if (microBolus >= maxBolus) {
+ rT.reason += "; maxBolus " + maxBolus;
+ }
+ if (durationReq > 0) {
+ rT.reason += "; setting " + durationReq + "m zero temp";
+ }
+ rT.reason += ". ";
+
+ //allow SMBs every 3 minutes
+ var nextBolusMins = round(3-lastBolusAge,1);
+ //console.error(naive_eventualBG, insulinReq, worstCaseInsulinReq, durationReq);
+ console.error("naive_eventualBG",naive_eventualBG+",",durationReq+"m zero temp needed; last bolus",lastBolusAge+"m ago; maxBolus: "+maxBolus);
+ if (lastBolusAge > 3) {
+ if (microBolus > 0) {
+ rT.units = microBolus;
+ rT.reason += "Microbolusing " + microBolus + "U. ";
+ }
+ } else {
+ rT.reason += "Waiting " + nextBolusMins + "m to microbolus again. ";
+ }
+ //rT.reason += ". ";
+
+ // if no zero temp is required, don't return yet; allow later code to set a high temp
+ if (durationReq > 0) {
+ rT.rate = 0;
+ rT.duration = durationReq;
+ return rT;
+ }
+
+ // if insulinReq is negative, snoozeBG > target_bg, and lastCOBpredBG > target_bg, set a neutral temp
+ if (insulinReq < 0 && snoozeBG > target_bg && lastCOBpredBG > target_bg) {
+ rT.reason += "; SMB bolus snooze: setting current basal of " + basal + " as temp. ";
+ return tempBasalFunctions.setTempBasal(basal, 30, profile, rT, currenttemp);
+ }
+ }
var maxSafeBasal = tempBasalFunctions.getMaxSafeBasal(profile);
@@ -476,25 +864,25 @@ var determine_basal = function determine_basal(glucose_status, currenttemp, iob_
var insulinScheduled = currenttemp.duration * (currenttemp.rate - basal) / 60;
if (insulinScheduled >= insulinReq * 2) { // if current temp would deliver >2x more than the required insulin, lower the rate
- rT.reason += currenttemp.duration + "m@" + (currenttemp.rate - basal).toFixed(3) + " = " + insulinScheduled.toFixed(3) + " > 2 * req " + insulinReq + ". Setting temp basal of " + rate + "U/hr";
+ rT.reason += currenttemp.duration + "m@" + (currenttemp.rate).toFixed(2) + " > 2 * insulinReq. Setting temp basal of " + rate + "U/hr. ";
return tempBasalFunctions.setTempBasal(rate, 30, profile, rT, currenttemp);
}
if (typeof currenttemp.duration == 'undefined' || currenttemp.duration == 0) { // no temp is set
- rT.reason += "no temp, setting " + rate + "U/hr";
+ rT.reason += "no temp, setting " + rate + "U/hr. ";
return tempBasalFunctions.setTempBasal(rate, 30, profile, rT, currenttemp);
}
if (currenttemp.duration > 5 && (round_basal(rate, profile) <= round_basal(currenttemp.rate, profile))) { // if required temp <~ existing temp basal
- rT.reason += "temp " + currenttemp.rate + " >~ req " + rate + "U/hr";
+ rT.reason += "temp " + currenttemp.rate + " >~ req " + rate + "U/hr. ";
return rT;
}
// required temp > existing temp basal
- rT.reason += "temp " + currenttemp.rate + "<" + rate + "U/hr";
+ rT.reason += "temp " + currenttemp.rate + "<" + rate + "U/hr. ";
return tempBasalFunctions.setTempBasal(rate, 30, profile, rT, currenttemp);
}
};
-module.exports = determine_basal;
+module.exports = determine_basal;
\ No newline at end of file
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment