function middleware(iob, currenttemp, glucose, profile, autosens, meal, reservoir, clock, pumphistory, preferences) { // This middleware only works if you have added pumphistory to middleware in FreeAPS X code (my pumphistory branch). const BG = glucose[0].glucose; // Change to false to turn off Chris Wilson's formula var chrisFormula = preferences.enableChris; const minLimitChris = profile.autosens_min; const maxLimitChris = profile.autosens_max; const adjustmentFactor = preferences.adjustmentFactor; const currentMinTarget = profile.min_bg; var exerciseSetting = false; var enoughData = false; var pumpData = 0; var log = ""; var logTDD = ""; var logBasal = ""; var logBolus = ""; var logTempBasal = ""; var current = 0; // If you have not set this to 0.05 in FAX settings (Omnipod), this will be set to 0.1 in code. var minimalDose = profile.bolus_increment; var TDD = 0; var insulin = 0; var tempInsulin = 0; var bolusInsulin = 0; var scheduledBasalInsulin = 0; var incrementsRaw = 0; var incrementsRounded = 0; var quota = 0; if (profile.high_temptarget_raises_sensitivity == true || profile.exercise_mode == true) { exerciseSetting = true; } // Turn off Chris' formula (and AutoISF) when using a temp target >= 118 (6.5 mol/l) and if an exercise setting is enabled. // If using AutoISF uncomment the profile.use_autoisf = false if (currentMinTarget >= 118 && exerciseSetting == true) { // profile.use_autoisf = false; chrisFormula = false; log = "Chris' formula is off due to a high temp target/exercising. Current min target: " + currentMinTarget; } // Check that there is enough pump history data (>23 hours) for TDD calculation, else end this middleware. if (chrisFormula == true) { let ph_length = pumphistory.length; let endDate = new Date(pumphistory[ph_length-1].timestamp); let startDate = new Date(pumphistory[0].timestamp); // > 23 hours pumpData = (startDate - endDate) / 36e5; if (pumpData >= 23) { enoughData = true; } else { chrisFormula = false; return "Chris' formula is temporarily off. 24 hours of data is required for a correct TDD calculation. Currently only " + pumpData.toPrecision(3) + " hours of pump history data available."; } } // Calculate TDD -------------------------------------- //Bolus: for (let i = 0; i < pumphistory.length; i++) { if (pumphistory[i]._type == "Bolus") { bolusInsulin += pumphistory[i].amount; } } // Temp basals: if (minimalDose != 0.05) { minimalDose = 0.1; } for (let j = 1; j < pumphistory.length; j++) { if (pumphistory[j]._type == "TempBasal" && pumphistory[j].rate > 0) { current = j; quota = pumphistory[j].rate; var duration = pumphistory[j-1]['duration (min)'] / 60; var origDur = duration; var pastTime = new Date(pumphistory[j-1].timestamp); // If temp basal hasn't yet ended, use now as end date for calculation do { j--; if (j <= 0) { morePresentTime = new Date(); break; } else if (pumphistory[j]._type == "TempBasal" || pumphistory[j]._type == "PumpSuspend") { morePresentTime = new Date(pumphistory[j].timestamp); break; } } while (j >= 0); var diff = (morePresentTime - pastTime) / 36e5; if (diff < origDur) { duration = diff; } insulin = quota * duration; // Account for smallest possible pump dosage incrementsRaw = insulin / minimalDose; if (incrementsRaw >= 1) { incrementsRounded = Math.floor(incrementsRaw); insulin = incrementsRounded * minimalDose; tempInsulin += insulin; } else { insulin = 0} j = current; } } // Check and count for when basals are delivered with a scheduled basal rate or an Autotuned basal rate. // 1. Check for 0 temp basals with 0 min duration. This is for when ending a manual temp basal and (perhaps) continuing in open loop for a while. // 2. Check for temp basals that completes. This is for when disconected from link/iphone, or when in open loop. // To do: need to check for more circumstances when scheduled basal rates are used. // for (let i = 0; i < pumphistory.length; i++) { // Check for 0 temp basals with 0 min duration. insulin = 0; if (pumphistory[i]['duration (min)'] == 0) { let time1 = new Date(pumphistory[i].timestamp); let time2 = time1; let j = i; do { --j; if (pumphistory[j]._type == "TempBasal" && j >= 0) { time2 = new Date(pumphistory[j].timestamp); break; } } while (j > 0); // duration of current scheduled basal in h let basDuration = (time2 - time1) / 36e5; if (basDuration > 0) { let hour = time1.getHours(); let minutes = time1.getMinutes(); let seconds = "00"; let string = "" + hour + ":" + minutes + ":" + seconds; let baseTime = new Date(string); let basalScheduledRate = profile.basalprofile[0].start; for (let k = 0; k < profile.basalprofile.length; k++) { if (profile.basalprofile[k].start == baseTime) { basalScheduledRate = profile.basalprofile[k].rate; insulin = basalScheduledRate * basDuration; break; } else if (k + 1 < profile.basalprofile.length) { if (profile.basalprofile[k].start < baseTime && profile.basalprofile[k+1].start > baseTime) { basalScheduledRate = profile.basalprofile[k].rate; insulin = basalScheduledRate * basDuration; break; } } else if (k == profile.basalprofile.length - 1) { basalScheduledRate = profile.basalprofile[k].rate; insulin = basalScheduledRate * basDuration; break; } } // Account for smallest possible pump dosage incrementsRaw = insulin / minimalDose; if (incrementsRaw >= 1) { incrementsRounded = Math.floor(incrementsRaw); insulin = incrementsRounded * minimalDose; scheduledBasalInsulin += insulin; } else { insulin = 0} } } } // Check for temp basals that completes for (let i = pumphistory.length -1; i > 0; i--) { if (pumphistory[i]._type == "TempBasalDuration" && pumphistory[i]['duration (min)'] > 0) { let time2Duration = pumphistory[i]['duration (min)'] / 60; let time2 = new Date(pumphistory[i].timestamp);; let time1 = time2; let m = i; do { --m; if (m >= 0) { if (pumphistory[m]._type == "TempBasal" && m >= 0) { // next (newer) temp basal let time1 = new Date(pumphistory[m].timestamp); break; } } } while (m > 0); // Time difference in hours let tempBasalTimeDifference = (time1 - time2) / 36e5; if (time2Duration < tempBasalTimeDifference) { let timeOfbasal = tempBasalTimeDifference - time2Duration; // let hour = time2.getHours(); let minutes = time2.getMinutes(); let seconds = "00"; let string = "" + hour + ":" + minutes + ":" + seconds; let baseTime = new Date(string); // Default if correct basal schedule rate not found let basalScheduledRate = profile.basalprofile[0].rate; for (let k = 0; k < profile.basalprofile.length; ++k) { if (profile.basalprofile[k].start == baseTime) { basalScheduledRate = profile.basalprofile[k].rate; break; } else if (k+1 < profile.basalprofile.length) { if (profile.basalprofile[k].start < baseTime && profile.basalprofile[k+1].start > baseTime) { basalScheduledRate = profile.basalprofile[k].rate; break; } } else if (k == (profile.basalprofile.length - 1)) { basalScheduledRate = profile.basalprofile[k].rate; break; } } insulin = basalScheduledRate * timeOfbasal; // Account for smallest possible pump dosage incrementsRaw = insulin / minimalDose; if (incrementsRaw >= 1) { incrementsRounded = Math.floor(incrementsRaw); scheduledBasalInsulin += incrementsRounded * minimalDose; } else { insulin = 0} } } } TDD = bolusInsulin + tempInsulin + scheduledBasalInsulin; logBolus = ". Bolus insulin: " + bolusInsulin.toPrecision(5) + " U"; logTempBasal = ". Temporary basal insulin: " + tempInsulin.toPrecision(5) + " U"; logBasal = ". Delivered scheduled basal rate insulin: " + scheduledBasalInsulin.toPrecision(5) + " U"; logTDD = ". TDD past 24h is: " + TDD.toPrecision(5) + " U"; // ---------------------------------------------------- // Chris' formula with added adjustmentFactor for tuning: if (chrisFormula == true && TDD > 0) { var newRatio = profile.sens / (277700 / (adjustmentFactor * TDD * BG)); log = "New ratio using Chris' formula is " + newRatio.toPrecision(3) + " with ISF: " + (profile.sens / newRatio).toPrecision(3) + " (" + ((profile.sens / newRatio) * 0.0555).toPrecision(3) + " mmol/l/U)"; // Respect autosens.max and autosens.min limits if (newRatio > maxLimitChris) { newRatio = maxLimitChris; log = "Chris' formula hit limit by autosens_max setting: " + maxLimitChris + ". ISF: " + (profile.sens / newRatio).toPrecision(3) + " (" + ((profile.sens / newRatio) * 0.0555).toPrecision(3) + " mmol/l/U)"; } else if (newRatio < minLimitChris) { newRatio = minLimitChris; log = "Chris' formula hit limit by autosens_min setting: " + minLimitChris + ". ISF: " + (profile.sens / newRatio).toPrecision(3) + " (" + ((profile.sens / newRatio) * 0.0555).toPrecision(3) + " mmol/l/U)"; } // Set the new ratio autosens.ratio = newRatio; // Print to log return log + logTDD + logBolus + logTempBasal + logBasal; // + " TimeDiff: " + timeOfbasal + " baseTime: " + baseTime + " insulin:" + insulin; } else { return "Chris' formula is off." } }