Mercurial > public > mk2
comparison code_part1/OSTC_code_c_part2/p2_deco_main - 090915b.c @ 0:96a35aeda5f2
Initial setup
author | heinrichsweikamp |
---|---|
date | Tue, 12 Jan 2010 15:05:59 +0100 |
parents | |
children |
comparison
equal
deleted
inserted
replaced
-1:000000000000 | 0:96a35aeda5f2 |
---|---|
1 /* | |
2 * OSTC - diving computer code | |
3 * =========================== | |
4 * PART 2 : C code | |
5 * | |
6 * p2_deco_main.c for OSTC Mk.2 | |
7 * Created on: 31.08.2009 | |
8 * Author: christian.w @ heinrichsweikamp.com | |
9 * | |
10 * #include <p2_tables.romdata> | |
11 * #include "ostc28.drx.txt" | |
12 * #include "ostc28.tbl.txt" | |
13 * #include "ostc48.tbl.txt" | |
14 * #include "ostc48.drx.txt" | |
15 * #include "ostc90.drx.txt" | |
16 * #include "ostc90.tbl.txt" | |
17 */ | |
18 | |
19 // Copyright (C) 2009 HeinrichsWeikamp GbR | |
20 | |
21 // This program is free software: you can redistribute it and/or modify | |
22 // it under the terms of the GNU General Public License as published by | |
23 // the Free Software Foundation, either version 3 of the License, or | |
24 // (at your option) any later version. | |
25 | |
26 // This program is distributed in the hope that it will be useful, | |
27 // but WITHOUT ANY WARRANTY; without even the implied warranty of | |
28 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | |
29 // GNU General Public License for more details. | |
30 | |
31 // You should have received a copy of the GNU General Public License | |
32 // along with this program. If not, see <http://www.gnu.org/licenses/>. | |
33 | |
34 | |
35 // ***************************** | |
36 // ** I N T R O D U C T I O N ** | |
37 // ***************************** | |
38 // | |
39 // OSTC | |
40 // | |
41 // code: | |
42 // p2_deco_main_c_v101.c | |
43 // part2 of the OSTC code | |
44 // code with constant O2 partial pressure routines | |
45 // experimental | |
46 // | |
47 // summary: | |
48 // decompression routines | |
49 // and display routines | |
50 // for the OSTC experimental project | |
51 // written by Christian Weikamp | |
52 // last revision __________ | |
53 // comments added _________ | |
54 // | |
55 // additional files: | |
56 // #include <p2_tables.romdata> | |
57 // #include "ostc28.drx.txt" | |
58 // #include "ostc28.tbl.txt" | |
59 // #include "ostc48.tbl.txt" | |
60 // #include "ostc48.drx.txt" | |
61 // #include "ostc90.drx.txt" | |
62 // #include "ostc90.tbl.txt" | |
63 // assembler code (PART 1) for working OSTC experimental plattform | |
64 // | |
65 // history: | |
66 // 14/09/09 v201: first alpha candidate for Mk. 2 based on v108 for OSTC | |
67 // | |
68 // literature: | |
69 // B"uhlmann, Albert: Tauchmedizin; 4. Auflage; | |
70 // Schr"oder, Kai & Reith, Steffen; 2000; S"attigungsvorg"ange beim Tauchen, das Modell ZH-L16, Funktionsweise von Tauchcomputern; http://www.achim-und-kai.de/kai/tausim/saett_faq | |
71 // Morrison, Stuart; 2000; DIY DECOMPRESSION; http://www.lizardland.co.uk/DIYDeco.html | |
72 // Balthasar, Steffen; Dekompressionstheorie I: Neo Haldane Modelle; http://www.txfreak.de/dekompressionstheorie_1.pdf | |
73 // Baker, Erik C.; Clearing Up The Confusion About "Deep Stops" | |
74 // Baker, Erik C.; Understanding M-values; http://www.txfreak.de/understanding_m-values.pdf | |
75 | |
76 | |
77 // ********************* | |
78 // ** I N C L U D E S ** | |
79 // ********************* | |
80 #include <p18f4685.h> | |
81 #include <math.h> | |
82 | |
83 // ******************************** | |
84 // ** C O N F I G U R A T I O N ** | |
85 // ** for simulation without asm ** | |
86 // ******************************** | |
87 #pragma config OSC = IRCIO67 | |
88 #pragma config FCMEN = OFF | |
89 #pragma config IESO = OFF | |
90 #pragma config PWRT = ON | |
91 #pragma config BOREN = OFF | |
92 #pragma config WDT = OFF | |
93 #pragma config WDTPS = 128 | |
94 #pragma config MCLRE = ON | |
95 #pragma config LPT1OSC = OFF | |
96 #pragma config PBADEN = OFF | |
97 #pragma config DEBUG = OFF | |
98 #pragma config XINST = OFF | |
99 #pragma config LVP = OFF | |
100 #pragma config STVREN = OFF | |
101 | |
102 // **************************** | |
103 // ** D E F I N E S ** | |
104 // ** missing in p18f4685.h ** | |
105 // **************************** | |
106 #define INT0IF 1 | |
107 #define INT1IF 0 | |
108 #define TMR1IF 0 | |
109 | |
110 | |
111 # define DBG_c_gas 0b0000000000000001 | |
112 # define DBG_c_ppO2 0b0000000000000010 | |
113 # define DBG_RUN 0b0000000000000100 | |
114 # define DBG_RESTART 0b0000000000001000 | |
115 | |
116 # define DBG_CdeSAT 0b0000000000010000 | |
117 # define DBG_C_MODE 0b0000000000100000 | |
118 # define DBG_C_SURF 0b0000000001000000 | |
119 # define DBG_HEwoHE 0b0000000010000000 | |
120 | |
121 # define DBG_C_DPPO2 0b0000000100000000 | |
122 # define DBG_C_DGAS 0b0000001000000000 | |
123 # define DBG_C_DIST 0b0000010000000000 | |
124 # define DBG_C_LAST 0b0000100000000000 | |
125 | |
126 # define DBG_C_GF 0b0001000000000000 | |
127 # define DBG_ZH16ERR 0b0010000000000000 | |
128 # define DBG_PHIGH 0b0100000000000000 | |
129 # define DBG_PLOW 0b1000000000000000 | |
130 | |
131 | |
132 # define DBS_mode 0b0000000000000001 | |
133 # define DBS_ppO2 0b0000000000000010 | |
134 # define DBS_HE_sat 0b0000000000000100 | |
135 # define DBS_ppO2chg 0b0000000000001000 | |
136 | |
137 # define DBS_SAT2l 0b0000000000010000 | |
138 # define DBS_SAT2h 0b0000000000100000 | |
139 # define DBS_GFLOW2l 0b0000000001000000 | |
140 # define DBS_GFLOW2h 0b0000000010000000 | |
141 | |
142 # define DBS_GFHGH2l 0b0000000100000000 | |
143 # define DBS_GFHGH2h 0b0000001000000000 | |
144 # define DBS_GASO22l 0b0000010000000000 | |
145 # define DBS_GASO22h 0b0000100000000000 | |
146 | |
147 # define DBS_DIST2h 0b0001000000000000 | |
148 # define DBS_LAST2h 0b0010000000000000 | |
149 # define DBS_DECOO2l 0b0100000000000000 | |
150 # define DBS_DECOO2h 0b1000000000000000 | |
151 | |
152 | |
153 # define DBS2_PRES2h 0b0000000000000001 | |
154 # define DBS2_PRES2l 0b0000000000000010 | |
155 # define DBS2_SURF2l 0b0000000000000100 | |
156 # define DBS2_SURF2h 0b0000000000001000 | |
157 | |
158 # define DBS2_DESAT2l 0b0000000000010000 | |
159 # define DBS2_DESAT2h 0b0000000000100000 | |
160 # define DBS2_GFDneg 0b0000000001000000 | |
161 # define DBS2_ 0b000000000000000 | |
162 | |
163 # define DBS2_ 0b000000000000000 | |
164 # define DBS2_ 0b000000000000000 | |
165 # define DBS2_ 0b000000000000000 | |
166 # define DBS2_ 0b000000000000000 | |
167 | |
168 // NDL_at_20mtr | |
169 | |
170 | |
171 | |
172 // *********************** | |
173 // ** V A R I A B L E S ** | |
174 // *********************** | |
175 // prefixes etc: | |
176 // _O_ = output for use in the assembler code | |
177 // _I_ = input from the assembler code for the c code | |
178 // char_ and int_ = used to identify output and input size | |
179 // var = variable (from b"uhlmann) | |
180 // pres = pressure | |
181 // gtissue = guiding tissue, the one limiting the ascent | |
182 // e2secs = exp of the b"uhlmann formula precalculated for a 2 second step | |
183 // e1min = same for 1 minute step | |
184 // sim = used in simulating the ascent to the surface | |
185 // nullzeit = remaining ground/bottom time for "no deco" | |
186 // hauptroutine = main | |
187 | |
188 #define WP_FONT_SMALL_HEIGHT 24 | |
189 #define WP_FONT_MEDIUM_HEIGHT 32 | |
190 #define WP_FONT_LARGE_HEIGHT 58 | |
191 | |
192 #define oled_rw PORTA,2,0 | |
193 #define oled_rs PORTE,0,0 | |
194 | |
195 #define U8 unsigned char | |
196 #define U16 unsigned int | |
197 | |
198 // IO assembler | |
199 #pragma udata bank0a=0x060 | |
200 volatile unsigned char wp_stringstore[26]; | |
201 volatile U8 wp_color1; | |
202 volatile U8 wp_color2; | |
203 volatile U8 wp_top; | |
204 volatile U8 wp_leftx2; | |
205 volatile U8 wp_font; | |
206 volatile U8 wp_invert; | |
207 volatile U8 wp_temp_U8; | |
208 // internal C | |
209 #pragma udata bank0b=0x081 | |
210 volatile U8 wp_txtptr; | |
211 volatile unsigned char wp_char; | |
212 volatile U8 wp_command; | |
213 volatile U16 wp_data_16bit; | |
214 volatile U8 wp_data_8bit_one; | |
215 volatile U8 wp_data_8bit_two; | |
216 volatile U16 wp_start; | |
217 volatile U16 wp_end; | |
218 volatile U16 wp_i; | |
219 volatile U8 wp_black; | |
220 // some spare | |
221 volatile U8 wp_debug_U8; | |
222 | |
223 // asm only | |
224 #pragma udata bank0c=0x0D0 | |
225 #define LENGTH_kf_bank0 48 | |
226 volatile unsigned char keep_free_bank0[LENGTH_kf_bank0]; | |
227 | |
228 | |
229 #pragma udata bank1=0x100 | |
230 #define LENGTH_kf_bank1 256 | |
231 volatile unsigned char keep_free_bank1[LENGTH_kf_bank1]; // used by the assembler code | |
232 | |
233 #pragma udata bank2a=0x200 | |
234 // output: | |
235 static unsigned int int_O_tissue_for_debug[32]; | |
236 static unsigned int int_O_GF_spare____; // 0x240 | |
237 static unsigned int int_O_GF_step; // 0x242 | |
238 static unsigned int int_O_gtissue_limit; // 0x244 | |
239 static unsigned int int_O_gtissue_press; // 0x246 | |
240 static unsigned int int_O_limit_GF_low; // 0x248 | |
241 static unsigned int int_O_gtissue_press_at_GF_low; // 0x24A | |
242 // 0x24C + 0x24D noch unbenutzt | |
243 | |
244 #pragma udata bank2b=0x24E | |
245 static unsigned char char_O_GF_low_pointer; // 0x24E | |
246 static unsigned char char_O_actual_pointer; // 0x24F | |
247 #pragma udata bank2c=0x250 | |
248 static unsigned char char_O_deco_table[32]; // 0x250 | |
249 #pragma udata bank2d=0x270 | |
250 static unsigned char char_I_table_deco_done[32]; | |
251 #pragma udata bank2e=0x290 | |
252 static unsigned int int_O_calc_tissue_call_counter; // 0x290 | |
253 // internal: | |
254 unsigned char lock_GF_depth_list; | |
255 static float temp_limit; | |
256 static float GF_low; | |
257 static float GF_high; | |
258 static float GF_delta; | |
259 static float GF_temp; | |
260 static float GF_step; | |
261 static float GF_step2; | |
262 static float temp_pres_gtissue; | |
263 static float temp_pres_gtissue_diff; | |
264 static float temp_pres_gtissue_limit_GF_low; | |
265 static float temp_pres_gtissue_limit_GF_low_below_surface; | |
266 static unsigned int temp_depth_limit; | |
267 static unsigned char temp_decotime; | |
268 static unsigned char temp_gtissue_no; | |
269 static unsigned int temp_depth_last_deco; // new in v.101 | |
270 | |
271 static unsigned char temp_depth_GF_low_meter; | |
272 static unsigned char temp_depth_GF_low_number; | |
273 static unsigned char internal_deco_pointer; | |
274 static unsigned char internal_deco_table[32]; | |
275 static float temp_pres_deco_GF_low; | |
276 | |
277 static unsigned int debug_temp; | |
278 | |
279 | |
280 #pragma udata bank3a=0x300 | |
281 static char output[32]; | |
282 | |
283 #pragma udata bank3b=0x380 | |
284 volatile float pres_tissue_vault[32]; | |
285 | |
286 #pragma udata bank4a=0x400 | |
287 // internal: | |
288 unsigned char ci ; // don't move - used in _asm routines - if moved then modify movlb commands | |
289 unsigned char x; | |
290 unsigned int main_i; | |
291 unsigned int int_temp; | |
292 unsigned int int_temp_decostatus; | |
293 static float pres_respiration; | |
294 static float pres_surface; | |
295 static float temp1; | |
296 static float temp2; | |
297 static float temp3; | |
298 static float temp4; | |
299 static float temp_deco; | |
300 static float temp_atem; | |
301 static float temp2_atem; | |
302 static float temp_tissue; | |
303 static float temp_surface; | |
304 static float N2_ratio; | |
305 static float He_ratio; | |
306 static float temp_ratio; | |
307 static float var_a; | |
308 static float var2_a; | |
309 static float var_b; | |
310 static float var2_b; | |
311 static float var_t05nc; | |
312 static float var2_t05nc; | |
313 static float var_e2secs; | |
314 static float var2_e2secs; | |
315 static float var_e1min; | |
316 static float var2_e1min; | |
317 static float var_halftimes; | |
318 static float var2_halftimes; | |
319 static float pres_gtissue_limit; | |
320 static float temp_pres_gtissue_limit; | |
321 static float actual_ppO2; // new in v.102 | |
322 | |
323 #pragma udata bank4b=0x480 | |
324 static float pres_tissue[32]; | |
325 | |
326 #pragma udata bank5=0x500 | |
327 // don't move positions in this bank, the registers are addressed directly from assembler code | |
328 // input: | |
329 static unsigned int int_I_pres_respiration; // 0x500 | |
330 static unsigned int int_I_pres_surface; // 0x502 | |
331 static unsigned int int_I_temp; // 0x504 new in v101 | |
332 static unsigned char char_I_temp; // 0x506 new in v101 | |
333 static unsigned char char_I_actual_ppO2; // 0x507 | |
334 static unsigned int int_I_spare_3; | |
335 static unsigned int int_I_spare_4; | |
336 static unsigned int int_I_spare_5; | |
337 static unsigned int int_I_spare_6; | |
338 static unsigned char char_I_N2_ratio; // 0x510 | |
339 static unsigned char char_I_He_ratio; // 0x511 | |
340 static unsigned char char_I_saturation_multiplier; // for conservatism/safety values 1.0 (no conservatism) to 1.5 (50% faster saturation | |
341 static unsigned char char_I_desaturation_multiplier; // for conservatism/safety values 0.66 (50% slower desaturation) to 1.0 (no conservatism)// consveratism used in calc_tissue(), calc_tissue_step_1_min() and sim_tissue_1min() | |
342 static unsigned char char_I_GF_High_percentage; // 0x514 new in v.102 | |
343 static unsigned char char_I_GF_Low_percentage; // 0x515 new in v.102 | |
344 static unsigned char char_I_spare; // 0x516 | |
345 static unsigned char char_I_deco_distance; // 0x517 | |
346 static unsigned char char_I_const_ppO2; // 0x518 new in v.101 | |
347 static unsigned char char_I_deco_ppO2_change; // 0x519 new in v.101 | |
348 static unsigned char char_I_deco_ppO2; // 0x51A new in v.101 | |
349 static unsigned char char_I_deco_gas_change; // 0x51B new in v.101 | |
350 static unsigned char char_I_deco_N2_ratio; // 0x51C new in v.101 | |
351 static unsigned char char_I_deco_He_ratio; // 0x51D new in v.101 | |
352 static unsigned char char_I_depth_last_deco; // 0x51E new in v.101 unit: [m] | |
353 static unsigned char char_I_deco_model; // 0x51F new in v.102 ( 1 = MultiGraF, sonst Std. mit (de-)saturation_multiplier) | |
354 // output: | |
355 static unsigned int int_O_desaturation_time; // 0x520 | |
356 static unsigned char char_O_nullzeit; // 0x522 | |
357 static unsigned char char_O_deco_status; // 0x523 | |
358 static unsigned char char_O_array_decotime[7]; // 0x524 | |
359 static unsigned char char_O_array_decodepth[6]; // 0x52B | |
360 static unsigned char char_O_ascenttime; // 0x531 | |
361 static unsigned char char_O_gradient_factor; // 0x532 | |
362 static unsigned char char_O_tissue_saturation[32]; // 0x533 | |
363 static unsigned char char_O_array_gradient_weighted[16]; // 0x553 | |
364 static unsigned char char_O_gtissue_no; // 0x563 | |
365 static unsigned char char_O_diluent; // 0x564 new in v.101 | |
366 static unsigned char char_O_CNS_fraction; // 0x565 new in v.101 | |
367 static unsigned char char_O_relative_gradient_GF; // 0x566 new in v.102 | |
368 | |
369 // internal: | |
370 static float pres_tissue_limit[16]; | |
371 static float sim_pres_tissue_limit[16]; | |
372 static float pres_diluent; // new in v.101 | |
373 static float deco_diluent; // new in v.101 | |
374 static float const_ppO2; // new in v.101 | |
375 static float deco_ppO2_change; // new in v.101 | |
376 static float deco_ppO2; // new in v.101 | |
377 | |
378 #pragma udata bank6=0x600 | |
379 // internal: | |
380 static float sim_pres_tissue[32]; | |
381 static float sim_pres_tissue_backup[32]; | |
382 | |
383 //#pragma udata bank7=0x700 | |
384 //const unsigned char keep_free_bank7[256]; // used by the assembler code | |
385 | |
386 #pragma udata bank8=0x800 | |
387 static char md_pi_subst[256]; | |
388 | |
389 #pragma udata bank9a=0x900 | |
390 // output: | |
391 static char md_state[48]; // DONT MOVE !! // has to be at the beginning of bank 9 for the asm code!!! | |
392 #pragma udata bank9b=0x930 | |
393 // output: | |
394 static unsigned int int_O_DBS_bitfield; // 0x930 new in v.108 | |
395 static unsigned int int_O_DBS2_bitfield; // 0x932 new in v.108 | |
396 static unsigned int int_O_DBG_pre_bitfield; // 0x934 new in v.108 | |
397 static unsigned int int_O_DBG_post_bitfield; // 0x936 new in v.108 | |
398 static char char_O_NDL_at_20mtr; // 0x938 new in v.108 // 0xFF == undefined, max. 254 | |
399 // internal: | |
400 static char md_t; | |
401 static char md_buffer[16]; | |
402 static char md_cksum[16]; | |
403 static char md_i; | |
404 static char md_j; | |
405 static char md_temp; | |
406 static unsigned int md_pointer; | |
407 static float deco_N2_ratio; // new in v.101 | |
408 static float deco_He_ratio; // new in v.101 | |
409 static float calc_N2_ratio; // new in v.101 | |
410 static float calc_He_ratio; // new in v.101 | |
411 static float deco_gas_change; // new in v.101 | |
412 static float CNS_fraction; // new in v.101 | |
413 static float float_saturation_multiplier; // new in v.101 | |
414 static float float_desaturation_multiplier; // new in v.101 | |
415 static float float_deco_distance; // new in v.101 | |
416 // internal, dbg: | |
417 static unsigned char DBG_char_I_deco_model; // new in v.108 | |
418 static unsigned char DBG_char_I_depth_last_deco; // new in v.108 | |
419 static float DBG_pres_surface; // new in v.108 | |
420 static float DBG_GF_low; // new in v.108 | |
421 static float DBG_GF_high; // new in v.108 | |
422 static float DBG_const_ppO2; // new in v.108 | |
423 static float DBG_deco_ppO2_change; // new in v.108 | |
424 static float DBG_deco_ppO2; // new in v.108 | |
425 static float DBG_deco_N2_ratio; // new in v.108 | |
426 static float DBG_deco_He_ratio; // new in v.108 | |
427 static float DBG_deco_gas_change; // new in v.108 | |
428 static float DBG_float_saturation_multiplier; // new in v.108 | |
429 static float DBG_float_desaturation_multiplier; // new in v.108 | |
430 static float DBG_float_deco_distance; // new in v.108 | |
431 static float DBG_deco_N2_ratio; // new in v.108 | |
432 static float DBG_deco_He_ratio; // new in v.108 | |
433 static float DBG_N2_ratio; // new in v.108 | |
434 static float DBG_He_ratio; // new in v.108 | |
435 static char flag_in_divemode; // new in v.108 | |
436 static int int_dbg_i; // new in v.108 | |
437 unsigned int temp_DBS; | |
438 | |
439 // ************************* | |
440 // ** P R O T O T Y P E S ** | |
441 // ************************* | |
442 void main_calc_hauptroutine(void); | |
443 void main_calc_without_deco(void); | |
444 void main_clear_tissue(void); | |
445 void main_calc_percentage(void); | |
446 void main_calc_wo_deco_step_1_min(void); | |
447 void main_debug(void); | |
448 void main_gradient_array(void); | |
449 void main_hash(void); | |
450 | |
451 void calc_hauptroutine(void); | |
452 void calc_tissue(void); | |
453 void calc_nullzeit(void); | |
454 void backup_sim_pres_tissue(void); | |
455 void restore_sim_pres_tissue(void); | |
456 | |
457 void calc_without_deco(void); | |
458 void clear_tissue(void); | |
459 void calc_ascenttime(void); | |
460 void update_startvalues(void); | |
461 void clear_decoarray(void); | |
462 void update_decoarray(void); | |
463 void sim_tissue_1min(void); | |
464 void sim_tissue_10min(void); | |
465 void calc_gradient_factor(void); | |
466 void calc_gradient_array_only(void); | |
467 void calc_desaturation_time(void); | |
468 void calc_wo_deco_step_1_min(void); | |
469 void calc_tissue_step_1_min(void); | |
470 void hash(void); | |
471 void clear_CNS_fraction(void); | |
472 void calc_CNS_fraction(void); | |
473 void calc_CNS_decrease_15min(void); | |
474 void calc_percentage(void); | |
475 void main(void); | |
476 void calc_hauptroutine_data_input(void); | |
477 void calc_hauptroutine_update_tissues(void); | |
478 void calc_hauptroutine_calc_deco(void); | |
479 void calc_hauptroutine_calc_ascend_to_deco(void); | |
480 void calc_nextdecodepth_GF(void); | |
481 void copy_deco_table_GF(void); | |
482 void clear_internal_deco_table_GF(void); | |
483 void update_internal_deco_table_GF(void); | |
484 void push_tissues_to_vault(void); | |
485 void pull_tissues_from_vault(void); | |
486 void main_push_tissues_to_vault(void); | |
487 void main_pull_tissues_from_vault(void); | |
488 void wordprocessor(void); | |
489 | |
490 // ******************************* | |
491 // ** start ** | |
492 // ** necessary for compilation ** | |
493 // ******************************* | |
494 #pragma romdata der_code = 0x0000 | |
495 #pragma code der_start = 0x0000 | |
496 void der_start(void) | |
497 { | |
498 _asm | |
499 goto main | |
500 _endasm | |
501 } | |
502 | |
503 // *********************************** | |
504 // ** main code for simulation / ** | |
505 // ** tests without assembler code ** | |
506 // ** is NOT a part of the OSTC ** | |
507 // *********************************** | |
508 #pragma code main = 0x9000 | |
509 void main(void) | |
510 { | |
511 for(wp_temp_U8=0;wp_temp_U8<LENGTH_kf_bank0 - 1;wp_temp_U8++) | |
512 keep_free_bank0[wp_temp_U8] = 7; | |
513 keep_free_bank0[LENGTH_kf_bank0 - 1] = 7; | |
514 | |
515 for(wp_temp_U8=0;wp_temp_U8<LENGTH_kf_bank1 - 1;wp_temp_U8++) | |
516 keep_free_bank1[wp_temp_U8] = 7; | |
517 keep_free_bank1[LENGTH_kf_bank1 - 1] = 7; | |
518 | |
519 #if 1 | |
520 // new main to test DR-5 | |
521 | |
522 wp_top = 10; | |
523 wp_leftx2 = 10; | |
524 wp_color1 = 255; | |
525 wp_color2 = 255; | |
526 wp_font = 0; | |
527 wp_invert = 0; | |
528 wp_stringstore[0] = ' '; | |
529 wp_stringstore[1] = '!'; | |
530 wp_stringstore[2] = '"'; | |
531 wp_stringstore[3] = ':'; | |
532 wp_stringstore[4] = 0; | |
533 wordprocessor(); | |
534 | |
535 GF_low = 1.0; | |
536 GF_high = 1.0; | |
537 | |
538 GF_temp = GF_low * GF_high; | |
539 | |
540 clear_CNS_fraction(); | |
541 //char_I_const_ppO2 = 100; | |
542 //for (main_i=0;main_i<255;main_i++) | |
543 //{ | |
544 //calc_CNS_fraction(); | |
545 //} //for | |
546 | |
547 | |
548 | |
549 | |
550 int_I_pres_respiration = 1000;//980; | |
551 int_I_pres_surface = 1000;//980; | |
552 char_I_N2_ratio = 39; //38; | |
553 char_I_He_ratio = 40; //50; | |
554 char_I_deco_distance = 0; // 10 = 1 meter | |
555 char_I_depth_last_deco = 3; // values below 3 (meter) are ignored | |
556 | |
557 char_I_const_ppO2 = 0; | |
558 char_I_deco_ppO2_change = 0; // [dm] 10 = 1 meter | |
559 char_I_deco_ppO2 = 0; | |
560 | |
561 char_I_deco_gas_change = 0; // [m] 1 = 1 meter | |
562 char_I_deco_N2_ratio = 0; | |
563 char_I_deco_He_ratio = 0; | |
564 | |
565 //char_I_actual_ppO2; // 0x507 | |
566 char_I_GF_High_percentage = 100; // 0x514 new in v.102 | |
567 char_I_GF_Low_percentage = 100; // 0x515 new in v.102 | |
568 | |
569 char_I_saturation_multiplier = 110; | |
570 char_I_desaturation_multiplier = 90; | |
571 | |
572 char_I_deco_model = 0; | |
573 | |
574 main_clear_tissue(); | |
575 | |
576 int_I_pres_respiration = 1000 + int_I_pres_surface; | |
577 main_calc_wo_deco_step_1_min(); | |
578 int_I_pres_respiration = 3000 + int_I_pres_surface; | |
579 main_calc_wo_deco_step_1_min(); | |
580 int_I_pres_respiration = 5000 + int_I_pres_surface; | |
581 main_calc_wo_deco_step_1_min(); | |
582 | |
583 /* | |
584 int_I_pres_respiration = 6000 + int_I_pres_surface; | |
585 for (main_i=0;main_i<27;main_i++) | |
586 main_calc_wo_deco_step_1_min(); | |
587 */ | |
588 | |
589 char_O_deco_status = 255; | |
590 while (char_O_deco_status) | |
591 main_calc_hauptroutine(); | |
592 _asm | |
593 nop | |
594 _endasm | |
595 | |
596 for (main_i=0;main_i<50;main_i++) | |
597 { | |
598 main_calc_hauptroutine(); | |
599 } | |
600 int_I_pres_respiration = 10000; | |
601 for (main_i=0;main_i<1500;main_i++) | |
602 { | |
603 main_calc_hauptroutine(); | |
604 } | |
605 | |
606 _asm | |
607 nop | |
608 _endasm | |
609 | |
610 | |
611 int_I_pres_respiration = 3000; | |
612 for (main_i=0;main_i<150;main_i++) | |
613 { | |
614 calc_hauptroutine_data_input(); | |
615 calc_hauptroutine_update_tissues(); | |
616 } //for | |
617 | |
618 update_startvalues(); | |
619 clear_decoarray(); | |
620 clear_internal_deco_table_GF(); | |
621 calc_hauptroutine_calc_ascend_to_deco(); | |
622 if (char_O_deco_status > 15) // can't go up to first deco, too deep to calculate in the given time slot | |
623 { | |
624 char_O_deco_status = 2; | |
625 // char_O_lock_depth_list = 255; | |
626 } | |
627 else | |
628 { | |
629 // char_O_lock_depth_list = lock_GF_depth_list; | |
630 calc_hauptroutine_calc_deco(); | |
631 } | |
632 // build_debug_output(); | |
633 | |
634 _asm | |
635 nop | |
636 _endasm | |
637 while (char_O_deco_status == 1) | |
638 { | |
639 char_O_deco_status = 0; | |
640 // char_O_lock_depth_list = 255; | |
641 calc_hauptroutine_calc_deco(); | |
642 // build_debug_output(); | |
643 _asm | |
644 nop | |
645 _endasm | |
646 }; | |
647 debug_temp = 60; // [mtr Aufstieg in 10 mtr/min (30steps'2sec/min] | |
648 int_I_pres_respiration = 9980; | |
649 for (main_i=0;main_i<debug_temp;main_i++) | |
650 { | |
651 int_I_pres_respiration = int_I_pres_respiration - 33; | |
652 calc_hauptroutine_data_input(); | |
653 calc_hauptroutine_update_tissues(); | |
654 int_I_pres_respiration = int_I_pres_respiration - 33; | |
655 calc_hauptroutine_data_input(); | |
656 calc_hauptroutine_update_tissues(); | |
657 int_I_pres_respiration = int_I_pres_respiration - 34; | |
658 calc_hauptroutine_data_input(); | |
659 calc_hauptroutine_update_tissues(); | |
660 } //for | |
661 _asm | |
662 nop | |
663 _endasm | |
664 | |
665 update_startvalues(); | |
666 clear_decoarray(); | |
667 clear_internal_deco_table_GF(); | |
668 calc_hauptroutine_calc_ascend_to_deco(); | |
669 if (char_O_deco_status > 15) // can't go up to first deco, too deep to calculate in the given time slot | |
670 { | |
671 char_O_deco_status = 2; | |
672 // char_O_lock_depth_list = 255; | |
673 } | |
674 else | |
675 { | |
676 // char_O_lock_depth_list = lock_GF_depth_list; | |
677 calc_hauptroutine_calc_deco(); | |
678 } | |
679 // build_debug_output(); | |
680 | |
681 _asm | |
682 nop | |
683 _endasm | |
684 while (char_O_deco_status == 1) | |
685 { | |
686 char_O_deco_status = 0; | |
687 // char_O_lock_depth_list = 255; | |
688 calc_hauptroutine_calc_deco(); | |
689 // build_debug_output(); | |
690 _asm | |
691 nop | |
692 _endasm | |
693 }; | |
694 _asm | |
695 nop | |
696 _endasm | |
697 debug_temp = 60; // [mtr Aufstieg in 10 mtr/min (30steps'2sec/min] | |
698 int_I_pres_respiration = 9980; | |
699 debug_temp = debug_temp * 3; | |
700 for (main_i=0;main_i<debug_temp;main_i++) | |
701 { | |
702 calc_hauptroutine_data_input(); | |
703 calc_hauptroutine_update_tissues(); | |
704 } //for | |
705 _asm | |
706 nop | |
707 _endasm | |
708 #endif | |
709 // ----------------------- | |
710 | |
711 } // main | |
712 | |
713 // ****************************************************** | |
714 // ****************************************************** | |
715 // ** THE FOLLOWING CODE HAS TO BE COPPIED TO THE OSTC ** | |
716 // ****************************************************** | |
717 // ****************************************************** | |
718 | |
719 // *************** | |
720 // *************** | |
721 // ** THE FONTS ** | |
722 // *************** | |
723 // *************** | |
724 // all new for bigscreen | |
725 | |
726 // a bit is left (C4) | |
727 #pragma romdata font_data_small = 0x0A0C4 | |
728 rom const rom U16 wp_small_data[] = | |
729 { | |
730 #include "ostc28.drx.txt" | |
731 }; | |
732 | |
733 #pragma romdata font_table_small = 0x0B548 | |
734 rom const rom U16 wp_small_table[] = | |
735 { | |
736 #include "ostc28.tbl.txt" // length 0xEC | |
737 }; | |
738 | |
739 #pragma romdata font_table_medium = 0x0B632 | |
740 rom const rom U16 wp_medium_table[] = | |
741 { | |
742 #include "ostc48.tbl.txt" // length 0x22 | |
743 }; | |
744 | |
745 #pragma romdata font_data_medium = 0x0B656 | |
746 rom const rom U16 wp_medium_data[] = | |
747 { | |
748 #include "ostc48.drx.txt" // length 0x390 | |
749 }; | |
750 | |
751 #pragma romdata font_data_large = 0x0BA46 | |
752 rom const rom U16 wp_large_data[] = | |
753 { | |
754 #include "ostc90.drx.txt" | |
755 }; | |
756 | |
757 #pragma romdata font_table_large = 0x0BFE6 | |
758 rom const rom U16 wp_large_table[] = | |
759 { | |
760 #include "ostc90.tbl.txt" | |
761 }; | |
762 | |
763 // *********************** | |
764 // *********************** | |
765 // ** THE SUBROUTINES 2 ** | |
766 // *********************** | |
767 // *********************** | |
768 // all new in v.102 | |
769 // moved from 0x0D000 to 0x0C000 in v.108 | |
770 | |
771 #pragma code subroutines2 = 0x0C000 // can be adapted to fit the romdata tables ahead | |
772 | |
773 // ------------------------------- | |
774 // DBS - debug on start of dive // | |
775 // ------------------------------- | |
776 void create_dbs_set_dbg_and_ndl20mtr(void) | |
777 { | |
778 int_O_DBS_bitfield = 0; | |
779 int_O_DBS2_bitfield = 0; | |
780 if(int_O_DBG_pre_bitfield & DBG_RUN) | |
781 int_O_DBG_pre_bitfield = DBG_RESTART; | |
782 else | |
783 int_O_DBG_pre_bitfield = DBG_RUN; | |
784 int_O_DBG_post_bitfield = 0; | |
785 char_O_NDL_at_20mtr = 255; | |
786 | |
787 DBG_N2_ratio = N2_ratio; | |
788 DBG_He_ratio = He_ratio; | |
789 DBG_char_I_deco_model = char_I_deco_model; | |
790 DBG_char_I_depth_last_deco = char_I_depth_last_deco; | |
791 DBG_pres_surface = pres_surface; | |
792 DBG_GF_low = GF_low; | |
793 DBG_GF_high = GF_high; | |
794 DBG_const_ppO2 = const_ppO2; | |
795 DBG_deco_ppO2_change = deco_ppO2_change; | |
796 DBG_deco_ppO2 = deco_ppO2; | |
797 DBG_deco_N2_ratio = deco_N2_ratio; | |
798 DBG_deco_He_ratio = deco_He_ratio; | |
799 DBG_deco_gas_change = deco_gas_change; | |
800 DBG_float_saturation_multiplier = float_saturation_multiplier; | |
801 DBG_float_desaturation_multiplier = float_desaturation_multiplier; | |
802 DBG_float_deco_distance = float_deco_distance; | |
803 | |
804 if(char_I_deco_model) | |
805 int_O_DBS_bitfield |= DBS_mode; | |
806 if(const_ppO2) | |
807 int_O_DBS_bitfield |= DBS_ppO2; | |
808 for(int_dbg_i = 16; int_dbg_i < 32; int_dbg_i++) | |
809 if(pres_tissue[int_dbg_i]) | |
810 int_O_DBS_bitfield |= DBS_HE_sat; | |
811 if(deco_ppO2_change) | |
812 int_O_DBS_bitfield |= DBS_ppO2chg; | |
813 if(float_saturation_multiplier < 0.99) | |
814 int_O_DBS_bitfield |= DBS_SAT2l; | |
815 if(float_saturation_multiplier > 1.3) | |
816 int_O_DBS_bitfield |= DBS_SAT2h; | |
817 if(GF_low < 0.19) | |
818 int_O_DBS_bitfield |= DBS_GFLOW2l; | |
819 if(GF_low > 1.01) | |
820 int_O_DBS_bitfield |= DBS_GFLOW2h; | |
821 if(GF_high < 0.6) | |
822 int_O_DBS_bitfield |= DBS_GFHGH2l; | |
823 if(GF_high > 1.01) | |
824 int_O_DBS_bitfield |= DBS_GFHGH2h; | |
825 if((N2_ratio + He_ratio) > 0.95) | |
826 int_O_DBS_bitfield |= DBS_GASO22l; | |
827 if((N2_ratio + He_ratio) < 0.05) | |
828 int_O_DBS_bitfield |= DBS_GASO22h; | |
829 if(float_deco_distance > 0.25) | |
830 int_O_DBS_bitfield |= DBS_DIST2h; | |
831 if(char_I_depth_last_deco > 8) | |
832 int_O_DBS_bitfield |= DBS_LAST2h; | |
833 if(DBG_deco_gas_change && ((deco_N2_ratio + deco_He_ratio) > 0.95)) | |
834 int_O_DBS_bitfield |= DBS_DECOO2l; | |
835 if(DBG_deco_gas_change && ((deco_N2_ratio + deco_He_ratio) < 0.05)) | |
836 int_O_DBS_bitfield |= DBS_DECOO2h; | |
837 if(pres_respiration > 3.0) | |
838 int_O_DBS2_bitfield |= DBS2_PRES2h; | |
839 if(pres_surface - pres_respiration > 0.2) | |
840 int_O_DBS2_bitfield |= DBS2_PRES2l; | |
841 if(pres_surface < 0.75) | |
842 int_O_DBS2_bitfield |= DBS2_SURF2l; | |
843 if(pres_surface > 1.11) | |
844 int_O_DBS2_bitfield |= DBS2_SURF2h; | |
845 if(float_desaturation_multiplier < 0.70) | |
846 int_O_DBS2_bitfield |= DBS2_DESAT2l; | |
847 if(float_desaturation_multiplier > 1.01) | |
848 int_O_DBS2_bitfield |= DBS2_DESAT2h; | |
849 if(GF_low > GF_high) | |
850 int_O_DBS2_bitfield |= DBS2_GFDneg; | |
851 } | |
852 | |
853 // ------------------------------- | |
854 // DBG - set DBG to end_of_dive // | |
855 // ------------------------------- | |
856 void set_dbg_end_of_dive(void) | |
857 { | |
858 int_O_DBG_pre_bitfield &= (~DBG_RUN); | |
859 int_O_DBG_post_bitfield &= (~DBG_RUN); | |
860 } | |
861 | |
862 // ------------------------------- | |
863 // DBG - NDL at first 20 m. hit // | |
864 // ------------------------------- | |
865 void check_ndl(void) | |
866 { | |
867 if((char_O_NDL_at_20mtr == -1) && (int_I_pres_respiration > 3000)) | |
868 { | |
869 char_O_NDL_at_20mtr = char_O_nullzeit; | |
870 if(char_O_NDL_at_20mtr == 255) | |
871 char_O_NDL_at_20mtr == 254; | |
872 } | |
873 } | |
874 | |
875 // ------------------------------- | |
876 // DBG - multi main during dive // | |
877 // ------------------------------- | |
878 void check_dbg(char is_post_check) | |
879 { | |
880 temp_DBS = 0; | |
881 if( (DBG_N2_ratio != N2_ratio) || (DBG_He_ratio != He_ratio) ) | |
882 temp_DBS |= DBG_c_gas; | |
883 if(DBG_const_ppO2 != const_ppO2) | |
884 temp_DBS |= DBG_c_ppO2; | |
885 if((DBG_float_saturation_multiplier != float_saturation_multiplier) || (DBG_float_desaturation_multiplier != float_desaturation_multiplier)) | |
886 temp_DBS |= DBG_CdeSAT; | |
887 if(DBG_char_I_deco_model != char_I_deco_model) | |
888 temp_DBS |= DBG_C_MODE; | |
889 if(DBG_pres_surface != pres_surface) | |
890 temp_DBS |= DBG_C_SURF; | |
891 if((!DBS_HE_sat) && (!He_ratio)) | |
892 for(int_dbg_i = 16; int_dbg_i < 32; int_dbg_i++) | |
893 if(pres_tissue[int_dbg_i]) | |
894 temp_DBS |= DBG_HEwoHE; | |
895 if(DBG_deco_ppO2 != deco_ppO2) | |
896 temp_DBS |= DBG_C_DPPO2; | |
897 if((DBG_deco_gas_change != deco_gas_change) || (DBG_deco_N2_ratio != deco_N2_ratio) || (DBG_deco_He_ratio != deco_He_ratio)) | |
898 temp_DBS |= DBG_C_DGAS; | |
899 if(DBG_float_deco_distance != float_deco_distance) | |
900 temp_DBS |= DBG_C_DIST; | |
901 if(DBG_char_I_depth_last_deco != char_I_depth_last_deco) | |
902 temp_DBS |= DBG_C_LAST; | |
903 if((DBG_GF_low != GF_low) || (DBG_GF_high != GF_high)) | |
904 temp_DBS |= DBG_C_GF; | |
905 if(pres_respiration > 13.0) | |
906 temp_DBS |= DBG_PHIGH; | |
907 if(pres_surface - pres_respiration > 0.2) | |
908 temp_DBS |= DBG_PLOW; | |
909 /* | |
910 if() | |
911 temp_DBS |= ; | |
912 if() | |
913 temp_DBS |= ; | |
914 */ | |
915 if(is_post_check) | |
916 int_O_DBG_post_bitfield |= temp_DBS; | |
917 else | |
918 int_O_DBG_pre_bitfield |= temp_DBS; | |
919 } | |
920 | |
921 // ------------------------------- | |
922 // DBG - prior to calc. of dive // | |
923 // ------------------------------- | |
924 void check_pre_dbg(void) | |
925 { | |
926 check_dbg(0); | |
927 } | |
928 | |
929 // ------------------------------- | |
930 // DBG - after decocalc of dive // | |
931 // ------------------------------- | |
932 void check_post_dbg(void) | |
933 { | |
934 check_dbg(1); | |
935 } | |
936 | |
937 | |
938 | |
939 // ------------------------- | |
940 // calc_next_decodepth_GF // | |
941 // ------------------------- | |
942 // new in v.102 | |
943 void calc_nextdecodepth_GF(void) | |
944 { | |
945 // INPUT, changing during dive: | |
946 // temp_pres_gtissue_limit_GF_low | |
947 // temp_pres_gtissue_limit_GF_low_below_surface | |
948 // temp_pres_gtissue | |
949 // temp_pres_gtissue_diff | |
950 // lock_GF_depth_list | |
951 | |
952 // INPUT, fixed during dive: | |
953 // pres_surface | |
954 // GF_delta | |
955 // GF_high | |
956 // GF_low | |
957 // temp_depth_last_deco | |
958 // float_deco_distance | |
959 | |
960 // OUTPUT | |
961 // GF_step | |
962 // temp_deco | |
963 // temp_depth_limt | |
964 // lock_GF_depth_list | |
965 | |
966 // USES | |
967 // temp1 | |
968 // temp2 | |
969 // int_temp | |
970 | |
971 char_I_table_deco_done[0] = 0; // safety if changed somewhere else. Needed for exit | |
972 if (char_I_deco_model == 1) | |
973 { | |
974 if (lock_GF_depth_list == 0) | |
975 { | |
976 temp2 = temp_pres_gtissue_limit_GF_low_below_surface / 0.29985; // = ... / 99.95 / 0.003; | |
977 int_temp = (int) (temp2 + 0.99); | |
978 if (int_temp > 31) | |
979 int_temp = 31; // deepest deco at 93 meter (31 deco stops) | |
980 if (int_temp < 0) | |
981 int_temp = 0; | |
982 temp_depth_GF_low_number = int_temp; | |
983 temp_depth_GF_low_meter = 3 * temp_depth_GF_low_number; | |
984 temp2 = (float)temp_depth_GF_low_meter * 0.09995; | |
985 temp_pres_deco_GF_low = temp2 + float_deco_distance + pres_surface; | |
986 if (temp_depth_GF_low_number == 0) | |
987 GF_step = 0; | |
988 else | |
989 GF_step = GF_delta / (float)temp_depth_GF_low_number; | |
990 if (GF_step < 0) | |
991 GF_step = 0; | |
992 if (GF_step > GF_delta) | |
993 GF_step = GF_delta; | |
994 int_O_GF_step = (int)(GF_step * 10000); | |
995 int_O_limit_GF_low = (int)(temp_pres_deco_GF_low * 1000); | |
996 int_O_gtissue_press_at_GF_low = (int)(temp_pres_gtissue * 1000); | |
997 char_O_GF_low_pointer = temp_depth_GF_low_number; | |
998 lock_GF_depth_list = 1; | |
999 internal_deco_pointer = 0; | |
1000 } | |
1001 if (internal_deco_pointer == 0) // new run | |
1002 { | |
1003 internal_deco_pointer = temp_depth_GF_low_number; | |
1004 GF_temp = GF_high - ((float)internal_deco_pointer * GF_step); | |
1005 int_temp = char_I_table_deco_done[internal_deco_pointer]; | |
1006 output[8] = int_temp; | |
1007 output[9] = 33; | |
1008 } | |
1009 else | |
1010 { | |
1011 int_temp = 1; | |
1012 } | |
1013 while (int_temp == 1) | |
1014 { | |
1015 int_temp = internal_deco_pointer - 1; | |
1016 if (int_temp == 1) // new in v104 | |
1017 { | |
1018 temp2 = (float)(temp_depth_last_deco * int_temp) * 0.09995; | |
1019 GF_step2 = GF_step/3.0 * ((float)(6 - temp_depth_last_deco)); | |
1020 } | |
1021 else | |
1022 if (int_temp == 0) | |
1023 { | |
1024 temp2 = 0.0; | |
1025 GF_step2 = GF_high - GF_temp; | |
1026 } | |
1027 else | |
1028 { | |
1029 temp2 = (float)(3 *int_temp) * 0.09995; | |
1030 GF_step2 = GF_step; | |
1031 } | |
1032 temp2 = temp2 + pres_surface; // next deco stop to be tested | |
1033 temp1 = ((GF_temp + GF_step2)* temp_pres_gtissue_diff) + temp_pres_gtissue; // upper limit (lowest pressure allowed) // changes GF_step2 in v104 | |
1034 if (temp1 > temp2) // check if ascent to next deco stop is ok | |
1035 { | |
1036 int_temp = 0; // no | |
1037 } | |
1038 else | |
1039 { | |
1040 internal_deco_pointer = int_temp; | |
1041 GF_temp = GF_temp + GF_step2; // changed in v104 | |
1042 int_temp = char_I_table_deco_done[internal_deco_pointer]; // yes and check for ascent to even next stop if deco_done is set | |
1043 } | |
1044 } // while | |
1045 if (internal_deco_pointer > 0) | |
1046 { | |
1047 temp2 = (float)(0.29985 * internal_deco_pointer); | |
1048 temp_deco = temp2 + float_deco_distance + pres_surface; | |
1049 if (internal_deco_pointer == 1) // new in v104 | |
1050 temp_depth_limit = temp_depth_last_deco; | |
1051 else | |
1052 temp_depth_limit = 3 * internal_deco_pointer; | |
1053 if (output[9] == 33) | |
1054 { | |
1055 output[9] = internal_deco_pointer; | |
1056 output[10] = char_I_table_deco_done[internal_deco_pointer]; | |
1057 output[12] = output[12] + 1; | |
1058 if (output[12] == 100) | |
1059 output[12] = 0; | |
1060 } | |
1061 } | |
1062 else // if (char_I_deco_model == 1) | |
1063 { | |
1064 temp_deco = pres_surface; | |
1065 temp_depth_limit = 0; | |
1066 } | |
1067 } | |
1068 else | |
1069 { | |
1070 // calc_nextdecodepth - original | |
1071 // optimized in v.101 | |
1072 // depth_last_deco included in v.101 | |
1073 | |
1074 temp1 = temp_pres_gtissue_limit - pres_surface; | |
1075 if (temp1 >= 0) | |
1076 { | |
1077 temp1 = temp1 / 0.29985; // = temp1 / 99.95 / 0.003; | |
1078 temp_depth_limit = (int) (temp1 + 0.99); | |
1079 temp_depth_limit = 3 * temp_depth_limit; // depth for deco [m] | |
1080 if (temp_depth_limit == 0) | |
1081 temp_deco = pres_surface; | |
1082 else | |
1083 { | |
1084 if (temp_depth_limit < temp_depth_last_deco) | |
1085 temp_depth_limit = temp_depth_last_deco; | |
1086 temp1 = (float)temp_depth_limit * 0.09995; | |
1087 temp_deco = temp1 + float_deco_distance + pres_surface; // depth for deco [bar] | |
1088 } // if (temp_depth_limit == 0) | |
1089 } // if (temp1 >= 0) | |
1090 else | |
1091 { | |
1092 temp_deco = pres_surface; | |
1093 temp_depth_limit = 0; | |
1094 } // if (temp1 >= 0) | |
1095 } // calc_nextdecodepth original | |
1096 } // calc_nextdecodepth_GF | |
1097 | |
1098 | |
1099 #if 0 | |
1100 void build_debug_output(void) | |
1101 { | |
1102 output[0] = 0; // not used in asm PLED output | |
1103 output[1] = (int) (GF_low * 100); | |
1104 output[2] = (int) (GF_high * 100); | |
1105 output[3] = (int) (GF_step * 100); | |
1106 output[4] = (int) temp_depth_GF_low_number; | |
1107 output[5] = (int) temp_depth_GF_low_meter; | |
1108 //output[6] | |
1109 output[7] = (int) internal_deco_pointer; | |
1110 //output[8] = char_I_table_deco_done[temp_depth_GF_low_number] | |
1111 //output[9] = internal_deco_pointer @ new run | |
1112 //output[10] = char_I_table_deco_done[internal_deco_pointer] @ new run | |
1113 output [11] = (int) (temp_pres_deco_GF_low * 10); | |
1114 } // build_debug_output | |
1115 #endif | |
1116 | |
1117 // --------------------- | |
1118 // copy_deco_table_GF // | |
1119 // --------------------- | |
1120 // new in v.102 | |
1121 void copy_deco_table_GF(void) | |
1122 { | |
1123 if (char_I_deco_model == 1) | |
1124 { | |
1125 int_temp = 32; | |
1126 for (ci=0;ci<int_temp;ci++) | |
1127 char_O_deco_table[ci] = internal_deco_table[ci]; | |
1128 } | |
1129 } // copy_deco_table_GF | |
1130 | |
1131 | |
1132 // ------------------------------ | |
1133 // clear_internal_deco_table_GF// | |
1134 // ------------------------------ | |
1135 // new in v.102 | |
1136 void clear_internal_deco_table_GF(void) | |
1137 { | |
1138 if (char_I_deco_model == 1) | |
1139 { | |
1140 for (ci=0;ci<32;ci++) // cycle through the 16 b"uhlmann tissues for Helium | |
1141 { | |
1142 internal_deco_table[ci] = 0; | |
1143 } | |
1144 } | |
1145 } // clear_internal_deco_table_GF | |
1146 | |
1147 | |
1148 // -------------------------------- | |
1149 // update_internal_deco_table_GF // | |
1150 // -------------------------------- | |
1151 // new in v.102 | |
1152 void update_internal_deco_table_GF(void) | |
1153 { | |
1154 if ((char_I_deco_model == 1) && (internal_deco_table[internal_deco_pointer] < 255)) | |
1155 internal_deco_table[internal_deco_pointer] = internal_deco_table[internal_deco_pointer] + 1; | |
1156 } // update_internal_deco_table_GF | |
1157 | |
1158 | |
1159 // --------------------- | |
1160 // temp_tissue_safety // | |
1161 // --------------------- | |
1162 // outsourced in v.102 | |
1163 void temp_tissue_safety(void) | |
1164 { | |
1165 if (char_I_deco_model == 1) | |
1166 { | |
1167 } | |
1168 else | |
1169 { | |
1170 if (temp_tissue < 0.0) | |
1171 temp_tissue = temp_tissue * float_desaturation_multiplier; | |
1172 else | |
1173 temp_tissue = temp_tissue * float_saturation_multiplier; | |
1174 } | |
1175 } // temp_tissue_safety | |
1176 | |
1177 | |
1178 // ********************** | |
1179 // ********************** | |
1180 // ** THE JUMP-IN CODE ** | |
1181 // ** for the asm code ** | |
1182 // ********************** | |
1183 // ********************** | |
1184 #pragma code main_calc_hauptroutine = 0x10000 | |
1185 void main_calc_hauptroutine(void) | |
1186 { | |
1187 calc_hauptroutine(); | |
1188 int_O_desaturation_time = 65535; | |
1189 } // divemode | |
1190 #pragma code main_without_deco = 0x10020 | |
1191 void main_calc_without_deco(void) | |
1192 { | |
1193 calc_without_deco(); | |
1194 calc_desaturation_time(); | |
1195 } | |
1196 | |
1197 #pragma code main_clear_CNS_fraction = 0x10030 | |
1198 void main_clear_CNS_fraction(void) | |
1199 { | |
1200 clear_CNS_fraction(); | |
1201 } | |
1202 | |
1203 #pragma code main_calc_CNS_decrease_15min = 0x10034 | |
1204 void main_calc_CNS_decrease_15min(void) | |
1205 { | |
1206 calc_CNS_decrease_15min(); | |
1207 } | |
1208 | |
1209 #pragma code main_calc_percentage = 0x10038 | |
1210 void main_calc_percentage (void) | |
1211 { | |
1212 calc_percentage(); | |
1213 } | |
1214 | |
1215 #pragma code main_clear_tissue = 0x10040 | |
1216 void main_clear_tissue(void) | |
1217 { | |
1218 clear_tissue(); | |
1219 char_I_depth_last_deco = 0; // for compatibility with v.101pre_no_last_deco | |
1220 } | |
1221 | |
1222 #pragma code main_calc_CNS_fraction = 0x10050 | |
1223 void main_calc_CNS_fraction(void) | |
1224 { | |
1225 calc_CNS_fraction(); | |
1226 } | |
1227 | |
1228 #pragma code main_calc_desaturation_time = 0x10060 | |
1229 void main_calc_desaturation_time(void) | |
1230 { | |
1231 calc_desaturation_time(); | |
1232 } | |
1233 | |
1234 #pragma code main_calc_wo_deco_step_1_min = 0x10080 | |
1235 void main_calc_wo_deco_step_1_min(void) | |
1236 { | |
1237 calc_wo_deco_step_1_min(); | |
1238 char_O_deco_status = 3; // surface new in v.102 overwrites value of calc_wo_deco_step_1_min | |
1239 calc_desaturation_time(); | |
1240 } // surface mode | |
1241 | |
1242 #pragma code main_debug = 0x100A0 | |
1243 void main_debug(void) | |
1244 { | |
1245 //debug(); | |
1246 } | |
1247 | |
1248 #pragma code main_DD2_write_incon42 = 0x100B0 | |
1249 void main_DD2_write_incon42(void) | |
1250 { | |
1251 return; | |
1252 } | |
1253 | |
1254 #pragma code main_DD2_write_incon24 = 0x100B4 | |
1255 void main_DD2_write_incon24(void) | |
1256 { | |
1257 return; | |
1258 } | |
1259 | |
1260 #pragma code main_wordprocessor = 0x100B8 | |
1261 void main_wordprocessor(void) | |
1262 { | |
1263 wordprocessor(); | |
1264 } | |
1265 | |
1266 #pragma code main_gradient_array = 0x100C0 | |
1267 void main_gradient_array(void) | |
1268 { | |
1269 calc_gradient_array_only(); | |
1270 } | |
1271 #pragma code main_push_tissues = 0x100C4 | |
1272 void main_push_tissues_to_vault(void) | |
1273 { | |
1274 push_tissues_to_vault(); | |
1275 } | |
1276 #pragma code main_pull_tissues = 0x100C8 | |
1277 void main_pull_tissues_from_vault(void) | |
1278 { | |
1279 pull_tissues_from_vault(); | |
1280 } | |
1281 | |
1282 #pragma code main_hash = 0x100E0 | |
1283 void main_hash(void) | |
1284 { | |
1285 hash(); | |
1286 } | |
1287 | |
1288 // *********************** | |
1289 // *********************** | |
1290 // ** THE LOOKUP TABLES ** | |
1291 // *********************** | |
1292 // *********************** | |
1293 | |
1294 #pragma romdata tables = 0x10200 | |
1295 #include <p2_tables.romdata> // new table for deco_main_v.101 (var_a modified) | |
1296 | |
1297 #pragma romdata tables2 = 0x10600 | |
1298 rom const rom unsigned int md_pi[] = | |
1299 { | |
1300 0x292E, 0x43C9, 0xA2D8, 0x7C01, 0x3D36, 0x54A1, 0xECF0, 0x0613 | |
1301 , 0x62A7, 0x05F3, 0xC0C7, 0x738C, 0x9893, 0x2BD9, 0xBC4C, 0x82CA | |
1302 , 0x1E9B, 0x573C, 0xFDD4, 0xE016, 0x6742, 0x6F18, 0x8A17, 0xE512 | |
1303 , 0xBE4E, 0xC4D6, 0xDA9E, 0xDE49, 0xA0FB, 0xF58E, 0xBB2F, 0xEE7A | |
1304 , 0xA968, 0x7991, 0x15B2, 0x073F, 0x94C2, 0x1089, 0x0B22, 0x5F21 | |
1305 , 0x807F, 0x5D9A, 0x5A90, 0x3227, 0x353E, 0xCCE7, 0xBFF7, 0x9703 | |
1306 , 0xFF19, 0x30B3, 0x48A5, 0xB5D1, 0xD75E, 0x922A, 0xAC56, 0xAAC6 | |
1307 , 0x4FB8, 0x38D2, 0x96A4, 0x7DB6, 0x76FC, 0x6BE2, 0x9C74, 0x04F1 | |
1308 , 0x459D, 0x7059, 0x6471, 0x8720, 0x865B, 0xCF65, 0xE62D, 0xA802 | |
1309 , 0x1B60, 0x25AD, 0xAEB0, 0xB9F6, 0x1C46, 0x6169, 0x3440, 0x7E0F | |
1310 , 0x5547, 0xA323, 0xDD51, 0xAF3A, 0xC35C, 0xF9CE, 0xBAC5, 0xEA26 | |
1311 , 0x2C53, 0x0D6E, 0x8528, 0x8409, 0xD3DF, 0xCDF4, 0x4181, 0x4D52 | |
1312 , 0x6ADC, 0x37C8, 0x6CC1, 0xABFA, 0x24E1, 0x7B08, 0x0CBD, 0xB14A | |
1313 , 0x7888, 0x958B, 0xE363, 0xE86D, 0xE9CB, 0xD5FE, 0x3B00, 0x1D39 | |
1314 , 0xF2EF, 0xB70E, 0x6658, 0xD0E4, 0xA677, 0x72F8, 0xEB75, 0x4B0A | |
1315 , 0x3144, 0x50B4, 0x8FED, 0x1F1A, 0xDB99, 0x8D33, 0x9F11, 0x8314 | |
1316 }; | |
1317 | |
1318 // ********************* | |
1319 // ********************* | |
1320 // ** THE SUBROUTINES ** | |
1321 // ********************* | |
1322 // ********************* | |
1323 | |
1324 #pragma code subroutines = 0x10700 // can be adapted to fit the romdata tables ahead | |
1325 | |
1326 | |
1327 // --------------- | |
1328 // CLEAR tissue // | |
1329 // --------------- | |
1330 // optimized in v.101 (var_a) | |
1331 | |
1332 void clear_tissue(void) // preload tissues with standard pressure for the given ambient pressure | |
1333 { | |
1334 | |
1335 flag_in_divemode = 0; | |
1336 int_O_DBS_bitfield = 0; | |
1337 int_O_DBS2_bitfield = 0; | |
1338 int_O_DBG_pre_bitfield = 0; | |
1339 int_O_DBG_post_bitfield = 0; | |
1340 char_O_NDL_at_20mtr = 255; | |
1341 | |
1342 _asm | |
1343 lfsr 1, 0x300 // C math routines shall use this variable bank | |
1344 movlw 0x01 | |
1345 movwf TBLPTRU,0 | |
1346 _endasm | |
1347 | |
1348 // N2_ratio = (float)char_I_N2_ratio; // the 0.0002 of 0.7902 are missing with standard air | |
1349 N2_ratio = 0.7902; // N2_ratio / 100.0; | |
1350 pres_respiration = (float)int_I_pres_respiration / 1000.0; | |
1351 for (ci=0;ci<16;ci++) // cycle through the 16 b"uhlmann tissues | |
1352 { | |
1353 pres_tissue[ci] = N2_ratio * (pres_respiration - 0.0627) ; | |
1354 _asm | |
1355 movlw 0x02 | |
1356 movwf TBLPTRH,0 | |
1357 movlb 4 // fuer ci | |
1358 movf ci,0,1 | |
1359 addwf ci,0,1 | |
1360 addwf ci,0,1 | |
1361 addwf ci,0,1 | |
1362 addlw 0x80 | |
1363 movwf TBLPTRL,0 | |
1364 TBLRDPOSTINC | |
1365 movff TABLAT,var_a+1 | |
1366 TBLRDPOSTINC | |
1367 movff TABLAT,var_a | |
1368 TBLRDPOSTINC | |
1369 movff TABLAT,var_a+3 | |
1370 TBLRD | |
1371 movff TABLAT,var_a+2 | |
1372 addlw 0x80 | |
1373 movwf TBLPTRL,0 | |
1374 incf TBLPTRH,1,0 | |
1375 TBLRDPOSTINC | |
1376 movff TABLAT,var_b+1 | |
1377 TBLRDPOSTINC | |
1378 movff TABLAT,var_b | |
1379 TBLRDPOSTINC | |
1380 movff TABLAT,var_b+3 | |
1381 TBLRD | |
1382 movff TABLAT,var_b+2 | |
1383 _endasm | |
1384 | |
1385 pres_tissue_limit[ci] = (pres_tissue[ci] - var_a) * var_b ; | |
1386 // now update the guiding tissue | |
1387 if (pres_tissue_limit[ci] < 0) | |
1388 pres_tissue_limit[ci] = 0; | |
1389 } // for 0 to 16 | |
1390 | |
1391 for (ci=16;ci<32;ci++) // cycle through the 16 b"uhlmann tissues for Helium | |
1392 { | |
1393 pres_tissue[ci] = 0.0; | |
1394 } // for | |
1395 | |
1396 clear_decoarray(); | |
1397 char_O_deco_status = 0; | |
1398 char_O_nullzeit = 0; | |
1399 char_O_ascenttime = 0; | |
1400 char_O_gradient_factor = 0; | |
1401 char_O_relative_gradient_GF = 0; | |
1402 } // clear_tissue(void) | |
1403 | |
1404 | |
1405 // -------------------- | |
1406 // calc_without_deco // | |
1407 // fixed N2_ratio ! // | |
1408 // -------------------- | |
1409 // optimized in v.101 (float_..saturation_multiplier) | |
1410 | |
1411 void calc_without_deco(void) | |
1412 { | |
1413 _asm | |
1414 lfsr 1, 0x300 | |
1415 _endasm | |
1416 N2_ratio = 0.7902; // FIXED RATIO !! sum as stated in b"uhlmann | |
1417 pres_respiration = (float)int_I_pres_respiration / 1000.0; // assembler code uses different digit system | |
1418 pres_surface = (float)int_I_pres_surface / 1000.0; | |
1419 temp_atem = N2_ratio * (pres_respiration - 0.0627); // 0.0627 is the extra pressure in the body | |
1420 temp2_atem = 0.0; | |
1421 temp_surface = pres_surface; // the b"uhlmann formula using temp_surface does apply to the pressure without any inert ratio | |
1422 float_desaturation_multiplier = char_I_desaturation_multiplier / 100.0; | |
1423 float_saturation_multiplier = char_I_saturation_multiplier / 100.0; | |
1424 | |
1425 calc_tissue(); // update the pressure in the 16 tissues in accordance with the new ambient pressure | |
1426 | |
1427 clear_decoarray(); | |
1428 char_O_deco_status = 0; | |
1429 char_O_nullzeit = 0; | |
1430 char_O_ascenttime = 0; | |
1431 calc_gradient_factor(); | |
1432 | |
1433 } // calc_without_deco | |
1434 | |
1435 | |
1436 // -------------------- | |
1437 // calc_hauptroutine // | |
1438 // -------------------- | |
1439 // this is the major code in dive mode | |
1440 // calculates: | |
1441 // the tissues, | |
1442 // the bottom time | |
1443 // and simulates the ascend with all deco stops | |
1444 | |
1445 void calc_hauptroutine(void) | |
1446 { | |
1447 calc_hauptroutine_data_input(); | |
1448 | |
1449 if(!flag_in_divemode) | |
1450 { | |
1451 flag_in_divemode = 1; | |
1452 create_dbs_set_dbg_and_ndl20mtr(); | |
1453 } | |
1454 else | |
1455 check_pre_dbg(); | |
1456 | |
1457 calc_hauptroutine_update_tissues(); | |
1458 calc_gradient_factor(); | |
1459 | |
1460 | |
1461 switch (char_O_deco_status) // toggle between calculation for nullzeit (bottom time), deco stops and more deco stops (continue) | |
1462 { | |
1463 case 0: | |
1464 update_startvalues(); | |
1465 calc_nullzeit(); | |
1466 check_ndl(); | |
1467 char_O_deco_status = 255; // calc deco next time | |
1468 break; | |
1469 case 1: | |
1470 if (char_O_deco_status == 3) | |
1471 break; | |
1472 char_O_deco_status = 0; | |
1473 // char_O_lock_depth_list = 255; | |
1474 calc_hauptroutine_calc_deco(); | |
1475 // build_debug_output(); | |
1476 break; | |
1477 case 3: // new dive | |
1478 clear_decoarray(); | |
1479 clear_internal_deco_table_GF(); | |
1480 copy_deco_table_GF(); | |
1481 internal_deco_pointer = 0; | |
1482 lock_GF_depth_list = 0; | |
1483 update_startvalues(); | |
1484 calc_nextdecodepth_GF(); | |
1485 char_O_deco_status = 0; | |
1486 break; | |
1487 default: | |
1488 update_startvalues(); | |
1489 clear_decoarray(); | |
1490 clear_internal_deco_table_GF(); | |
1491 output[6] = 1; | |
1492 calc_hauptroutine_calc_ascend_to_deco(); | |
1493 if (char_O_deco_status > 15) // can't go up to first deco, too deep to calculate in the given time slot | |
1494 { | |
1495 char_O_deco_status = 2; | |
1496 // char_O_lock_depth_list = 255; | |
1497 } | |
1498 else | |
1499 { | |
1500 // char_O_lock_depth_list = lock_GF_depth_list; | |
1501 calc_hauptroutine_calc_deco(); | |
1502 } | |
1503 // build_debug_output(); | |
1504 break; | |
1505 } | |
1506 calc_ascenttime(); | |
1507 check_post_dbg(); | |
1508 } | |
1509 | |
1510 void calc_hauptroutine_data_input(void) | |
1511 { | |
1512 pres_respiration = (float)int_I_pres_respiration / 1000.0; | |
1513 pres_surface = (float)int_I_pres_surface / 1000.0; | |
1514 | |
1515 N2_ratio = (float)char_I_N2_ratio / 100.0;; // the 0.0002 of 0.7902 are missing with standard air | |
1516 He_ratio = (float)char_I_He_ratio / 100.0;; | |
1517 deco_N2_ratio = (float)char_I_deco_N2_ratio / 100.0; | |
1518 deco_He_ratio = (float)char_I_deco_He_ratio / 100.0; | |
1519 float_deco_distance = (float)char_I_deco_distance / 100.0; | |
1520 if(char_I_deco_gas_change) | |
1521 { | |
1522 deco_gas_change = (float)char_I_deco_gas_change / 9.995 + pres_surface; | |
1523 deco_gas_change = deco_gas_change + float_deco_distance; | |
1524 } | |
1525 else | |
1526 deco_gas_change = 0; | |
1527 const_ppO2 = (float)char_I_const_ppO2 / 100.0; | |
1528 deco_ppO2_change = (float)char_I_deco_ppO2_change / 99.95 + pres_surface; | |
1529 deco_ppO2_change = deco_ppO2_change + float_deco_distance; | |
1530 deco_ppO2 = (float)char_I_deco_ppO2 / 100.0; | |
1531 float_desaturation_multiplier = char_I_desaturation_multiplier / 100.0; | |
1532 float_saturation_multiplier = char_I_saturation_multiplier / 100.0; | |
1533 GF_low = (float)char_I_GF_Low_percentage / 100.0; | |
1534 GF_high = (float)char_I_GF_High_percentage / 100.0; | |
1535 GF_delta = GF_high - GF_low; | |
1536 | |
1537 temp2 = (pres_respiration - pres_surface) / 0.29985; | |
1538 int_temp = (int)(temp2); | |
1539 if (int_temp < 0) | |
1540 int_temp = 0; | |
1541 if (int_temp > 255) | |
1542 int_temp = 255; | |
1543 char_O_actual_pointer = int_temp; | |
1544 | |
1545 temp_depth_last_deco = (int)char_I_depth_last_deco; | |
1546 } | |
1547 | |
1548 void calc_hauptroutine_update_tissues(void) | |
1549 { | |
1550 int_O_calc_tissue_call_counter = int_O_calc_tissue_call_counter + 1; | |
1551 if (char_I_const_ppO2 == 0) // new in v.101 | |
1552 pres_diluent = pres_respiration; // new in v.101 | |
1553 else // new in v.101 | |
1554 pres_diluent = ((pres_respiration - const_ppO2)/(N2_ratio + He_ratio)); // new in v.101 | |
1555 if (pres_diluent > pres_respiration) // new in v.101 | |
1556 pres_diluent = pres_respiration; // new in v.101 | |
1557 if (pres_diluent > 0.0627) // new in v.101 | |
1558 { | |
1559 temp_atem = N2_ratio * (pres_diluent - 0.0627); // changed in v.101 | |
1560 temp2_atem = He_ratio * (pres_diluent - 0.0627); // changed in v.101 | |
1561 char_O_diluent = (char)(pres_diluent/pres_respiration*100.0); | |
1562 } | |
1563 else // new in v.101 | |
1564 { | |
1565 temp_atem = 0.0; // new in v.101 | |
1566 temp2_atem = 0.0; // new in v.101 | |
1567 char_O_diluent = 0; | |
1568 } | |
1569 temp_surface = pres_surface; | |
1570 calc_tissue(); | |
1571 int_O_gtissue_limit = (int)(pres_tissue_limit[char_O_gtissue_no] * 1000); | |
1572 int_O_gtissue_press = (int)((pres_tissue[char_O_gtissue_no] + pres_tissue[char_O_gtissue_no+16]) * 1000); | |
1573 if (char_I_deco_model == 1) | |
1574 { | |
1575 temp1 = temp1 * GF_high; | |
1576 } | |
1577 else | |
1578 { | |
1579 temp1 = temp_surface; | |
1580 } | |
1581 if (pres_gtissue_limit > temp1 && char_O_deco_status == 0) // if guiding tissue can not be exposed to surface pressure immediately | |
1582 { | |
1583 char_O_nullzeit = 0; // deco necessary | |
1584 char_O_deco_status = 255; // calculate deco skip nullzeit calculation | |
1585 } | |
1586 } // calc_hauptroutine_update_tissues | |
1587 void calc_hauptroutine_calc_deco(void) | |
1588 { | |
1589 do | |
1590 { | |
1591 int_temp_decostatus = 0; | |
1592 calc_nextdecodepth_GF(); | |
1593 if (temp_depth_limit > 0) | |
1594 { | |
1595 if (char_I_const_ppO2 == 0) // new in v.101 | |
1596 { | |
1597 deco_diluent = temp_deco; // new in v.101 | |
1598 if (temp_deco > deco_gas_change) | |
1599 { | |
1600 calc_N2_ratio = N2_ratio; | |
1601 calc_He_ratio = He_ratio; | |
1602 } | |
1603 else | |
1604 { | |
1605 calc_N2_ratio = deco_N2_ratio; | |
1606 calc_He_ratio = deco_He_ratio; | |
1607 } | |
1608 } | |
1609 else // new in v.101 | |
1610 { | |
1611 calc_N2_ratio = N2_ratio; | |
1612 calc_He_ratio = He_ratio; | |
1613 if (temp_deco > deco_ppO2_change) | |
1614 { | |
1615 deco_diluent = ((temp_deco - const_ppO2)/(N2_ratio + He_ratio)); // new in v.101 | |
1616 } | |
1617 else | |
1618 { | |
1619 deco_diluent = ((temp_deco - deco_ppO2)/(N2_ratio + He_ratio)); // new in v.101 | |
1620 } | |
1621 } | |
1622 if (deco_diluent > temp_deco) // new in v.101 | |
1623 deco_diluent = temp_deco; // new in v.101 | |
1624 if (deco_diluent > 0.0627) // new in v.101 | |
1625 { | |
1626 temp_atem = calc_N2_ratio * (deco_diluent - 0.0627); // changed in v.101 | |
1627 temp2_atem = calc_He_ratio * (deco_diluent - 0.0627); // changed in v.101 | |
1628 } | |
1629 else // new in v.101 | |
1630 { | |
1631 temp_atem = 0.0; // new in v.101 | |
1632 temp2_atem = 0.0; // new in v.101 | |
1633 } | |
1634 sim_tissue_1min(); | |
1635 update_internal_deco_table_GF(); | |
1636 temp_decotime = 1; | |
1637 update_decoarray(); | |
1638 char_O_deco_status = char_O_deco_status + 1; | |
1639 if (char_O_deco_status < 16) | |
1640 int_temp_decostatus = 1; | |
1641 } | |
1642 else // if (temp_depth_limit > 0) | |
1643 { | |
1644 char_O_deco_status = 0; | |
1645 } | |
1646 } while (int_temp_decostatus == 1); | |
1647 if (char_O_deco_status > 15) | |
1648 { | |
1649 char_O_deco_status = 1; | |
1650 } | |
1651 else | |
1652 { | |
1653 copy_deco_table_GF(); | |
1654 char_O_deco_status = 0; | |
1655 } | |
1656 } | |
1657 | |
1658 void calc_hauptroutine_calc_ascend_to_deco(void) | |
1659 { | |
1660 update_startvalues(); | |
1661 char_O_deco_status = 0; | |
1662 temp_deco = pres_respiration; | |
1663 lock_GF_depth_list = 1; // new in v.102 | |
1664 do // go up to first deco | |
1665 { | |
1666 int_temp_decostatus = 0; | |
1667 temp_deco = temp_deco - 1.0; | |
1668 if ( char_I_deco_model == 1) // new in v.102 , 4 = deep stops | |
1669 temp_limit = temp_pres_gtissue_limit_GF_low; | |
1670 else | |
1671 temp_limit = temp_pres_gtissue_limit; | |
1672 if ((temp_deco > temp_limit) && (temp_deco > pres_surface)) // changes in v.102 | |
1673 { | |
1674 lock_GF_depth_list = 0; // new in v.102, distance to first stop > 10 mtr. | |
1675 output[6] = 0; | |
1676 if (char_I_const_ppO2 == 0) // new in v.101 // calculate at half of the ascent | |
1677 { | |
1678 deco_diluent = temp_deco + 0.5; // new in v.101 | |
1679 if (temp_deco + 0.5 > deco_gas_change) | |
1680 { | |
1681 calc_N2_ratio = N2_ratio; | |
1682 calc_He_ratio = He_ratio; | |
1683 } | |
1684 else | |
1685 { | |
1686 calc_N2_ratio = deco_N2_ratio; | |
1687 calc_He_ratio = deco_He_ratio; | |
1688 } | |
1689 } | |
1690 else // new in v.101 | |
1691 { | |
1692 calc_N2_ratio = N2_ratio; | |
1693 calc_He_ratio = He_ratio; | |
1694 if (temp_deco + 0.5 > deco_ppO2_change) | |
1695 deco_diluent = ((temp_deco + 0.5 - const_ppO2)/(N2_ratio + He_ratio)); // new in v.101 // calculate at half of the ascent | |
1696 else | |
1697 deco_diluent = ((temp_deco + 0.5 - deco_ppO2)/(N2_ratio + He_ratio)); // new in v.101 // calculate at half of the ascent | |
1698 if (deco_diluent > (temp_deco +0.5)) // new in v.101 | |
1699 deco_diluent = temp_deco + 0.5; // new in v.101 // calculate at half of the ascent | |
1700 } | |
1701 if (deco_diluent > 0.0627) // new in v.101 | |
1702 { | |
1703 temp_atem = calc_N2_ratio * (deco_diluent - 0.0627); // changed in v.101 | |
1704 temp2_atem = calc_He_ratio * (deco_diluent - 0.0627); // changed in v.101 | |
1705 } | |
1706 else // new in v.101 | |
1707 { | |
1708 temp_atem = 0.0; // new in v.101 | |
1709 temp2_atem = 0.0; // new in v.101 | |
1710 } | |
1711 sim_tissue_1min(); | |
1712 char_O_deco_status = char_O_deco_status + 1; | |
1713 if (char_O_deco_status < 16) // 16 is the limit of calculations for one time slot | |
1714 int_temp_decostatus = 1; | |
1715 } | |
1716 } while (int_temp_decostatus == 1); | |
1717 } // calc_hauptroutine_calc_ascend_to_deco | |
1718 | |
1719 // -------------- | |
1720 // calc_tissue // | |
1721 // -------------- | |
1722 // optimized in v.101 | |
1723 | |
1724 void calc_tissue(void) | |
1725 { | |
1726 _asm | |
1727 lfsr 1, 0x300 | |
1728 movlw 0x01 | |
1729 movwf TBLPTRU,0 | |
1730 _endasm | |
1731 | |
1732 char_O_gtissue_no = 255; | |
1733 pres_gtissue_limit = 0.0; | |
1734 | |
1735 for (ci=0;ci<16;ci++) | |
1736 { | |
1737 _asm | |
1738 movlw 0x02 | |
1739 movwf TBLPTRH,0 | |
1740 movlb 4 // fuer ci | |
1741 movf ci,0,1 | |
1742 addwf ci,0,1 | |
1743 addwf ci,0,1 | |
1744 addwf ci,0,1 | |
1745 movwf TBLPTRL,0 | |
1746 TBLRDPOSTINC | |
1747 movff TABLAT,var_e2secs+1 // the order is confussing | |
1748 TBLRDPOSTINC | |
1749 movff TABLAT,var_e2secs // low byte first, high afterwards | |
1750 TBLRDPOSTINC | |
1751 movff TABLAT,var_e2secs+3 | |
1752 TBLRD | |
1753 movff TABLAT,var_e2secs+2 | |
1754 addlw 0x40 | |
1755 movwf TBLPTRL,0 | |
1756 TBLRDPOSTINC | |
1757 movff TABLAT,var2_e2secs+1 | |
1758 TBLRDPOSTINC | |
1759 movff TABLAT,var2_e2secs | |
1760 TBLRDPOSTINC | |
1761 movff TABLAT,var2_e2secs+3 | |
1762 TBLRD | |
1763 movff TABLAT,var2_e2secs+2 | |
1764 addlw 0x40 | |
1765 movwf TBLPTRL,0 | |
1766 TBLRDPOSTINC | |
1767 movff TABLAT,var_a+1 | |
1768 TBLRDPOSTINC | |
1769 movff TABLAT,var_a | |
1770 TBLRDPOSTINC | |
1771 movff TABLAT,var_a+3 | |
1772 TBLRD | |
1773 movff TABLAT,var_a+2 | |
1774 addlw 0x40 | |
1775 movwf TBLPTRL,0 | |
1776 TBLRDPOSTINC | |
1777 movff TABLAT,var2_a+1 | |
1778 TBLRDPOSTINC | |
1779 movff TABLAT,var2_a | |
1780 TBLRDPOSTINC | |
1781 movff TABLAT,var2_a+3 | |
1782 TBLRD | |
1783 movff TABLAT,var2_a+2 | |
1784 addlw 0x40 | |
1785 movwf TBLPTRL,0 | |
1786 incf TBLPTRH,1,0 | |
1787 TBLRDPOSTINC | |
1788 movff TABLAT,var_b+1 | |
1789 TBLRDPOSTINC | |
1790 movff TABLAT,var_b | |
1791 TBLRDPOSTINC | |
1792 movff TABLAT,var_b+3 | |
1793 TBLRD | |
1794 movff TABLAT,var_b+2 | |
1795 addlw 0x40 | |
1796 movwf TBLPTRL,0 | |
1797 TBLRDPOSTINC | |
1798 movff TABLAT,var2_b+1 | |
1799 TBLRDPOSTINC | |
1800 movff TABLAT,var2_b | |
1801 TBLRDPOSTINC | |
1802 movff TABLAT,var2_b+3 | |
1803 TBLRD | |
1804 movff TABLAT,var2_b+2 | |
1805 _endasm | |
1806 // the start values are the previous end values // write new values in temp | |
1807 | |
1808 if( (var_e2secs < 0.0000363) | |
1809 || (var_e2secs > 0.00577) | |
1810 || (var2_e2secs < 0.0000961) | |
1811 || (var2_e2secs > 0.150) | |
1812 || (var_a < 0.231) | |
1813 || (var_a > 1.27) | |
1814 || (var_b < 0.504) | |
1815 || (var_b > 0.966) | |
1816 || (var2_a < 0.510) | |
1817 || (var2_a > 1.75) | |
1818 || (var2_b < 0.423) | |
1819 || (var2_b > 0.927) | |
1820 ) | |
1821 int_O_DBG_pre_bitfield |= DBG_ZH16ERR; | |
1822 | |
1823 // N2 | |
1824 temp_tissue = (temp_atem - pres_tissue[ci]) * var_e2secs; | |
1825 temp_tissue_safety(); | |
1826 pres_tissue[ci] = pres_tissue[ci] + temp_tissue; | |
1827 | |
1828 // He | |
1829 temp_tissue = (temp2_atem - pres_tissue[ci+16]) * var2_e2secs; | |
1830 temp_tissue_safety(); | |
1831 pres_tissue[ci+16] = pres_tissue[ci+16] + temp_tissue; | |
1832 | |
1833 temp_tissue = pres_tissue[ci] + pres_tissue[ci+16]; | |
1834 | |
1835 var_a = (var_a * pres_tissue[ci] + var2_a * pres_tissue[ci+16]) / temp_tissue; | |
1836 var_b = (var_b * pres_tissue[ci] + var2_b * pres_tissue[ci+16]) / temp_tissue; | |
1837 pres_tissue_limit[ci] = (temp_tissue - var_a) * var_b; | |
1838 if (pres_tissue_limit[ci] < 0) | |
1839 pres_tissue_limit[ci] = 0; | |
1840 if (pres_tissue_limit[ci] > pres_gtissue_limit) | |
1841 { | |
1842 pres_gtissue_limit = pres_tissue_limit[ci]; | |
1843 char_O_gtissue_no = ci; | |
1844 }//if | |
1845 } // for | |
1846 }//calc_tissue(void) | |
1847 | |
1848 // ---------------- | |
1849 // calc_nullzeit // | |
1850 // ---------------- | |
1851 // calculates the remaining bottom time | |
1852 | |
1853 // unchanged in v.101 | |
1854 | |
1855 void calc_nullzeit(void) | |
1856 { | |
1857 char_O_nullzeit = 0; | |
1858 int_temp = 1; | |
1859 do | |
1860 { | |
1861 backup_sim_pres_tissue(); | |
1862 sim_tissue_10min(); | |
1863 char_O_nullzeit = char_O_nullzeit + 10; | |
1864 int_temp = int_temp + 1; | |
1865 if (char_I_deco_model == 1) | |
1866 temp1 = GF_high * temp_pres_gtissue_diff + temp_pres_gtissue; | |
1867 else | |
1868 temp1 = temp_pres_gtissue_limit; | |
1869 if (temp1 > temp_surface) // changed in v.102 , if guiding tissue can not be exposed to surface pressure immediately | |
1870 int_temp = 255; | |
1871 } while (int_temp < 17); | |
1872 if (int_temp == 255) | |
1873 { | |
1874 restore_sim_pres_tissue(); | |
1875 char_O_nullzeit = char_O_nullzeit - 10; | |
1876 } //if int_temp == 255] | |
1877 int_temp = 1; | |
1878 if (char_O_nullzeit < 60) | |
1879 { | |
1880 do | |
1881 { | |
1882 sim_tissue_1min(); | |
1883 char_O_nullzeit = char_O_nullzeit + 1; | |
1884 int_temp = int_temp + 1; // new in v.102a | |
1885 if (char_I_deco_model == 1) | |
1886 temp1 = GF_high * temp_pres_gtissue_diff + temp_pres_gtissue; | |
1887 else | |
1888 temp1 = temp_pres_gtissue_limit; | |
1889 if (temp1 > temp_surface) // changed in v.102 , if guiding tissue can not be exposed to surface pressure immediately | |
1890 int_temp = 255; | |
1891 } while (int_temp < 10); | |
1892 if (int_temp == 255) | |
1893 char_O_nullzeit = char_O_nullzeit - 1; | |
1894 } // if char_O_nullzeit < 60 | |
1895 } //calc_nullzeit | |
1896 | |
1897 // ------------------------- | |
1898 // backup_sim_pres_tissue // | |
1899 // ------------------------- | |
1900 void backup_sim_pres_tissue(void) | |
1901 { | |
1902 for (x = 0;x<16;x++) | |
1903 { | |
1904 sim_pres_tissue_backup[x] = sim_pres_tissue[x]; | |
1905 sim_pres_tissue_backup[x+16] = sim_pres_tissue[x+16]; | |
1906 } | |
1907 } // backup_sim | |
1908 | |
1909 // -------------------------- | |
1910 // restore_sim_pres_tissue // | |
1911 // -------------------------- | |
1912 void restore_sim_pres_tissue(void) | |
1913 { | |
1914 for (x = 0;x<16;x++) | |
1915 { | |
1916 sim_pres_tissue[x] = sim_pres_tissue_backup[x]; | |
1917 sim_pres_tissue[x+16] = sim_pres_tissue_backup[x+16]; | |
1918 } | |
1919 } // restore_sim | |
1920 | |
1921 // ------------------ | |
1922 // calc_ascenttime // | |
1923 // ------------------ | |
1924 | |
1925 void calc_ascenttime(void) | |
1926 { | |
1927 if (pres_respiration > pres_surface) | |
1928 { | |
1929 switch (char_O_deco_status) | |
1930 { | |
1931 case 2: | |
1932 char_O_ascenttime = 255; | |
1933 break; | |
1934 case 1: | |
1935 break; | |
1936 default: | |
1937 temp1 = pres_respiration - pres_surface + 0.6; // + 0.6 hence 1 minute ascent time from a depth of 4 meter on | |
1938 if (temp1 < 0) | |
1939 temp1 = 0; | |
1940 if (temp1 > 255) | |
1941 temp1 = 255; | |
1942 char_O_ascenttime = (char)temp1; | |
1943 | |
1944 for(ci=0;ci<7;ci++) | |
1945 { | |
1946 x = char_O_ascenttime + char_O_array_decotime[ci]; | |
1947 if (x < char_O_ascenttime) | |
1948 char_O_ascenttime = 255; | |
1949 else | |
1950 char_O_ascenttime = x; | |
1951 } | |
1952 } | |
1953 } | |
1954 else | |
1955 char_O_ascenttime = 0; | |
1956 } // calc_ascenttime() | |
1957 | |
1958 | |
1959 // --------------------- | |
1960 // update_startvalues // | |
1961 // --------------------- | |
1962 // updated in v.102 | |
1963 | |
1964 void update_startvalues(void) | |
1965 { | |
1966 temp_pres_gtissue_limit = pres_gtissue_limit; | |
1967 temp_pres_gtissue = pres_tissue[char_O_gtissue_no] + pres_tissue[char_O_gtissue_no+16]; | |
1968 temp_pres_gtissue_diff = temp_pres_gtissue_limit - temp_pres_gtissue; // negative number | |
1969 temp_pres_gtissue_limit_GF_low = GF_low * temp_pres_gtissue_diff + temp_pres_gtissue; | |
1970 temp_pres_gtissue_limit_GF_low_below_surface = temp_pres_gtissue_limit_GF_low - pres_surface; | |
1971 if (temp_pres_gtissue_limit_GF_low_below_surface < 0) | |
1972 temp_pres_gtissue_limit_GF_low_below_surface = 0; | |
1973 | |
1974 temp_gtissue_no = char_O_gtissue_no; | |
1975 for (x = 0;x<16;x++) | |
1976 { | |
1977 sim_pres_tissue[x] = pres_tissue[x]; | |
1978 sim_pres_tissue[x+16] = pres_tissue[x+16]; | |
1979 sim_pres_tissue_limit[x] = pres_tissue_limit[x]; | |
1980 } | |
1981 } // update_startvalues | |
1982 | |
1983 | |
1984 // ------------------ | |
1985 // sim_tissue_1min // | |
1986 // ------------------ | |
1987 // optimized in v.101 | |
1988 | |
1989 void sim_tissue_1min(void) | |
1990 { | |
1991 temp_pres_gtissue_limit = 0.0; | |
1992 temp_gtissue_no = 255; | |
1993 | |
1994 _asm | |
1995 lfsr 1, 0x300 | |
1996 movlw 0x01 | |
1997 movwf TBLPTRU,0 | |
1998 _endasm | |
1999 | |
2000 | |
2001 for (ci=0;ci<16;ci++) | |
2002 { | |
2003 _asm | |
2004 movlw 0x02 | |
2005 movwf TBLPTRH,0 | |
2006 movlb 4 // fuer ci | |
2007 movf ci,0,1 | |
2008 addwf ci,0,1 | |
2009 addwf ci,0,1 | |
2010 addwf ci,0,1 | |
2011 addlw 0x80 | |
2012 movwf TBLPTRL,0 | |
2013 TBLRDPOSTINC | |
2014 movff TABLAT,var_a+1 | |
2015 TBLRDPOSTINC | |
2016 movff TABLAT,var_a | |
2017 TBLRDPOSTINC | |
2018 movff TABLAT,var_a+3 | |
2019 TBLRD | |
2020 movff TABLAT,var_a+2 | |
2021 addlw 0x40 | |
2022 movwf TBLPTRL,0 | |
2023 TBLRDPOSTINC | |
2024 movff TABLAT,var2_a+1 | |
2025 TBLRDPOSTINC | |
2026 movff TABLAT,var2_a | |
2027 TBLRDPOSTINC | |
2028 movff TABLAT,var2_a+3 | |
2029 TBLRD | |
2030 movff TABLAT,var2_a+2 | |
2031 addlw 0x40 | |
2032 movwf TBLPTRL,0 | |
2033 incf TBLPTRH,1,0 | |
2034 TBLRDPOSTINC | |
2035 movff TABLAT,var_b+1 | |
2036 TBLRDPOSTINC | |
2037 movff TABLAT,var_b | |
2038 TBLRDPOSTINC | |
2039 movff TABLAT,var_b+3 | |
2040 TBLRD | |
2041 movff TABLAT,var_b+2 | |
2042 addlw 0x40 | |
2043 movwf TBLPTRL,0 | |
2044 TBLRDPOSTINC | |
2045 movff TABLAT,var2_b+1 | |
2046 TBLRDPOSTINC | |
2047 movff TABLAT,var2_b | |
2048 TBLRDPOSTINC | |
2049 movff TABLAT,var2_b+3 | |
2050 TBLRD | |
2051 movff TABLAT,var2_b+2 | |
2052 addlw 0xC0 | |
2053 movwf TBLPTRL,0 | |
2054 incf TBLPTRH,1,0 | |
2055 TBLRDPOSTINC | |
2056 movff TABLAT,var_e1min+1 | |
2057 TBLRDPOSTINC | |
2058 movff TABLAT,var_e1min | |
2059 TBLRDPOSTINC | |
2060 movff TABLAT,var_e1min+3 | |
2061 TBLRD | |
2062 movff TABLAT,var_e1min+2 | |
2063 addlw 0x40 | |
2064 movwf TBLPTRL,0 | |
2065 TBLRDPOSTINC | |
2066 movff TABLAT,var2_e1min+1 | |
2067 TBLRDPOSTINC | |
2068 movff TABLAT,var2_e1min | |
2069 TBLRDPOSTINC | |
2070 movff TABLAT,var2_e1min+3 | |
2071 TBLRD | |
2072 movff TABLAT,var2_e1min+2 | |
2073 _endasm | |
2074 // N2 | |
2075 temp_tissue = (temp_atem - sim_pres_tissue[ci]) * var_e1min; | |
2076 temp_tissue_safety(); | |
2077 sim_pres_tissue[ci] = sim_pres_tissue[ci] + temp_tissue; | |
2078 // He | |
2079 temp_tissue = (temp2_atem - sim_pres_tissue[ci+16]) * var2_e1min; | |
2080 temp_tissue_safety(); | |
2081 sim_pres_tissue[ci+16] = sim_pres_tissue[ci+16] + temp_tissue; | |
2082 // pressure limit | |
2083 temp_tissue = sim_pres_tissue[ci] + sim_pres_tissue[ci+16]; | |
2084 var_a = (var_a * sim_pres_tissue[ci] + var2_a * sim_pres_tissue[ci+16]) / temp_tissue; | |
2085 var_b = (var_b * sim_pres_tissue[ci] + var2_b * sim_pres_tissue[ci+16]) / temp_tissue; | |
2086 sim_pres_tissue_limit[ci] = (temp_tissue - var_a) * var_b; | |
2087 | |
2088 if (sim_pres_tissue_limit[ci] < 0) | |
2089 sim_pres_tissue_limit[ci] = 0; | |
2090 if (sim_pres_tissue_limit[ci] > temp_pres_gtissue_limit) | |
2091 { | |
2092 temp_pres_gtissue = temp_tissue; | |
2093 temp_pres_gtissue_limit = sim_pres_tissue_limit[ci]; | |
2094 temp_gtissue_no = ci; | |
2095 } | |
2096 } // for | |
2097 temp_pres_gtissue_diff = temp_pres_gtissue_limit - temp_pres_gtissue; | |
2098 temp_pres_gtissue_limit_GF_low = GF_low * temp_pres_gtissue_diff + temp_pres_gtissue; | |
2099 temp_pres_gtissue_limit_GF_low_below_surface = temp_pres_gtissue_limit_GF_low - pres_surface; | |
2100 if (temp_pres_gtissue_limit_GF_low_below_surface < 0) | |
2101 temp_pres_gtissue_limit_GF_low_below_surface = 0; | |
2102 } //sim_tissue_1min() | |
2103 | |
2104 //-------------------- | |
2105 // sim_tissue_10min // | |
2106 //-------------------- | |
2107 | |
2108 // Attention!! uses var_e1min und var2_e1min to load 10min data !!! | |
2109 // is identical to sim_tissue_1min routine except for the different load of those variables | |
2110 | |
2111 // optimized in v.101 | |
2112 | |
2113 void sim_tissue_10min(void) | |
2114 { | |
2115 temp_pres_gtissue_limit = 0.0; | |
2116 temp_gtissue_no = 255; | |
2117 | |
2118 _asm | |
2119 lfsr 1, 0x300 | |
2120 movlw 0x01 | |
2121 movwf TBLPTRU,0 | |
2122 _endasm | |
2123 | |
2124 for (ci=0;ci<16;ci++) | |
2125 { | |
2126 _asm | |
2127 movlw 0x02 | |
2128 movwf TBLPTRH,0 | |
2129 movlb 4 // fuer ci | |
2130 movf ci,0,1 | |
2131 addwf ci,0,1 | |
2132 addwf ci,0,1 | |
2133 addwf ci,0,1 | |
2134 addlw 0x80 | |
2135 movwf TBLPTRL,0 | |
2136 TBLRDPOSTINC | |
2137 movff TABLAT,var_a+1 | |
2138 TBLRDPOSTINC | |
2139 movff TABLAT,var_a | |
2140 TBLRDPOSTINC | |
2141 movff TABLAT,var_a+3 | |
2142 TBLRD | |
2143 movff TABLAT,var_a+2 | |
2144 addlw 0x40 | |
2145 movwf TBLPTRL,0 | |
2146 TBLRDPOSTINC | |
2147 movff TABLAT,var2_a+1 | |
2148 TBLRDPOSTINC | |
2149 movff TABLAT,var2_a | |
2150 TBLRDPOSTINC | |
2151 movff TABLAT,var2_a+3 | |
2152 TBLRD | |
2153 movff TABLAT,var2_a+2 | |
2154 addlw 0x40 | |
2155 movwf TBLPTRL,0 | |
2156 incf TBLPTRH,1,0 | |
2157 TBLRDPOSTINC | |
2158 movff TABLAT,var_b+1 | |
2159 TBLRDPOSTINC | |
2160 movff TABLAT,var_b | |
2161 TBLRDPOSTINC | |
2162 movff TABLAT,var_b+3 | |
2163 TBLRD | |
2164 movff TABLAT,var_b+2 | |
2165 addlw 0x40 | |
2166 movwf TBLPTRL,0 | |
2167 TBLRDPOSTINC | |
2168 movff TABLAT,var2_b+1 | |
2169 TBLRDPOSTINC | |
2170 movff TABLAT,var2_b | |
2171 TBLRDPOSTINC | |
2172 movff TABLAT,var2_b+3 | |
2173 TBLRD | |
2174 movff TABLAT,var2_b+2 | |
2175 addlw 0xC0 // different to 1 min | |
2176 movwf TBLPTRL,0 | |
2177 incf TBLPTRH,1,0 | |
2178 incf TBLPTRH,1,0 // different to 1 min | |
2179 TBLRDPOSTINC | |
2180 movff TABLAT,var_e1min+1 | |
2181 TBLRDPOSTINC | |
2182 movff TABLAT,var_e1min | |
2183 TBLRDPOSTINC | |
2184 movff TABLAT,var_e1min+3 | |
2185 TBLRD | |
2186 movff TABLAT,var_e1min+2 | |
2187 addlw 0x40 | |
2188 movwf TBLPTRL,0 | |
2189 //incf TBLPTRH,1,0 // different to 1 min | |
2190 TBLRDPOSTINC | |
2191 movff TABLAT,var2_e1min+1 | |
2192 TBLRDPOSTINC | |
2193 movff TABLAT,var2_e1min | |
2194 TBLRDPOSTINC | |
2195 movff TABLAT,var2_e1min+3 | |
2196 TBLRD | |
2197 movff TABLAT,var2_e1min+2 | |
2198 _endasm | |
2199 // N2 | |
2200 temp_tissue = (temp_atem - sim_pres_tissue[ci]) * var_e1min; | |
2201 temp_tissue_safety(); | |
2202 sim_pres_tissue[ci] = sim_pres_tissue[ci] + temp_tissue; | |
2203 // He | |
2204 temp_tissue = (temp2_atem - sim_pres_tissue[ci+16]) * var2_e1min; | |
2205 temp_tissue_safety(); | |
2206 sim_pres_tissue[ci+16] = sim_pres_tissue[ci+16] + temp_tissue; | |
2207 // pressure limit | |
2208 temp_tissue = sim_pres_tissue[ci] + sim_pres_tissue[ci+16]; | |
2209 var_a = (var_a * sim_pres_tissue[ci] + var2_a * sim_pres_tissue[ci+16]) / temp_tissue; | |
2210 var_b = (var_b * sim_pres_tissue[ci] + var2_b * sim_pres_tissue[ci+16]) / temp_tissue; | |
2211 | |
2212 sim_pres_tissue_limit[ci] = (temp_tissue - var_a) * var_b; | |
2213 if (sim_pres_tissue_limit[ci] < 0) | |
2214 sim_pres_tissue_limit[ci] = 0; | |
2215 if (sim_pres_tissue_limit[ci] > temp_pres_gtissue_limit) | |
2216 { | |
2217 temp_pres_gtissue = temp_tissue; | |
2218 temp_pres_gtissue_limit = sim_pres_tissue_limit[ci]; | |
2219 temp_gtissue_no = ci; | |
2220 } | |
2221 } // for | |
2222 temp_pres_gtissue_diff = temp_pres_gtissue_limit - temp_pres_gtissue; // negative number | |
2223 temp_pres_gtissue_limit_GF_low = GF_low * temp_pres_gtissue_diff + temp_pres_gtissue; | |
2224 temp_pres_gtissue_limit_GF_low_below_surface = temp_pres_gtissue_limit_GF_low - pres_surface; | |
2225 if (temp_pres_gtissue_limit_GF_low_below_surface < 0) | |
2226 temp_pres_gtissue_limit_GF_low_below_surface = 0; | |
2227 } //sim_tissue_10min() | |
2228 | |
2229 | |
2230 // ------------------ | |
2231 // clear_decoarray // | |
2232 // ------------------ | |
2233 // unchanged in v.101 | |
2234 | |
2235 void clear_decoarray(void) | |
2236 { | |
2237 char_O_array_decodepth[0] = 0; | |
2238 char_O_array_decodepth[1] = 0; | |
2239 char_O_array_decodepth[2] = 0; | |
2240 char_O_array_decodepth[3] = 0; | |
2241 char_O_array_decodepth[4] = 0; | |
2242 char_O_array_decodepth[5] = 0; | |
2243 char_O_array_decotime[0] = 0; | |
2244 char_O_array_decotime[1] = 0; | |
2245 char_O_array_decotime[2] = 0; | |
2246 char_O_array_decotime[3] = 0; | |
2247 char_O_array_decotime[4] = 0; | |
2248 char_O_array_decotime[5] = 0; | |
2249 char_O_array_decotime[6] = 0; | |
2250 } // clear_decoarray | |
2251 | |
2252 | |
2253 // ------------------- | |
2254 // update_decoarray // | |
2255 // ------------------- | |
2256 // unchanged in v.101 | |
2257 | |
2258 void update_decoarray() | |
2259 { | |
2260 x = 0; | |
2261 do | |
2262 { | |
2263 if (char_O_array_decodepth[x] == temp_depth_limit) | |
2264 { | |
2265 int_temp = char_O_array_decotime[x] + temp_decotime; | |
2266 if (int_temp < 0) | |
2267 int_temp = 0; | |
2268 if (int_temp > 240) | |
2269 int_temp = 240; | |
2270 char_O_array_decotime[x] = int_temp; | |
2271 x = 10; // exit | |
2272 } // if | |
2273 else | |
2274 { | |
2275 if (char_O_array_decodepth[x] == 0) | |
2276 { | |
2277 if (temp_depth_limit > 255) | |
2278 char_O_array_decodepth[x] = 255; | |
2279 else | |
2280 char_O_array_decodepth[x] = (char)temp_depth_limit; | |
2281 int_temp = char_O_array_decotime[x] + temp_decotime; | |
2282 if (int_temp > 240) | |
2283 char_O_array_decotime[x] = 240; | |
2284 else | |
2285 char_O_array_decotime[x] = (char)int_temp; | |
2286 x = 10; // exit | |
2287 } // if | |
2288 else | |
2289 x++; | |
2290 } // else | |
2291 } while (x<6); | |
2292 if (x == 6) | |
2293 { | |
2294 int_temp = char_O_array_decotime[6] + temp_decotime; | |
2295 if (int_temp > 220) | |
2296 char_O_array_decotime[6] = 220; | |
2297 else | |
2298 char_O_array_decotime[6] = (char)int_temp; | |
2299 } // if x == 6 | |
2300 } // update_decoarray | |
2301 | |
2302 | |
2303 // ----------------------- | |
2304 // calc_gradient_factor // | |
2305 // ----------------------- | |
2306 // optimized in v.101 (var_a) | |
2307 // new code in v.102 | |
2308 | |
2309 void calc_gradient_factor(void) | |
2310 { | |
2311 // tissue > respiration (entsaettigungsvorgang) | |
2312 // gradient ist wieviel prozent an limit mit basis tissue | |
2313 // dh. 0% = respiration == tissue | |
2314 // dh. 100% = respiration == limit | |
2315 temp_tissue = pres_tissue[char_O_gtissue_no] + pres_tissue[char_O_gtissue_no+16]; | |
2316 temp1 = temp_tissue - pres_respiration; | |
2317 temp2 = temp_tissue - pres_tissue_limit[char_O_gtissue_no]; // changed in v.102 | |
2318 temp2 = temp1/temp2; | |
2319 temp2 = temp2 * 100; // displayed in percent | |
2320 if (temp2 < 0) | |
2321 temp2 = 0; | |
2322 if (temp2 > 255) | |
2323 temp2 = 255; | |
2324 if (temp1 < 0) | |
2325 char_O_gradient_factor = 0; | |
2326 else | |
2327 char_O_gradient_factor = (char)temp2; | |
2328 | |
2329 temp3 = temp2; | |
2330 | |
2331 if (char_I_deco_model == 1) // calculate relative gradient factor | |
2332 { | |
2333 temp1 = (float)temp_depth_GF_low_meter * 0.09995; | |
2334 temp2 = pres_respiration - pres_surface; | |
2335 if (temp2 <= 0) | |
2336 temp1 = GF_high; | |
2337 else | |
2338 if (temp2 >= temp1) | |
2339 temp1 = GF_low; | |
2340 else | |
2341 temp1 = GF_low + (temp1 - temp2)/temp1*GF_delta; | |
2342 if (temp_depth_GF_low_meter == 0) | |
2343 temp1 = GF_high; | |
2344 temp2 = temp3 / temp1; // temp3 is already in percent | |
2345 if (temp2 < 0) | |
2346 temp2 = 0; | |
2347 if (temp2 > 255) | |
2348 temp2 = 255; | |
2349 char_O_relative_gradient_GF = (char)temp2; | |
2350 } // calc relative gradient factor | |
2351 else | |
2352 { | |
2353 char_O_relative_gradient_GF = char_O_gradient_factor; | |
2354 } | |
2355 } // calc_gradient | |
2356 | |
2357 // --------------------------- | |
2358 // calc_gradient_array_only // | |
2359 // --------------------------- | |
2360 // optimized in v.101 (var_a) | |
2361 // new code in v.102 | |
2362 | |
2363 void calc_gradient_array_only() | |
2364 { | |
2365 pres_respiration = (float)int_I_pres_respiration / 1000.0; // assembler code uses different digit system | |
2366 for (ci=0;ci<16;ci++) | |
2367 { | |
2368 temp_tissue = pres_tissue[ci] + pres_tissue[ci+16]; | |
2369 temp1 = temp_tissue - pres_respiration; | |
2370 temp2 = temp_tissue - pres_tissue_limit[ci]; | |
2371 temp2 = temp1/temp2; | |
2372 temp2 = temp2 * 200; // because of output in (Double-)percentage | |
2373 if (temp2 < 0) | |
2374 temp2 = 0; | |
2375 if (temp2 > 255) | |
2376 temp2 = 255; | |
2377 if (temp1 < 0) | |
2378 char_O_array_gradient_weighted[ci] = 0; | |
2379 else | |
2380 char_O_array_gradient_weighted[ci] = (char)temp2; | |
2381 } // for | |
2382 } // calc_gradient_array_only | |
2383 | |
2384 | |
2385 // ------------------------- | |
2386 // calc_desaturation_time // | |
2387 // ------------------------- | |
2388 // FIXED N2_ratio | |
2389 // unchanged in v.101 | |
2390 | |
2391 void calc_desaturation_time(void) | |
2392 { | |
2393 _asm | |
2394 lfsr 1, 0x300 | |
2395 movlw 0x01 | |
2396 movwf TBLPTRU,0 | |
2397 _endasm | |
2398 N2_ratio = 0.7902; // FIXED sum as stated in b"uhlmann | |
2399 pres_surface = (float)int_I_pres_surface / 1000.0; | |
2400 temp_atem = N2_ratio * (pres_surface - 0.0627); | |
2401 int_O_desaturation_time = 0; | |
2402 float_desaturation_multiplier = char_I_desaturation_multiplier / 142.0; // new in v.101 (70,42%/100.=142) | |
2403 | |
2404 for (ci=0;ci<16;ci++) | |
2405 { | |
2406 _asm | |
2407 movlw 0x04 | |
2408 movwf TBLPTRH,0 | |
2409 movlb 4 // fuer ci | |
2410 movf ci,0,1 | |
2411 addwf ci,0,1 | |
2412 addwf ci,0,1 | |
2413 addwf ci,0,1 | |
2414 addlw 0x80 | |
2415 movwf TBLPTRL,0 | |
2416 TBLRDPOSTINC | |
2417 movff TABLAT,var_halftimes+1 | |
2418 TBLRDPOSTINC | |
2419 movff TABLAT,var_halftimes | |
2420 TBLRDPOSTINC | |
2421 movff TABLAT,var_halftimes+3 | |
2422 TBLRD | |
2423 movff TABLAT,var_halftimes+2 | |
2424 addlw 0x40 | |
2425 movwf TBLPTRL,0 | |
2426 TBLRDPOSTINC | |
2427 movff TABLAT,var2_halftimes+1 | |
2428 TBLRDPOSTINC | |
2429 movff TABLAT,var2_halftimes | |
2430 TBLRDPOSTINC | |
2431 movff TABLAT,var2_halftimes+3 | |
2432 TBLRD | |
2433 movff TABLAT,var2_halftimes+2 | |
2434 _endasm | |
2435 | |
2436 // saturation_time (for flight) and N2_saturation in multiples of halftime | |
2437 // version v.100: 1.1 = 10 percent distance to totally clean (totally clean is not possible, would take infinite time ) | |
2438 // new in version v.101: 1.07 = 7 percent distance to totally clean (totally clean is not possible, would take infinite time ) | |
2439 // changes in v.101: 1.05 = 5 percent dist to totally clean is new desaturation point for display and noFly calculations | |
2440 // N2 | |
2441 temp1 = 1.05 * temp_atem; | |
2442 temp1 = temp1 - pres_tissue[ci]; | |
2443 temp2 = temp_atem - pres_tissue[ci]; | |
2444 if (temp2 >= 0.0) | |
2445 { | |
2446 temp1 = 0; | |
2447 temp2 = 0; | |
2448 } | |
2449 else | |
2450 temp1 = temp1 / temp2; | |
2451 if (temp1 > 0.0) | |
2452 { | |
2453 temp1 = log(1.0 - temp1); | |
2454 temp1 = temp1 / -0.6931; // temp1 is the multiples of half times necessary. | |
2455 // 0.6931 is ln(2), because the math function log() calculates with a base of e not 2 as requested. | |
2456 // minus because log is negative | |
2457 temp2 = var_halftimes * temp1 / float_desaturation_multiplier; // time necessary (in minutes ) for complete desaturation (see comment about 10 percent) , new in v.101: float_desaturation_multiplier | |
2458 } | |
2459 else | |
2460 { | |
2461 temp1 = 0; | |
2462 temp2 = 0; | |
2463 } | |
2464 | |
2465 // He | |
2466 temp3 = 0.1 - pres_tissue[ci+16]; | |
2467 if (temp3 >= 0.0) | |
2468 { | |
2469 temp3 = 0; | |
2470 temp4 = 0; | |
2471 } | |
2472 else | |
2473 temp3 = -1.0 * temp3 / pres_tissue[ci+16]; | |
2474 if (temp3 > 0.0) | |
2475 { | |
2476 temp3 = log(1.0 - temp3); | |
2477 temp3 = temp3 / -0.6931; // temp1 is the multiples of half times necessary. | |
2478 // 0.6931 is ln(2), because the math function log() calculates with a base of e not 2 as requested. | |
2479 // minus because log is negative | |
2480 temp4 = var2_halftimes * temp3 / float_desaturation_multiplier; // time necessary (in minutes ) for "complete" desaturation, new in v.101 float_desaturation_multiplier | |
2481 } | |
2482 else | |
2483 { | |
2484 temp3 = 0; | |
2485 temp4 = 0; | |
2486 } | |
2487 | |
2488 // saturation_time (for flight) | |
2489 if (temp4 > temp2) | |
2490 int_temp = (int)temp4; | |
2491 else | |
2492 int_temp = (int)temp2; | |
2493 if(int_temp > int_O_desaturation_time) | |
2494 int_O_desaturation_time = int_temp; | |
2495 | |
2496 // N2 saturation in multiples of halftime for display purposes | |
2497 temp2 = temp1 * 20.0; // 0 = 1/8, 120 = 0, 249 = 8 | |
2498 temp2 = temp2 + 80.0; // set center | |
2499 if (temp2 < 0.0) | |
2500 temp2 = 0.0; | |
2501 if (temp2 > 255.0) | |
2502 temp2 = 255.0; | |
2503 char_O_tissue_saturation[ci] = (char)temp2; | |
2504 // He saturation in multiples of halftime for display purposes | |
2505 temp4 = temp3 * 20.0; // 0 = 1/8, 120 = 0, 249 = 8 | |
2506 temp4 = temp4 + 80.0; // set center | |
2507 if (temp4 < 0.0) | |
2508 temp4 = 0.0; | |
2509 if (temp4 > 255.0) | |
2510 temp4 = 255.0; | |
2511 char_O_tissue_saturation[ci+16] = (char)temp4; | |
2512 } // for | |
2513 } // calc_desaturation_time | |
2514 | |
2515 | |
2516 // -------------------------- | |
2517 // calc_wo_deco_step_1_min // | |
2518 // -------------------------- | |
2519 // FIXED N2 Ratio | |
2520 // optimized in v.101 (...saturation_multiplier) | |
2521 // desaturation slowed down to 70,42% | |
2522 | |
2523 void calc_wo_deco_step_1_min(void) | |
2524 { | |
2525 if(flag_in_divemode) | |
2526 { | |
2527 flag_in_divemode = 0; | |
2528 set_dbg_end_of_dive(); | |
2529 } | |
2530 _asm | |
2531 lfsr 1, 0x300 | |
2532 _endasm | |
2533 N2_ratio = 0.7902; // FIXED, sum lt. buehlmann | |
2534 pres_respiration = (float)int_I_pres_respiration / 1000.0; // assembler code uses different digit system | |
2535 pres_surface = (float)int_I_pres_surface / 1000.0; | |
2536 temp_atem = N2_ratio * (pres_respiration - 0.0627); // 0.0627 is the extra pressure in the body | |
2537 temp2_atem = 0.0; | |
2538 temp_surface = pres_surface; // the b"uhlmann formula using temp_surface does not use the N2_ratio | |
2539 float_desaturation_multiplier = char_I_desaturation_multiplier / 142.0; // new in v.101 (70,42%/100.=142) | |
2540 float_saturation_multiplier = char_I_saturation_multiplier / 100.0; | |
2541 | |
2542 calc_tissue_step_1_min(); // update the pressure in the 16 tissues in accordance with the new ambient pressure | |
2543 clear_decoarray(); | |
2544 char_O_deco_status = 0; | |
2545 char_O_nullzeit = 0; | |
2546 char_O_ascenttime = 0; | |
2547 calc_gradient_factor(); | |
2548 | |
2549 } // calc_wo_deco_step_1_min(void) | |
2550 | |
2551 | |
2552 // ------------------------- | |
2553 // calc_tissue_step_1_min // | |
2554 // ------------------------- | |
2555 // optimized in v.101 | |
2556 | |
2557 void calc_tissue_step_1_min(void) | |
2558 { | |
2559 _asm | |
2560 lfsr 1, 0x300 | |
2561 movlw 0x01 | |
2562 movwf TBLPTRU,0 | |
2563 _endasm | |
2564 | |
2565 char_O_gtissue_no = 255; | |
2566 pres_gtissue_limit = 0.0; | |
2567 | |
2568 for (ci=0;ci<16;ci++) | |
2569 { | |
2570 _asm | |
2571 movlw 0x02 | |
2572 movwf TBLPTRH,0 | |
2573 movlb 4 // fuer ci | |
2574 movf ci,0,1 | |
2575 addwf ci,0,1 | |
2576 addwf ci,0,1 | |
2577 addwf ci,0,1 | |
2578 addlw 0x80 | |
2579 movwf TBLPTRL,0 | |
2580 TBLRDPOSTINC | |
2581 movff TABLAT,var_a+1 | |
2582 TBLRDPOSTINC | |
2583 movff TABLAT,var_a | |
2584 TBLRDPOSTINC | |
2585 movff TABLAT,var_a+3 | |
2586 TBLRD | |
2587 movff TABLAT,var_a+2 | |
2588 addlw 0x40 | |
2589 movwf TBLPTRL,0 | |
2590 TBLRDPOSTINC | |
2591 movff TABLAT,var2_a+1 | |
2592 TBLRDPOSTINC | |
2593 movff TABLAT,var2_a | |
2594 TBLRDPOSTINC | |
2595 movff TABLAT,var2_a+3 | |
2596 TBLRD | |
2597 movff TABLAT,var2_a+2 | |
2598 addlw 0x40 | |
2599 movwf TBLPTRL,0 | |
2600 incf TBLPTRH,1,0 | |
2601 TBLRDPOSTINC | |
2602 movff TABLAT,var_b+1 | |
2603 TBLRDPOSTINC | |
2604 movff TABLAT,var_b | |
2605 TBLRDPOSTINC | |
2606 movff TABLAT,var_b+3 | |
2607 TBLRD | |
2608 movff TABLAT,var_b+2 | |
2609 addlw 0x40 | |
2610 movwf TBLPTRL,0 | |
2611 TBLRDPOSTINC | |
2612 movff TABLAT,var2_b+1 | |
2613 TBLRDPOSTINC | |
2614 movff TABLAT,var2_b | |
2615 TBLRDPOSTINC | |
2616 movff TABLAT,var2_b+3 | |
2617 TBLRD | |
2618 movff TABLAT,var2_b+2 | |
2619 addlw 0xC0 | |
2620 movwf TBLPTRL,0 | |
2621 incf TBLPTRH,1,0 | |
2622 TBLRDPOSTINC | |
2623 movff TABLAT,var_e1min+1 | |
2624 TBLRDPOSTINC | |
2625 movff TABLAT,var_e1min | |
2626 TBLRDPOSTINC | |
2627 movff TABLAT,var_e1min+3 | |
2628 TBLRD | |
2629 movff TABLAT,var_e1min+2 | |
2630 addlw 0x40 | |
2631 movwf TBLPTRL,0 | |
2632 TBLRDPOSTINC | |
2633 movff TABLAT,var2_e1min+1 | |
2634 TBLRDPOSTINC | |
2635 movff TABLAT,var2_e1min | |
2636 TBLRDPOSTINC | |
2637 movff TABLAT,var2_e1min+3 | |
2638 TBLRD | |
2639 movff TABLAT,var2_e1min+2 | |
2640 _endasm | |
2641 | |
2642 // N2 1 min | |
2643 temp_tissue = (temp_atem - pres_tissue[ci]) * var_e1min; | |
2644 temp_tissue_safety(); | |
2645 pres_tissue[ci] = pres_tissue[ci] + temp_tissue; | |
2646 | |
2647 // He 1 min | |
2648 temp_tissue = (temp2_atem - pres_tissue[ci+16]) * var2_e1min; | |
2649 temp_tissue_safety(); | |
2650 pres_tissue[ci+16] = pres_tissue[ci+16] + temp_tissue; | |
2651 | |
2652 temp_tissue = pres_tissue[ci] + pres_tissue[ci+16]; | |
2653 var_a = (var_a * pres_tissue[ci] + var2_a * pres_tissue[ci+16]) / temp_tissue; | |
2654 var_b = (var_b * pres_tissue[ci] + var2_b * pres_tissue[ci+16]) / temp_tissue; | |
2655 pres_tissue_limit[ci] = (temp_tissue - var_a) * var_b; | |
2656 if (pres_tissue_limit[ci] < 0) | |
2657 pres_tissue_limit[ci] = 0; | |
2658 if (pres_tissue_limit[ci] > pres_gtissue_limit) | |
2659 { | |
2660 pres_gtissue_limit = pres_tissue_limit[ci]; | |
2661 char_O_gtissue_no = ci; | |
2662 }//if | |
2663 | |
2664 // gradient factor array for graphical display | |
2665 // display range is 0 to 250! in steps of 5 for 1 pixel | |
2666 // the display is divided in 6 blocks | |
2667 // -> double the gradient 100% = 200 | |
2668 // tissue > respiration (entsaettigungsvorgang) | |
2669 // gradient ist wieviel prozent an limit von tissue aus | |
2670 // dh. 0% = respiration == tissue | |
2671 // dh. 100% = respiration == limit | |
2672 temp1 = temp_tissue - pres_respiration; | |
2673 temp2 = temp_tissue - pres_tissue_limit[ci]; // changed in v.102 | |
2674 temp2 = temp1/temp2; | |
2675 temp2 = temp2 * 200; // because of output in (Double-)percentage | |
2676 if (temp2 < 0) | |
2677 temp2 = 0; | |
2678 if (temp2 > 255) | |
2679 temp2 = 255; | |
2680 if (temp1 < 0) | |
2681 char_O_array_gradient_weighted[ci] = 0; | |
2682 else | |
2683 char_O_array_gradient_weighted[ci] = (char)temp2; | |
2684 | |
2685 } // for | |
2686 } // calc wo deco 1min | |
2687 | |
2688 #if 0 | |
2689 // -------- | |
2690 // debug // | |
2691 // -------- | |
2692 void debug(void) | |
2693 { | |
2694 for (ci=0;ci<32;ci++) | |
2695 { | |
2696 int_O_tissue_for_debug[ci] = (unsigned int)(pres_tissue[ci] *1000); | |
2697 } | |
2698 } // void debug(void) | |
2699 #endif | |
2700 | |
2701 // ---------- | |
2702 // md hash // | |
2703 // ---------- | |
2704 void hash(void) | |
2705 { | |
2706 // init | |
2707 for (md_i=0;md_i<16;md_i++) | |
2708 { | |
2709 md_state[md_i] = 0; | |
2710 md_cksum[md_i] = 0; | |
2711 } // for md_i 16 | |
2712 | |
2713 _asm | |
2714 movlw 0x01 | |
2715 movwf TBLPTRU,0 | |
2716 movlw 0x06 | |
2717 movwf TBLPTRH,0 | |
2718 movlw 0x00 | |
2719 movwf TBLPTRL,0 | |
2720 _endasm | |
2721 for (md_i=0;md_i<127;md_i++) | |
2722 { | |
2723 _asm | |
2724 TBLRDPOSTINC | |
2725 movff TABLAT,md_temp | |
2726 _endasm | |
2727 md_pi_subst[md_i] = md_temp; | |
2728 } // for md_i 256 | |
2729 _asm | |
2730 TBLRDPOSTINC | |
2731 movff TABLAT,md_temp | |
2732 _endasm | |
2733 md_pi_subst[127] = md_temp; | |
2734 for (md_i=0;md_i<127;md_i++) | |
2735 { | |
2736 _asm | |
2737 TBLRDPOSTINC | |
2738 movff TABLAT,md_temp | |
2739 _endasm | |
2740 md_pi_subst[md_i+128] = md_temp; | |
2741 } // for md_i 256 | |
2742 _asm | |
2743 TBLRD | |
2744 movff TABLAT,md_temp | |
2745 _endasm | |
2746 md_pi_subst[255] = md_temp; | |
2747 | |
2748 _asm | |
2749 movlw 0x00 | |
2750 movwf TBLPTRU,0 | |
2751 movlw 0x00 | |
2752 movwf TBLPTRH,0 | |
2753 movlw 0x00 | |
2754 movwf TBLPTRL,0 | |
2755 _endasm | |
2756 // cycle buffers | |
2757 for (md_pointer=0x0000;md_pointer<0x17f3;md_pointer++) | |
2758 { | |
2759 md_t = 0; | |
2760 for (md_i=0;md_i<16;md_i++) | |
2761 { | |
2762 if(md_pointer == 9) | |
2763 md_temp = md_cksum[md_i]; | |
2764 else | |
2765 { | |
2766 _asm | |
2767 TBLRDPOSTINC | |
2768 movff TABLAT,md_temp | |
2769 _endasm | |
2770 } // else | |
2771 md_buffer[md_i] = md_temp; | |
2772 md_state[md_i+16] = md_buffer[md_i]; | |
2773 md_state[md_i+32] = (unsigned char)(md_buffer[md_i] ^ md_state[md_i]); | |
2774 } // for md_i 16 | |
2775 | |
2776 for (md_i=0;md_i<18;md_i++) | |
2777 { | |
2778 for (md_j=0;md_j<48;md_j++) | |
2779 { | |
2780 md_state[md_j] = (unsigned char)(md_state[md_j] ^ md_pi_subst[md_t]); | |
2781 md_t = md_state[md_j]; | |
2782 } // for md_j 48 | |
2783 md_t = (unsigned char)(md_t+1); | |
2784 } // for md_i 18 | |
2785 md_t = md_cksum[15]; | |
2786 | |
2787 for (md_i=0;md_i<16;md_i++) | |
2788 { | |
2789 md_cksum[md_i] = (unsigned char)(md_cksum[md_i] ^ md_pi_subst[(md_buffer[md_i] ^ md_t)]); | |
2790 md_t = md_cksum[md_i]; | |
2791 } // for md_i 16 | |
2792 } // for md_pointer | |
2793 } // void hash(void) | |
2794 | |
2795 // --------------------- | |
2796 // clear_CNS_fraction // | |
2797 // --------------------- | |
2798 // new in v.101 | |
2799 | |
2800 void clear_CNS_fraction(void) | |
2801 { | |
2802 CNS_fraction = 0.0; | |
2803 char_O_CNS_fraction = 0; | |
2804 } // void clear_CNS_fraction(void) | |
2805 | |
2806 | |
2807 // -------------------- | |
2808 // calc_CNS_fraction // | |
2809 // -------------------- | |
2810 // new in v.101 | |
2811 // optimized in v.102 : with new variables char_I_actual_ppO2 and actual_ppO2 | |
2812 | |
2813 // Input: char_I_actual_ppO2 | |
2814 // Output: char_O_CNS_fraction | |
2815 // Uses and Updates: CNS_fraction | |
2816 // Uses: acutal_ppO2 | |
2817 | |
2818 void calc_CNS_fraction(void) | |
2819 { | |
2820 actual_ppO2 = (float)char_I_actual_ppO2 / 100.0; | |
2821 | |
2822 if (char_I_actual_ppO2 < 50) | |
2823 CNS_fraction = CNS_fraction;// no changes | |
2824 else if (char_I_actual_ppO2 < 60) | |
2825 CNS_fraction = 1/(-54000.0 * actual_ppO2 + 54000.0) + CNS_fraction; | |
2826 else if (char_I_actual_ppO2 < 70) | |
2827 CNS_fraction = 1/(-45000.0 * actual_ppO2 + 48600.0) + CNS_fraction; | |
2828 else if (char_I_actual_ppO2 < 80) | |
2829 CNS_fraction = 1/(-36000.0 * actual_ppO2 + 42300.0) + CNS_fraction; | |
2830 else if (char_I_actual_ppO2 < 90) | |
2831 CNS_fraction = 1/(-27000.0 * actual_ppO2 + 35100.0) + CNS_fraction; | |
2832 else if (char_I_actual_ppO2 < 110) | |
2833 CNS_fraction = 1/(-18000.0 * actual_ppO2 + 27000.0) + CNS_fraction; | |
2834 else if (char_I_actual_ppO2 < 150) | |
2835 CNS_fraction = 1/(-9000.0 * actual_ppO2 + 17100.0) + CNS_fraction; | |
2836 else if (char_I_actual_ppO2 < 160) | |
2837 CNS_fraction = 1/(-22500.0 * actual_ppO2 + 37350.0) + CNS_fraction; | |
2838 else if (char_I_actual_ppO2 < 165) | |
2839 CNS_fraction = 0.000755 + CNS_fraction; // Arieli et all.(2002): Modeling pulmonary and CNS O2 toxicity... Formula (A1) based on value for 1.55 and c=20 | |
2840 else if (char_I_actual_ppO2 < 170) | |
2841 CNS_fraction = 0.00102 + CNS_fraction; // example calculation: Sqrt((1.7/1.55)^20)*0.000404 | |
2842 else if (char_I_actual_ppO2 < 175) | |
2843 CNS_fraction = 0.00136 + CNS_fraction; | |
2844 else if (char_I_actual_ppO2 < 180) | |
2845 CNS_fraction = 0.00180 + CNS_fraction; | |
2846 else if (char_I_actual_ppO2 < 185) | |
2847 CNS_fraction = 0.00237 + CNS_fraction; | |
2848 else if (char_I_actual_ppO2 < 190) | |
2849 CNS_fraction = 0.00310 + CNS_fraction; | |
2850 else if (char_I_actual_ppO2 < 195) | |
2851 CNS_fraction = 0.00401 + CNS_fraction; | |
2852 else if (char_I_actual_ppO2 < 200) | |
2853 CNS_fraction = 0.00517 + CNS_fraction; | |
2854 else if (char_I_actual_ppO2 < 230) | |
2855 CNS_fraction = 0.0209 + CNS_fraction; | |
2856 else | |
2857 CNS_fraction = 0.0482 + CNS_fraction; // value for 2.5 | |
2858 | |
2859 if (CNS_fraction > 2.5) | |
2860 CNS_fraction = 2.5; | |
2861 if (CNS_fraction < 0.0) | |
2862 CNS_fraction = 0.0; | |
2863 char_O_CNS_fraction = (char)((CNS_fraction + 0.005)* 100.0); | |
2864 } // void calc_CNS_fraction(void) | |
2865 | |
2866 // -------------------------- | |
2867 // calc_CNS_decrease_15min // | |
2868 // -------------------------- | |
2869 // new in v.101 | |
2870 | |
2871 // calculates the half time of 90 minutes in 6 steps of 15 min | |
2872 | |
2873 // Output: char_O_CNS_fraction | |
2874 // Uses and Updates: CNS_fraction | |
2875 | |
2876 void calc_CNS_decrease_15min(void) | |
2877 { | |
2878 CNS_fraction = 0.890899 * CNS_fraction; | |
2879 char_O_CNS_fraction = (char)((CNS_fraction + 0.005)* 100.0); | |
2880 }// calc_CNS_decrease_15min(void) | |
2881 | |
2882 | |
2883 // ------------------ | |
2884 // calc_percentage // | |
2885 // ------------------ | |
2886 // new in v.101 | |
2887 | |
2888 // calculates int_I_temp * char_I_temp / 100 | |
2889 // output is int_I_temp | |
2890 | |
2891 void calc_percentage(void) | |
2892 { | |
2893 temp1 = (float)int_I_temp; | |
2894 temp2 = (float)char_I_temp / 100.0; | |
2895 temp3 = temp1 * temp2; | |
2896 int_I_temp = (int)temp3; | |
2897 } | |
2898 void push_tissues_to_vault(void) | |
2899 { | |
2900 for (ci=0;ci<32;ci++) | |
2901 pres_tissue_vault[ci] = pres_tissue[ci]; | |
2902 } | |
2903 void pull_tissues_from_vault(void) | |
2904 { | |
2905 for (ci=0;ci<32;ci++) | |
2906 pres_tissue[ci] = pres_tissue_vault[ci]; | |
2907 } | |
2908 | |
2909 void wp_write_command(void) | |
2910 { | |
2911 _asm | |
2912 bcf oled_rs | |
2913 movff wp_command,PORTD | |
2914 bcf oled_rw | |
2915 bsf oled_rw | |
2916 _endasm | |
2917 } | |
2918 | |
2919 void wp_write_data(void) | |
2920 { | |
2921 wp_data_8bit_one = wp_data_16bit >> 8; | |
2922 wp_data_8bit_two = wp_data_16bit; | |
2923 _asm | |
2924 bsf oled_rs | |
2925 movff wp_data_8bit_one,PORTD | |
2926 bcf oled_rw | |
2927 bsf oled_rw | |
2928 movff wp_data_8bit_two,PORTD | |
2929 bcf oled_rw | |
2930 bsf oled_rw | |
2931 _endasm | |
2932 } | |
2933 | |
2934 void wp_write_black(void) | |
2935 { | |
2936 _asm | |
2937 movff wp_black,PORTD | |
2938 bcf oled_rw | |
2939 bsf oled_rw | |
2940 bcf oled_rw | |
2941 bsf oled_rw | |
2942 _endasm | |
2943 } | |
2944 | |
2945 void wp_write_color(void) | |
2946 { | |
2947 _asm | |
2948 movff wp_color1,PORTD | |
2949 bcf oled_rw | |
2950 bsf oled_rw | |
2951 movff wp_color2,PORTD | |
2952 bcf oled_rw | |
2953 bsf oled_rw | |
2954 _endasm | |
2955 } | |
2956 | |
2957 void wp_set_window(void) | |
2958 { | |
2959 // x axis start ( 0 - 319) | |
2960 wp_command = 0x35; | |
2961 wp_write_command(); | |
2962 wp_data_16bit = ((U16)wp_leftx2) << 1; | |
2963 wp_write_data(); | |
2964 // x axis end ( 0 - 319) | |
2965 wp_command = 0x36; | |
2966 wp_write_command(); | |
2967 wp_data_16bit = 319; | |
2968 wp_write_data(); | |
2969 // y axis start + end ( 0 - 239 ) | |
2970 wp_command = 0x37; | |
2971 wp_write_command(); | |
2972 // the bottom part | |
2973 wp_data_16bit = wp_top; | |
2974 if(wp_font == 2) | |
2975 wp_data_16bit += WP_FONT_LARGE_HEIGHT; | |
2976 else if(wp_font == 1) | |
2977 wp_data_16bit += WP_FONT_MEDIUM_HEIGHT; | |
2978 else | |
2979 wp_data_16bit += WP_FONT_SMALL_HEIGHT; | |
2980 wp_data_16bit--; | |
2981 if(wp_data_16bit > 239) | |
2982 wp_data_16bit = 239; | |
2983 // the top part | |
2984 wp_data_16bit |= ((U16)wp_top) << 8; | |
2985 // all together in one 16bit transfer | |
2986 wp_write_data(); | |
2987 | |
2988 // start | |
2989 wp_command = 0x20; | |
2990 wp_write_command(); | |
2991 wp_data_16bit = wp_top; | |
2992 wp_write_data(); | |
2993 | |
2994 wp_command = 0x21; | |
2995 wp_write_command(); | |
2996 wp_data_16bit = ((U16)wp_leftx2) << 1; | |
2997 wp_write_data(); | |
2998 } | |
2999 | |
3000 void wp_set_char_font_small(void) | |
3001 { | |
3002 if (wp_char == ' ') | |
3003 wp_char = '¶'; | |
3004 | |
3005 if((wp_char < '!') || (wp_char > '¶')) // font has 34 chars after ~ // ¾ + 4 chars limit to end of battery at the moment | |
3006 wp_char = '¤'; // 0x82; // ¤ | |
3007 | |
3008 if (wp_char > 0x7E) // skip space between ~ and ¡ | |
3009 wp_char -= 34; | |
3010 | |
3011 wp_start = wp_small_table[wp_char - '!']; | |
3012 wp_end = wp_small_table[1 + wp_char - '!']; | |
3013 } | |
3014 | |
3015 void wp_set_char_font_medium(void) | |
3016 { | |
3017 // space is 3E | |
3018 if (wp_char == 0x27) // 0x27 == ' | |
3019 wp_char = 0x3B; | |
3020 if (wp_char == '"') | |
3021 wp_char = 0x3C; | |
3022 if (wp_char == 'm') | |
3023 wp_char = 0x3D; | |
3024 if (wp_char == ' ') | |
3025 wp_char = 0x3E; | |
3026 | |
3027 if((wp_char < '.') || (wp_char > 0x3E)) | |
3028 wp_char = 0x3E; | |
3029 wp_start = wp_medium_table[wp_char - '.']; | |
3030 wp_end = wp_medium_table[1 + wp_char - '.']; | |
3031 } | |
3032 | |
3033 void wp_set_char_font_large(void) | |
3034 { | |
3035 // space is / = 0x2F | |
3036 if (wp_char == ' ') | |
3037 wp_char = 0x2F; | |
3038 | |
3039 if((wp_char < '.') || (wp_char > '9')) | |
3040 wp_char = 0x2F; | |
3041 wp_start = wp_large_table[wp_char - '.']; | |
3042 wp_end = wp_large_table[1 + wp_char - '.']; | |
3043 } | |
3044 | |
3045 void wordprocessor(void) | |
3046 { | |
3047 wp_set_window(); | |
3048 | |
3049 // access to GRAM | |
3050 wp_command = 0x22; | |
3051 wp_write_command(); | |
3052 _asm | |
3053 bsf oled_rs | |
3054 _endasm | |
3055 | |
3056 wp_txtptr = 0; | |
3057 wp_char = wp_stringstore[wp_txtptr]; | |
3058 | |
3059 while(wp_char) | |
3060 { | |
3061 if(wp_font == 2) | |
3062 wp_set_char_font_large(); | |
3063 else if(wp_font == 1) | |
3064 wp_set_char_font_medium(); | |
3065 else | |
3066 wp_set_char_font_small(); | |
3067 | |
3068 wp_black = 0; | |
3069 | |
3070 for(wp_i = wp_start; wp_i<wp_end;wp_i++) | |
3071 { | |
3072 if(wp_font == 2) | |
3073 wp_data_16bit = wp_large_data[wp_i / 2]; | |
3074 else if(wp_font == 1) | |
3075 wp_data_16bit = wp_medium_data[wp_i / 2]; | |
3076 else | |
3077 wp_data_16bit = wp_small_data[wp_i / 2]; | |
3078 if(wp_i & 1) | |
3079 wp_temp_U8 = wp_data_16bit & 0xFF; | |
3080 else | |
3081 wp_temp_U8 = wp_data_16bit >> 8; | |
3082 if((wp_temp_U8 & 128)) | |
3083 { | |
3084 wp_temp_U8 -= 127; | |
3085 if(wp_invert) | |
3086 { | |
3087 while(wp_temp_U8 > 0) | |
3088 { | |
3089 wp_temp_U8--; | |
3090 wp_write_color(); | |
3091 } | |
3092 } | |
3093 else | |
3094 { | |
3095 _asm | |
3096 movff wp_black,PORTD | |
3097 _endasm | |
3098 while(wp_temp_U8 > 0) | |
3099 { | |
3100 wp_temp_U8--; | |
3101 _asm | |
3102 bcf oled_rw | |
3103 bsf oled_rw | |
3104 bcf oled_rw | |
3105 bsf oled_rw | |
3106 _endasm | |
3107 } | |
3108 } | |
3109 } | |
3110 else | |
3111 { | |
3112 wp_temp_U8++; | |
3113 if(wp_invert) | |
3114 { | |
3115 _asm | |
3116 movff wp_black,PORTD | |
3117 _endasm | |
3118 while(wp_temp_U8 > 0) | |
3119 { | |
3120 wp_temp_U8--; | |
3121 _asm | |
3122 bcf oled_rw | |
3123 bsf oled_rw | |
3124 bcf oled_rw | |
3125 bsf oled_rw | |
3126 _endasm | |
3127 } | |
3128 } | |
3129 else | |
3130 { | |
3131 while(wp_temp_U8 > 0) | |
3132 { | |
3133 wp_temp_U8--; | |
3134 wp_write_color(); | |
3135 } | |
3136 } | |
3137 } | |
3138 } | |
3139 wp_txtptr++; | |
3140 wp_char = wp_stringstore[wp_txtptr]; | |
3141 } | |
3142 wp_command = 0x00; | |
3143 wp_write_command(); | |
3144 } |