TweetFollow Us on Twitter

68881 Access
Volume Number:4
Issue Number:4
Column Tag:Forth Forum

Direct 68881 Floating Point Access

By Jörg Langowski, MacTutor Editorial Board, Grenoble, France

Being the proud owner of a MacII since the beginning of this year, I was - like so many others - disappointed at the relative slowness of the SANE package. Even with the floating point coprocessor, the speedup is only a factor of 5-10 compared to the Mac SE. The coprocessor itself allows you to go much faster than that, and it is Apple’s SANE implementation that slows down the operations.

The reason for this is, of course, that Apple is trying to guarantee the ‘compatibility’ of the SANE on the MacII with the old SANE implementation. That means that the results are supposed to be the same, down to the least significant bit (or even to the last guard bit?). The 68881 uses different algorithms for calculating transcendental functions than SANE does, therefore some of its built-in operations could not be used and had to be replaced by software. Of course, that slows down things a lot.

As a side note, I’m not at all in favor of such a strategy. Artificially restraining a high performance chip just because its results, although accurate enough, don’t match the old - also accurate - results to the last bit seems a little exaggerated. Stable numerical algorithms should take into account the possibility that the hardware changes slightly and the machine errors are different, and they should be immune against such changes; or at least the new SANE should have had an option built in that uses the 68881 directly! [Amen to that statement! -Ed]

Many of the development systems for the Macintosh now come with the possibility to generate code that supports the floating point coprocessor directly, many others don’t; working with Forth, we don’t have a real problem since we can easily redefine our floating point operators. This month I’ll show you how to do that.

The 68881

Let’s first have a look at the floating point coprocessor itself.

The 68881 is accessed from the 68020 in a special address area. When the 68020 encounters an instruction of the form $Fxxx (previously the F line trap), it will set its function code lines (pins FC0-FC2) all to high, indicating ‘CPU space’. It will then exchange information with the coprocessor’s internal registers to perform the floating point operation requested.

This information exchange occurs automatically, it is part of the 68020’s design; when using the floating point instructions, you don’t notice the communication between the two processors at all. The 68881 appears as an extension to the 68020, just as if we had a new set of registers available, with special instructions operating on them.

The registers FP0-FP7 can each hold one extended precision (80-bit) floating point number, in the format given by the IEEE standard: bits 0-63 contain the mantissa, 64-78 the 15-bit exponent offset by 16383, and bit 79 the sign. The instruction set of the 68881 allows you to do floating point operations on any or between any two of these registers, and transfer data between them, the 68020’s registers and memory.

Mach2’s assembler fully supports all 68881 instructions. The redefinition of the words in the SANE vocabulary is therefore quite straightforward. Register D7 is used by Mach2 as the floating point stack pointer, all we have to do is to transfer floating point numbers from the FP stack to the 68881’s registers, do the FP operation and transfer back the result.

The 68881 expects floating point numbers in a different format than SANE does. In order to keep FP numbers aligned to the boundaries of a 32-bit long word, which speeds up operations, an extended precision number in memory will be 96 bits long instead of 80; there is a 16 bit gap between the 64-bit mantissa and the 16-bit sign and exponent field. Since SANE and the floating point stack in Mach2 use 80 bit numbers, we have to convert their format when we use the coprocessor.

Listing 1 provides the macros f>2, f>1 and 1>f for this purpose. The general format for a binary floating point operation is then

 f>2  
 fop.x  FPn, FPm 
 1>f

where fop.x is an extended format floating point instruction that operates on floating point registers FPn and FPm. The example always uses registers FP0 and FP1. Unary operations are encoded the same way, using only FP0.

In the example, I provided a new vocabulary f68881 that contains redefinitions of the most important SANE operations for the 68881. Of course, the concept can be far extended. We have eight registers at our disposition and can use them to optimize more complicated numerical algorithms, using only assembly language. I might give some examples how to do this in a later column.

For the moment, let’s content ourselves with the speed improvement that we have achieved so far, which is already remarkable. Some simple benchmarks are listed at the end of listing 1, and the results are given here:

( MacII, direct 68881 access, 100000 loops each )

bmark1 2 Secs 18 Ticks ok <0> [0]

bmark2 4 Secs 48 Ticks ok <0> [0]

bmark3 4 Secs 14 Ticks ok <0> [0]

bmark4 4 Secs 17 Ticks ok <0> [0]

bmark5 4 Secs 14 Ticks ok <0> [0]

( MacII, SANE w/68881, 10000 loops each )

smark1 15 Ticks ok <0> [0]

smark2 37 Secs 12 Ticks ok <0> [0]

smark3 24 Ticks ok <0> [0]

smark4 1 Secs 33 Ticks ok <0> [0]

smark5 1 Secs 33 Ticks ok <0> [0]

( Mac +, SANE, 10000 loops each )

ok <0> [0]

smark1 24 Ticks ok <0> [0]

smark2 181 Secs 54 Ticks ok <0> [0]

smark3 49 Ticks ok <0> [0]

smark4 4 Secs 47 Ticks ok <0> [0]

smark5 5 Secs 14 Ticks ok <0> [0]

For the Mac+, the fnull1 and fnull2 operations had been replaced by simple fdrops. As you see, the speedup going from Mac+ to MacII’s SANE is not so breathtaking: a factor of 6 for the exponential, 4 for addition and subtraction; but when we access the 68881 directly, we gain another factor of 3 for simple addition and multiplication and 78 for the exponential. It is in the calculation of the transcendentals where the 68881 really shines.

Pop up menus

Someone approached me lately on the question of how to do pop up menu selection. Since the technical notes contain only sketchy references to popup menus at the time I write this (e.g. TN156, TN172), I’d like to give you a practical example how to use pop up menus from Mach2 in a simple way.

Listing 2 explains the process. The PopUpMenuSelect trap takes four parameters:

- a handle to the menu to be displayed (32 bits),

- the top and left global coordinate of the point at which to display the menu (2*16 bits),

- the menu item which should be positioned at that point for the default selection.

Although Mach2 knows the trap name, the interface to this routine is not (yet) correct, so we have to redefine it in assembler. Note that the point returned by the @mouse function is in local coordinates, while PopUpMenuSelect expects it in global coordinates.

The example defines a menu using the Mach2 interface; the menu is created with -1 as the insertion parameter (-1 150 mymenu BOUNDS) so that after insertion into the tasks’s menu bar the menu will stay invisible (just as we did for the hierarchical menus). Note that a pop up menu has to be inserted into the menu bar before using it.

The content handler for the default Mach2 window is then rewritten so that on a mouse down event it will select the example pop up menu. The menu handler will just beep a number of times depending on the item selected. dopop activates the new content handler while nopop deactivates it.

Feedback dept.

This letter comes from Vassili Dzuba, Paris:

“In January’s issue ‘Mousehole Report’, Alan Dall put in his wish list the ability to define ‘ghost copies’ of applications. Even without Unix’ capability of defining links, this can be done with a small program using the _Launch trap. This program takes only 1K on the disk. The path name of the application to launch is stored in ‘STR ‘ resource 1000. It’s possible to set the creator and the bundle bit of the ghost to have it share the same icon as the real thing. Of course, double-clicking on a document can then launch the ghost instead of the application, but the slowing down is only marginal.

The program is the following

(using MPW’s assembler):

 INCLUDE‘traps.a’
ghost MAIN
 MOVE.W #0,-(SP) 
 ; the context data (_launch parameter)
 MOVE.L #’STR ‘,-(SP)
 ; 1st parameter of GetResource
 MOVE.W #1000,-(SP)
 ; 2nd parameter of GetResource
 _GetResource
 ; handle to string in (SP)
 TST.L  (SP)
 ; test if null handle (no resource available)
 BEQ.B  exit
 ; if null, go to exit
 MOVE.L (SP),A0
 MOVE.L (A0),(SP)
 ; handle dereferenced
 MOVE.L SP,A0
 ; stack pointer in A0
 _Launch
exit  _ExitToShell
 END

The resource file is something like this (using Rez format):

#include “types.r”
resource ‘STR ‘ (1000) {
 “Sys:myDirectory:MyApp”
};

A ghost can be easily created using a small shell script (assuming the original ghost’s directory being {MPW}dev) which sets up the string resource:

duplicate {MPW}’dev:ghost’ {2}
echo ‘#include “types.r”’n
‘resource ‘’’STR ‘’’ (1000) { “‘{1}’” };’
| rez -a -o {2}

Assuming this script is named ‘summon’, the creation in the current directory of a ghost of MacPaint would be something like:

summon ‘sys:appli ƒ:mac paint’ ‘macpaint.Ghost’

Sincerely yours”

Thank you, Vassili, for this helpful little utility.


Listing 1: direct access 68881 floating point words for Mach2
\ 68881 access, © J. Langowski/MacTutor Jan 1988

only forth also assembler also sane
vocabulary f68881
also f68881 definitions

code f>2
 add.l  #20,d7
 move.l d7,a0
 move.l -(a0),-(a7); move mantissa
 move.l -(a0),-(a7); in two 32-bit chunks
 subq.l #2,a7    ; 16-bit gap
 move.w -(a0),-(a7); move exponent + sign
 fmove.x(a7)+,fp0; transfer from stack to fp0
 move.l -(a0),-(a7); same for fp1...
 move.l -(a0),-(a7)
 subq.l #2,a7
 move.w -(a0),-(a7)
 fmove.x(a7)+,fp1
 add.l  #10,a0 ; a0 points to 2nd fp stack item
 rts
end-code mach

code f>1
 add.l  #10,d7
 move.l d7,a0
 move.l -(a0),-(a7)
 move.l -(a0),-(a7)
 subq.l #2,a7
 move.w -(a0),-(a7)
 fmove.x(a7)+,fp0
 rts
end-code mach

code 1>f
 fmove.xfp0,-(a7)
 move.w (a7)+,(a0)+; transfer exponent + sign
 addq.l #2,a7    ; skip16 bit gap
 move.l (a7)+,(a0)+; transfer mantissa 
 move.l (a7)+,(a0)+; in 2 steps
 sub.l  #10,d7   ; adjust FP stack pointer
 rts
end-code mach

( note: f>1 or f>2 and 1>f should always occur in pairs since 1>f expects 
A0 to point to second floating point stack position; this is assured 
by f>1 and f>2 )

code f+
 f>2
 fadd.x fp1,fp0
 1>f
 rts
end-code

code f-
 f>2
 fsub.x fp1,fp0
 1>f
 rts
end-code

code f/
 f>2
 fdiv.x fp1,fp0
 1>f
 rts
end-code

code f*
 f>2
 fmul.x fp1,fp0
 1>f
 rts
end-code

code fmod
 f>2
 fmod.x fp1,fp0
 1>f
 rts
end-code

code frem
 f>2
 frem.x fp1,fp0
 1>f
 rts
end-code

code fabs
 f>1
 fabs.x fp0
 1>f
 rts
end-code

code facos
 f>1
 facos.xfp0
 1>f
 rts
end-code

code fasin
 f>1
 fasin.xfp0
 1>f
 rts
end-code

code fatan
 f>1
 fatan.xfp0
 1>f
 rts
end-code

code fatanh
 f>1
 fatanh.x fp0
 1>f
 rts
end-code

code fcos
 f>1
 fcos.x fp0
 1>f
 rts
end-code

code fcosh
 f>1
 fcosh.xfp0
 1>f
 rts
end-code

code fe^x
 f>1
 fetox.xfp0
 1>f
 rts
end-code
 
code fe^x-1
 f>1
 fetoxm1.x fp0
 1>f
 rts
end-code

code fgetexp
 f>1
 fgetexp.x fp0
 1>f
 rts
end-code

code fgetman
 f>1
 fgetexp.x fp0
 1>f
 rts
end-code

code fint
 f>1
 fint.x fp0
 1>f
 rts
end-code

code fintrz
 f>1
 fintrz.x fp0
 1>f
 rts
end-code

code fln
 f>1
 flogn.x  fp0
 1>f
 rts
end-code

code fln+1
 f>1
 flognp1.x fp0
 1>f
 rts
end-code

code flog10
 f>1
 flog10.x fp0
 1>f
 rts
end-code

code flog2
 f>1
 flog2.x  fp0
 1>f
 rts
end-code

code fneg
 f>1
 fneg.x fp0
 1>f
 rts
end-code

code fsin
 f>1
 fsin.x fp0
 1>f
 rts
end-code

code fsinh
 f>1
 fsinh.x  fp0
 1>f
 rts
end-code

code fsqrt
 f>1
 fsqrt.x  fp0
 1>f
 rts
end-code

code ftan
 f>1
 ftan.x fp0
 1>f
 rts
end-code

code ftanh
 f>1
 ftanh.x  fp0
 1>f
 rts
end-code

code f10^x
 f>1
 ftentox.x fp0
 1>f
 rts
end-code

code f2^x
 f>1
 ftwotox.x fp0
 1>f
 rts
end-code

fp

also forth definitions
code fnull1
 f>1
 1>f
 rts
end-code

code fnull2
 f>2
 1>f
 rts
end-code

: bmark1 counter 1.0 
 100000 0 do fdup fnull1 fdrop loop 
timer fdrop ;

: bmark2 counter 4.3352 
 100000 0 do fdup fe^x fdrop loop 
timer fdrop ;

: bmark3 counter 
 3.5 4.5 100000 0 do fover fover fnull2 fdrop loop 
timer fdrop fdrop ;

: bmark4 counter 
 3.5 4.5 100000 0 do fover fover f+ fdrop loop 
timer fdrop fdrop ;

: bmark5 counter 
 3.5 4.5 100000 0 do fover fover f* fdrop loop 
timer fdrop fdrop ;

also sane

: smark1 counter 
 1.0 10000 0 do fdup fnull1 fdrop loop 
timer fdrop ;

: smark2 counter 
 4.3352 10000 0 do fdup fe^x fdrop loop 
timer fdrop ;

: smark3 counter 
 3.5 4.5 10000 0 do fover fover fnull2 fdrop loop 
timer fdrop fdrop ;

: smark4 counter 
 3.5 4.5 10000 0 do fover fover f+ fdrop loop 
timer fdrop fdrop ;

: smark5 counter 
 3.5 4.5 10000 0 do fover fover f* fdrop loop 
timer fdrop fdrop ;
Listing 2: Pop up menus
\ popup menu interface from Forth
\ © J. Langowski/MacTutor Jan 1988

\ FUNCTION PopUpMenuSelect 
\ (menu: MenuHandle; top, left, popUpItem: INTEGER):
\LONGINT;
\  INLINE $A80B;

code popup ( hMenu top left item# | menuID item  -- )
 exg  d4,a7
 clr.l  -(a7)
 move.l 12(a6),-(a7) \ handle
 move.w 10(a6),-(a7) \ top
 move.w 6(a6),-(a7)\ left
 move.w 2(a6),-(a7)\ item
 adda.l #16,A6
 _popupmenuselect
 clr.l  d0
 clr.l  d1
 move.w (a7)+,d0
 move.w (a7)+,d1
 exg  d4,a7
 move.l d0,-(a6)
 move.l d1,-(a6)
 rts
end-code

CODE  @MOUSE
    SUBQ.L    #4,A6
    MOVE.L    A6,-(A7)
    _GETMOUSE
    RTS
END-CODE

CODE unpack
 MOVE.L (A6),D0
 CLR.L  D1
 MOVE.W D0,D1
 CLR.W  D0
 SWAP.W D0
 MOVE.L D0,(A6)
 MOVE.L D1,-(A6)
 RTS
END-CODE MACH

108 user taskmenubar
152 user content-hook

NEW.MENU myMenu
“ TestMenu” myMenu TITLE
-1 150 myMenu BOUNDS \ invisible menu
“ Item 1;Item 2;Item 3” myMenu ITEMS

taskmenubar @ mymenu add
call drawmenubar

: beep 
 ?dup if 
 0 do 5 call sysbeep loop 
 then
;

: content-handler { | pt -- }
 @mouse -> pt
 ^ pt call localtoglobal
 mymenu @
 pt unpack
 1
 popup
 beep drop
 run-content
;

: dopop [‘] content-handler content-hook ! ;
: nopop [‘] run-content content-hook ! ;
 

Community Search:
MacTech Search:

Software Updates via MacUpdate

Adobe Animate CC 2017 16.2.0 - Advanced...
Animate CC 2017 is available as part of Adobe Creative Cloud for as little as $19.99/month (or $9.99/month if you're a previous Flash Professional customer). Animate CC 2017 (was Flash CC) lets you... Read more
Tunnelblick 3.7.0 - GUI for OpenVPN.
Tunnelblick is a free, open source graphic user interface for OpenVPN on OS X. It provides easy control of OpenVPN client and/or server connections. It comes as a ready-to-use application with all... Read more
DEVONthink Pro 2.9.11 - Knowledge base,...
DEVONthink Pro is your essential assistant for today's world, where almost everything is digital. From shopping receipts to important research papers, your life often fills your hard drive in the... Read more
DiskCatalogMaker 6.8.1 - Catalog your di...
DiskCatalogMaker is a simple disk management tool which catalogs disks. Simple, light-weight, and fast Finder-like intuitive look and feel Super-fast search algorithm Can compress catalog data for... Read more
OmniGraffle 7.3 - Create diagrams, flow...
OmniGraffle helps you draw beautiful diagrams, family trees, flow charts, org charts, layouts, and (mathematically speaking) any other directed or non-directed graphs. We've had people use Graffle to... Read more
OmniGraffle Pro 7.3 - Create diagrams, f...
OmniGraffle Pro helps you draw beautiful diagrams, family trees, flow charts, org charts, layouts, and (mathematically speaking) any other directed or non-directed graphs. We've had people use... Read more
A Better Finder Rename 10.15 - File, pho...
A Better Finder Rename is the most complete renaming solution available on the market today. That's why, since 1996, tens of thousands of hobbyists, professionals and businesses depend on A Better... Read more
Adobe Dreamweaver CC 2017 17.0.2.9391 -...
Dreamweaver CC 2017 is available as part of Adobe Creative Cloud for as little as $19.99/month (or $9.99/month if you're a previous Dreamweaver customer). Adobe Dreamweaver CC 2017 allows you to... Read more
BBEdit 11.6.5 - Powerful text and HTML e...
BBEdit is the leading professional HTML and text editor for the Mac. Specifically crafted in response to the needs of Web authors and software developers, this award-winning product provides a... Read more
BBEdit 11.6.5 - Powerful text and HTML e...
BBEdit is the leading professional HTML and text editor for the Mac. Specifically crafted in response to the needs of Web authors and software developers, this award-winning product provides a... Read more

Meganoid(2017) (Games)
Meganoid(2017) 1.0 Device: iOS Universal Category: Games Price: $3.99, Version: 1.0 (iTunes) Description: LAUNCH DISCOUNT 20% UNTIL APRIL 2nd! Support, tip and tricks: http://www.orangepixel.net/forum/ Subscribe to our newsletter... | Read more »
Telltale's Guardians of the Galaxy...
Telltale will be releasing their rendition of Guardians of the Galaxy later this month. The first episode, Tangled Up in Blue, features familiar faces including Star-Lord, Groot, Rocket, Gamora, and Drax. If the first episode's title is any... | Read more »
Royal Dungeon (Games)
Royal Dungeon 1.0 Device: iOS Universal Category: Games Price: $2.99, Version: 1.0 (iTunes) Description: The king and his queen are trapped in their castle which suddenly turned out as a very dangerous place. The goal is to escape... | Read more »
Tom Clancy's ShadowBreak is a real-...
Ubisoft is treating Tom Clancy fans to the series' very first mobile-exclusive game in ShadowBreak, a real-time, multiplayer shooter in which players snipe at enemies in fast-paced tactics-driven combat. [Read more] | Read more »
Power Rangers: Legacy Wars beginner...
Rita Repulsa is back, but this time she's invading your mobile phone in Power Rangers: Legacy Wars. What looks to be a straightforward beat 'em up is actually a tough-as-nails multiplayer strategy game that requires some deft tactical maneuvering.... | Read more »
Hearthstone celebrates the upcoming Jour...
Hearthstone gets a new expansion, Journey to Un'Goro, in a little over a week, and they'll be welcoming the Year of the Mammoth, the next season, at the same time. There's a lot to be excited about, so Blizzard is celebrating in kind. Players will... | Read more »
4 smart and stylish puzzle games like Ty...
TypeShift launched a little over a week ago, offering some puzzling new challenges for word nerds equipped with an iOS device. Created by Zach Gage, the mind behind Spelltower, TypeShift boasts, like its predecessor, a sleak design and some very... | Read more »
The best deals on the App Store this wee...
Deals, deals, deals. We're all about a good bargain here on 148Apps, and luckily this was another fine week in App Store discounts. There's a big board game sale happening right now, and a few fine indies are still discounted through the weekend.... | Read more »
The best new games we played this week
It's been quite the week, but now that all of that business is out of the way, it's time to hunker down with some of the excellent games that were released over the past few days. There's a fair few to help you relax in your down time or if you're... | Read more »
Orphan Black: The Game (Games)
Orphan Black: The Game 1.0 Device: iOS Universal Category: Games Price: $4.99, Version: 1.0 (iTunes) Description: Dive into a dark and twisted puzzle-adventure that retells the pivotal events of Orphan Black. | Read more »

Price Scanner via MacPrices.net

Is Apple Planning An iPhone Based Modular Doc...
Today’s more powerful and larger-screened smartphones and phablets are becoming the default anchor computing device for more and more users computing devices, but even a five or six inch panel is not... Read more
Razer Launches New Razer Blade Pro World’s Fi...
Razer, the gaming and high performance hardware specialists, have announced the new Razer Blade Pro laptop — the first laptop to be qualified for THX Mobile Certification, an accreditation reserved... Read more
Gro CRM’s Apple Small Business Mac And iOS CR...
Gro Software, developers of the Mac CRM software for small business and enterprise, are included in FinancesOnline 2017 CRM Rising Stars and Great User Experience lists by business software review... Read more
Deal alert! 15-inch and 13-inch MacBook Pros...
B&H Photo has the new 2016 15″ and 13″ Apple MacBook Pros in stock today and on sale for up to $200 off MSRP. Shipping is free, and B&H charges NY sales tax only: - 15″ 2.7GHz Touch Bar... Read more
Save up to $420 on a new MacBook Pro with App...
Apple is offering Certified Refurbished 2016 15″ and 13″ MacBook Pros, including some Touch Bar models, for up to $420 off original MSRP. An Apple one-year warranty is included with each model, and... Read more
12-inch 1.2GHz Retina MacBooks on sale for $1...
B&H has 12″ 1.2GHz Retina MacBooks on sale for up to $200 off MSRP. Shipping is free, and B&H charges NY sales tax only: - 12″ 1.2GHz Space Gray Retina MacBook: $1449 $150 off MSRP - 12″ 1.... Read more
Is A New 10.5-inch iPad Still Coming In April...
There was no sign or mention of a long-rumored and much anticipated 10.5-inch iPad Pro in Apple’s product announcements last week. The exciting iPad news was release of an upgraded iPad Air with a... Read more
T-Mobile’s Premium Device Protection Now Incl...
Good news for T-Mobile customers who love their iPhones and iPads. The “Un-carrier” has become the first national wireless company to give customers AppleCare Services at zero additional cost as part... Read more
FileWave Ensures Support for Latest Apple OS...
FileWave multi-platform device management providers announced support for Apple’s release yesterday of iOS 10.3, macOS Sierra 10.12.4, and tvOS 11.2. FileWave has a history of providing zero-day... Read more
Use Apple’s Education discount to save up to...
Purchase a new Mac or iPad using Apple’s Education Store and take up to $300 off MSRP. All teachers, students, and staff of any educational institution qualify for the discount. Shipping is free: -... Read more

Jobs Board

Fulltime aan de slag als shopmanager in een h...
Ben jij helemaal gek van Apple -producten en vind je het helemaal super om fulltime shopmanager te zijn in een jonge en hippe elektronicazaak? Wil jij werken in Read more
Fulltime aan de slag als shopmanager in een h...
Ben jij helemaal gek van Apple -producten en vind je het helemaal super om fulltime shopmanager te zijn in een jonge en hippe elektronicazaak? Wil jij werken in Read more
Desktop Analyst - *Apple* Products - Montef...
…technology to improve patient care. JOB RESPONSIBILITIES: Provide day-to-day support for Apple Hardware and Software in the environment based on the team's support Read more
*Apple* Mobile Master - Best Buy (United Sta...
**493168BR** **Job Title:** Apple Mobile Master **Location Number:** 000827-Denton-Store **Job Description:** **What does a Best Buy Apple Mobile Master do?** At Read more
Fulltime aan de slag als shopmanager in een h...
Ben jij helemaal gek van Apple -producten en vind je het helemaal super om fulltime shopmanager te zijn in een jonge en hippe elektronicazaak? Wil jij werken in Read more
All contents are Copyright 1984-2011 by Xplain Corporation. All rights reserved. Theme designed by Icreon.