TweetFollow Us on Twitter

Background TCL
Volume Number:9
Issue Number:12
Column Tag:TCL Workshop

A Piece of Apple Pi

Background processing with TCL

By Malcolm H. Teas, Rye, New Hampshire

Note: Source code files accompanying article are located on MacTech CD-ROM or source code disks.

About the author

Malcolm Teas has been programming the Macintosh for five years. He’s active with object programming in TCL and also uses C with Think and MPW. He lives, works, and consults from his house on the seacoast of New Hampshire.

I needed a project to test the new Think C version, since I have an interest in number theory, the project I picked was a Π calculation program. Most similar programs take over the machine while they’re at work. I wanted this one (just to keep life interesting) to optionally calculate Π in the background and let other applications continue. As the Mac doesn’t have preemptive multitasking, background applications only get time as permitted from the foreground application and vice versa. So, when we program time-intensive tasks, we’ve got to be nice to other software that may be running. The approach I took here is to let the user decide if this application should hog the system or not by setting an option in a menu. The default is to not hog.

To do this, I break the calculation up into a number of steps. The steps are implemented by a subclass of a TCL Chore object. This work is done at idle event time indirectly by the Dawdle() method of the bureaurocrats. The amount of CPU time that the calculation takes is the same if it’s done in the background or if the application hogs the system. However, the amount of time to the user is significantly longer if the work’s done in the background. This is due to the overhead of the time through the toolbox getting the idle event to run the calculation code. In addition, the calculation code is done in small steps and interspersed with other applications, etc. In the hog mode, it never leaves the calculation code at all but other apps don’t get CPU time.

Backgrounding is best for an application that’s got pre-calculation to do or can intersperse work with user input. As humans are so slow compared to the computer, the imbalance of speed can be put to good use this way.

Designing a chore for background processing

There are a couple of issues for chore-based background work. First, the chore has to be handled or know where to get all the information it needs for it’s processing. Second, if it’s work that gets completed that something else depends on, the chore has to be able to signal work complete. Last, chores need to be disposed of. One thing at a time.

In my case, my CPiChore class is a subclass of the CChore class. You’ll have to subclass the CChore class too as the standard class does nothing. The only method in CChore is the Perform() method (and of course those methods it inherits from CObject). My CPiChore subclass overrides the Perform() and Dispose() methods. It also has two new methods InitChore() and GetTicks() for setting the necessary data and getting information.

The InitChore() method takes all the necessary data as parameters. If I’d had a lot of data I could’ve passed a structure with the data loaded into it. However, as I only had four items, that seemed overkill for this project. Another alternative would’ve been for the chore to ask the necessary objects for the relevant data.

InitChore() takes the number of digits of Π to calculate, a pointer to the document object, and two flags to determine if the calculation should be timed and if hog mode is on. This data is stored locally. Note that the flags are stored in the object, but the document and the number of digits are stored in global variables that are defined and declared in the CPiChore.cp file. This makes them local to that file. The document object CAppPiDoc declares the CPiChore in it’s include file (CAppPiDoc.h). As a short-cut, I included the CAppPiDoc.h file in the CPiChore.cp file, and declared a global “document”. This means that I don’t have to type-cast as often when I use the “document” variable. The “Digits” variable is where the numDigits value is put. As this Π calculation code is ported, I didn’t convert all the calculation code to methods, instead, I left them and their global data in this CPiChore.cp file as straight C functions and data.

Perform() is the function in a chore that actually does the work. It’s called by the Idle() function in CApplication which is called every time an idle event happens. TCL considers an idle event to occur when we either receive a null event or a mouse-moved event. Note that if you’re doing background processing, you want null events while you’re app is in the background. To get them, the “Background null events” flag in the Size resource must be set. Use the SIZE flags popup in the Set Project Type dialog to do this.

The Idle() function also does some other things like check memory and do the Dawdle() method in the CBureaucrat object. Idle() runs through the chain of command to call all the Dawdle() methods there. It then runs through the list of chores to call Perform() on all the chores in the list.

When it gets to the Perform() method in CPiChore, CPiChore uses a stage counter to figure out what’s been done and what to do next. Basically, this determines and remembers what step it’s on. CPiChore uses the Machin algorithm to calculate Π (see the sidebar). The last step in the calculation stops the count of the time and tells the CAppPiDoc object that the chore is done by calling PiDone() in the document with the result of the calculation as it’s parameter.

The document gets the time the calculation took, sets flags for the state of the calculation and sends the result to the pane for display. The next time and idle event happens and the document’s Dawdle() method is called, the chore will be disposed of. This approach to deallocating the chore is known as lazy management of memory. It has the advantage that, as the Dawdle() method is called all the time, it’s more sure in deallocating the chore. The completion and the deallocation are in independent routines.

The pane to display this information isn’t a standard text edit pane. Text edit can only handle information up to around 30K of text characters. It’s easy to generate a value of Π that takes more space that this. I subclassed CPanorama to make a pane that just displayed a lot of text for this application. In addition, this CPiPane object also knows how to translate the digit numbers to text characters and how to display the time information too.

Implementation details

All the key user actions are handled through the DoCommand() method of the CAppPiDoc object. The one exception is the about box which is handled in the application’s DoCommand() method.

The number of digits is set in a dialog box that I use the regular TCL dialog classes to run. The DoCommand() method extracts the digit count from the dialog before it’s disposed. Although Symantec’s documentation states that the dialog objects must be subclassed to be really useful, I didn’t find the need to do that with the number of dialog view objects that know about the format and range of information. In this case, I used an integer text dialog view object that was told, via the template, what range of numbers to accept. Then, I use the FindViewByID() method to get that object and extract the now-conditioned value.

The time and hog (or “faster”) commands just toggle their values. In the UpdateMenus() method, these values are used to determine if those menu items are checked or not.

The start command is the important one. If the Π calculation chore is running, then this item is labeled “Stop”. If selected, it cancels and disposes of the chore. If the chore isn’t running, it starts it using the currently set number of digits and the timing and faster flags. It then assigns the chore to the application, this puts the chore in the chore list to be processed on an idle event. The chore will stay in the list until removed.

Once the chore is in the list, then the Perform() method is called with every idle event. If the faster or hog mode is on, then the Perform() method doesn’t return until the calculation is complete. Otherwise, it does one step and returns.

When the calculation is complete, the Perform() method calls the PiDone() method in the document which hands the result data off to the pane for display with the pane’s SetContent() method. This result data looks like the other data used in the calculation. It’s an array of character, but I use the characters as small integers. Actually, there’s a fair amount of wasted space as the integers range from 0-9. The calculation is carried out in decimal and there are custom routines to handle multiple place arithmetic.

The SetContent() method disposes any previous pointer to data (the data is held as pointers due to the ported code) and converts the data to characters for display. It does this by simply adding the integer value of the character zero to each number in the array.

When the Draw() method for the pane is called (it’s called automatically when the pane needs to be refreshed or when Refresh() is called), it sets the font and size of the text to draw. Draw() then creates a header line, draws it, skips down a little and draws the characters for Π. To do the drawing of the text, Draw() calls a routine called DrawWrapText() to wrap the text at the margins correctly.

DrawWrapText() takes parameters to the text, it’s size, and the margin, line information and whether or not to force a margin. This latter is similar to typing the return key to force a new line.

Listing:  CPiChore.cp
/*
 CPiChore
 Written and copyright by Malcolm H. Teas 1993.
 All rights reserved.
 
*/

#include "CAppPiDoc.h"

/*
 * Global variables
 */
int   Sign, 
 Zero, 
 Digits, 
 Pass,  /* two passes */
 Exp, /* exponent for divide */
 Divide;
static int constant[3]={0,25,239};

CAppPiDoc *document;
 
/*
 * These character arrays represent the digits for extended
 * precision arithmatic. (These may need to be int on 
 * some machines).
 */
#define DIGIT char
DIGIT *Power, *Term, *Result;

void copy (void);
void init(void);
void divi(register DIGIT *array);
void sub (void);
void add (void);

void CPiChore::Perform (long *maxSleep)
{
 long start;
 
 if (timeCalc)
 start = TickCount ();
 
 do  {
 switch (stage)  {
 case 0:
 init ();
 stage++;
 break;
 
 case 1:
 copy();
 Divide = Exp;
 divi(Term);
 if ( Sign > 0 ) 
 add();
 if ( Sign < 0 ) 
 sub();
 Exp = Exp + 2;
 Sign *= -1;
 Divide = constant[Pass];
 divi(Power);
 if (Pass == 2) 
 divi(Power);
 if (Zero == 0)
 stage++;
 break;
 
 case 2:
 Pass++;
 init ();
 stage++;
 break;
 
 case 3:
 copy();
 Divide = Exp;
 divi(Term);
 if ( Sign > 0 ) 
 add();
 if ( Sign < 0 ) 
 sub();
 Exp = Exp + 2;
 Sign *= -1;
 Divide = constant[Pass];
 divi(Power);
 if (Pass == 2) 
 divi(Power);
 if (Zero == 0)
 stage++;
 break;
 
 case 4:
 if (timeCalc)  {
 ticks += TickCount() - start;
 timeCalc = false;
 }
 document->PiDone (Result);
 stage++;
 break;
 
 case 5:
 break;
 }
 }  while (faster && stage != 5);
 
 if (timeCalc)
 ticks += TickCount() - start;
}

void CPiChore::InitChore (long numDigits, CDocument *doc, 
 Boolean timeIt, Boolean faster)
{
 stage = 0;
 Pass = 1;
 document = (CAppPiDoc *) doc;
 timeCalc = timeIt;
 ticks = 0;
 this->faster = faster;
 
 Digits = numDigits;
 FailNIL (Power = (char *) NewPtr (numDigits + 1));
 FailNIL (Term = (char *) NewPtr (numDigits + 1));
 FailNIL (Result = (char *) NewPtr (numDigits + 1));
}

long CPiChore::GetTicks (void)
{
 return ticks;
}

void CPiChore::Dispose (void)
{
 DisposPtr (Power);
 DisposPtr (Term);

 inherited::Dispose (); 
}

void add (void)
{
 register DIGIT *r, *t;
 register int sum, carry = 0;
    
 r = Result + Digits;
 t = Term + Digits;
 
 while( r >= Result )
 {
 sum = *r + *t + carry;
 carry = 0;
 if( sum >= 10 )
 {
 sum -= 10;
 carry++;
 }
 *r-- = sum;
 --t;
 }
}
 
void sub (void)
{
 register DIGIT *r, *t;
 register int diff, loan = 0; 
 
 r = Result + Digits;
 t = Term + Digits;
 
 while( r >= Result )
 {
 diff = *r - *t - loan;
 loan = 0;
 if( diff < 0 )
 {
 diff += 10;
 loan++;
 }
 *r-- = diff;
 --t;
 }
}
 
void divi(register DIGIT *array)
{
 register DIGIT *end;
 register int quotient, residue, digit = 0;
 
 Zero = 0;
 
 for( end = array + Digits; array <= end; )
 {
 digit    += *array;
 quotient =  digit / Divide;
 residue  =  digit % Divide;
 
 if((Zero !=0) || ((quotient+residue) != 0)) 
 Zero = 1;
 else 
 Zero = 0;
 
 *array++ = quotient;
 digit = 10 * residue;
 }
}
 
void init(void)
{
 register DIGIT *p, *t, *r, *end;
 
 p = Power;
 t = Term;
 r = Result;
 end = Power+Digits;
 
 while( p <= end )
 {
 *p++ = 0;
 *t++ = 0;
 
 if (Pass == 1) 
 *r++ = 0;
 }
 
 *Power = 16 / (Pass * Pass);
 
 if( Pass == 1 ) 
 {
 Divide = 5;
 }
 else
 {
 Divide = 239;
 }
 
 divi(Power);
 Exp = 1;
 Sign = 3 - (2 * Pass);
}
 
void copy (void)
{
 register DIGIT *t, *p, *end;
 
 t = Term;
 p = Power;
 end = Term + Digits;
 
 while (t <= end)
 {
 *t++ = *p++;
 }
}

How to calculate Π

The number Π in it’s simplest form, is the ratio of the circumference of a circle to that circle’s diameter. It’s an irrational number so it can’t be represented exactly by a ratio, or fraction, between two integers. As decimal numbering systems are essentially another way of expressing a fraction (i.e. the number 2.5 is actually the fraction 25/10), then Π, expressed as a decimal number, has an infinite number of places. A good reference for Π calculations is the book “History of Π” by Petr Beckmann.

The Π calculation algorithms work by interating over an approximation to Π. Each iteration, or time through the algorithm, the accuracy of the Π estimate increases. We say that we’re calculating Π, but we’re actually calculating an approximation to it. John Wallis, a sixteenth-century mathemetician, discovered a series of fractions, that when multiplied together result in Π. These are two ways to write Wallis’ product:

While this works fine, it’s also quite slow to converge or increase in accuracy, each fraction added doesn’t increase the accuracy very much. The Apple Π program in the article uses the Machin algorithm which converges much faster. It sums more complex fractions to speed up the process. More recently, mathemeticians using elliptical integrals discovered the holy grail of Π calculations: an algorithm, that while quite complex, doubles the number of accurate digits of Π for each iteration.

I ported the Machin algorithm from a version by Bryan Costales, who in turn ported it from a version by Robert Bishop. The algorithm is based on the equation:

[This is a Taylor expansion of 16 * Arctan(1/5) - 4 * Arctan(1/239). Tech. Ed.]

I use arrays of digits to implement the numbers and have written my own routines to do the arithmetic with these numbers. The calculation above is broken down into a number of steps. A step is done each time through the idle loop in the application.

Some have claimed that calculating Π to ever increasing numbers of digits (the record is around 40 million digits currently) is unecessary and not useful. However, when companies install a new supercomputer, this kind of calculation can be used for testing as it repeats intensive CPU work.

Number theorists use the results of long Π calculations to analyze the characteristics of the number, and through it, other similar numbers. Since number theory is a fundamental branch of mathematics that’s used in cryptography and computer theory this work may shed light on some aspects of mathematics that can be of practical use. Besides, it’s fun.

 

Community Search:
MacTech Search:

Software Updates via MacUpdate

ExpanDrive 6.1.8 - Access cloud storage...
ExpanDrive builds cloud storage in every application, acts just like a USB drive plugged into your Mac. With ExpanDrive, you can securely access any remote file server directly from the Finder or... Read more
DiskCatalogMaker 7.2.7 - 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
DiskCatalogMaker 7.2.7 - 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
Iridient Developer 3.2.1 - Powerful imag...
Iridient Developer (was RAW Developer) is a powerful image-conversion application designed specifically for OS X. Iridient Developer gives advanced photographers total control over every aspect of... Read more
BusyContacts 1.2.7 - Fast, efficient con...
BusyContacts is a contact manager for OS X that makes creating, finding, and managing contacts faster and more efficient. It brings to contact management the same power, flexibility, and sharing... Read more
MegaSeg 6.0.5 - Professional DJ and radi...
MegaSeg is a complete solution for pro audio/video DJ mixing, radio automation, and music scheduling with rock-solid performance and an easy-to-use design. Mix with visual waveforms and Magic... Read more
BusyContacts 1.2.7 - Fast, efficient con...
BusyContacts is a contact manager for OS X that makes creating, finding, and managing contacts faster and more efficient. It brings to contact management the same power, flexibility, and sharing... Read more
MegaSeg 6.0.5 - Professional DJ and radi...
MegaSeg is a complete solution for pro audio/video DJ mixing, radio automation, and music scheduling with rock-solid performance and an easy-to-use design. Mix with visual waveforms and Magic... Read more
Iridient Developer 3.2.1 - Powerful imag...
Iridient Developer (was RAW Developer) is a powerful image-conversion application designed specifically for OS X. Iridient Developer gives advanced photographers total control over every aspect of... Read more
iFFmpeg 6.6.1 - Convert multimedia files...
iFFmpeg is a comprehensive media tool to convert movie, audio and media files between formats. The FFmpeg command line instructions can be very hard to master/understand, so iFFmpeg does all the hard... Read more

Latest Forum Discussions

See All

The best deals on the App Store this wee...
A new week means new discounts on the App Store. This week's deals run the gamut of action-adventure titles, puzzle games, and one of the best narrative adventure series out there. If you're looking to fill out your mobile gaming library on a... | Read more »
What you need to know about Animal Cross...
We hope you've been hard at work on collecting all of those holiday items in Animal Crossing: Pocket Camp, because you're about to get a whole new list of fun things to do as the game receives its first big update sometime soon. There are a lot of... | Read more »
Reigns: Her Majesty guide - how to use e...
Ruling a kingdom isn't easy--doubly so for a queen whose every decision is questioned by the other factions seeking a slice of power. Reigns: Her Majesty builds on the original game's swipey tactics, adding items that you can use to move the story... | Read more »
The best new games we played this week -...
Friday has crept up on us once again, so it's time to honor the best new games we've played over the past few days. This past week was a pretty exciting one, with the debut of lots of beautiful new indies and some familiar faces returning to the... | Read more »
Portal Knights guide- beginner tips and...
Portal Knights is finally making the jump to iOS and Android, and it's already climbing the ranks to become the next big MMO experience on mobile. This sprawling sandbox game will let you pursue any adventure you wish, whether you want to sling... | Read more »
Reigns: Her Majesty guide - how to swipe...
Reigns: Her Majesty is storming the App Store this week, bringing more tinder-esque kingdom building to eager players everywhere. If you've played the original Reigns, you'll know that leading a kingdom is never easy. It's a careful balancing act... | Read more »
Getting Over It (Games)
Getting Over It 1.0 Device: iOS Universal Category: Games Price: $4.99, Version: 1.0 (iTunes) Description: A game I madeFor a certain kind of person To hurt them. • Climb up an enormous mountain with nothing but a hammer and a pot.•... | Read more »
Reigns: Her Majesty (Games)
Reigns: Her Majesty 1.0 Device: iOS Universal Category: Games Price: $2.99, Version: 1.0 (iTunes) Description: | Read more »
Pocket Legends Adventures guide - how to...
Pocket Legends Adventures is a fun action adventure RPG that takes control when you want it to, but also opens itself for player input, too, if you're looking to tkae a more active role in combat. Regardless of play style, the game can be quite... | Read more »
Portal Knights (Games)
Portal Knights 1.2.4 Device: iOS Universal Category: Games Price: $4.99, Version: 1.2.4 (iTunes) Description: Craft your adventure. Forge your hero. Become the ultimate Portal Knight! | Read more »

Price Scanner via MacPrices.net

Green Monday deal: 15″ 2.8GHz MacBook Pro on...
B&H Photo has the 15″ 2.8GHz Space Gray MacBook Pro on sale for $250 off MSRP for today only as part of their Green Monday/Holiday sale. Shipping is free, and B&H charges sales tax for NY... Read more
Green Monday sale: B&H offers 12″ Apple i...
B&H Photo has 12″ iPad Pros on sale for up to $150 off MSRP as part of their Green Monday/Holiday sale. Shipping is free, and B&H charges sales tax in NY & NJ only: – 12″ 64GB WiFi iPad... Read more
Holiday deal: 21″ and 27″ Apple iMacs on sale...
MacMall has 2017 21″ and 27″ Apple iMacs on sale for up to $200 off MSRP. Shipping is free: – 21″ 2.3GHz iMac: $999 $100 off MSRP – 21″ 3.0GHz iMac: $1199 $100 off MSRP – 21″ 3.4GHz iMac: $1379 $120... Read more
Holiday deal: Apple Mac minis for up to $150...
MacMall has Mac minis on sale for up to $100 off MSRP, each including free shipping: – 1.4GHz Mac mini: $399 $100 off MSRP – 2.6GHz Mac mini: $599 $100 off MSRP – 2.8GHz Mac mini: $949 $50 off MSRP... Read more
Beats by Dr. Dre – BeatsX Earphones on sale f...
Best Buy has BeatsX Earphones on sale for $109, $40 off, on their online store. Sale price for online orders only. Choose free store pickup, if available, or choose free shipping. Read more
10″ 64GB WiFi Apple iPad Pros on sale for $59...
MacMall has 10.5″ 64GB Apple iPad Pros on sale for $599 including free shipping. That’s $50 off MSRP and among the lowest prices available for these iPads from any Apple reseller. Read more
15″ 2.2GHz MacBook Pros on sale for $200-$300...
B&H Photo has the 15″ 2.2GHz MacBook Pro available for $200 off MSRP including free shipping plus NY & NJ sales tax only: – 15″ 2.2GHz MacBook Pro (MJLQ2LL/A): $1799 $200 off MSRP Apple has... Read more
Holiday sale: 15″ MacBook Pros for $200-$420...
MacMall has 15″ MacBook Pros on sale for $220-$300 off MSRP, each including free shipping: – 15″ 2.8GHz MacBook Pro Space Gray (MPTR2LL/A): $2179, $220 off MSRP – 15″ 2.8GHz MacBook Pro Silver (... Read more
Holiday sale: 13″ MacBook Airs for up to $150...
B&H Photo has 13″ MacBook Airs on sale for $100 off MSRP as part of their Holiday sale. Shipping is free, and B&H charges sales tax for NY & NJ residents only: – 13″ 1.8GHz/128GB MacBook... Read more
The best Holiday sale prices on 13″ MacBook P...
B&H Photo has 13″ MacBook Pros on sale this weekend, with models available for $100-$150 off MSRP. Shipping is free, and B&H charges sales tax for NY & NJ residents only: – 13-inch 2.3GHz... Read more

Jobs Board

*Apple* Retail - Multiple Positions - Apple,...
Job Description:SalesSpecialist - Retail Customer Service and SalesTransform Apple Store visitors into loyal Apple customers. When customers enter the store, Read more
*Apple* Retail - Multiple Positions - Apple,...
Job Description: Sales Specialist - Retail Customer Service and Sales Transform Apple Store visitors into loyal Apple customers. When customers enter the store, Read more
*Apple* Retail - Multiple Positions - Apple,...
Job Description:SalesSpecialist - Retail Customer Service and SalesTransform Apple Store visitors into loyal Apple customers. When customers enter the store, Read more
*Apple* Information Security - Security Data...
# Apple Information Security - Security Data Analyst Job Number: 113119545 Austin, Texas, United States Posted: 10-Nov-2017 Weekly Hours: 40.00 **Job Summary** This Read more
*Apple* Retail - Multiple Positions - Apple,...
Job Description: Sales Specialist - Retail Customer Service and Sales Transform Apple Store visitors into loyal Apple customers. When customers enter the store, Read more
All contents are Copyright 1984-2011 by Xplain Corporation. All rights reserved. Theme designed by Icreon.