Mercurial > public > ostc4
annotate Discovery/Src/test_vpm.c @ 726:8f3a8c85a6c4
Bugfix data synchronization after RTE start:
In previous version invalid data may be transferred in case the communication is corrupted directly after RTE startup. Root cause was that the data was evaluated without caring if a communication error was reported. In the new version the startup condition and evaluation of data is only done if valid data has been received by the RTE => SPI communication is up and runnig.
author | Ideenmodellierer |
---|---|
date | Sat, 14 Jan 2023 20:35:40 +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 } |