mk2: code_part1/OSTC_code_c_part2/p2

comparison code_part1/OSTC_code_c_part2/p2_deco.c @ 165:8b5b848db4cd

+ Clean Buhlmann tables addressing + OSTC 1.80 uses ZH-L16a coefficient, version 1995 (?), with compartiment 1a.

author	JeanDo
date	Fri, 21 Jan 2011 01:44:38 +0100
parents	999abb01c78f
children	cb055a7d75f3

comparison

equal deleted inserted replaced

-:999abb01c78f
+:8b5b848db4cd
 static float			temp_tissue;
 static float			temp_surface;
 static float			N2_ratio;
 static float			He_ratio;
 static float			temp_ratio;
-static float 			var_a;
+static float 			var_N2_a;
-static float 			var2_a;
+static float 			var_N2_b;
-static float 			var_b;
+static float 			var_He_a;
-static float 			var2_b;
+static float 			var_He_b;
-static float 			var_t05nc;
+static float  			var_N2_e;
-static float 			var2_t05nc;
+static float  			var_He_e;
-static float  			var_e2secs;
+static float  			var_N2_halftime;
-static float  			var2_e2secs;
+static float  			var_He_halftime;
-static float  			var_e1min;
-static float  			var2_e1min;
-static float  			var_halftimes;
-static float  			var2_halftimes;
 static float			pres_gtissue_limit;
 static float			temp_pres_gtissue_limit;
 static float			actual_ppO2;                // new in v.102
 static float            pres_diluent;               // new in v.101
 // ***********************
 // ** THE LOOKUP TABLES **
 // ***********************
 // ***********************
-#pragma romdata tables = 0x10200
+#pragma romdata buhlmann_tables = 0x010200  // Needs to be in UPPER bank.
-#include	"p2_tables.romdata" 		// new table for deco_main_v.101 (var_a modified)
+#include	"p2_tables.romdata" 		// new table for deco_main_v.101 (var_N2_a modified)
-#pragma romdata tables2 = 0x10600
+#pragma romdata hash_tables = 0x010600  // Address fixed by ASM access...
 rom const rom unsigned int md_pi[] =
 {
 0x292E, 0x43C9, 0xA2D8, 0x7C01, 0x3D36, 0x54A1, 0xECF0, 0x0613
 , 0x62A7, 0x05F3, 0xC0C7, 0x738C, 0x9893, 0x2BD9, 0xBC4C, 0x82CA
 , 0x1E9B, 0x573C, 0xFDD4, 0xE016, 0x6742, 0x6F18, 0x8A17, 0xE512
 void check_post_dbg(void)
 {
 	check_dbg(1);
 }
-// -------------------------
+//////////////////////////////////////////////////////////////////////////////
+// read buhlmann tables for compatriment ci
+// If mode == 0 : 2sec interval
+//            1 : 1 min interval
+//            2 : 10 min interval.
+// TODO: an assembleur version of the table reads...
+static void read_buhlmann_compartment(static char mode)
+{
+var_N2_a = buhlmann_a[ci];
+var_N2_b = buhlmann_b[ci];
+var_He_a = buhlmann_a[ci + 16];
+var_He_b = buhlmann_b[ci + 16];
+// Check reading consistency:
+	if(	(var_N2_a < 0.231)
+	 || (var_N2_a > 1.27)
+	 || (var_N2_b < 0.504)
+	 || (var_N2_b > 0.966)
+	 || (var_He_a < 0.510)
+	 || (var_He_a > 1.75)
+	 || (var_He_b < 0.423)
+	 || (var_He_b > 0.927)
+)
+int_O_DBG_pre_bitfield |= DBG_ZH16ERR;
+// Integration intervals.
+switch(mode)
+{
+case 0: //---- 2 sec -----------------------------------------------------
+var_N2_e = e2secs[ci];
+var_He_e = e2secs[ci + 16];
+// Check reading consistency:
+	if(	(var_N2_e < 0.0000363)
+	 || (var_N2_e > 0.00577)
+	 || (var_He_e < 0.0000961)
+	 || (var_He_e > 0.150)
+)
+int_O_DBG_pre_bitfield |= DBG_ZH16ERR;
+break;
+case 1: //---- 1 min -----------------------------------------------------
+var_N2_e = e1min[ci];
+var_He_e = e1min[ci + 16];
+// Check reading consistency:
+	if(	(var_N2_e < 1.09E-3)
+	 || (var_N2_e > 0.1592)
+	 || (var_He_e < 0.00288)
+	 || (var_He_e > 0.3682)
+)
+int_O_DBG_pre_bitfield |= DBG_ZH16ERR;
+break;
+case 2: //---- 10 min ----------------------------------------------------
+var_N2_e = e10min[ci];
+var_He_e = e10min[ci + 16];
+// Check reading consistency:
+	if(	(var_N2_e < 0.01085)
+	 || (var_N2_e > 0.82323)
+	 || (var_He_e < 0.02846)
+	 || (var_He_e > 0.98986)
+)
+int_O_DBG_pre_bitfield |= DBG_ZH16ERR;
+break;
+}
+}
+//////////////////////////////////////////////////////////////////////////////
 // calc_next_decodepth_GF //
-// -------------------------
+//
 // new in v.102
+//
 void calc_nextdecodepth_GF(void)
 {
 // INPUT, changing during dive:
 // temp_pres_gtissue_limit_GF_low
 // temp_pres_gtissue_limit_GF_low_below_surface
 }
 // ---------------
 // CLEAR tissue //
 // ---------------
-// optimized in v.101 (var_a)
+// optimized in v.101 (var_N2_a)
 #pragma code p2_deco_suite = 0x10700
 void clear_tissue(void)    // preload tissues with standard pressure for the given ambient pressure
 {
 movlw 1
 movwf TBLPTRU,0
 _endasm
 for(ci=0;ci<16;ci++)
 {
-overlay float p;
 // cycle through the 16 b"uhlmann tissues
-p = N2_ratio * (pres_respiration -  0.0627) ;
+overlay float p = N2_ratio * (pres_respiration -  0.0627);
 pres_tissue[ci] = p;
-var_a = ((rom float*)a_N2)[ci];
+read_buhlmann_compartment(0);
-var_b = ((rom float*)b_N2)[ci];
-p = (p - var_a) * var_b ;
+p = (p - var_N2_a) * var_N2_b ;
 if( p < 0.0 )
 p = 0.0;
 pres_tissue_limit[ci] = p;
 // cycle through the 16 b"uhlmann tissues for Helium
 char_O_gtissue_no = 255;
 pres_gtissue_limit = 0.0;
 for (ci=0;ci<16;ci++)
 {
-_asm
+read_buhlmann_compartment(0);   // 2 sec mode.
-movlw	0x02
-movwf	TBLPTRH,0
-movlb	5 // fuer ci
-movf    ci,0,1
-addwf	ci,0,1
-addwf	ci,0,1
-addwf	ci,0,1
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var_e2secs+1 // the order is confussing
-TBLRDPOSTINC
-movff	TABLAT,var_e2secs	// low byte first, high afterwards
-TBLRDPOSTINC
-movff	TABLAT,var_e2secs+3
-TBLRD
-movff	TABLAT,var_e2secs+2
-addlw	0x40
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var2_e2secs+1
-TBLRDPOSTINC
-movff	TABLAT,var2_e2secs
-TBLRDPOSTINC
-movff	TABLAT,var2_e2secs+3
-TBLRD
-movff	TABLAT,var2_e2secs+2
-addlw	0x40
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var_a+1
-TBLRDPOSTINC
-movff	TABLAT,var_a
-TBLRDPOSTINC
-movff	TABLAT,var_a+3
-TBLRD
-movff	TABLAT,var_a+2
-addlw	0x40
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var2_a+1
-TBLRDPOSTINC
-movff	TABLAT,var2_a
-TBLRDPOSTINC
-movff	TABLAT,var2_a+3
-TBLRD
-movff	TABLAT,var2_a+2
-addlw	0x40
-movwf	TBLPTRL,0
-incf	TBLPTRH,1,0
-TBLRDPOSTINC
-movff	TABLAT,var_b+1
-TBLRDPOSTINC
-movff	TABLAT,var_b
-TBLRDPOSTINC
-movff	TABLAT,var_b+3
-TBLRD
-movff	TABLAT,var_b+2
-addlw	0x40
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var2_b+1
-TBLRDPOSTINC
-movff	TABLAT,var2_b
-TBLRDPOSTINC
-movff	TABLAT,var2_b+3
-TBLRD
-movff	TABLAT,var2_b+2
-_endasm
-// the start values are the previous end values // write new values in temp
-	if(	(var_e2secs < 0.0000363)
-		|| (var_e2secs > 0.00577)
-		|| (var2_e2secs < 0.0000961)
-		|| (var2_e2secs > 0.150)
-		|| (var_a < 0.231)
-		|| (var_a > 1.27)
-		|| (var_b < 0.504)
-		|| (var_b > 0.966)
-		|| (var2_a < 0.510)
-		|| (var2_a > 1.75)
-		|| (var2_b < 0.423)
-		|| (var2_b > 0.927)
-		)
-		int_O_DBG_pre_bitfield |= DBG_ZH16ERR;
 // N2
-temp_tissue = (temp_atem - pres_tissue[ci]) * var_e2secs;
+temp_tissue = (temp_atem - pres_tissue[ci]) * var_N2_e;
 temp_tissue_safety();
 pres_tissue[ci] = pres_tissue[ci] + temp_tissue;
 // He
-temp_tissue = (temp2_atem - pres_tissue[ci+16]) * var2_e2secs;
+temp_tissue = (temp2_atem - pres_tissue[ci+16]) * var_He_e;
 temp_tissue_safety();
 pres_tissue[ci+16] = pres_tissue[ci+16] + temp_tissue;
 temp_tissue = pres_tissue[ci] + pres_tissue[ci+16];
-var_a = (var_a * pres_tissue[ci] + var2_a * pres_tissue[ci+16]) / temp_tissue;
+var_N2_a = (var_N2_a * pres_tissue[ci] + var_He_a * pres_tissue[ci+16]) / temp_tissue;
-var_b = (var_b * pres_tissue[ci] + var2_b * pres_tissue[ci+16]) / temp_tissue;
+var_N2_b = (var_N2_b * pres_tissue[ci] + var_He_b * pres_tissue[ci+16]) / temp_tissue;
-pres_tissue_limit[ci] = (temp_tissue - var_a) * var_b;
+pres_tissue_limit[ci] = (temp_tissue - var_N2_a) * var_N2_b;
 if (pres_tissue_limit[ci] < 0)
 pres_tissue_limit[ci] = 0;
 if (pres_tissue_limit[ci] > pres_gtissue_limit)
 {
 pres_gtissue_limit = pres_tissue_limit[ci];
 temp_pres_gtissue_limit = 0.0;
 temp_gtissue_no = 255;
 for (ci=0;ci<16;ci++)
 {
-_asm
+read_buhlmann_compartment(1);       // 1 minute interval
-movlw	0x02
-movwf	TBLPTRH,0
-movlb	5 // fuer ci
-movf    ci,0,1
-addwf	ci,0,1
-addwf	ci,0,1
-addwf	ci,0,1
-addlw	0x80
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var_a+1
-TBLRDPOSTINC
-movff	TABLAT,var_a
-TBLRDPOSTINC
-movff	TABLAT,var_a+3
-TBLRD
-movff	TABLAT,var_a+2
-addlw	0x40
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var2_a+1
-TBLRDPOSTINC
-movff	TABLAT,var2_a
-TBLRDPOSTINC
-movff	TABLAT,var2_a+3
-TBLRD
-movff	TABLAT,var2_a+2
-addlw	0x40
-movwf	TBLPTRL,0
-incf	TBLPTRH,1,0
-TBLRDPOSTINC
-movff	TABLAT,var_b+1
-TBLRDPOSTINC
-movff	TABLAT,var_b
-TBLRDPOSTINC
-movff	TABLAT,var_b+3
-TBLRD
-movff	TABLAT,var_b+2
-addlw	0x40
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var2_b+1
-TBLRDPOSTINC
-movff	TABLAT,var2_b
-TBLRDPOSTINC
-movff	TABLAT,var2_b+3
-TBLRD
-movff	TABLAT,var2_b+2
-addlw	0xC0
-movwf	TBLPTRL,0
-incf	TBLPTRH,1,0
-TBLRDPOSTINC
-movff	TABLAT,var_e1min+1
-TBLRDPOSTINC
-movff	TABLAT,var_e1min
-TBLRDPOSTINC
-movff	TABLAT,var_e1min+3
-TBLRD
-movff	TABLAT,var_e1min+2
-addlw	0x40
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var2_e1min+1
-TBLRDPOSTINC
-movff	TABLAT,var2_e1min
-TBLRDPOSTINC
-movff	TABLAT,var2_e1min+3
-TBLRD
-movff	TABLAT,var2_e1min+2
-_endasm
 // N2
-temp_tissue = (temp_atem - sim_pres_tissue[ci]) * var_e1min;
+temp_tissue = (temp_atem - sim_pres_tissue[ci]) * var_N2_e;
 temp_tissue_safety();
 sim_pres_tissue[ci] = sim_pres_tissue[ci] + temp_tissue;
 // He
-temp_tissue = (temp2_atem - sim_pres_tissue[ci+16]) * var2_e1min;
+temp_tissue = (temp2_atem - sim_pres_tissue[ci+16]) * var_He_e;
 temp_tissue_safety();
 sim_pres_tissue[ci+16] = sim_pres_tissue[ci+16] + temp_tissue;
 // pressure limit
 temp_tissue = sim_pres_tissue[ci] + sim_pres_tissue[ci+16];
-var_a = (var_a * sim_pres_tissue[ci] + var2_a * sim_pres_tissue[ci+16]) / temp_tissue;
+var_N2_a = (var_N2_a * sim_pres_tissue[ci] + var_He_a * sim_pres_tissue[ci+16]) / temp_tissue;
-var_b = (var_b * sim_pres_tissue[ci] + var2_b * sim_pres_tissue[ci+16]) / temp_tissue;
+var_N2_b = (var_N2_b * sim_pres_tissue[ci] + var_He_b * sim_pres_tissue[ci+16]) / temp_tissue;
-sim_pres_tissue_limit[ci] = (temp_tissue - var_a) * var_b;
+sim_pres_tissue_limit[ci] = (temp_tissue - var_N2_a) * var_N2_b;
 if (sim_pres_tissue_limit[ci] < 0)
 sim_pres_tissue_limit[ci] = 0;
 if (sim_pres_tissue_limit[ci] > temp_pres_gtissue_limit)
 {
 //--------------------
 // sim_tissue_10min //
 //--------------------
-// Attention!! uses var_e1min und var2_e1min to load 10min data !!!
+// Attention!! uses var_N2_e und var_He_e to load 10min data !!!
 // is identical to sim_tissue_1min routine except for the different load of those variables
 // optimized in v.101
 void sim_tissue_10min(void)
 temp_pres_gtissue_limit = 0.0;
 temp_gtissue_no = 255;
 for (ci=0;ci<16;ci++)
 {
-_asm
+read_buhlmann_compartment(2);       // 10 minute interval
-movlw	0x02
-movwf	TBLPTRH,0
-movlb	5 // fuer ci
-movf    ci,0,1
-addwf	ci,0,1
-addwf	ci,0,1
-addwf	ci,0,1
-addlw	0x80
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var_a+1
-TBLRDPOSTINC
-movff	TABLAT,var_a
-TBLRDPOSTINC
-movff	TABLAT,var_a+3
-TBLRD
-movff	TABLAT,var_a+2
-addlw	0x40
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var2_a+1
-TBLRDPOSTINC
-movff	TABLAT,var2_a
-TBLRDPOSTINC
-movff	TABLAT,var2_a+3
-TBLRD
-movff	TABLAT,var2_a+2
-addlw	0x40
-movwf	TBLPTRL,0
-incf	TBLPTRH,1,0
-TBLRDPOSTINC
-movff	TABLAT,var_b+1
-TBLRDPOSTINC
-movff	TABLAT,var_b
-TBLRDPOSTINC
-movff	TABLAT,var_b+3
-TBLRD
-movff	TABLAT,var_b+2
-addlw	0x40
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var2_b+1
-TBLRDPOSTINC
-movff	TABLAT,var2_b
-TBLRDPOSTINC
-movff	TABLAT,var2_b+3
-TBLRD
-movff	TABLAT,var2_b+2
-addlw	0xC0				// different to 1 min
-movwf	TBLPTRL,0
-incf	TBLPTRH,1,0
-incf	TBLPTRH,1,0			// different to 1 min
-TBLRDPOSTINC
-movff	TABLAT,var_e1min+1
-TBLRDPOSTINC
-movff	TABLAT,var_e1min
-TBLRDPOSTINC
-movff	TABLAT,var_e1min+3
-TBLRD
-movff	TABLAT,var_e1min+2
-addlw	0x40
-movwf	TBLPTRL,0
-//incf	TBLPTRH,1,0			// different to 1 min
-TBLRDPOSTINC
-movff	TABLAT,var2_e1min+1
-TBLRDPOSTINC
-movff	TABLAT,var2_e1min
-TBLRDPOSTINC
-movff	TABLAT,var2_e1min+3
-TBLRD
-movff	TABLAT,var2_e1min+2
-_endasm
 // N2
-temp_tissue = (temp_atem - sim_pres_tissue[ci]) * var_e1min;
+temp_tissue = (temp_atem - sim_pres_tissue[ci]) * var_N2_e;
 temp_tissue_safety();
 sim_pres_tissue[ci] = sim_pres_tissue[ci] + temp_tissue;
 // He
-temp_tissue = (temp2_atem - sim_pres_tissue[ci+16]) * var2_e1min;
+temp_tissue = (temp2_atem - sim_pres_tissue[ci+16]) * var_He_e;
 temp_tissue_safety();
 sim_pres_tissue[ci+16] = sim_pres_tissue[ci+16] + temp_tissue;
 // pressure limit
 temp_tissue = sim_pres_tissue[ci] + sim_pres_tissue[ci+16];
-var_a = (var_a * sim_pres_tissue[ci] + var2_a * sim_pres_tissue[ci+16]) / temp_tissue;
+var_N2_a = (var_N2_a * sim_pres_tissue[ci] + var_He_a * sim_pres_tissue[ci+16]) / temp_tissue;
-var_b = (var_b * sim_pres_tissue[ci] + var2_b * sim_pres_tissue[ci+16]) / temp_tissue;
+var_N2_b = (var_N2_b * sim_pres_tissue[ci] + var_He_b * sim_pres_tissue[ci+16]) / temp_tissue;
-sim_pres_tissue_limit[ci] = (temp_tissue - var_a) * var_b;
+sim_pres_tissue_limit[ci] = (temp_tissue - var_N2_a) * var_N2_b;
 if (sim_pres_tissue_limit[ci] < 0)
 sim_pres_tissue_limit[ci] = 0;
 if (sim_pres_tissue_limit[ci] > temp_pres_gtissue_limit)
 {
 temp_pres_gtissue = temp_tissue;
 // -----------------------
 // calc_gradient_factor //
 // -----------------------
-// optimized in v.101 (var_a)
+// optimized in v.101 (var_N2_a)
 // new code in v.102
 void calc_gradient_factor(void)
 {
 	// tissue > respiration (entsaettigungsvorgang)
 int_O_desaturation_time = 0;
 float_desaturation_multiplier = char_I_desaturation_multiplier / 142.0; // new in v.101	(70,42%/100.=142)
 for (ci=0;ci<16;ci++)
 {
-_asm
+var_N2_halftime = buhlmann_ht[ci];
-movlw	0x04
+var_He_halftime = buhlmann_ht[ci + 16];
-movwf	TBLPTRH,0
-movlb	5 // fuer ci
-movf    ci,0,1
-addwf	ci,0,1
-addwf	ci,0,1
-addwf	ci,0,1
-addlw	0x80
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var_halftimes+1
-TBLRDPOSTINC
-movff	TABLAT,var_halftimes
-TBLRDPOSTINC
-movff	TABLAT,var_halftimes+3
-TBLRD
-movff	TABLAT,var_halftimes+2
-addlw	0x40
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var2_halftimes+1
-TBLRDPOSTINC
-movff	TABLAT,var2_halftimes
-TBLRDPOSTINC
-movff	TABLAT,var2_halftimes+3
-TBLRD
-movff	TABLAT,var2_halftimes+2
-_endasm
 // saturation_time (for flight) and N2_saturation in multiples of halftime
 // version v.100: 1.1 = 10 percent distance to totally clean (totally clean is not possible, would take infinite time )
 // new in version v.101: 1.07 = 7 percent distance to totally clean (totally clean is not possible, would take infinite time )
 // changes in v.101: 1.05 = 5 percent dist to totally clean is new desaturation point for display and noFly calculations
 {
 temp1 = log(1.0 - temp1);
 temp1 = temp1 / -0.6931; // temp1 is the multiples of half times necessary.
 						 // 0.6931 is ln(2), because the math function log() calculates with a base of e not 2 as requested.
 						 // minus because log is negative
-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
+temp2 = var_N2_halftime * temp1 / float_desaturation_multiplier; // time necessary (in minutes ) for complete desaturation (see comment about 10 percent) , new in v.101: float_desaturation_multiplier
 }
 else
 {
 temp1 = 0;
 temp2 = 0;
 	{
 	temp3 = log(1.0 - temp3);
 	temp3 = temp3 / -0.6931; // temp1 is the multiples of half times necessary.
 							 // 0.6931 is ln(2), because the math function log() calculates with a base of e  not 2 as requested.
 							 // minus because log is negative
-	temp4 = var2_halftimes * temp3 / float_desaturation_multiplier; // time necessary (in minutes ) for "complete" desaturation, new in v.101 float_desaturation_multiplier
+	temp4 = var_He_halftime * temp3 / float_desaturation_multiplier; // time necessary (in minutes ) for "complete" desaturation, new in v.101 float_desaturation_multiplier
 	}
 else
 	{
 	temp3 = 0;
 	temp4 = 0;
 char_O_gtissue_no = 255;
 pres_gtissue_limit = 0.0;
 for (ci=0;ci<16;ci++)
 {
-_asm
+read_buhlmann_compartment(1);       // 1 minute interval
-movlw	0x02
-movwf	TBLPTRH,0
-movlb	5 // fuer ci
-movf    ci,0,1
-addwf	ci,0,1
-addwf	ci,0,1
-addwf	ci,0,1
-addlw	0x80
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var_a+1
-TBLRDPOSTINC
-movff	TABLAT,var_a
-TBLRDPOSTINC
-movff	TABLAT,var_a+3
-TBLRD
-movff	TABLAT,var_a+2
-addlw	0x40
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var2_a+1
-TBLRDPOSTINC
-movff	TABLAT,var2_a
-TBLRDPOSTINC
-movff	TABLAT,var2_a+3
-TBLRD
-movff	TABLAT,var2_a+2
-addlw	0x40
-movwf	TBLPTRL,0
-incf	TBLPTRH,1,0
-TBLRDPOSTINC
-movff	TABLAT,var_b+1
-TBLRDPOSTINC
-movff	TABLAT,var_b
-TBLRDPOSTINC
-movff	TABLAT,var_b+3
-TBLRD
-movff	TABLAT,var_b+2
-addlw	0x40
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var2_b+1
-TBLRDPOSTINC
-movff	TABLAT,var2_b
-TBLRDPOSTINC
-movff	TABLAT,var2_b+3
-TBLRD
-movff	TABLAT,var2_b+2
-addlw	0xC0
-movwf	TBLPTRL,0
-incf	TBLPTRH,1,0
-TBLRDPOSTINC
-movff	TABLAT,var_e1min+1
-TBLRDPOSTINC
-movff	TABLAT,var_e1min
-TBLRDPOSTINC
-movff	TABLAT,var_e1min+3
-TBLRD
-movff	TABLAT,var_e1min+2
-addlw	0x40
-movwf	TBLPTRL,0
-TBLRDPOSTINC
-movff	TABLAT,var2_e1min+1
-TBLRDPOSTINC
-movff	TABLAT,var2_e1min
-TBLRDPOSTINC
-movff	TABLAT,var2_e1min+3
-TBLRD
-movff	TABLAT,var2_e1min+2
-_endasm
 // N2 1 min
-temp_tissue = (temp_atem - pres_tissue[ci]) * var_e1min;
+temp_tissue = (temp_atem - pres_tissue[ci]) * var_N2_e;
 temp_tissue_safety();
 pres_tissue[ci] = pres_tissue[ci] + temp_tissue;
 // He 1 min
-temp_tissue = (temp2_atem - pres_tissue[ci+16]) * var2_e1min;
+temp_tissue = (temp2_atem - pres_tissue[ci+16]) * var_He_e;
 temp_tissue_safety();
 pres_tissue[ci+16] = pres_tissue[ci+16] + temp_tissue;
 temp_tissue = pres_tissue[ci] + pres_tissue[ci+16];
-var_a = (var_a * pres_tissue[ci] + var2_a * pres_tissue[ci+16]) / temp_tissue;
+var_N2_a = (var_N2_a * pres_tissue[ci] + var_He_a * pres_tissue[ci+16]) / temp_tissue;
-var_b = (var_b * pres_tissue[ci] + var2_b * pres_tissue[ci+16]) / temp_tissue;
+var_N2_b = (var_N2_b * pres_tissue[ci] + var_He_b * pres_tissue[ci+16]) / temp_tissue;
-pres_tissue_limit[ci] = (temp_tissue - var_a) * var_b;
+pres_tissue_limit[ci] = (temp_tissue - var_N2_a) * var_N2_b;
 if (pres_tissue_limit[ci] < 0)
 pres_tissue_limit[ci] = 0;
 if (pres_tissue_limit[ci] > pres_gtissue_limit)
 {
 pres_gtissue_limit = pres_tissue_limit[ci];
 md_state[md_i] = 0;
 char_O_hash[md_i] = 0;
 } // for md_i 16
 _asm
-movlw	0x01
+movlw	0x01            // md_pi address.
 movwf	TBLPTRU,0
 movlw	0x06
 movwf	TBLPTRH,0
 movlw	0x00
 movwf	TBLPTRL,0
 _endasm;
-for(md_i=0;md_i<=255;md_i++)
+md_i = 0;
-{
+do {
 _asm
 TBLRDPOSTINC
 movff	TABLAT,md_temp
 _endasm
-md_pi_subst[md_i] = md_temp;
+md_pi_subst[md_i++] = md_temp;
-} // for md_i 256
+} while( md_i != 0 );
 _asm
 movlw	0x00
 movwf	TBLPTRU,0
 movlw	0x00
 movwf	TBLPTRH,0

Mercurial > public > mk2

comparison code_part1/OSTC_code_c_part2/p2_deco.c @ 165:8b5b848db4cd