+- HP Forums (https://www.hpmuseum.org/forum)
+-- Forum: HP Software Libraries (/forum-10.html)
+--- Forum: HP-41C Software Library (/forum-11.html)
+--- Thread: (41) Jet fuel Density calculator (/thread-18068.html)
(41) Jet fuel Density calculator - WrongWay - 02-23-2022 05:38 AM
The attached program computes the fuel density of JET A fuel and compute the needed volumetric fuel uplift. This is my first time back to the 41 after a fair bit so the code I'm sure could easily be optimized and constructive feedback is always welcome. The compressed raw file is attached.
XEQ FDENSTY to run.
and follow on screen instructions. The initial prompt asks for an entry of 0 or 1 to set the desired units.
size[ 00E4 ] ( 228 bytes )
Code:
LBL "FDENSTY"
CLRG
CLST
CLA
CF 00
FIX 3
"LB: 0, KG:1"
PROMPT
X=0?
SF 00
1.013879343
STO 01
-2.987077809E-07
STO 02
-9.200271003E-04
STO 03
"OAT C?"
PROMPT
STO 00
*
X<>Y
RCL 00
X^2
*
+
+
FS? 00
GTO "L"
GTO "K"
LBL "L"
100
15
/
*
STO 01
GTO 00
LBL "K"
0.8
*
STO 01
GTO 00
LBL 00
"DEN: "
ARCL X
FS? 00
GTO 01
GTO 02
LBL 01
>"LB/GAL"
GTO 03
LBL 02
>"KG/M3"
GTO 03
LBL 03
STO 01
AVIEW
STOP
FIX 0
"FPL FUEL?"
PROMPT
"FOB?"
PROMPT
-
RCL 01
/
"NEED: "
ARCL X
FS?C 00
GTO 04
GTO 05
LBL 04
>"G"
GTO 06
LBL 05
>"L"
GTO 06
LBL 06
AVIEW
FIX 9
ENDRE: (41) Jet fuel Density calculator - Thomas Klemm - 06-25-2022 03:47 AM
(02-23-2022 05:38 AM)WrongWay Wrote: … the code I'm sure could easily be optimized and constructive feedback is always welcome.
if-else Statements
In general an if-else statement like this:
Code:
if predicate:
# true-case
...
else:
# false-case
...… can be translated into:
Code:
predicate?
GTO 00
# false-case
...
GTO 01
LBL 00
# true-case
...
LBL 01This uses only two labels instead of three:
Code:
33 FS? 00
34 GTO 00
35 ├"KG/M3"
36 GTO 01
37▸LBL 00
38 ├"LB/GAL"
39▸LBL 01Since we use only forward jumps, both labels can be reused in other occurrences as well:
Code:
52 FS?C 00
53 GTO 00
54 ├"L"
55 GTO 01
56 ├"G"
57▸LBL 01This also avoids using the expensive global labels "K" and "L".
Polynomial Evaluation
A quadratic polynomial can be rewritten using Horner's method:
\(
\begin{align}
ax^2 + bx + c = (ax + b)x + c
\end{align}
\)
This avoids calculating powers of \(x\).
If we first fill the stack with \(x\) it will be copied from register T down the stack with each operation.
Since the coefficients are not used elsewhere we don't have to store them in registers:
Code:
12 ENTER
13 ENTER
14 ENTER
15 -2.987077809ᴇ-07
16 ×
17 -9.200271003ᴇ-04
18 +
19 ×
20 1.013879343
21 +Saving Results
I decided not to save the density in a register but leave the result on the stack for further computations:
Code:
43 "FPL FUEL?"
44 PROMPT
45 "FOB?"
46 PROMPT
47 -
48 X<>Y
49 ÷
50 "NEED: "
51 ARCL ST XOf course the user might mess with the stack which may lead to wrong results.
This leaves us with the following program for the HP-42S:
Code:
00 { 191-Byte Prgm }
01▸LBL "FDENSTY"
02 CLST
03 CLA
04 CF 00
05 FIX 03
06 "LB: 0, KG:1"
07 PROMPT
08 X=0?
09 SF 00
10 "OAT C?"
11 PROMPT
12 ENTER
13 ENTER
14 ENTER
15 -2.987077809ᴇ-07
16 ×
17 -9.200271003ᴇ-04
18 +
19 ×
20 1.013879343
21 +
22 FS? 00
23 GTO 00
24 .8
25 ×
26 GTO 01
27▸LBL 00
28 .15
29 ÷
30▸LBL 01
31 "DEN: "
32 ARCL ST X
33 FS? 00
34 GTO 00
35 ├"KG/M3"
36 GTO 01
37▸LBL 00
38 ├"LB/GAL"
39▸LBL 01
40 AVIEW
41 STOP
42 FIX 00
43 "FPL FUEL?"
44 PROMPT
45 "FOB?"
46 PROMPT
47 -
48 X<>Y
49 ÷
50 "NEED: "
51 ARCL ST X
52 FS?C 00
53 GTO 00
54 ├"L"
55 GTO 01
56 ├"G"
57▸LBL 01
58 AVIEW
59 FIX 09
60 ENDIt should also work with the HP-41C after the obvious transformations.
We end up with a program with 20 fewer lines, using 191 bytes instead of 230, and using no registers.
RE: (41) Jet fuel Density calculator - WrongWay - 05-12-2024 09:36 PM
Thank you for the awesome feedback and excellent explanation.