Mercurial > public > ostc4

view Discovery/Src/test_vpm.c @ 196:2885628ab3ba div-fixes-cleaup-2

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bugfix, minor: color the overview customview correctly When using a custom color (from the SYS2 menu, layout), the overview customview was the only one presented in the wrong color. While the fix for this is rather simple (in hindsight), it was a non trivial search trough rather obfuscated code. There is a specific function to set the text color index on this page, but this only worked for a single line text string, that has the color index byte as the very first character. The original author already recognized that this function "could be extended", but left this extension as an exercise to the reader. So, the coloring is fixed by extending the function, and actually using it for the overview customview. Signed-off-by: Jan Mulder <jlmulder@xs4all.nl>
author Jan Mulder <jlmulder@xs4all.nl>
date Wed, 20 Mar 2019 16:24:10 +0100
parents 5f11787b4f42
children f9b17e898a7a
line wrap: on
line source

///////////////////////////////////////////////////////////////////////////////
/// -*- coding: UTF-8 -*-
///
/// \file   Discovery/Src/test_vpm.c
/// \brief  test 101
/// \author Heinrichs Weikamp
/// \date   26-Oct-2014
///
/// \details
///
/// $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 <stdio.h>
#include <stdint.h>
//#include "LED.h"
//#include "Keyboard.h"
//#include "stm32f4xx_hal.h"
#include "buehlmann.h"
#include "calc_crush.h"
#include "vpm.h"
#include "display.h"
#include "test_vpm.h"
#include "math.h"
#include "data_central.h"
#include "decom.h"
#include "logbook.h"
#include "tInfoLog.h"

#define true 1
#define false 0
//#define uint8_t unsigned char

extern SSettings Settings;

_Bool simulate_descent(SDiveState* pInput, float ending_depth_meter, float rate_meter_per_minutes);
void init_buehlmann(SDiveState* pInput);

_Bool test1(void);
uint8_t test2_unapproved(void);
uint8_t test3_unapproved(void);

_Bool simulate_descent(SDiveState* pInput,
                       float ending_depth_meter,
                       float rate_meter_per_minutes)
{
    int i =0;
    static	float initial_helium_pressure[16];
    static	float initial_nitrogen_pressure[16];
    static	float initial_inspired_he_pressure;
    static	float initial_inspired_n2_pressure;
    static	float fraction_nitrogen_begin;
    static	float fraction_nitrogen_end;
    static	float fraction_helium_begin;
    static	float fraction_helium_end;
    static	float nitrogen_rate;
    static	float helium_rate;
    static	float time;

    extern const float WATER_VAPOR_PRESSURE;
    extern const float HELIUM_TIME_CONSTANT[];
    extern const float NITROGEN_TIME_CONSTANT[];

    float starting_ambient_pressure = pInput->lifeData.pressure_ambient_bar * 10;
    float ending_ambient_pressure = ending_depth_meter + pInput->lifeData.pressure_surface_bar * 10;

    if((rate_meter_per_minutes <= 0) || (starting_ambient_pressure >= ending_ambient_pressure))
        return 0;

    for(i=0; i<16; i++)
    {
        initial_helium_pressure[i] = pInput->lifeData.tissue_helium_bar[i] * 10.0f;
        initial_nitrogen_pressure[i] = pInput->lifeData.tissue_nitrogen_bar[i] * 10.0f;
    }

    //New
    time = (ending_ambient_pressure - starting_ambient_pressure) / rate_meter_per_minutes;
    decom_get_inert_gases(starting_ambient_pressure / 10, &pInput->lifeData.actualGas, &fraction_nitrogen_begin, &fraction_helium_begin );
    decom_get_inert_gases(ending_ambient_pressure   / 10, &pInput->lifeData.actualGas, &fraction_nitrogen_end, &fraction_helium_end );
    initial_inspired_he_pressure =	(starting_ambient_pressure - WATER_VAPOR_PRESSURE) * fraction_helium_begin;
    initial_inspired_n2_pressure =	(starting_ambient_pressure - WATER_VAPOR_PRESSURE) * fraction_nitrogen_begin;
    helium_rate = ((ending_ambient_pressure  - WATER_VAPOR_PRESSURE)* fraction_helium_end - initial_inspired_he_pressure)/time;
    nitrogen_rate = ((ending_ambient_pressure  - WATER_VAPOR_PRESSURE)* fraction_nitrogen_end - initial_inspired_n2_pressure)/time;
    pInput->lifeData.pressure_ambient_bar = ending_ambient_pressure/10;

    for( i = 0; i < 16; i++)
    {
        pInput->lifeData.tissue_helium_bar[i] = schreiner_equation__2(&initial_inspired_he_pressure,
        &helium_rate,
        &time,
        &HELIUM_TIME_CONSTANT[i],
        &initial_helium_pressure[i])/10.0f;

         pInput->lifeData.tissue_nitrogen_bar[i] = schreiner_equation__2(&initial_inspired_n2_pressure,
        &nitrogen_rate,
        &time,
        &NITROGEN_TIME_CONSTANT[i],
        &initial_nitrogen_pressure[i]) / 10.0f;
    }


    calc_crushing_pressure(&pInput->lifeData, &pInput->vpm,initial_helium_pressure,initial_nitrogen_pressure,starting_ambient_pressure, rate_meter_per_minutes);

    pInput->lifeData.dive_time_seconds += ((ending_ambient_pressure - starting_ambient_pressure)/rate_meter_per_minutes) * 60;
    return 1;
}
void init_buehlmann(SDiveState* pInput)
{
    pInput->diveSettings.internal__pressure_first_stop_ambient_bar_as_upper_limit_for_gf_low_otherwise_zero = 0;
    pInput->lifeData.dive_time_seconds = 0;
    for(int i=0;i<BUEHLMANN_STRUCT_MAX_GASES;i++)
    {
        pInput->diveSettings.decogaslist[i].change_during_ascent_depth_meter_otherwise_zero = 0;
        pInput->diveSettings.decogaslist[i].nitrogen_percentage = 79 - i;
        pInput->diveSettings.decogaslist[i].helium_percentage = i;
        pInput->diveSettings.decogaslist[i].setPoint_cbar = 0;
    }
    pInput->lifeData.actualGas = pInput->diveSettings.decogaslist[0];
    pInput->diveSettings.last_stop_depth_bar = 0.3f;
    pInput->diveSettings.input_next_stop_increment_depth_bar = 0.3f;

    pInput->decolistVPM.output_time_to_surface_seconds = 0;
    pInput->decolistFutureVPM.output_time_to_surface_seconds = 0;
    pInput->decolistBuehlmann.output_time_to_surface_seconds = 0;
    pInput->decolistFutureBuehlmann.output_time_to_surface_seconds = 0;
    for(int i=0;i<DECOINFO_STRUCT_MAX_STOPS;i++)
    {
        pInput->decolistVPM.output_stop_length_seconds[i] = 0;
        pInput->decolistFutureVPM.output_stop_length_seconds[i] = 0;

        pInput->decolistBuehlmann.output_stop_length_seconds[i] = 0;
        pInput->decolistFutureBuehlmann.output_stop_length_seconds[i] = 0;
    }
    for(int i=0;i<16;i++)
    {
        pInput->lifeData.tissue_nitrogen_bar[i] = 0.750927f;
        pInput->lifeData.tissue_helium_bar[i] = 0;
    }
    pInput->diveSettings.gf_high = 80;
    pInput->diveSettings.gf_low = 20;
    pInput->diveSettings.vpm_conservatism = 2;

    pInput->lifeData.pressure_surface_bar = 1.0f;
    pInput->lifeData.pressure_ambient_bar = 1.0f;

    pInput->warnings.decoMissed = 0;
    pInput->events.gasChange = 0;
    pInput->events.info_GasChange = 0;
    pInput->events.info_manuelGasSetO2 = 0;
    pInput->events.info_manuelGasSetHe = 0;
    pInput->events.manuelGasSet = 0;
    pInput->warnings.ppO2High = 0;
    pInput->warnings.ppO2Low = 0;
    pInput->warnings.slowWarning = 0;


    //pInput->decolistVPM.UNUSED_input_necessary_stop_length_seconds_otherwise_zero[i] = 0;
    /*for(i=0;i<BUEHLMANN_STRUCT_MAX_ASCENDRATES;i++)
    {
        pInput->lifeData.ascentrate[i].rate_bar_per_minute = 1.2f;
        pInput->lifeData.ascentrate[i].use_from_depth_bar = 0; // only one ascendrate at the moment
    }*/
    //pInput->diveSettings.input_second_stop_depth_bar = 0.6f;
    //pInput->lifeData.actual_gas_id = 0;
    //pInput->lifeData.actual_setpoint_bar_if_rebreather_otherwise_zero = 0;
    //pInput->lifeData.distance_used_below_stop_levels_bar = 0;
//	pInput->lifeData.pressure_first_stop_ambient_bar_as_upper_limit_for_gf_low_otherwise_zero = 0;

}


void init_buehlmann2(SDiveState* pInput)
{
    pInput->diveSettings.internal__pressure_first_stop_ambient_bar_as_upper_limit_for_gf_low_otherwise_zero = 0;
    pInput->lifeData.dive_time_seconds = 0;
    int i=0;
    for(i=0;i<BUEHLMANN_STRUCT_MAX_GASES;i++)
    {
        pInput->diveSettings.decogaslist[i].change_during_ascent_depth_meter_otherwise_zero = 0;
        pInput->diveSettings.decogaslist[i].nitrogen_percentage = 20;
        pInput->diveSettings.decogaslist[i].helium_percentage = 70;
        pInput->diveSettings.decogaslist[i].setPoint_cbar = 0;
    }
    pInput->lifeData.actualGas = pInput->diveSettings.decogaslist[0];
    /*for(i=0;i<BUEHLMANN_STRUCT_MAX_ASCENDRATES;i++)
    {
        pInput->lifeData.ascentrate[i].rate_bar_per_minute = 1.2f;
        pInput->lifeData.ascentrate[i].use_from_depth_bar = 0; // only one ascendrate at the moment
    }*/
    pInput->diveSettings.last_stop_depth_bar = 0.3f;
    //pInput->diveSettings.input_second_stop_depth_bar = 0.6f;
    pInput->diveSettings.input_next_stop_increment_depth_bar = 0.3f;
    pInput->decolistVPM.output_time_to_surface_seconds = 0;
    pInput->decolistFutureVPM.output_time_to_surface_seconds = 0;
    pInput->decolistBuehlmann.output_time_to_surface_seconds = 0;
    pInput->decolistFutureBuehlmann.output_time_to_surface_seconds = 0;
    for(int i=0;i<DECOINFO_STRUCT_MAX_STOPS;i++)
    {
        pInput->decolistVPM.output_stop_length_seconds[i] = 0;
        pInput->decolistFutureVPM.output_stop_length_seconds[i] = 0;

        pInput->decolistBuehlmann.output_stop_length_seconds[i] = 0;
        pInput->decolistFutureBuehlmann.output_stop_length_seconds[i] = 0;
    }
    for(i=0;i<16;i++)
    {
        pInput->lifeData.tissue_nitrogen_bar[i] = 0.750927f;
        pInput->lifeData.tissue_helium_bar[i] = 0;
    }

//	pInput->lifeData.distance_used_below_stop_levels_bar = 0;
    pInput->diveSettings.gf_high = 80;
    pInput->diveSettings.gf_low = 20;
    pInput->diveSettings.vpm_conservatism = 2;
//	pInput->lifeData.pressure_first_stop_ambient_bar_as_upper_limit_for_gf_low_otherwise_zero = 0;
    pInput->lifeData.pressure_surface_bar = 1.0f;
    pInput->lifeData.pressure_ambient_bar = 1.0f;
}


_Bool test1()
{
    /* debug code with watch */
    static int32_t output_time_to_surface_minutes;
    static int32_t counter = 0;
    //static float decotable_minutes[DECOINFO_STRUCT_MAX_STOPS];

    /* all the rest */
    SDiveState input;

    init_buehlmann(&input);
    //vpm conservatism = 0, repetitive = false,
    vpm_init(&input.vpm,0,false,0);

    //runter auf 70 meter mit 26 meter/minute
    simulate_descent(&input, 70.0f, 26.0f);
    //10 minuten settigung
    //buehlmann__test__saturate_tissues(&input,  10 * 60);
     decom_tissues_exposure(10 * 60, &input.lifeData);
    //buehlmann_calc_deco(&input);

    vpm_calc(&(input.lifeData),&(input.diveSettings),&(input.vpm),&(input.decolistVPM), DECOSTOPS);

    //Check time to surface 46 min +- 0.6
    // MultiDeco hw: 42 min

    output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60;
    if (output_time_to_surface_minutes != 46)
        counter = 0;
    else
        counter++;

    if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 46.0f) >= 0.6f)
        return false;
/*	for(i=0;i<DECOINFO_STRUCT_MAX_STOPS;i++)
    {

        if(decotable_minutes[i] != ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f)
        {
            counter2++;
            decotable_minutes[i] = ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f;
        }
    }
    i = i;*/

    vpm_saturation_after_ascent(&input.lifeData);
    input.vpm.decomode_vpm_plus_conservatism_last_dive = input.diveSettings.vpm_conservatism;

    //Pause 60 min
     decom_tissues_exposure(60 * 60, &input.lifeData );
    //buehlmann__test__saturate_tissues(&input,  60 * 60);
    vpm_init(&input.vpm,0,true, 60 * 60);
    //runter auf 70 meter mit 26 meter/minute
    simulate_descent(&input, 70.0f, 26.0f);
    //10 minuten settigung
    //buehlmann__test__saturate_tissues(&input,  10 * 60);
    decom_tissues_exposure(10 * 60, &input.lifeData);
    vpm_calc(&(input.lifeData),&(input.diveSettings),&(input.vpm),&(input.decolistVPM), DECOSTOPS);
    //Check time to surface 46 min +- 0.6
    // MultiDeco hw: 42 min

    output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60;

    if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 57.0f) >= 0.6f)
        return false;

    return true;
}


uint8_t test2_unapproved(void)
{
    /* debug code with watch */
    static int32_t output_time_to_surface_minutes;
    static int32_t counter = 0;
    static float decotable_minutes[DECOINFO_STRUCT_MAX_STOPS];
    static int32_t counter2 = 0;

    /* all the rest */
    SDiveState input;
    int i;

    init_buehlmann(&input);
    //vpm conservatism = 3, repetitive = false,
    vpm_init(&(input.vpm),3,false,0);

    //runter auf 70 meter mit 26 meter/minute
    simulate_descent(&input, 70.0f, 26.0f);
    //30 minuten saetigung
    //buehlmann__test__saturate_tissues(&input,  30 * 60);
     decom_tissues_exposure(30 * 60, &input.lifeData );
    //buehlmann_calc_deco(&input);

        vpm_calc(&(input.lifeData),&(input.diveSettings),&(input.vpm),&(input.decolistVPM), DECOSTOPS);
    //Check time to surface 179.833 min (Peter Version 140415) +- 0.6, MultiDeco is 195 min

    output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60;
    if (output_time_to_surface_minutes != 180)
        counter = 0;
    else
        counter++;

    if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 180.0f) >= 0.6f)
        return false;
    for(i=0;i<DECOINFO_STRUCT_MAX_STOPS;i++)
    {
        if(decotable_minutes[i] != ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f)
        {
            counter2++;
            decotable_minutes[i] = ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f;
        }
    }
    i = i;
    return true;
}

/**
  ******************************************************************************
  * @brief		test 3
  *						Trimix 10/70
    *						everything else identical to test1 by Peter Ryser
  * @version 	V0.0.1
  * @date   	19-April-2014
  * @retval 	1 for result is similar to DRx code, 0 otherwise
  ******************************************************************************
  */
uint8_t test3_unapproved(void)
{

    /* debug code with watch */
    static int32_t output_time_to_surface_minutes;
    static int32_t counter = 0;
    static float decotable_minutes[DECOINFO_STRUCT_MAX_STOPS];
    static int32_t counter2 = 0;

    /* all the rest */
    SDiveState input;
    int i;

    init_buehlmann2(&input);
    //vpm conservatism = 0, repetitive = false,
    vpm_init(&(input.vpm),0,false,0);

    //runter auf 70 meter mit 26 meter/minute
    simulate_descent(&input, 70.0f, 26.0f);
    //10 minuten settigung
     decom_tissues_exposure(10 * 60, &input.lifeData);
    //buehlmann__test__saturate_tissues(&input,  10 * 60);
    //buehlmann_calc_deco(&input);

       vpm_calc(&(input.lifeData),&(input.diveSettings),&(input.vpm),&(input.decolistVPM), DECOSTOPS);

    //Check time to surface 46 min +- 0.6

    output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60;
    if (output_time_to_surface_minutes != 46)
        counter = 0;
    else
        counter++;

    if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 46.0f) >= 0.6f)
        return false;
    for(i=0;i<DECOINFO_STRUCT_MAX_STOPS;i++)
    {
        if(decotable_minutes[i] != ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f)
        {
            counter2++;
            decotable_minutes[i] = ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f;
        }
    }
    i = i;
    return true;
}

/**
  ******************************************************************************
  * @brief		test 4 - find the limit
  *						Trimix 10/70
    *						200 Meter, 30 Minuten
  * @version 	V0.0.1
  * @date   	19-April-2014
  * @retval 	1 for result is similar to DRx code, 0 otherwise
  ******************************************************************************
  */
uint8_t test4_unapproved(void)
{

    /* debug code with watch */
    static int32_t output_time_to_surface_minutes;
    static int32_t counter = 0;
    static float decotable_minutes[DECOINFO_STRUCT_MAX_STOPS];
    static int32_t counter2 = 0;

    /* all the rest */
    SDiveState input;
    int i;

    init_buehlmann2(&input);
    //vpm conservatism = 0, repetitive = false,
    vpm_init(&input.vpm,0,false,0);

    //runter auf 70 meter mit 26 meter/minute
    simulate_descent(&input, 200.0f, 26.0f);
    //10 minuten settigung
     decom_tissues_exposure(10 * 60, &input.lifeData );
    //buehlmann__test__saturate_tissues(&input,  30 * 60);
    //buehlmann_calc_deco(&input);

    vpm_calc(&(input.lifeData),&(input.diveSettings),&(input.vpm),&(input.decolistVPM), DECOSTOPS);

    //Check time to surface 1270 min
    // Multi Deco 1270 Minuten

    output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60;
    if (output_time_to_surface_minutes != 1270)
        counter = 0;
    else
        counter++;

    if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 1270.0f) >= 0.6f)
        return false;
    for(i=0;i<DECOINFO_STRUCT_MAX_STOPS;i++)
    {
        if(decotable_minutes[i] != ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f)
        {
            counter2++;
            decotable_minutes[i] = ((float)input.decolistVPM.output_stop_length_seconds[i]) / 60.0f;
        }
    }
    i = i;
    return true;
}

/*uint8_t test5_unapproved(uint32_t frame1, uint32_t frame2, uint32_t frame3, uint32_t frame4)*/
uint8_t test5_unapproved(void)
{
    /* debug code with watch */
    static int32_t output_time_to_surface_minutes;
    static int32_t counter = 0;

//	static int32_t counter2 = 0;
    /* all the rest */
    SDiveState input;

    //uint32_t frame[5];

    uint8_t vpm_count;

    /*
    frame[0] = frame1;
    frame[1] = frame2;
    frame[2] = frame3;
    frame[3] = frame4;
    frame[4] = frame[0];
*/
    init_buehlmann(&input);
    vpm_init(&input.vpm,0,false,0);
    logbook_initNewdiveProfile(&input,&Settings);
    setSimulationValues(12, 26 , 70, 30);

    long time = 60 * 70 / 26 + 10 *60;

    vpm_count = 0;
    while(input.lifeData.dive_time_seconds < time )
    {
/*				frame[4] = frame[0];
                frame[0] = frame[1];
                frame[1] = frame[2];
                frame[2] = frame[3];
                frame[3] = frame[4];*/
        UpdateLifeDataTest(&input);

                vpm_count++;
                if(vpm_count > 20)
                {
                    vpm_calc(&input.lifeData, &(input.diveSettings),&input.vpm, &input.decolistVPM, DECOSTOPS);
                    vpm_count = 0;
                }

/*
#ifdef VGAOUT
                tVGA_refresh(frame[1], &input);
                GFX_VGA_transform(frame[1],frame[0]);
                GFX_SetFrameBuffer(frame[0], TOP_LAYER);
                GFX_clear_buffer(frame[3]); // frame[3] is the previous frame[0]
#endif
*/
            if(input.lifeData.dive_time_seconds == 60 *5)
            {
                input.events.gasChange = 1;
                input.events.info_GasChange = 2;
            }
            else
            {
                input.events.gasChange = 0;
                input.events.info_GasChange = 0;
            }
            logbook_writeSample(input);
    }
        volatile SLogbookHeader* logbookHeader = logbook_getCurrentHeader();

        logbookHeader->total_diveTime_seconds = input.lifeData.dive_time_seconds;
        logbookHeader->maxDepth = input.lifeData.max_depth_meter * 100;
         logbook_EndDive();

    output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60;
    if (output_time_to_surface_minutes != 46)
        counter = 0;
    else
        counter++;

    if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 46.0f) >= 0.6f)
        return false;

    return true;
}


uint8_t test6_unapproved(void)
{
    /* debug code with watch */
    static int32_t output_time_to_surface_minutes;
    static int32_t counter = 0;

//	static int32_t counter2 = 0;
    /* all the rest */
    SDiveState input;

    //uint32_t frame[5];

    uint8_t vpm_count;

    init_buehlmann(&input);
    vpm_init(&input.vpm,0,false,0);
    logbook_initNewdiveProfile(&input,&Settings);
    setSimulationValues(12, 26 , 65, 15);

    long time = 60 * 70 / 26 + 10 *60;

    vpm_count = 0;
    while(input.lifeData.dive_time_seconds < time )
    {
        UpdateLifeDataTest(&input);

                vpm_count++;
                if(vpm_count > 20)
                {
                    vpm_calc(&input.lifeData, &(input.diveSettings),&input.vpm, &input.decolistVPM, DECOSTOPS);
                    vpm_count = 0;
                }
                if(input.lifeData.dive_time_seconds == 60 *5)
                {
                        input.events.gasChange = 1;
                        input.events.info_GasChange = 2;
                }
                else
                {
                        input.events.gasChange = 0;
                        input.events.info_GasChange = 0;
                }
            logbook_writeSample(input);
    }
    volatile SLogbookHeader* logbookHeader = logbook_getCurrentHeader();

    logbookHeader->total_diveTime_seconds = input.lifeData.dive_time_seconds;
    logbookHeader->maxDepth = input.lifeData.max_depth_meter * 100;
    logbook_EndDive();

    output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60;
    if (output_time_to_surface_minutes != 46)
        counter = 0;
    else
        counter++;

    if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 46.0f) >= 0.6f)
        return false;

    return true;
}

uint8_t test7_unapproved(void)
{
    /* debug code with watch */
    static int32_t output_time_to_surface_minutes;
    static int32_t counter = 0;

//	static int32_t counter2 = 0;
    /* all the rest */
    SDiveState input;

    //uint32_t frame[5];

    uint8_t vpm_count;

    init_buehlmann(&input);
    vpm_init(&input.vpm,0,false,0);
    logbook_initNewdiveProfile(&input,&Settings);
    setSimulationValues(12, 26 , 40, 45);

    long time = 60 * 70 / 26 + 10 *60;

    vpm_count = 0;
    while(input.lifeData.dive_time_seconds < time )
    {
        UpdateLifeDataTest(&input);

                vpm_count++;
                if(vpm_count > 20)
                {
                    vpm_calc(&input.lifeData,&input.diveSettings, &input.vpm, &input.decolistVPM, DECOSTOPS);
                    vpm_count = 0;
                }
                if(input.lifeData.dive_time_seconds == 60 *5)
                {
                        input.events.gasChange = 1;
                        input.events.info_GasChange = 2;
                }
                else
                {
                        input.events.gasChange = 0;
                        input.events.info_GasChange = 0;
                }
            logbook_writeSample(input);
    }
    volatile SLogbookHeader* logbookHeader = logbook_getCurrentHeader();

    logbookHeader->total_diveTime_seconds = input.lifeData.dive_time_seconds;
    logbookHeader->maxDepth = input.lifeData.max_depth_meter * 100;
    logbook_EndDive();

    output_time_to_surface_minutes = input.decolistVPM.output_time_to_surface_seconds / 60;
    if (output_time_to_surface_minutes != 46)
        counter = 0;
    else
        counter++;

    if(fabsf( ((float)input.decolistVPM.output_time_to_surface_seconds / 60.0f) - 46.0f) >= 0.6f)
        return false;

    return true;
}

void test_log_only(uint8_t max_depth_meter, uint16_t divetime_minutes)
{
    SDiveState input;
    float ascendrate_seconds;
    float descendrate_seconds;
    uint32_t divetime_seconds;
    uint32_t divetime_start_ascend;

    init_buehlmann(&input);
    input.lifeData.max_depth_meter = 0.0;
    input.lifeData.depth_meter = 0.0;
    input.lifeData.temperature_celsius = 22.7;

    ascendrate_seconds = 12.0 / 60.0;
    descendrate_seconds = 20.0 / 60.0;
    divetime_seconds = divetime_minutes * 60;
    divetime_start_ascend = divetime_seconds - (uint32_t)(max_depth_meter / ascendrate_seconds);

    logbook_initNewdiveProfile(&input,&Settings);

    while(input.lifeData.dive_time_seconds < divetime_seconds )
    {
            input.lifeData.dive_time_seconds += 1;

            if(input.lifeData.max_depth_meter < (float)max_depth_meter)
            {
                input.lifeData.depth_meter += descendrate_seconds;
                input.lifeData.max_depth_meter = input.lifeData.depth_meter;
            }
            else if((input.lifeData.dive_time_seconds >= divetime_start_ascend) && (input.lifeData.depth_meter > 0))
            {
                input.lifeData.depth_meter -= ascendrate_seconds;
                if(input.lifeData.depth_meter < 0)
                    input.lifeData.depth_meter = 0;
            }

            logbook_writeSample(input);
    }
    volatile SLogbookHeader* logbookHeader = logbook_getCurrentHeader();
    logbookHeader->total_diveTime_seconds = input.lifeData.dive_time_seconds;
    logbookHeader->maxDepth = input.lifeData.max_depth_meter * 100;
    logbook_EndDive();
}


/**
  ******************************************************************************
  * @brief		test 101
  *						a) for air
    *						b) air + oxygen
    * 					c) Trimix 10/70 + oxygen
    *						65 Meter, 42 Minuten with descent
  * @version 	V0.0.1
  * @date   	26-Oct-2014
    * @retval 	ToDo: 1 for result is similar to MultiDeco
  ******************************************************************************
  */

uint8_t test101_buehlmann_unapproved(void)
{
    /* all the rest */
    SDiveState input;

    init_buehlmann(&input);

    //Gas Change at 6 meter to oxygin
    //input.diveSettings.decogaslist[1].change_during_ascent_depth_bar_otherwise_zero = 0.6f;
    //input.diveSettings.decogaslist[1].nitrogen_percentage = 0;
    //input.diveSettings.decogaslist[1].helium_percentage = 0;
    input.diveSettings.gf_high = 100;
    input.diveSettings.gf_low = 100;
    input.diveSettings.ascentRate_meterperminute = 10.0f;
    input.diveSettings.last_stop_depth_bar = 0.3f; //	input.diveSettings.last_stop_depth_bar = 0.6f; /* ist egal bei oxygen */

    //runter auf 65 meter mit 20 meter/minute
    simulate_descent(&input, 65.0f, 20.0f);
    //38min 45sec saettigung == 2325 sec
     decom_tissues_exposure(30*60, &input.lifeData );
//	 decom_tissues_exposure(2325, &input.lifeData );

  //vpm_calc(&(input.lifeData),&(input.diveSettings),&(input.vpm),&(input.decolistVPM));
    buehlmann_calc_deco(&input.lifeData,&input.diveSettings,&input.decolistBuehlmann);


    //Check time to surface MultiDeco 4.04
    // 308 min with Air
    // 190,5 min with Air + 6m last stop with oxygen
    // 538 min with Trimix 10/70 and oxygen at 6m

    // ...

    return true;
}