Mercurial > public > ostc4
annotate Discovery/Src/test_vpm.c @ 747:df0d43da1614
Added pressure compensation to CO2 detection:
A pressure compensation is needed if the ExplorIR shall be used under extended pressure conditions. The procedure recommended by the application note has been integrated. To keep things simple the focus of the CO2 measurement is not the precision, as it is needed for the decompression calculation, but the indication of a critical increase of CO2 in the breathing loop. That's why only a lookup table with 1000ppm steps has been implemented instead of calculating the polynom for every measurement.
| author | Ideenmodellierer |
|---|---|
| date | Sun, 05 Mar 2023 22:06:47 +0100 |
| parents | 54d14bc2083c |
| children | aa6006975e76 |
| rev | line source |
|---|---|
| 38 | 1 /////////////////////////////////////////////////////////////////////////////// |
| 2 /// -*- coding: UTF-8 -*- | |
| 3 /// | |
| 4 /// \file Discovery/Src/test_vpm.c | |
| 5 /// \brief test 101 | |
| 6 /// \author Heinrichs Weikamp | |
| 7 /// \date 26-Oct-2014 | |
| 8 /// | |
| 9 /// \details | |
| 10 /// | |
| 11 /// $Id$ | |
| 12 /////////////////////////////////////////////////////////////////////////////// | |
| 13 /// \par Copyright (c) 2014-2018 Heinrichs Weikamp gmbh | |
| 14 /// | |
| 15 /// This program is free software: you can redistribute it and/or modify | |
| 16 /// it under the terms of the GNU General Public License as published by | |
| 17 /// the Free Software Foundation, either version 3 of the License, or | |
| 18 /// (at your option) any later version. | |
| 19 /// | |
| 20 /// This program is distributed in the hope that it will be useful, | |
| 21 /// but WITHOUT ANY WARRANTY; without even the implied warranty of | |
| 22 /// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | |
| 23 /// GNU General Public License for more details. | |
| 24 /// | |
| 25 /// You should have received a copy of the GNU General Public License | |
| 26 /// along with this program. If not, see <http://www.gnu.org/licenses/>. | |
| 27 ////////////////////////////////////////////////////////////////////////////// | |
| 28 | |
| 29 #include <stdio.h> | |
| 30 #include <stdint.h> | |
| 31 //#include "LED.h" | |
| 32 //#include "Keyboard.h" | |
| 33 //#include "stm32f4xx_hal.h" | |
| 34 #include "buehlmann.h" | |
| 35 #include "calc_crush.h" | |
| 36 #include "vpm.h" | |
| 37 #include "display.h" | |
| 38 #include "test_vpm.h" | |
| 39 #include "math.h" | |
| 40 #include "data_central.h" | |
| 41 #include "decom.h" | |
| 42 #include "logbook.h" | |
| 43 #include "tInfoLog.h" | |
| 44 | |
| 45 #define true 1 | |
| 46 #define false 0 | |
| 47 //#define uint8_t unsigned char | |
| 48 | |
| 49 extern SSettings Settings; | |
| 50 | |
| 51 _Bool simulate_descent(SDiveState* pInput, float ending_depth_meter, float rate_meter_per_minutes); | |
| 52 void init_buehlmann(SDiveState* pInput); | |
| 53 | |
| 54 _Bool test1(void); | |
| 55 uint8_t test2_unapproved(void); | |
| 56 uint8_t test3_unapproved(void); | |
| 57 | |
| 58 _Bool simulate_descent(SDiveState* pInput, | |
| 59 float ending_depth_meter, | |
| 60 float rate_meter_per_minutes) | |
| 61 { | |
| 62 int i =0; | |
| 63 static float initial_helium_pressure[16]; | |
| 64 static float initial_nitrogen_pressure[16]; | |
| 65 static float initial_inspired_he_pressure; | |
| 66 static float initial_inspired_n2_pressure; | |
| 67 static float fraction_nitrogen_begin; | |
| 68 static float fraction_nitrogen_end; | |
| 69 static float fraction_helium_begin; | |
| 70 static float fraction_helium_end; | |
| 71 static float nitrogen_rate; | |
| 72 static float helium_rate; | |
| 73 static float time; | |
| 74 | |
| 75 extern const float WATER_VAPOR_PRESSURE; | |
| 76 extern const float HELIUM_TIME_CONSTANT[]; | |
| 77 extern const float NITROGEN_TIME_CONSTANT[]; | |
| 78 | |
| 79 float starting_ambient_pressure = pInput->lifeData.pressure_ambient_bar * 10; | |
| 80 float ending_ambient_pressure = ending_depth_meter + pInput->lifeData.pressure_surface_bar * 10; | |
| 81 | |
| 82 if((rate_meter_per_minutes <= 0) || (starting_ambient_pressure >= ending_ambient_pressure)) | |
| 83 return 0; | |
| 84 | |
| 85 for(i=0; i<16; i++) | |
| 86 { | |
| 87 initial_helium_pressure[i] = pInput->lifeData.tissue_helium_bar[i] * 10.0f; | |
| 88 initial_nitrogen_pressure[i] = pInput->lifeData.tissue_nitrogen_bar[i] * 10.0f; | |
| 89 } | |
| 90 | |
| 91 //New | |
| 92 time = (ending_ambient_pressure - starting_ambient_pressure) / rate_meter_per_minutes; | |
| 93 decom_get_inert_gases(starting_ambient_pressure / 10, &pInput->lifeData.actualGas, &fraction_nitrogen_begin, &fraction_helium_begin ); | |
| 94 decom_get_inert_gases(ending_ambient_pressure / 10, &pInput->lifeData.actualGas, &fraction_nitrogen_end, &fraction_helium_end ); | |
| 95 initial_inspired_he_pressure = (starting_ambient_pressure - WATER_VAPOR_PRESSURE) * fraction_helium_begin; | |
| 96 initial_inspired_n2_pressure = (starting_ambient_pressure - WATER_VAPOR_PRESSURE) * fraction_nitrogen_begin; | |
| 97 helium_rate = ((ending_ambient_pressure - WATER_VAPOR_PRESSURE)* fraction_helium_end - initial_inspired_he_pressure)/time; | |
| 98 nitrogen_rate = ((ending_ambient_pressure - WATER_VAPOR_PRESSURE)* fraction_nitrogen_end - initial_inspired_n2_pressure)/time; | |
| 99 pInput->lifeData.pressure_ambient_bar = ending_ambient_pressure/10; | |
| 100 | |
| 101 for( i = 0; i < 16; i++) | |
| 102 { | |
| 103 pInput->lifeData.tissue_helium_bar[i] = schreiner_equation__2(&initial_inspired_he_pressure, | |
| 104 &helium_rate, | |
| 105 &time, | |
| 106 &HELIUM_TIME_CONSTANT[i], | |
| 107 &initial_helium_pressure[i])/10.0f; | |
| 108 | |
| 109 pInput->lifeData.tissue_nitrogen_bar[i] = schreiner_equation__2(&initial_inspired_n2_pressure, | |
| 110 &nitrogen_rate, | |
| 111 &time, | |
| 112 &NITROGEN_TIME_CONSTANT[i], | |
| 113 &initial_nitrogen_pressure[i]) / 10.0f; | |
| 114 } | |
| 115 | |
| 116 | |
| 117 calc_crushing_pressure(&pInput->lifeData, &pInput->vpm,initial_helium_pressure,initial_nitrogen_pressure,starting_ambient_pressure, rate_meter_per_minutes); | |
| 118 | |
| 119 pInput->lifeData.dive_time_seconds += ((ending_ambient_pressure - starting_ambient_pressure)/rate_meter_per_minutes) * 60; | |
| 120 return 1; | |
| 121 } | |
| 122 void init_buehlmann(SDiveState* pInput) | |
| 123 { | |
| 124 pInput->diveSettings.internal__pressure_first_stop_ambient_bar_as_upper_limit_for_gf_low_otherwise_zero = 0; | |
| 125 pInput->lifeData.dive_time_seconds = 0; | |
| 126 for(int i=0;i<BUEHLMANN_STRUCT_MAX_GASES;i++) | |
| 127 { | |
| 128 pInput->diveSettings.decogaslist[i].change_during_ascent_depth_meter_otherwise_zero = 0; | |
| 129 pInput->diveSettings.decogaslist[i].nitrogen_percentage = 79 - i; | |
| 130 pInput->diveSettings.decogaslist[i].helium_percentage = i; | |
| 131 pInput->diveSettings.decogaslist[i].setPoint_cbar = 0; | |
| 132 } | |
| 133 pInput->lifeData.actualGas = pInput->diveSettings.decogaslist[0]; | |
| 134 pInput->diveSettings.last_stop_depth_bar = 0.3f; | |
| 135 pInput->diveSettings.input_next_stop_increment_depth_bar = 0.3f; | |
| 136 | |
| 137 pInput->decolistVPM.output_time_to_surface_seconds = 0; | |
| 138 pInput->decolistFutureVPM.output_time_to_surface_seconds = 0; | |
| 139 pInput->decolistBuehlmann.output_time_to_surface_seconds = 0; | |
| 140 pInput->decolistFutureBuehlmann.output_time_to_surface_seconds = 0; | |
| 141 for(int i=0;i<DECOINFO_STRUCT_MAX_STOPS;i++) | |
| 142 { | |
| 143 pInput->decolistVPM.output_stop_length_seconds[i] = 0; | |
| 144 pInput->decolistFutureVPM.output_stop_length_seconds[i] = 0; | |
| 145 | |
| 146 pInput->decolistBuehlmann.output_stop_length_seconds[i] = 0; | |
| 147 pInput->decolistFutureBuehlmann.output_stop_length_seconds[i] = 0; | |
| 148 } | |
| 149 for(int i=0;i<16;i++) | |
| 150 { | |
| 151 pInput->lifeData.tissue_nitrogen_bar[i] = 0.750927f; | |
| 152 pInput->lifeData.tissue_helium_bar[i] = 0; | |
| 153 } | |
| 154 pInput->diveSettings.gf_high = 80; | |
| 155 pInput->diveSettings.gf_low = 20; | |
| 156 pInput->diveSettings.vpm_conservatism = 2; | |
| 157 | |
| 158 pInput->lifeData.pressure_surface_bar = 1.0f; | |
| 159 pInput->lifeData.pressure_ambient_bar = 1.0f; | |
| 160 | |
| 161 pInput->warnings.decoMissed = 0; | |
| 162 pInput->events.gasChange = 0; | |
| 163 pInput->events.info_GasChange = 0; | |
| 281 | 164 pInput->events.info_manualGasSetO2 = 0; |
| 165 pInput->events.info_manualGasSetHe = 0; | |
| 166 pInput->events.manualGasSet = 0; | |
| 38 | 167 pInput->warnings.ppO2High = 0; |
| 168 pInput->warnings.ppO2Low = 0; | |
| 169 pInput->warnings.slowWarning = 0; | |
| 170 | |
| 171 | |
| 172 //pInput->decolistVPM.UNUSED_input_necessary_stop_length_seconds_otherwise_zero[i] = 0; | |
| 173 /*for(i=0;i<BUEHLMANN_STRUCT_MAX_ASCENDRATES;i++) | |
| 174 { | |
| 175 pInput->lifeData.ascentrate[i].rate_bar_per_minute = 1.2f; | |
| 176 pInput->lifeData.ascentrate[i].use_from_depth_bar = 0; // only one ascendrate at the moment | |
| 177 }*/ | |
| 178 //pInput->diveSettings.input_second_stop_depth_bar = 0.6f; | |
| 179 //pInput->lifeData.actual_gas_id = 0; | |
| 180 //pInput->lifeData.actual_setpoint_bar_if_rebreather_otherwise_zero = 0; | |
| 181 //pInput->lifeData.distance_used_below_stop_levels_bar = 0; | |
| 182 // pInput->lifeData.pressure_first_stop_ambient_bar_as_upper_limit_for_gf_low_otherwise_zero = 0; | |
| 183 | |
| 184 } | |
| 185 | |
| 186 | |
| 187 void init_buehlmann2(SDiveState* pInput) | |
| 188 { | |
| 189 pInput->diveSettings.internal__pressure_first_stop_ambient_bar_as_upper_limit_for_gf_low_otherwise_zero = 0; | |
| 190 pInput->lifeData.dive_time_seconds = 0; | |
| 191 int i=0; | |
| 192 for(i=0;i<BUEHLMANN_STRUCT_MAX_GASES;i++) | |
| 193 { | |
| 194 pInput->diveSettings.decogaslist[i].change_during_ascent_depth_meter_otherwise_zero = 0; | |
| 195 pInput->diveSettings.decogaslist[i].nitrogen_percentage = 20; | |
| 196 pInput->diveSettings.decogaslist[i].helium_percentage = 70; | |
| 197 pInput->diveSettings.decogaslist[i].setPoint_cbar = 0; | |
| 198 } | |
| 199 pInput->lifeData.actualGas = pInput->diveSettings.decogaslist[0]; | |
| 200 /*for(i=0;i<BUEHLMANN_STRUCT_MAX_ASCENDRATES;i++) | |
| 201 { | |
| 202 pInput->lifeData.ascentrate[i].rate_bar_per_minute = 1.2f; | |
| 203 pInput->lifeData.ascentrate[i].use_from_depth_bar = 0; // only one ascendrate at the moment | |
| 204 }*/ | |
| 205 pInput->diveSettings.last_stop_depth_bar = 0.3f; | |
| 206 //pInput->diveSettings.input_second_stop_depth_bar = 0.6f; | |
| 207 pInput->diveSettings.input_next_stop_increment_depth_bar = 0.3f; | |
| 208 pInput->decolistVPM.output_time_to_surface_seconds = 0; | |
| 209 pInput->decolistFutureVPM.output_time_to_surface_seconds = 0; | |
| 210 pInput->decolistBuehlmann.output_time_to_surface_seconds = 0; | |
| 211 pInput->decolistFutureBuehlmann.output_time_to_surface_seconds = 0; | |
| 212 for(int i=0;i<DECOINFO_STRUCT_MAX_STOPS;i++) | |
| 213 { | |
| 214 pInput->decolistVPM.output_stop_length_seconds[i] = 0; | |
| 215 pInput->decolistFutureVPM.output_stop_length_seconds[i] = 0; | |
| 216 | |
| 217 pInput->decolistBuehlmann.output_stop_length_seconds[i] = 0; | |
| 218 pInput->decolistFutureBuehlmann.output_stop_length_seconds[i] = 0; | |
| 219 } | |
| 220 for(i=0;i<16;i++) | |
| 221 { | |
| 222 pInput->lifeData.tissue_nitrogen_bar[i] = 0.750927f; | |
| 223 pInput->lifeData.tissue_helium_bar[i] = 0; | |
| 224 } | |
| 225 | |
| 226 // pInput->lifeData.distance_used_below_stop_levels_bar = 0; | |
| 227 pInput->diveSettings.gf_high = 80; | |
| 228 pInput->diveSettings.gf_low = 20; | |
| 229 pInput->diveSettings.vpm_conservatism = 2; | |
| 230 // pInput->lifeData.pressure_first_stop_ambient_bar_as_upper_limit_for_gf_low_otherwise_zero = 0; | |
| 231 pInput->lifeData.pressure_surface_bar = 1.0f; | |
| 232 pInput->lifeData.pressure_ambient_bar = 1.0f; | |
| 233 } | |
| 234 | |
| 235 | |
| 236 _Bool test1() | |
| 237 { | |
| 238 /* debug code with watch */ | |
| 239 static int32_t output_time_to_surface_minutes; | |
| 240 static int32_t counter = 0; | |
| 241 //static float decotable_minutes[DECOINFO_STRUCT_MAX_STOPS]; | |
| 242 | |
| 243 /* all the rest */ | |
| 244 SDiveState input; | |
| 245 | |
| 246 init_buehlmann(&input); | |
| 247 //vpm conservatism = 0, repetitive = false, | |
| 248 vpm_init(&input.vpm,0,false,0); | |
| 249 | |
| 250 //runter auf 70 meter mit 26 meter/minute | |
| 251 simulate_descent(&input, 70.0f, 26.0f); | |
| 252 //10 minuten settigung | |
| 253 //buehlmann__test__saturate_tissues(&input, 10 * 60); | |
| 254 decom_tissues_exposure(10 * 60, &input.lifeData); | |
| 255 //buehlmann_calc_deco(&input); | |
| 256 | |
| 257 vpm_calc(&(input.lifeData),&(input.diveSettings),&(input.vpm),&(input.decolistVPM), DECOSTOPS); | |
| 258 | |
| 259 //Check time to surface 46 min +- 0.6 | |
| 260 // MultiDeco hw: 42 min | |
| 261 | |
| 262 output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60; | |
| 263 if (output_time_to_surface_minutes != 46) | |
| 264 counter = 0; | |
| 265 else | |
| 266 counter++; | |
| 267 | |
| 268 if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 46.0f) >= 0.6f) | |
| 269 return false; | |
| 270 /* for(i=0;i<DECOINFO_STRUCT_MAX_STOPS;i++) | |
| 271 { | |
| 272 | |
| 273 if(decotable_minutes[i] != ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f) | |
| 274 { | |
| 275 counter2++; | |
| 276 decotable_minutes[i] = ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f; | |
| 277 } | |
| 278 } | |
| 279 i = i;*/ | |
| 280 | |
| 281 vpm_saturation_after_ascent(&input.lifeData); | |
| 282 input.vpm.decomode_vpm_plus_conservatism_last_dive = input.diveSettings.vpm_conservatism; | |
| 283 | |
| 284 //Pause 60 min | |
| 285 decom_tissues_exposure(60 * 60, &input.lifeData ); | |
| 286 //buehlmann__test__saturate_tissues(&input, 60 * 60); | |
| 287 vpm_init(&input.vpm,0,true, 60 * 60); | |
| 288 //runter auf 70 meter mit 26 meter/minute | |
| 289 simulate_descent(&input, 70.0f, 26.0f); | |
| 290 //10 minuten settigung | |
| 291 //buehlmann__test__saturate_tissues(&input, 10 * 60); | |
| 292 decom_tissues_exposure(10 * 60, &input.lifeData); | |
| 293 vpm_calc(&(input.lifeData),&(input.diveSettings),&(input.vpm),&(input.decolistVPM), DECOSTOPS); | |
| 294 //Check time to surface 46 min +- 0.6 | |
| 295 // MultiDeco hw: 42 min | |
| 296 | |
| 297 output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60; | |
| 298 | |
| 299 if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 57.0f) >= 0.6f) | |
| 300 return false; | |
| 301 | |
| 302 return true; | |
| 303 } | |
| 304 | |
| 305 | |
| 306 uint8_t test2_unapproved(void) | |
| 307 { | |
| 308 /* debug code with watch */ | |
| 309 static int32_t output_time_to_surface_minutes; | |
| 310 static int32_t counter = 0; | |
| 311 static float decotable_minutes[DECOINFO_STRUCT_MAX_STOPS]; | |
| 312 static int32_t counter2 = 0; | |
| 313 | |
| 314 /* all the rest */ | |
| 315 SDiveState input; | |
| 316 int i; | |
| 317 | |
| 318 init_buehlmann(&input); | |
| 319 //vpm conservatism = 3, repetitive = false, | |
| 320 vpm_init(&(input.vpm),3,false,0); | |
| 321 | |
| 322 //runter auf 70 meter mit 26 meter/minute | |
| 323 simulate_descent(&input, 70.0f, 26.0f); | |
| 324 //30 minuten saetigung | |
| 325 //buehlmann__test__saturate_tissues(&input, 30 * 60); | |
| 326 decom_tissues_exposure(30 * 60, &input.lifeData ); | |
| 327 //buehlmann_calc_deco(&input); | |
| 328 | |
| 329 vpm_calc(&(input.lifeData),&(input.diveSettings),&(input.vpm),&(input.decolistVPM), DECOSTOPS); | |
| 330 //Check time to surface 179.833 min (Peter Version 140415) +- 0.6, MultiDeco is 195 min | |
| 331 | |
| 332 output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60; | |
| 333 if (output_time_to_surface_minutes != 180) | |
| 334 counter = 0; | |
| 335 else | |
| 336 counter++; | |
| 337 | |
| 338 if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 180.0f) >= 0.6f) | |
| 339 return false; | |
| 340 for(i=0;i<DECOINFO_STRUCT_MAX_STOPS;i++) | |
| 341 { | |
| 342 if(decotable_minutes[i] != ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f) | |
| 343 { | |
| 344 counter2++; | |
| 345 decotable_minutes[i] = ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f; | |
| 346 } | |
| 347 } | |
| 348 return true; | |
| 349 } | |
| 350 | |
| 351 /** | |
| 352 ****************************************************************************** | |
| 353 * @brief test 3 | |
| 354 * Trimix 10/70 | |
| 355 * everything else identical to test1 by Peter Ryser | |
| 356 * @version V0.0.1 | |
| 357 * @date 19-April-2014 | |
| 358 * @retval 1 for result is similar to DRx code, 0 otherwise | |
| 359 ****************************************************************************** | |
| 360 */ | |
| 361 uint8_t test3_unapproved(void) | |
| 362 { | |
| 363 | |
| 364 /* debug code with watch */ | |
| 365 static int32_t output_time_to_surface_minutes; | |
| 366 static int32_t counter = 0; | |
| 367 static float decotable_minutes[DECOINFO_STRUCT_MAX_STOPS]; | |
| 368 static int32_t counter2 = 0; | |
| 369 | |
| 370 /* all the rest */ | |
| 371 SDiveState input; | |
| 372 int i; | |
| 373 | |
| 374 init_buehlmann2(&input); | |
| 375 //vpm conservatism = 0, repetitive = false, | |
| 376 vpm_init(&(input.vpm),0,false,0); | |
| 377 | |
| 378 //runter auf 70 meter mit 26 meter/minute | |
| 379 simulate_descent(&input, 70.0f, 26.0f); | |
| 380 //10 minuten settigung | |
| 381 decom_tissues_exposure(10 * 60, &input.lifeData); | |
| 382 //buehlmann__test__saturate_tissues(&input, 10 * 60); | |
| 383 //buehlmann_calc_deco(&input); | |
| 384 | |
| 385 vpm_calc(&(input.lifeData),&(input.diveSettings),&(input.vpm),&(input.decolistVPM), DECOSTOPS); | |
| 386 | |
| 387 //Check time to surface 46 min +- 0.6 | |
| 388 | |
| 389 output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60; | |
| 390 if (output_time_to_surface_minutes != 46) | |
| 391 counter = 0; | |
| 392 else | |
| 393 counter++; | |
| 394 | |
| 395 if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 46.0f) >= 0.6f) | |
| 396 return false; | |
| 397 for(i=0;i<DECOINFO_STRUCT_MAX_STOPS;i++) | |
| 398 { | |
| 399 if(decotable_minutes[i] != ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f) | |
| 400 { | |
| 401 counter2++; | |
| 402 decotable_minutes[i] = ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f; | |
| 403 } | |
| 404 } | |
| 405 return true; | |
| 406 } | |
| 407 | |
| 408 /** | |
| 409 ****************************************************************************** | |
| 410 * @brief test 4 - find the limit | |
| 411 * Trimix 10/70 | |
| 412 * 200 Meter, 30 Minuten | |
| 413 * @version V0.0.1 | |
| 414 * @date 19-April-2014 | |
| 415 * @retval 1 for result is similar to DRx code, 0 otherwise | |
| 416 ****************************************************************************** | |
| 417 */ | |
| 418 uint8_t test4_unapproved(void) | |
| 419 { | |
| 420 | |
| 421 /* debug code with watch */ | |
| 422 static int32_t output_time_to_surface_minutes; | |
| 423 static int32_t counter = 0; | |
| 424 static float decotable_minutes[DECOINFO_STRUCT_MAX_STOPS]; | |
| 425 static int32_t counter2 = 0; | |
| 426 | |
| 427 /* all the rest */ | |
| 428 SDiveState input; | |
| 429 int i; | |
| 430 | |
| 431 init_buehlmann2(&input); | |
| 432 //vpm conservatism = 0, repetitive = false, | |
| 433 vpm_init(&input.vpm,0,false,0); | |
| 434 | |
| 435 //runter auf 70 meter mit 26 meter/minute | |
| 436 simulate_descent(&input, 200.0f, 26.0f); | |
| 437 //10 minuten settigung | |
| 438 decom_tissues_exposure(10 * 60, &input.lifeData ); | |
| 439 //buehlmann__test__saturate_tissues(&input, 30 * 60); | |
| 440 //buehlmann_calc_deco(&input); | |
| 441 | |
| 442 vpm_calc(&(input.lifeData),&(input.diveSettings),&(input.vpm),&(input.decolistVPM), DECOSTOPS); | |
| 443 | |
| 444 //Check time to surface 1270 min | |
| 445 // Multi Deco 1270 Minuten | |
| 446 | |
| 447 output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60; | |
| 448 if (output_time_to_surface_minutes != 1270) | |
| 449 counter = 0; | |
| 450 else | |
| 451 counter++; | |
| 452 | |
| 453 if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 1270.0f) >= 0.6f) | |
| 454 return false; | |
| 455 for(i=0;i<DECOINFO_STRUCT_MAX_STOPS;i++) | |
| 456 { | |
| 457 if(decotable_minutes[i] != ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f) | |
| 458 { | |
| 459 counter2++; | |
| 460 decotable_minutes[i] = ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f; | |
| 461 } | |
| 462 } | |
| 463 return true; | |
| 464 } | |
| 465 | |
| 466 /*uint8_t test5_unapproved(uint32_t frame1, uint32_t frame2, uint32_t frame3, uint32_t frame4)*/ | |
| 467 uint8_t test5_unapproved(void) | |
| 468 { | |
| 469 /* debug code with watch */ | |
| 470 static int32_t output_time_to_surface_minutes; | |
| 471 static int32_t counter = 0; | |
| 472 | |
| 473 // static int32_t counter2 = 0; | |
| 474 /* all the rest */ | |
| 475 SDiveState input; | |
| 476 | |
| 477 //uint32_t frame[5]; | |
| 478 | |
| 479 uint8_t vpm_count; | |
| 480 | |
| 481 /* | |
| 482 frame[0] = frame1; | |
| 483 frame[1] = frame2; | |
| 484 frame[2] = frame3; | |
| 485 frame[3] = frame4; | |
| 486 frame[4] = frame[0]; | |
| 487 */ | |
| 488 init_buehlmann(&input); | |
| 489 vpm_init(&input.vpm,0,false,0); | |
| 490 logbook_initNewdiveProfile(&input,&Settings); | |
| 491 setSimulationValues(12, 26 , 70, 30); | |
| 492 | |
| 493 long time = 60 * 70 / 26 + 10 *60; | |
| 494 | |
| 495 vpm_count = 0; | |
| 496 while(input.lifeData.dive_time_seconds < time ) | |
| 497 { | |
| 498 /* frame[4] = frame[0]; | |
| 499 frame[0] = frame[1]; | |
| 500 frame[1] = frame[2]; | |
| 501 frame[2] = frame[3]; | |
| 502 frame[3] = frame[4];*/ | |
| 503 UpdateLifeDataTest(&input); | |
| 504 | |
| 505 vpm_count++; | |
| 506 if(vpm_count > 20) | |
| 507 { | |
| 508 vpm_calc(&input.lifeData, &(input.diveSettings),&input.vpm, &input.decolistVPM, DECOSTOPS); | |
| 509 vpm_count = 0; | |
| 510 } | |
| 511 | |
| 512 /* | |
| 513 #ifdef VGAOUT | |
| 514 tVGA_refresh(frame[1], &input); | |
| 515 GFX_VGA_transform(frame[1],frame[0]); | |
| 516 GFX_SetFrameBuffer(frame[0], TOP_LAYER); | |
| 517 GFX_clear_buffer(frame[3]); // frame[3] is the previous frame[0] | |
| 518 #endif | |
| 519 */ | |
| 520 if(input.lifeData.dive_time_seconds == 60 *5) | |
| 521 { | |
| 522 input.events.gasChange = 1; | |
| 523 input.events.info_GasChange = 2; | |
| 524 } | |
| 525 else | |
| 526 { | |
| 527 input.events.gasChange = 0; | |
| 528 input.events.info_GasChange = 0; | |
| 529 } | |
|
269
6e78137952af
cleanup: do not pass large struct by value
Jan Mulder <jlmulder@xs4all.nl>
parents:
219
diff
changeset
|
530 logbook_writeSample(&input); |
| 38 | 531 } |
| 532 volatile SLogbookHeader* logbookHeader = logbook_getCurrentHeader(); | |
| 533 | |
| 534 logbookHeader->total_diveTime_seconds = input.lifeData.dive_time_seconds; | |
| 535 logbookHeader->maxDepth = input.lifeData.max_depth_meter * 100; | |
| 536 logbook_EndDive(); | |
| 537 | |
| 538 output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60; | |
| 539 if (output_time_to_surface_minutes != 46) | |
| 540 counter = 0; | |
| 541 else | |
| 542 counter++; | |
| 543 | |
| 544 if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 46.0f) >= 0.6f) | |
| 545 return false; | |
| 546 | |
| 547 return true; | |
| 548 } | |
| 549 | |
| 550 | |
| 551 uint8_t test6_unapproved(void) | |
| 552 { | |
| 553 /* debug code with watch */ | |
| 554 static int32_t output_time_to_surface_minutes; | |
| 555 static int32_t counter = 0; | |
| 556 | |
| 557 // static int32_t counter2 = 0; | |
| 558 /* all the rest */ | |
| 559 SDiveState input; | |
| 560 | |
| 561 //uint32_t frame[5]; | |
| 562 | |
| 563 uint8_t vpm_count; | |
| 564 | |
| 565 init_buehlmann(&input); | |
| 566 vpm_init(&input.vpm,0,false,0); | |
| 567 logbook_initNewdiveProfile(&input,&Settings); | |
| 568 setSimulationValues(12, 26 , 65, 15); | |
| 569 | |
| 570 long time = 60 * 70 / 26 + 10 *60; | |
| 571 | |
| 572 vpm_count = 0; | |
| 573 while(input.lifeData.dive_time_seconds < time ) | |
| 574 { | |
| 575 UpdateLifeDataTest(&input); | |
| 576 | |
| 577 vpm_count++; | |
| 578 if(vpm_count > 20) | |
| 579 { | |
| 580 vpm_calc(&input.lifeData, &(input.diveSettings),&input.vpm, &input.decolistVPM, DECOSTOPS); | |
| 581 vpm_count = 0; | |
| 582 } | |
| 583 if(input.lifeData.dive_time_seconds == 60 *5) | |
| 584 { | |
| 585 input.events.gasChange = 1; | |
| 586 input.events.info_GasChange = 2; | |
| 587 } | |
| 588 else | |
| 589 { | |
| 590 input.events.gasChange = 0; | |
| 591 input.events.info_GasChange = 0; | |
| 592 } | |
|
269
6e78137952af
cleanup: do not pass large struct by value
Jan Mulder <jlmulder@xs4all.nl>
parents:
219
diff
changeset
|
593 logbook_writeSample(&input); |
| 38 | 594 } |
| 595 volatile SLogbookHeader* logbookHeader = logbook_getCurrentHeader(); | |
| 596 | |
| 597 logbookHeader->total_diveTime_seconds = input.lifeData.dive_time_seconds; | |
| 598 logbookHeader->maxDepth = input.lifeData.max_depth_meter * 100; | |
| 599 logbook_EndDive(); | |
| 600 | |
| 601 output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60; | |
| 602 if (output_time_to_surface_minutes != 46) | |
| 603 counter = 0; | |
| 604 else | |
| 605 counter++; | |
| 606 | |
| 607 if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 46.0f) >= 0.6f) | |
| 608 return false; | |
| 609 | |
| 610 return true; | |
| 611 } | |
| 612 | |
| 613 uint8_t test7_unapproved(void) | |
| 614 { | |
| 615 /* debug code with watch */ | |
| 616 static int32_t output_time_to_surface_minutes; | |
| 617 static int32_t counter = 0; | |
| 618 | |
| 619 // static int32_t counter2 = 0; | |
| 620 /* all the rest */ | |
| 621 SDiveState input; | |
| 622 | |
| 623 //uint32_t frame[5]; | |
| 624 | |
| 625 uint8_t vpm_count; | |
| 626 | |
| 627 init_buehlmann(&input); | |
| 628 vpm_init(&input.vpm,0,false,0); | |
| 629 logbook_initNewdiveProfile(&input,&Settings); | |
| 630 setSimulationValues(12, 26 , 40, 45); | |
| 631 | |
| 632 long time = 60 * 70 / 26 + 10 *60; | |
| 633 | |
| 634 vpm_count = 0; | |
| 635 while(input.lifeData.dive_time_seconds < time ) | |
| 636 { | |
| 637 UpdateLifeDataTest(&input); | |
| 638 | |
| 639 vpm_count++; | |
| 640 if(vpm_count > 20) | |
| 641 { | |
| 642 vpm_calc(&input.lifeData,&input.diveSettings, &input.vpm, &input.decolistVPM, DECOSTOPS); | |
| 643 vpm_count = 0; | |
| 644 } | |
| 645 if(input.lifeData.dive_time_seconds == 60 *5) | |
| 646 { | |
| 647 input.events.gasChange = 1; | |
| 648 input.events.info_GasChange = 2; | |
| 649 } | |
| 650 else | |
| 651 { | |
| 652 input.events.gasChange = 0; | |
| 653 input.events.info_GasChange = 0; | |
| 654 } | |
|
269
6e78137952af
cleanup: do not pass large struct by value
Jan Mulder <jlmulder@xs4all.nl>
parents:
219
diff
changeset
|
655 logbook_writeSample(&input); |
| 38 | 656 } |
| 657 volatile SLogbookHeader* logbookHeader = logbook_getCurrentHeader(); | |
| 658 | |
| 659 logbookHeader->total_diveTime_seconds = input.lifeData.dive_time_seconds; | |
| 660 logbookHeader->maxDepth = input.lifeData.max_depth_meter * 100; | |
| 661 logbook_EndDive(); | |
| 662 | |
| 663 output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60; | |
| 664 if (output_time_to_surface_minutes != 46) | |
| 665 counter = 0; | |
| 666 else | |
| 667 counter++; | |
| 668 | |
| 669 if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 46.0f) >= 0.6f) | |
| 670 return false; | |
| 671 | |
| 672 return true; | |
| 673 } | |
| 674 | |
| 675 void test_log_only(uint8_t max_depth_meter, uint16_t divetime_minutes) | |
| 676 { | |
| 677 SDiveState input; | |
| 678 float ascendrate_seconds; | |
| 679 float descendrate_seconds; | |
| 680 uint32_t divetime_seconds; | |
| 681 uint32_t divetime_start_ascend; | |
| 682 | |
| 683 init_buehlmann(&input); | |
| 684 input.lifeData.max_depth_meter = 0.0; | |
| 685 input.lifeData.depth_meter = 0.0; | |
| 686 input.lifeData.temperature_celsius = 22.7; | |
| 687 | |
| 688 ascendrate_seconds = 12.0 / 60.0; | |
| 689 descendrate_seconds = 20.0 / 60.0; | |
| 690 divetime_seconds = divetime_minutes * 60; | |
| 691 divetime_start_ascend = divetime_seconds - (uint32_t)(max_depth_meter / ascendrate_seconds); | |
| 692 | |
| 693 logbook_initNewdiveProfile(&input,&Settings); | |
| 694 | |
| 695 while(input.lifeData.dive_time_seconds < divetime_seconds ) | |
| 696 { | |
| 697 input.lifeData.dive_time_seconds += 1; | |
| 698 | |
| 699 if(input.lifeData.max_depth_meter < (float)max_depth_meter) | |
| 700 { | |
| 701 input.lifeData.depth_meter += descendrate_seconds; | |
| 702 input.lifeData.max_depth_meter = input.lifeData.depth_meter; | |
| 703 } | |
| 704 else if((input.lifeData.dive_time_seconds >= divetime_start_ascend) && (input.lifeData.depth_meter > 0)) | |
| 705 { | |
| 706 input.lifeData.depth_meter -= ascendrate_seconds; | |
| 707 if(input.lifeData.depth_meter < 0) | |
| 708 input.lifeData.depth_meter = 0; | |
| 709 } | |
| 710 | |
|
269
6e78137952af
cleanup: do not pass large struct by value
Jan Mulder <jlmulder@xs4all.nl>
parents:
219
diff
changeset
|
711 logbook_writeSample(&input); |
| 38 | 712 } |
| 713 volatile SLogbookHeader* logbookHeader = logbook_getCurrentHeader(); | |
| 714 logbookHeader->total_diveTime_seconds = input.lifeData.dive_time_seconds; | |
| 715 logbookHeader->maxDepth = input.lifeData.max_depth_meter * 100; | |
| 716 logbook_EndDive(); | |
| 717 } | |
| 718 | |
| 719 | |
| 720 /** | |
| 721 ****************************************************************************** | |
| 722 * @brief test 101 | |
| 723 * a) for air | |
| 724 * b) air + oxygen | |
| 725 * c) Trimix 10/70 + oxygen | |
| 726 * 65 Meter, 42 Minuten with descent | |
| 727 * @version V0.0.1 | |
| 728 * @date 26-Oct-2014 | |
| 729 * @retval ToDo: 1 for result is similar to MultiDeco | |
| 730 ****************************************************************************** | |
| 731 */ | |
| 732 | |
| 733 uint8_t test101_buehlmann_unapproved(void) | |
| 734 { | |
| 735 /* all the rest */ | |
| 736 SDiveState input; | |
| 737 | |
| 738 init_buehlmann(&input); | |
| 739 | |
| 740 //Gas Change at 6 meter to oxygin | |
| 741 //input.diveSettings.decogaslist[1].change_during_ascent_depth_bar_otherwise_zero = 0.6f; | |
| 742 //input.diveSettings.decogaslist[1].nitrogen_percentage = 0; | |
| 743 //input.diveSettings.decogaslist[1].helium_percentage = 0; | |
| 744 input.diveSettings.gf_high = 100; | |
| 745 input.diveSettings.gf_low = 100; | |
| 746 input.diveSettings.ascentRate_meterperminute = 10.0f; | |
| 747 input.diveSettings.last_stop_depth_bar = 0.3f; // input.diveSettings.last_stop_depth_bar = 0.6f; /* ist egal bei oxygen */ | |
| 748 | |
| 749 //runter auf 65 meter mit 20 meter/minute | |
| 750 simulate_descent(&input, 65.0f, 20.0f); | |
| 751 //38min 45sec saettigung == 2325 sec | |
| 752 decom_tissues_exposure(30*60, &input.lifeData ); | |
| 753 // decom_tissues_exposure(2325, &input.lifeData ); | |
| 754 | |
| 755 //vpm_calc(&(input.lifeData),&(input.diveSettings),&(input.vpm),&(input.decolistVPM)); | |
| 756 buehlmann_calc_deco(&input.lifeData,&input.diveSettings,&input.decolistBuehlmann); | |
| 757 | |
| 758 | |
| 759 //Check time to surface MultiDeco 4.04 | |
| 760 // 308 min with Air | |
| 761 // 190,5 min with Air + 6m last stop with oxygen | |
| 762 // 538 min with Trimix 10/70 and oxygen at 6m | |
| 763 | |
| 764 // ... | |
| 765 | |
| 766 return true; | |
| 767 } |