Add a 'Reverse' action to the compass heading dive menu.
This allows the compass heading to be reversed for the return leg of the
dive - as a bonus the reversal will be logged, establishing at which
point in the log the dive was turned (mikeller)
line source
///////////////////////////////////////////////////////////////////////////////+ −
/// -*- coding: UTF-8 -*-+ −
///+ −
/// \file Discovery/Src/simulation.c+ −
/// \brief Contains dive simulation functionality+ −
/// \author Heinrichs Weikamp gmbh+ −
/// \date 13-Oct-2014+ −
///+ −
/// \details+ −
/// The simulation uses "extern SDiveState stateSim" defined in dataCentral.h"+ −
///+ −
/// simulation_start(void) sets stateUsed to stateSim and initializes simulation+ −
/// simulation_UpdateLifeData should be called at least once per second+ −
/// simulation_end() sets stateUsed back to stateReal+ −
///+ −
/// $Id$+ −
///////////////////////////////////////////////////////////////////////////////+ −
/// \par Copyright (c) 2014-2018 Heinrichs Weikamp gmbh+ −
///+ −
/// This program is free software: you can redistribute it and/or modify+ −
/// it under the terms of the GNU General Public License as published by+ −
/// the Free Software Foundation, either version 3 of the License, or+ −
/// (at your option) any later version.+ −
///+ −
/// This program is distributed in the hope that it will be useful,+ −
/// but WITHOUT ANY WARRANTY; without even the implied warranty of+ −
/// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the+ −
/// GNU General Public License for more details.+ −
///+ −
/// You should have received a copy of the GNU General Public License+ −
/// along with this program. If not, see <http://www.gnu.org/licenses/>.+ −
//////////////////////////////////////////////////////////////////////////////+ −
+ −
#include <string.h>+ −
#include "simulation.h"+ −
+ −
#include "decom.h"+ −
#include "calc_crush.h"+ −
#include "data_exchange.h"+ −
#include "data_exchange_main.h"+ −
#include "timer.h"+ −
#include "check_warning.h"+ −
#include "vpm.h"+ −
#include "buehlmann.h"+ −
#include "logbook_miniLive.h"+ −
+ −
#include "configuration.h"+ −
+ −
//Private state variables+ −
static float sim_aim_depth_meter;+ −
static float sim_aim_time_minutes;+ −
static _Bool sim_heed_decostops = 1;+ −
+ −
static float sim_descent_rate_meter_per_min = 20;+ −
+ −
static uint16_t* pReplayData; /* pointer to source dive data */+ −
static uint8_t simReplayActive = 0;+ −
+ −
static uint16_t simScrubberTimeoutCount = 0;+ −
+ −
+ −
//Private functions+ −
static float sim_get_ambient_pressure(SDiveState * pDiveState);+ −
static void sim_reduce_deco_time_one_second(SDiveState* pDiveState);+ −
static void simulation_set_aim_depth(int depth_meter);+ −
+ −
#define NUM_OF_SENSORS (3u)+ −
#define SIM_PPO2_STEP (1.1f)+ −
static float simSensmVOffset[NUM_OF_SENSORS];+ −
+ −
/**+ −
******************************************************************************+ −
* @brief sets heed_decostops_while_ascending+ −
******************************************************************************+ −
* @param heed_decostops_while_ascending : true -> deco_stops are considered while ascending+ −
* @return void+ −
*/+ −
void simulation_set_heed_decostops(_Bool heed_decostops_while_ascending)+ −
{+ −
sim_heed_decostops = heed_decostops_while_ascending;+ −
}+ −
+ −
/**+ −
******************************************************************************+ −
* @brief start of simulation+ −
******************************************************************************+ −
* @return void+ −
*/+ −
void simulation_start(int aim_depth, uint16_t aim_time_minutes)+ −
{+ −
uint16_t replayDataLength = 0;+ −
uint8_t* pReplayMarker;+ −
uint16_t max_depth = 10;+ −
uint16_t diveMinutes = 0;+ −
+ −
copyDiveSettingsToSim();+ −
copyVpmRepetetiveDataToSim();+ −
+ −
//vpm_init(&stateSimGetPointerWrite()->vpm, stateSimGetPointerWrite()->diveSettings.vpm_conservatism, 0, 0);+ −
stateSimGetPointerWrite()->lifeData.counterSecondsShallowDepth = 0;+ −
stateSimGetPointerWrite()->mode = MODE_DIVE;+ −
if(aim_depth <= 0)+ −
aim_depth = 20;+ −
sim_descent_rate_meter_per_min = 20;+ −
simulation_set_aim_depth(aim_depth);+ −
sim_aim_time_minutes = aim_time_minutes;+ −
timer_init();+ −
set_stateUsedToSim();+ −
stateSim.lifeData.boolResetAverageDepth = 1;+ −
decoLock = DECO_CALC_init_as_is_start_of_dive;+ −
+ −
stateSim.lifeData.apnea_total_max_depth_meter = 0;+ −
+ −
memcpy(stateSim.scrubberDataDive, settingsGetPointer()->scrubberData, sizeof(stateSim.scrubberDataDive));+ −
memset(simSensmVOffset,0,sizeof(simSensmVOffset));+ −
if(getReplayOffset() != 0xFFFF)+ −
{+ −
simReplayActive = 1;+ −
getReplayInfo(&pReplayData, &pReplayMarker, &replayDataLength, &max_depth, &diveMinutes);+ −
}+ −
}+ −
+ −
/**+ −
******************************************************************************+ −
* @brief end of simulation+ −
******************************************************************************+ −
*+ −
* @return void+ −
*/+ −
void simulation_exit(void)+ −
{+ −
timer_Stopwatch_Stop();+ −
+ −
disableTimer();+ −
+ −
set_stateUsedToReal();+ −
}+ −
+ −
/**+ −
******************************************************************************+ −
* @brief simulates change of Lifedata (saturation, depth change, etc.) within one second+ −
******************************************************************************+ −
*+ −
* @param checkOncePerSecond : true -> simulation in real time (function is evaluated only once per second)+ −
* and copy of parts of LifeData from SmallCPU with each call from HAL_TIM_PeriodElapsedCallback()+ −
* : false -> fast simulation (many simulation cycles per second are possible)+ −
* @return void+ −
*/+ −
void simulation_UpdateLifeData( _Bool checkOncePerSecond)+ −
{+ −
SDiveState * pDiveState = &stateSim;+ −
const SDiveState * pRealState = stateRealGetPointer();+ −
SSettings *pSettings;+ −
+ −
static int last_second = -1;+ −
static _Bool two_second = 0;+ −
static float lastPressure_bar = 0;+ −
+ −
pSettings = settingsGetPointer();+ −
+ −
if ((sim_aim_time_minutes && sim_aim_time_minutes * 60 <= pDiveState->lifeData.dive_time_seconds)+ −
&& (!simReplayActive))+ −
{+ −
simulation_set_aim_depth(0);+ −
}+ −
+ −
float localCalibCoeff[3] = { 0.0, 0.0, 0.0 };+ −
uint8_t index, index2;+ −
+ −
if(checkOncePerSecond)+ −
{+ −
int now = current_second();+ −
if( last_second == now)+ −
return;+ −
last_second = now;+ −
+ −
if(!two_second)+ −
two_second = 1;+ −
else+ −
{+ −
two_second = 0;+ −
}+ −
+ −
for(index = 0; index < 3; index++)+ −
{+ −
localCalibCoeff[index] = pSettings->ppo2sensors_calibCoeff[index];+ −
if(localCalibCoeff[index] < 0.01)+ −
{+ −
for(index2 = 0; index2 < 3; index2++) /* no valid coeff => check other entries */+ −
{+ −
if(pSettings->ppo2sensors_calibCoeff[index2] > 0.01)+ −
{+ −
localCalibCoeff[index] = pSettings->ppo2sensors_calibCoeff[index2];+ −
break;+ −
}+ −
if(index2 == 3) /* no coeff at all => use default */+ −
{+ −
localCalibCoeff[index] = 0.02;+ −
}+ −
}+ −
}+ −
}+ −
+ −
pDiveState->lifeData.temperature_celsius = pRealState->lifeData.temperature_celsius;+ −
pDiveState->lifeData.battery_charge = pRealState->lifeData.battery_charge;+ −
pDiveState->lifeData.compass_heading = pRealState->lifeData.compass_heading;+ −
pDiveState->lifeData.compass_roll = pRealState->lifeData.compass_roll;+ −
pDiveState->lifeData.compass_pitch = pRealState->lifeData.compass_pitch;+ −
+ −
for(index = 0; index < 3; index++)+ −
{+ −
memcpy(&pDiveState->lifeData.extIf_sensor_data[index], &pRealState->lifeData.extIf_sensor_data[index], 32);+ −
}+ −
+ −
#ifdef ENABLE_BOTTLE_SENSOR+ −
pDiveState->lifeData.bottle_bar[pDiveState->lifeData.actualGas.GasIdInSettings] = pRealState->lifeData.bottle_bar[pRealState->lifeData.actualGas.GasIdInSettings];+ −
pDiveState->lifeData.bottle_bar_age_MilliSeconds[pDiveState->lifeData.actualGas.GasIdInSettings] = pRealState->lifeData.bottle_bar_age_MilliSeconds[pRealState->lifeData.actualGas.GasIdInSettings];+ −
#endif+ −
}+ −
else if(pDiveState->lifeData.depth_meter <= (float)(decom_get_actual_deco_stop(pDiveState) + 0.001))+ −
{+ −
if(decoLock == DECO_CALC_FINSHED_vpm)+ −
{+ −
sim_reduce_deco_time_one_second(&stateDeco);+ −
}+ −
else+ −
{+ −
sim_reduce_deco_time_one_second(pDiveState);+ −
}+ −
}+ −
+ −
pDiveState->lifeData.dive_time_seconds += 1;+ −
pDiveState->lifeData.pressure_ambient_bar = sim_get_ambient_pressure(pDiveState);+ −
if(pDiveState->lifeData.depth_meter < 1.5)+ −
{+ −
lastPressure_bar = 0;+ −
pDiveState->lifeData.ascent_rate_meter_per_min = 0;+ −
}+ −
+ −
if((pSettings->scrubTimerMode != SCRUB_TIMER_OFF) && (isLoopMode(pSettings->dive_mode)) && (pDiveState->mode == MODE_DIVE) && isLoopMode(pDiveState->diveSettings.diveMode))+ −
{+ −
simScrubberTimeoutCount++;+ −
if(simScrubberTimeoutCount >= 60) /* resolution is minutes */+ −
{+ −
simScrubberTimeoutCount = 0;+ −
if(pDiveState->scrubberDataDive[pSettings->scubberActiveId].TimerCur > MIN_SCRUBBER_TIME)+ −
{+ −
pDiveState->scrubberDataDive[pSettings->scubberActiveId].TimerCur--;+ −
}+ −
translateDate(stateUsed->lifeData.dateBinaryFormat, &stateUsedWrite->scrubberDataDive[pSettings->scubberActiveId].lastDive);+ −
}+ −
}+ −
+ −
+ −
if(lastPressure_bar > 0)+ −
{+ −
//1 second * 60 == 1 minute, bar * 10 = meter+ −
pDiveState->lifeData.ascent_rate_meter_per_min = (lastPressure_bar - pDiveState->lifeData.pressure_ambient_bar) * 600.0;+ −
}+ −
lastPressure_bar = pDiveState->lifeData.pressure_ambient_bar;+ −
+ −
pDiveState->lifeData.sensorVoltage_mV[0] = pRealState->lifeData.sensorVoltage_mV[0] + simSensmVOffset[0];+ −
if(pDiveState->lifeData.sensorVoltage_mV[0] < 0.0) { pDiveState->lifeData.sensorVoltage_mV[0] = 0.0; }+ −
pDiveState->lifeData.sensorVoltage_mV[1] = pRealState->lifeData.sensorVoltage_mV[1] + simSensmVOffset[1];+ −
if(pDiveState->lifeData.sensorVoltage_mV[1] < 0.0) { pDiveState->lifeData.sensorVoltage_mV[1] = 0.0; }+ −
pDiveState->lifeData.sensorVoltage_mV[2] = pRealState->lifeData.sensorVoltage_mV[2] + simSensmVOffset[2];+ −
if(pDiveState->lifeData.sensorVoltage_mV[2] < 0.0) { pDiveState->lifeData.sensorVoltage_mV[2] = 0.0; }+ −
+ −
pDiveState->lifeData.ppO2Sensor_bar[0] = pDiveState->lifeData.sensorVoltage_mV[0] * localCalibCoeff[0] * pDiveState->lifeData.pressure_ambient_bar;+ −
pDiveState->lifeData.ppO2Sensor_bar[1] = pDiveState->lifeData.sensorVoltage_mV[1] * localCalibCoeff[1] * pDiveState->lifeData.pressure_ambient_bar;+ −
pDiveState->lifeData.ppO2Sensor_bar[2] = pDiveState->lifeData.sensorVoltage_mV[2] * localCalibCoeff[2] * pDiveState->lifeData.pressure_ambient_bar;+ −
+ −
pDiveState->lifeData.CO2_data.CO2_ppm = pRealState->lifeData.CO2_data.CO2_ppm;+ −
+ −
if(is_ambient_pressure_close_to_surface(&pDiveState->lifeData)) // new hw 170214+ −
{+ −
if(!(stateSimGetPointer()->lifeData.counterSecondsShallowDepth))+ −
{+ −
if(pDiveState->diveSettings.diveMode != DIVEMODE_Apnea)+ −
pDiveState->lifeData.counterSecondsShallowDepth = settingsGetPointer()->timeoutDiveReachedZeroDepth - 15;+ −
else+ −
{+ −
pDiveState->lifeData.apnea_last_dive_time_seconds = pDiveState->lifeData.dive_time_seconds;+ −
if(pDiveState->lifeData.apnea_last_dive_time_seconds > pDiveState->lifeData.dive_time_seconds_without_surface_time)+ −
pDiveState->lifeData.apnea_last_dive_time_seconds = pDiveState->lifeData.dive_time_seconds_without_surface_time;+ −
pDiveState->lifeData.apnea_last_max_depth_meter = pDiveState->lifeData.max_depth_meter;+ −
pDiveState->lifeData.counterSecondsShallowDepth = 1;+ −
}+ −
}+ −
}+ −
else+ −
{+ −
pDiveState->lifeData.counterSecondsShallowDepth = 0;+ −
}+ −
+ −
if(!is_ambient_pressure_close_to_surface(&pDiveState->lifeData) && !(stateSimGetPointer()->lifeData.counterSecondsShallowDepth) )+ −
{+ −
pDiveState->lifeData.dive_time_seconds_without_surface_time += 1;+ −
}+ −
+ −
pDiveState->lifeData.depth_meter = (pDiveState->lifeData.pressure_ambient_bar - pDiveState->lifeData.pressure_surface_bar) * 10.0f;+ −
if(pDiveState->lifeData.max_depth_meter < pDiveState->lifeData.depth_meter)+ −
pDiveState->lifeData.max_depth_meter = pDiveState->lifeData.depth_meter;+ −
+ −
/* apnoe specials+ −
*/+ −
if(pDiveState->diveSettings.diveMode == DIVEMODE_Apnea)+ −
{+ −
if(pDiveState->lifeData.max_depth_meter > pDiveState->lifeData.apnea_total_max_depth_meter)+ −
pDiveState->lifeData.apnea_total_max_depth_meter = pDiveState->lifeData.max_depth_meter;+ −
+ −
if(pDiveState->lifeData.counterSecondsShallowDepth)+ −
{+ −
pDiveState->lifeData.dive_time_seconds = 0;+ −
pDiveState->lifeData.max_depth_meter = 0;+ −
pDiveState->lifeData.boolResetAverageDepth = 1;+ −
}+ −
}+ −
+ −
setAvgDepth(pDiveState);+ −
+ −
/* Exposure Tissues+ −
*/+ −
decom_tissues_exposure(1, &pDiveState->lifeData);+ −
decom_oxygen_calculate_cns_exposure(1, &pDiveState->lifeData.actualGas, pDiveState->lifeData.pressure_ambient_bar, &pDiveState->lifeData.cns);+ −
+ −
if(stateSimGetPointer()->lifeData.counterSecondsShallowDepth)+ −
{+ −
stateSimGetPointerWrite()->lifeData.counterSecondsShallowDepth += 1;+ −
if(stateSimGetPointer()->lifeData.counterSecondsShallowDepth >= settingsGetPointer()->timeoutDiveReachedZeroDepth)+ −
simulation_exit();+ −
}+ −
vpm_crush(pDiveState);+ −
}+ −
+ −
/**+ −
******************************************************************************+ −
* @brief adds extra time for fast simulation+ −
******************************************************************************+ −
*@param minutes+ −
* @return float : new pressure+ −
*/+ −
static void simulation_add_time(int minutes)+ −
{+ −
for(int i = 0; i < 60 * minutes; i++)+ −
{+ −
simulation_UpdateLifeData(0);+ −
updateMiniLiveLogbook(0);+ −
timer_UpdateSecond(0);+ −
}+ −
}+ −
+ −
/**+ −
******************************************************************************+ −
* @brief get aim_depth+ −
******************************************************************************+ −
* @return sim_aim_depth_meter;+ −
*/+ −
+ −
uint16_t simulation_get_aim_depth(void)+ −
{+ −
return (uint16_t)sim_aim_depth_meter;+ −
}+ −
+ −
/**+ −
******************************************************************************+ −
* @brief get heed decostops+ −
******************************************************************************+ −
* @return true if ascend follows decostops;+ −
*/+ −
+ −
_Bool simulation_get_heed_decostops(void)+ −
{+ −
return sim_heed_decostops;+ −
}+ −
+ −
/**+ −
******************************************************************************+ −
* @brief sets aim_depth+ −
******************************************************************************+ −
*@param depth_meter+ −
* @return float : new pressure+ −
*/+ −
static void simulation_set_aim_depth(int depth_meter)+ −
{+ −
sim_aim_depth_meter = depth_meter;+ −
}+ −
+ −
/**+ −
******************************************************************************+ −
* @brief simulates ambient pressure depending on aim depth+ −
******************************************************************************+ −
* @note if aim_depth != actual depth, the depth change within one second+ −
* (depending on descent or ascent) rate is calculated+ −
* @param SDiveState* pDiveState:+ −
* @return float : new ambient pressure+ −
*/+ −
static float sim_get_ambient_pressure(SDiveState * pDiveState)+ −
{+ −
//Calc next depth+ −
uint8_t actual_deco_stop = decom_get_actual_deco_stop(pDiveState);+ −
float depth_meter = pDiveState->lifeData.depth_meter;+ −
float surface_pressure_bar = pDiveState->lifeData.pressure_surface_bar;+ −
static uint8_t sampleToggle = 0;+ −
static float sim_ascent_rate_meter_per_min_local = 0;+ −
uint8_t sampleTime = getReplayDataResolution();+ −
+ −
if(simReplayActive) /* precondition: function is called once per second, sample rate is a multiple of second */+ −
{+ −
if(sampleToggle == 0)+ −
{+ −
sampleToggle = sampleTime - 1;+ −
sim_aim_depth_meter = (float)(*pReplayData++/100.0);+ −
if(sim_aim_depth_meter > depth_meter)+ −
{+ −
sim_descent_rate_meter_per_min = (sim_aim_depth_meter - depth_meter) * (60 / sampleTime);+ −
}+ −
else+ −
{+ −
sim_ascent_rate_meter_per_min_local = (depth_meter - sim_aim_depth_meter) * (60 / sampleTime);+ −
}+ −
}+ −
else+ −
{+ −
sampleToggle--;+ −
}+ −
}+ −
else+ −
{+ −
sim_ascent_rate_meter_per_min_local = pDiveState->diveSettings.ascentRate_meterperminute;+ −
}+ −
+ −
if(depth_meter < sim_aim_depth_meter)+ −
{+ −
depth_meter = depth_meter + sim_descent_rate_meter_per_min / 60;+ −
if(depth_meter > sim_aim_depth_meter)+ −
depth_meter = sim_aim_depth_meter;+ −
}+ −
else if(depth_meter > sim_aim_depth_meter)+ −
{+ −
+ −
depth_meter -= sim_ascent_rate_meter_per_min_local / 60;+ −
if(depth_meter < sim_aim_depth_meter)+ −
depth_meter = sim_aim_depth_meter;+ −
+ −
if(sim_heed_decostops && depth_meter < actual_deco_stop)+ −
{+ −
if(actual_deco_stop < (depth_meter + sim_ascent_rate_meter_per_min_local / 60))+ −
depth_meter = actual_deco_stop;+ −
else+ −
depth_meter += sim_ascent_rate_meter_per_min_local / 60;+ −
}+ −
+ −
}+ −
+ −
return surface_pressure_bar + depth_meter / 10;+ −
}+ −
+ −
+ −
/**+ −
******************************************************************************+ −
* @brief Reduces deco time of deepest stop by one second+ −
******************************************************************************+ −
* @note called during fast simulation+ −
* @param SDiveState* pDiveState:+ −
* @return void+ −
*/+ −
static void sim_reduce_deco_time_one_second(SDiveState* pDiveState)+ −
{+ −
SDecoinfo* pDecoinfo;+ −
int8_t index = 0;+ −
+ −
+ −
if(pDiveState->diveSettings.deco_type.ub.standard == GF_MODE)+ −
pDecoinfo = &pDiveState->decolistBuehlmann;+ −
else+ −
pDecoinfo = &pDiveState->decolistVPM;+ −
+ −
//Reduce deco time of deepest stop by one second+ −
for(index = DECOINFO_STRUCT_MAX_STOPS -1 ;index >= 0; index--)+ −
{+ −
if(pDecoinfo->output_stop_length_seconds[index] > 0)+ −
{+ −
pDecoinfo->output_stop_length_seconds[index]--;+ −
break;+ −
}+ −
}+ −
/* update TTS */+ −
if(pDecoinfo->output_time_to_surface_seconds)+ −
{+ −
pDecoinfo->output_time_to_surface_seconds--;+ −
}+ −
}+ −
+ −
SDecoinfo* simulation_decoplaner(uint16_t depth_meter, uint16_t intervall_time_minutes, uint16_t dive_time_minutes, uint8_t *gasChangeListDepthGas20x2)+ −
{+ −
uint8_t ptrGasChangeList = 0; // new hw 160704+ −
uint8_t index = 0;+ −
for (int i = 0; i < 40; i++)+ −
gasChangeListDepthGas20x2[i] = 0;+ −
+ −
SDiveState * pDiveState = &stateSim;+ −
copyDiveSettingsToSim();+ −
+ −
/* activate deco calculation for all deco gases */+ −
for(index = 0; index < 1 + (2*NUM_GASES); index++)+ −
{+ −
if(pDiveState->diveSettings.gas[index].note.ub.deco)+ −
{+ −
pDiveState->diveSettings.gas[index].note.ub.decocalc = 1;+ −
}+ −
}+ −
+ −
vpm_init(&pDiveState->vpm, pDiveState->diveSettings.vpm_conservatism, 0, 0);+ −
//buehlmann_init();+ −
//timer_init();+ −
memset(&pDiveState->events,0, sizeof(SEvents));+ −
pDiveState->diveSettings.internal__pressure_first_stop_ambient_bar_as_upper_limit_for_gf_low_otherwise_zero = 0;+ −
//Calc desaturation during intervall (with Air)+ −
setActualGasAir(&pDiveState->lifeData);+ −
if(intervall_time_minutes > 0)+ −
{+ −
decom_tissues_exposure(intervall_time_minutes * 60, &pDiveState->lifeData);+ −
decom_oxygen_calculate_cns_degrade(&pDiveState->lifeData.cns, intervall_time_minutes * 60);+ −
}+ −
+ −
//Switch to first Gas+ −
setActualGasFirst(&pDiveState->lifeData);+ −
+ −
// new hw 160704+ −
if(gasChangeListDepthGas20x2)+ −
{+ −
gasChangeListDepthGas20x2[ptrGasChangeList++] = 0;+ −
gasChangeListDepthGas20x2[ptrGasChangeList++] = pDiveState->lifeData.actualGas.GasIdInSettings;+ −
gasChangeListDepthGas20x2[0] =0; // depth zero+ −
}+ −
+ −
//Going down / descent+ −
simulation_set_aim_depth(depth_meter);+ −
sim_aim_time_minutes = 0;+ −
for(int i = 0; i < 60 * dive_time_minutes; i++)+ −
{+ −
simulation_UpdateLifeData(0);+ −
check_warning2(pDiveState);+ −
if(pDiveState->warnings.betterGas)+ −
{+ −
setActualGas(&pDiveState->lifeData,actualBetterGasId(),pDiveState->lifeData.actualGas.setPoint_cbar);+ −
if(gasChangeListDepthGas20x2 && (pDiveState->diveSettings.diveMode == DIVEMODE_OC))+ −
{+ −
gasChangeListDepthGas20x2[ptrGasChangeList++] = pDiveState->lifeData.depth_meter;+ −
gasChangeListDepthGas20x2[ptrGasChangeList++] = actualBetterGasId();+ −
}+ −
}+ −
}+ −
+ −
decom_CreateGasChangeList(&pDiveState->diveSettings, &pDiveState->lifeData); // was there before and needed for buehlmann_calc_deco and vpm_calc+ −
+ −
// new hw 160704+ −
if(gasChangeListDepthGas20x2 && (pDiveState->diveSettings.diveMode == DIVEMODE_OC))+ −
{+ −
// change direction from better gas to deco gas+ −
gasChangeListDepthGas20x2[ptrGasChangeList++] = 255;+ −
gasChangeListDepthGas20x2[ptrGasChangeList++] = 255;+ −
+ −
// ascend (deco) gases+ −
for(int i=1; i<=5;i++)+ −
{+ −
if((pDiveState->diveSettings.decogaslist[i].change_during_ascent_depth_meter_otherwise_zero == 0)+ −
|| (pDiveState->diveSettings.gas[pDiveState->diveSettings.decogaslist[i].GasIdInSettings].note.ub.decocalc == 0))+ −
break;+ −
gasChangeListDepthGas20x2[ptrGasChangeList++] = pDiveState->diveSettings.decogaslist[i].change_during_ascent_depth_meter_otherwise_zero;+ −
gasChangeListDepthGas20x2[ptrGasChangeList++] = pDiveState->diveSettings.decogaslist[i].GasIdInSettings;+ −
}+ −
gasChangeListDepthGas20x2[0] = 0;+ −
}+ −
+ −
// deco and ascend calc+ −
if(pDiveState->diveSettings.deco_type.ub.standard == GF_MODE)+ −
{+ −
/* this does modify the cns now 11.06.2015 */+ −
buehlmann_calc_deco(&pDiveState->lifeData,&pDiveState->diveSettings,&pDiveState->decolistBuehlmann);+ −
pDiveState->lifeData.cns += buehlmann_get_gCNS();+ −
return &pDiveState->decolistBuehlmann;+ −
}+ −
else+ −
{+ −
/* this does modify the cns now 11.06.2015 */+ −
vpm_calc(&pDiveState->lifeData,&pDiveState->diveSettings,&pDiveState->vpm,&pDiveState->decolistVPM, DECOSTOPS);+ −
pDiveState->lifeData.cns += vpm_get_CNS();+ −
+ −
while(decoLock == DECO_CALC_FINSHED_vpm)+ −
{+ −
HAL_Delay(2); /* The deco data is copied during the timer ISR => wait till this has happened */+ −
}+ −
return &pDiveState->decolistVPM;+ −
}+ −
}+ −
+ −
static float sGChelper_bar(uint16_t depth_meter)+ −
{+ −
SDiveState * pDiveState = &stateSim;+ −
float ambient, surface, density, meter;+ −
+ −
surface = pDiveState->lifeData.pressure_surface_bar;+ −
+ −
if(!depth_meter)+ −
return surface;+ −
+ −
density = ((float)( 100 + settingsGetPointer()->salinity)) / 100.0f;+ −
meter = depth_meter * (0.09807f * density);+ −
ambient = (meter + surface);+ −
+ −
return ambient;+ −
}+ −
+ −
+ −
/**+ −
******************************************************************************+ −
* @brief simulation_helper_change_points+ −
******************************************************************************+ −
* @param+ −
* @return void+ −
*/+ −
void simulation_helper_change_points(SSimDataSummary *outputSummary, uint16_t depth_meter, uint16_t dive_time_minutes, SDecoinfo *decoInfoInput, const uint8_t *gasChangeListDepthGas20x2)+ −
{+ −
uint8_t ptrDecoInfo = 0;+ −
uint16_t actualDepthPoint = 0;+ −
uint16_t nextDepthPoint = 0;+ −
uint8_t actualConsumGasId = 0;+ −
uint8_t nextGasChangeMeter = 0;+ −
uint8_t ptrChangeList = 0;+ −
+ −
float timeThis = 0;+ −
float timeSummary = 0;+ −
float sim_descent_rate_meter_per_min_local = 10;+ −
float sim_ascent_rate_meter_per_min_local = 10;+ −
+ −
SDiveState * pDiveState = &stateSim;+ −
+ −
uint8_t depthDecoNext, depthLast, depthSecond, depthInc;+ −
+ −
if(pDiveState->diveSettings.deco_type.ub.standard == GF_MODE)+ −
{+ −
sim_descent_rate_meter_per_min_local = sim_descent_rate_meter_per_min; // const float+ −
sim_ascent_rate_meter_per_min_local = pDiveState->diveSettings.ascentRate_meterperminute;+ −
}+ −
else+ −
{+ −
sim_descent_rate_meter_per_min_local = sim_descent_rate_meter_per_min; // const float+ −
sim_ascent_rate_meter_per_min_local = 10;// fix in vpm_calc_deco();+ −
}+ −
+ −
outputSummary->descentRateMeterPerMinute = sim_descent_rate_meter_per_min_local;+ −
outputSummary->ascentRateMeterPerMinute = sim_ascent_rate_meter_per_min_local;+ −
+ −
// bottom gas ppO2+ −
if(gasChangeListDepthGas20x2)+ −
{+ −
nextGasChangeMeter = gasChangeListDepthGas20x2[ptrChangeList++];+ −
actualConsumGasId = gasChangeListDepthGas20x2[ptrChangeList++];+ −
nextGasChangeMeter = gasChangeListDepthGas20x2[ptrChangeList++];+ −
+ −
while(actualDepthPoint < depth_meter)+ −
{+ −
if(nextGasChangeMeter && (nextGasChangeMeter < depth_meter) && (gasChangeListDepthGas20x2[ptrChangeList] != 255)) // list has 255,255 for turn from travel to deco+ −
{+ −
nextDepthPoint = nextGasChangeMeter;+ −
}+ −
else+ −
{+ −
nextDepthPoint = depth_meter;+ −
}+ −
+ −
if(actualConsumGasId > 5) // safety first+ −
actualConsumGasId = 0;+ −
+ −
actualDepthPoint = nextDepthPoint;+ −
+ −
if(actualDepthPoint != depth_meter)+ −
{+ −
actualConsumGasId = gasChangeListDepthGas20x2[ptrChangeList++];+ −
nextGasChangeMeter = gasChangeListDepthGas20x2[ptrChangeList++];+ −
}+ −
}+ −
}+ −
else+ −
{+ −
actualConsumGasId = pDiveState->lifeData.actualGas.GasIdInSettings;+ −
nextGasChangeMeter = 0;+ −
}+ −
outputSummary->ppO2AtBottom = (sGChelper_bar(depth_meter) - WATER_VAPOUR_PRESSURE) * pDiveState->diveSettings.gas[actualConsumGasId].oxygen_percentage / 100.0f;+ −
+ −
+ −
// going down+ −
actualDepthPoint = 0;+ −
nextDepthPoint = depth_meter;+ −
+ −
timeThis = ((float)(nextDepthPoint - actualDepthPoint)) / sim_descent_rate_meter_per_min_local;+ −
timeSummary += timeThis;+ −
outputSummary->timeToBottom = (uint16_t)timeThis;+ −
+ −
// bottom time+ −
timeThis = ((float)dive_time_minutes) - timeSummary;+ −
timeSummary += timeThis;+ −
outputSummary->timeAtBottom = (uint16_t)timeSummary;+ −
+ −
+ −
// ascend to first deco stop+ −
actualDepthPoint = depth_meter; // that is where we are+ −
timeThis = 0;+ −
+ −
if(!decoInfoInput->output_stop_length_seconds[0]) // NDL dive+ −
{+ −
depthLast = 0;+ −
ptrDecoInfo = 0;+ −
depthDecoNext = 0;+ −
}+ −
else+ −
{+ −
// prepare deco stop list+ −
depthLast = (uint8_t)(stateUsed->diveSettings.last_stop_depth_bar * 10);+ −
depthSecond = (uint8_t)(stateUsed->diveSettings.input_second_to_last_stop_depth_bar * 10);+ −
depthInc = (uint8_t)(stateUsed->diveSettings.input_next_stop_increment_depth_bar * 10);+ −
+ −
for(ptrDecoInfo=DECOINFO_STRUCT_MAX_STOPS-1; ptrDecoInfo>0; ptrDecoInfo--)+ −
if(decoInfoInput->output_stop_length_seconds[ptrDecoInfo]) break;+ −
+ −
if(ptrDecoInfo == 0)+ −
{+ −
depthDecoNext = depthLast;+ −
}+ −
else+ −
depthDecoNext = depthSecond + (( ptrDecoInfo - 1 )* depthInc);+ −
}+ −
+ −
nextDepthPoint = depthDecoNext;+ −
if(actualDepthPoint > nextDepthPoint)+ −
{+ −
// flip signs! It's going up+ −
timeThis = ((float)(actualDepthPoint - nextDepthPoint)) / sim_ascent_rate_meter_per_min_local;+ −
actualDepthPoint = nextDepthPoint; // that is where we are+ −
}+ −
timeSummary += timeThis;+ −
outputSummary->timeToFirstStop = (uint16_t)timeSummary;+ −
outputSummary->depthMeterFirstStop = actualDepthPoint;+ −
+ −
if(decoInfoInput->output_time_to_surface_seconds)+ −
{+ −
outputSummary->timeToSurface = outputSummary->timeAtBottom + (decoInfoInput->output_time_to_surface_seconds / 60);+ −
}+ −
else+ −
{+ −
outputSummary->timeToSurface = outputSummary->timeToFirstStop;+ −
}+ −
}+ −
+ −
+ −
/**+ −
******************************************************************************+ −
* @brief simulation_gas_consumption+ −
******************************************************************************+ −
* @note called by openEdit_PlanResult() in tMenuEditPlanner.c+ −
* @note the ascend and descend time is taken from pDiveState->lifeData.ascent_rate_meter_per_min and const float sim_descent_rate_meter_per_min+ −
* @param outputConsumptionList list from 1 to 5 for gas 1 to 5+ −
* @param depth_meter for descend+ −
* @param dive_time_minutes for descend and bottom time+ −
* @param the calculated deco list+ −
* @param gasConsumTravelInput: how many l/min for all but deco stops+ −
* @param gasConsumDecoInput: how many l/min for deco stops only+ −
* @return void+ −
*/+ −
+ −
void simulation_gas_consumption(uint16_t *outputConsumptionList, uint16_t depth_meter, uint16_t dive_time_minutes, SDecoinfo *decoInfoInput, uint8_t gasConsumTravelInput, uint8_t gasConsumDecoInput, const uint8_t *gasChangeListDepthGas20x2)+ −
{+ −
uint8_t ptrDecoInfo = 0;+ −
uint8_t ptrChangeList = 0;+ −
uint8_t actualConsumGasId = 0;+ −
uint8_t nextGasChangeMeter = 0;+ −
uint16_t actualDepthPoint = 0;+ −
uint16_t nextDepthPoint = 0;+ −
uint16_t inBetweenDepthPoint = 0;+ −
float timeThis = 0;+ −
float consumThis = 0;+ −
float timeSummary = 0;+ −
float outputConsumptionTempFloat[6];+ −
float sim_descent_rate_meter_per_min_local = 10;+ −
float sim_ascent_rate_meter_per_min_local = 10;+ −
+ −
SDiveState * pDiveState = &stateSim;+ −
+ −
uint8_t depthDecoNext = 0;+ −
uint8_t depthLast = 0;+ −
uint8_t depthSecond = 0;+ −
uint8_t depthInc = 0;+ −
+ −
for(int i = 1; i < 6; i++)+ −
outputConsumptionTempFloat[i] = 0;+ −
+ −
if(gasChangeListDepthGas20x2)+ −
{+ −
nextGasChangeMeter = gasChangeListDepthGas20x2[ptrChangeList++];+ −
actualConsumGasId = gasChangeListDepthGas20x2[ptrChangeList++];+ −
nextGasChangeMeter = gasChangeListDepthGas20x2[ptrChangeList++];+ −
}+ −
else+ −
{+ −
actualConsumGasId = pDiveState->lifeData.actualGas.GasIdInSettings;+ −
nextGasChangeMeter = 0;+ −
}+ −
+ −
if(pDiveState->diveSettings.deco_type.ub.standard == GF_MODE)+ −
{+ −
sim_descent_rate_meter_per_min_local = sim_descent_rate_meter_per_min; // const float+ −
sim_ascent_rate_meter_per_min_local = pDiveState->diveSettings.ascentRate_meterperminute;+ −
}+ −
else+ −
{+ −
sim_descent_rate_meter_per_min_local = sim_descent_rate_meter_per_min; // const float+ −
sim_ascent_rate_meter_per_min_local = 10;// fix in vpm_calc_deco();+ −
}+ −
+ −
// while((nextGasChangeMeter < depth_meter) && (actualDepthPoint < depth_meter))+ −
while(actualDepthPoint < depth_meter)+ −
{+ −
if(nextGasChangeMeter && (nextGasChangeMeter < depth_meter) && (gasChangeListDepthGas20x2[ptrChangeList] != 255)) // list has 255,255 for turn from travel to deco+ −
{+ −
nextDepthPoint = nextGasChangeMeter;+ −
}+ −
else+ −
{+ −
nextDepthPoint = depth_meter;+ −
}+ −
+ −
if(actualConsumGasId > 5) // safety first+ −
actualConsumGasId = 0;+ −
+ −
timeThis = ((float)(nextDepthPoint - actualDepthPoint)) / sim_descent_rate_meter_per_min_local;+ −
if(actualDepthPoint) // not if on surface+ −
{+ −
consumThis = ((float)gasConsumTravelInput) * sGChelper_bar(actualDepthPoint) * timeThis;+ −
}+ −
consumThis += ((float)gasConsumTravelInput) * sGChelper_bar(nextDepthPoint -actualDepthPoint) * timeThis / 2;+ −
outputConsumptionTempFloat[actualConsumGasId] += consumThis;+ −
timeSummary += timeThis;+ −
+ −
actualDepthPoint = nextDepthPoint;+ −
+ −
if(actualDepthPoint != depth_meter)+ −
{+ −
actualConsumGasId = gasChangeListDepthGas20x2[ptrChangeList++];+ −
nextGasChangeMeter = gasChangeListDepthGas20x2[ptrChangeList++];+ −
}+ −
}+ −
+ −
// bottom Time+ −
timeThis = ((float)dive_time_minutes) - timeSummary;+ −
+ −
if(timeThis > 0)+ −
{+ −
consumThis = ((float)gasConsumTravelInput) * sGChelper_bar(depth_meter) * timeThis;+ −
outputConsumptionTempFloat[actualConsumGasId] += consumThis;+ −
}+ −
+ −
// ascend with deco stops prepare+ −
if(gasChangeListDepthGas20x2)+ −
{+ −
ptrChangeList++;// gasChangeListDepthGas20x2[ptrChangeList++]; // should be the 255+ −
nextGasChangeMeter = gasChangeListDepthGas20x2[ptrChangeList++];+ −
}+ −
else+ −
{+ −
nextGasChangeMeter = 0;+ −
}+ −
+ −
+ −
if(!decoInfoInput->output_stop_length_seconds[0]) // NDL dive+ −
{+ −
depthLast = 0;+ −
ptrDecoInfo = 0;+ −
}+ −
else+ −
{+ −
// prepare deco stop list+ −
depthLast = (uint8_t)(stateUsed->diveSettings.last_stop_depth_bar * 10);+ −
depthSecond = (uint8_t)(stateUsed->diveSettings.input_second_to_last_stop_depth_bar * 10);+ −
depthInc = (uint8_t)(stateUsed->diveSettings.input_next_stop_increment_depth_bar * 10);+ −
+ −
for(ptrDecoInfo=DECOINFO_STRUCT_MAX_STOPS-1; ptrDecoInfo>0; ptrDecoInfo--)+ −
if(decoInfoInput->output_stop_length_seconds[ptrDecoInfo]) break;+ −
}+ −
+ −
actualDepthPoint = depth_meter; // that is where we are+ −
+ −
// ascend with deco stops+ −
while(actualDepthPoint)+ −
{+ −
if(ptrDecoInfo == 0)+ −
{+ −
depthDecoNext = depthLast;+ −
}+ −
else+ −
depthDecoNext = depthSecond + (( ptrDecoInfo - 1 )* depthInc);+ −
+ −
if(nextGasChangeMeter && (nextGasChangeMeter > depthDecoNext))+ −
{+ −
nextDepthPoint = nextGasChangeMeter;+ −
}+ −
else+ −
{+ −
nextDepthPoint = depthDecoNext;+ −
}+ −
+ −
if(actualConsumGasId > 5) // safety first+ −
actualConsumGasId = 0;+ −
+ −
if(actualDepthPoint > nextDepthPoint)+ −
{+ −
// flip signs! It's going up+ −
timeThis = ((float)(actualDepthPoint - nextDepthPoint)) / sim_ascent_rate_meter_per_min_local;+ −
inBetweenDepthPoint = nextDepthPoint + ((actualDepthPoint - nextDepthPoint)/2);+ −
consumThis = ((float)gasConsumDecoInput) * sGChelper_bar(inBetweenDepthPoint) * timeThis;+ −
/*+ −
if(nextDepthPoint)+ −
{+ −
consumThis = ((float)gasConsumDecoInput) * sGChelper_bar(nextDepthPoint) * timeThis;+ −
}+ −
else+ −
{+ −
consumThis = 0;+ −
}+ −
consumThis += ((float)gasConsumDecoInput) * sGChelper_bar(actualDepthPoint - nextDepthPoint) * timeThis / 2;+ −
*/+ −
outputConsumptionTempFloat[actualConsumGasId] += consumThis;+ −
}+ −
+ −
if(nextGasChangeMeter && (nextDepthPoint == nextGasChangeMeter))+ −
{+ −
actualConsumGasId = gasChangeListDepthGas20x2[ptrChangeList++];+ −
nextGasChangeMeter = gasChangeListDepthGas20x2[ptrChangeList++];+ −
}+ −
+ −
if(actualConsumGasId > 5) // safety first+ −
actualConsumGasId = 0;+ −
+ −
if(nextDepthPoint && (nextDepthPoint == depthDecoNext))+ −
{+ −
if(decoInfoInput->output_stop_length_seconds[ptrDecoInfo])+ −
{+ −
timeThis = ((float)(decoInfoInput->output_stop_length_seconds[ptrDecoInfo])) / 60.0f;+ −
consumThis = ((float)gasConsumDecoInput) * sGChelper_bar(nextDepthPoint) * timeThis;+ −
outputConsumptionTempFloat[actualConsumGasId] += consumThis;+ −
}+ −
if(ptrDecoInfo != 0)+ −
{+ −
ptrDecoInfo--;+ −
}+ −
else+ −
{+ −
depthLast = 0;+ −
}+ −
}+ −
actualDepthPoint = nextDepthPoint;+ −
}+ −
+ −
// copy and return+ −
for(int i = 1; i < 6; i++)+ −
outputConsumptionList[i] = (uint16_t)(outputConsumptionTempFloat[i]);+ −
}+ −
+ −
/**+ −
******************************************************************************+ −
* @brief Simulator control during simulated dive+ −
******************************************************************************+ −
* @note called by user via tHomeDiveMenuControl()+ −
* @param void+ −
* @return void+ −
*/+ −
+ −
+ −
void Sim_Descend (void)+ −
{+ −
stateSimGetPointerWrite()->lifeData.counterSecondsShallowDepth = 0;+ −
if(simulation_get_aim_depth() < 200)+ −
simulation_set_aim_depth(simulation_get_aim_depth() + 1);+ −
}+ −
+ −
+ −
void Sim_Ascend (void)+ −
{+ −
if(simulation_get_aim_depth() > 0)+ −
simulation_set_aim_depth(simulation_get_aim_depth() - 1);+ −
}+ −
+ −
+ −
void Sim_Divetime (void)+ −
{+ −
simulation_add_time(5);+ −
}+ −
+ −
+ −
void Sim_Quit (void)+ −
{+ −
if(stateSimGetPointer()->lifeData.counterSecondsShallowDepth)+ −
{+ −
simulation_exit();+ −
return;+ −
}+ −
+ −
if(simulation_get_aim_depth() > 0)+ −
{+ −
simulation_set_aim_depth(0);+ −
}+ −
else+ −
{+ −
stateSimGetPointerWrite()->lifeData.depth_meter = 0;+ −
if(stateSimGetPointer()->diveSettings.diveMode == DIVEMODE_Apnea)+ −
{+ −
stateSimGetPointerWrite()->lifeData.counterSecondsShallowDepth = 1;+ −
}+ −
else+ −
{+ −
stateSimGetPointerWrite()->lifeData.counterSecondsShallowDepth = settingsGetPointer()->timeoutDiveReachedZeroDepth - 15;+ −
}+ −
}+ −
}+ −
+ −
void Sim_IncreasePPO(uint8_t sensorIdx)+ −
{+ −
if((sensorIdx < NUM_OF_SENSORS) && (simSensmVOffset[sensorIdx] + SIM_PPO2_STEP < 100.0) && ((stateUsed->diveSettings.ppo2sensors_deactivated & (1 << sensorIdx)) == 0))+ −
{+ −
simSensmVOffset[sensorIdx] += SIM_PPO2_STEP;+ −
}+ −
}+ −
void Sim_DecreasePPO(uint8_t sensorIdx)+ −
{+ −
if((sensorIdx < NUM_OF_SENSORS) && (simSensmVOffset[sensorIdx] - SIM_PPO2_STEP >= -100.0))+ −
{+ −
simSensmVOffset[sensorIdx] -= SIM_PPO2_STEP;+ −
}+ −
}+ −