Newsgroups: sci.fractals,bit.listserv.frac-l
Path: unixg.ubc.ca!nntp.cs.ubc.ca!utcsri!utnut!torn!news2.uunet.ca!math.ohio-state.edu!usc!news.service.uci.edu!unogate!mvb.saic.com!ast.saic.com!jupiter!hilljr
From: hilljr@jupiter.saic.com (Jay R. Hill)
Subject: Re: Have we seen it all?
Message-ID: <1993Oct18.202157.2691@ast.saic.com>
Followup-To: sci.fractals
Sender: news@ast.saic.com
Organization: SAIC
References: <1993Oct14.170930.2402@ast.saic.com> <29mbqa$5t4@crcnis1.unl.edu>
Date: Mon, 18 Oct 1993 20:21:57 GMT
Lines: 47
Xref: unixg.ubc.ca sci.fractals:2204 bit.listserv.frac-l:2114
 
In article <29mbqa$5t4@crcnis1.unl.edu>, jepler@herbie.unl.edu (Jeff Epler) writes:
|> hilljr@jupiter.saic.com (Jay R. Hill) writes:
|> 
|> >>Distance-Estimator method, and setting it to view the area of
|> >>the complex plane centered at
|> 
|> >  -1.769 110 375 463 767 386 0  +  0.009 020 388 228 023 439 7i
|> 
|> >>with an edge length (height and width) of 6.4e-17 .  Make sure
|> >>the number of iterations is at least 10,000.
|> 
|> Anybody get this to go with Xfractint?  I'm using it on a Linux 386
|> machine, and it always rudely dumps me back to the shell prompt after
|> I enter those center coords and a magnification of 1.5625e16.  (I've
|> had other problems with xfractint, my binary may just have problems.)
 
To get this view you need 80 bit floats such as that provided by INTEL.
I rolled my own code in C++ (also you can use Turbo Pascal) using
long double.  If you are stuck with double, you only get 54 bits mantissa
or 10^-16 resolution.  Robert P Munafo <mrob@WORLD.STD.COM> did mention
this is his original post.  Unfortunately, this leaves you about 10^4
short of zomming into his 'view'.  No matter how many digits we have there
is always still a deeper zoom  ;-)  
 
For those on sci.fractals, to which I'm cross posting this, Robert did
provide an interesting view to try with normal double-precision:
 
   -1.749 925 438  +  0.000 000 000 i, size 2.98e-7 and
   -1.749 925 361 593 013 94  +  0.000 000 004 193 282 92 i, size 1.11e-15
 
These are near the tree strucrutes I mentioned in 
<1993Oct14.170930.2402@ast.saic.com>.  The tree structures I noted are
the 'large scale' filements along which are very small island msets.  It
is near one of these into which Robert suggests zooming.
 
|> --
|> Jeff Epler jepler@herbie.unl.edu (Preferred) or jepler@nyx.cs.du.edu
|> ____ "Nuke the unborn gay whales" -- Never seen on a protest sign
|> \bi/ |umop apisdn| First year comp sci major.  CRPG addict. 
|>  \/  1.5<kinsey<2.5 IRC Synger  Running Linux 0.99.13. DOS is evil
 
Jay
-- 
{ hilljr@jupiter.saic.com } begin writeln(3*ln(640320)/sqrt(163):17:15) end.
void main(){double sqrt(), y=1/sqrt(2.), a=.5, m=1,z; int n=0;
for(;m*=2,z=sqrt(1-y*y),y=(1-z)/(1+z),a=a*(1+y)*(1+y)-m*y,n<4;n++);
printf("%17.15lf\n",1/a);}