/COPYRIGHT 1970, DIGITAL EQUIPMENT CORP., MAYNARD, MASS. /EDIT 3 9-25-70 .TITLE .BF RAISE REAL BASE(A) TO DOUBLE POWER (B) A**B .GLOBL .BF,.CB,.DE,.AS,.DF,.ER,.AB / CALLING SEQUENCE (A**B = EXP(B * LOGE(A))) / JMS* (.BF) SUBR CALL (A IN FLOAT ACC) / CAL/XCT ADDR ADDR OF B (XCT IF INDIRECT / NEXT INSTRUCTION SUBR RETURN (RESULT IN FLOAT ACC) / .BF CAL 0 /ENTRY-EXIT JMS* .CB /SHORT GET ARG (.CB) BF02 CAL 0 /ADDR OF ARG B LAC* .AB /GET BASE SZA JMP BF06 ISZ BF02 /PASS BY EXP. LAC* BF02 SPA!SNA!CLA JMS* .ER /OTS 15 AND 15 /DUMMY INSTR. JMP* .BF /ANS=BASE BF06 JMS* .DE /COMPUTE LOG2(A) JMS* .AS /MULTIPLY DOUBLE (.AS) .DSA BF05 / (LOG2(A) * LOGE(2) = LOGE(A)) JMS* .AS /MULTIPLY DOUBLE (.AS) .DSA BF02+400000 / (B X LOGE(A)) JMS* .DF /COMPUTE EXP (.DF) (EXP(B * LOGE(A))) JMP* .BF /EXIT BF05 0 /LOGE(2) (1) (0.6931471806) 261344 / (2) 137700 / (3) .END