Mercurial > public > ostc4

view Small_CPU/Inc/vpm.h @ 250:822416168585 bm-2

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Buelmann: new implementation for ceiling Since my first functional fix in the ceiling computation in commit ceecabfddb57, I noticed that the computation used a linear search, that became rather computational expensive after that commit. The simple question is: why not a binary search? So, this commit implements the binary search. But there is a long story attached to this. Comparing ceiling results from hwOS and this OSTC4 code were very different. Basically, the original OSTC4 algorithm computed the ceiling using the same GFlow to GFhigh slope, in such a way, that the ceiling was in sync with the presented deco stops, where the hwOS code presents a GFhigh based ceiling. This said, it is more logical when the OSTC4 and hwOS code give similar results. This new recursive algorithm gives very similar results for the ceiling compared to hwOS. To be complete here, the Buelmann ceiling is the depth to which you can ascend, so that the leading tissue reaches GFhigh. This also explains why the deepest deco stop is normally deeper than the ceiling (unless one dives with GF like 80/80). The code implemented here is rather straightforward recursion. Signed-off-by: Jan Mulder <jlmulder@xs4all.nl>
author Jan Mulder <jlmulder@xs4all.nl>
date Thu, 11 Apr 2019 17:48:48 +0200
parents 5f11787b4f42
children
line wrap: on
line source

/*
 * vpm.h
 *
 * Created: 23.04.2013 11:38:19
 *  Author: produktion04
 */


#ifndef VPM_H
#define VPM_H

#include "data_central.h"
extern long vpm_time_calc_begin;
enum DECOLIST{DECOSTOPS,FUTURESTOPS,BAILOUTSTOPS,NULLZEIT,OFF = -1}; // order is important!!
enum VPM_CALC_STATUS{CALC_END, CALC_BEGIN, CALC_CRITICAL, CALC_FINAL_DECO, CALC_NULLZEIT };
//int calc_crushing_pressure(float *starting_depth, float *ending_depth, float *rate);
//void vpm_check_calc(unsigned short* stoplist);
//void vpm_init(void);
//_Bool vpm_crush(void);

void vpm_reset_variables(void);
_Bool vpm_build_variables_from_file(unsigned long in_sdram_start);
int vpm_store_variables_in_sdram_for_transfer(unsigned long in_sdram_start);
void vpm_calc_crushing_pressure(float starting_ambient_pressure, float ending_ambient_pressure, float rate);
void vpm_deco_plan(unsigned short divetime,unsigned short * divetime_first_stop, float* first_stop_depth);

float schreiner_equation__2(float *initial_inspired_gas_pressure,float *rate_change_insp_gas_pressure,float *interval_time_minutes,  const float *gas_time_constant,float *initial_gas_pressure);

int  vpm_calc(SLifeData* pINPUT, SVpm* pVPM, SDecoinfo*  pDECOINFO);
void  vpm_saturation_after_ascent(SLifeData* input);
extern  const float helium_time_constant[16];
extern  const float nitrogen_time_constant[16];

#endif /* VPM_H */