home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
ftp.xmission.com
/
2014.06.ftp.xmission.com.tar
/
ftp.xmission.com
/
pub
/
lists
/
fractint
/
archive
/
v01.n440
< prev
next >
Wrap
Internet Message Format
|
2000-01-20
|
56KB
From: owner-fractint-digest@lists.xmission.com (fractint-digest)
To: fractint-digest@lists.xmission.com
Subject: fractint-digest V1 #440
Reply-To: fractint-digest
Sender: owner-fractint-digest@lists.xmission.com
Errors-To: owner-fractint-digest@lists.xmission.com
Precedence: bulk
fractint-digest Thursday, January 20 2000 Volume 01 : Number 440
----------------------------------------------------------------------
Date: Mon, 17 Jan 2000 09:56:52 -0700
From: Phil McRevis <legalize@xmission.com>
Subject: Re: (fractint) Determining the M-Set
First, please turn off your HTML mail setting. The fractint mailing
list specifically requests that you post in plain text and not use
attachments.
In article <005101bf605f$e56319e0$b0fa343e@thomaspe>,
lochfrass@friendfactory.com writes:
> I know that FRACTINT stops the calculation of an orbit if |z|>2. I heard =
> that it is certain that a point is outside the M-Set if this happens. =
> Can anyone explain this to me?
The M set is defined as the set of points for which the iterated
sequence remains bounded. M is computed from the sequence
znew = zold^2 + c
Once the magnitude of zold becomes greater than 2, this sequence will
diverge to infinity.
- --
<http://www.xmission.com/~legalize/> Legalize Adulthood!
``Ain't it funny that they all fire the pistol,
at the wrong end of the race?''--PDBT
legalize@xmission.com <http://www.thewho.net>
- --------------------------------------------------------------
Thanks for using Fractint, The Fractals and Fractint Discussion List
Post Message: fractint@lists.xmission.com
Get Commands: majordomo@lists.xmission.com "help"
Administrator: twegner@swbell.net
Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
------------------------------
Date: Tue, 18 Jan 2000 01:09:47 -0500 (EST)
From: Jim Muth <jamth@mindspring.com>
Subject: (fractint) FOTD, 18-01-00 (Bits and Pieces) (c)
FOTD -- January 18, 2000
Fractal enthusiasts and visionaries:
I named today's fractal "Bits and Pieces". I gave it this name
because of the disconnected bits and pieces of fractal stuff
scattered throughout the area. These broken-off pieces are
there because the parent fractal is not totally connected.
The Mandelbrot set is connected. That is it has no parts that
are disconnected from its main body. Every midget is connected
to the main bay or one of the buds by an infinitely thin
filament of trapped points. But not all fractals are connected
in such a manner.
The parent fractal of today's image is not connected. It
contains islands of fractal material cut off from the main body.
In fact, it has no main body, appearing instead as several large
areas of fractal chaos filled with scattered small quadratic
M-sets of various distorted shapes and sizes. Around the larger
chaotic areas lie smaller disconnected areas of chaos filled
with tiny M-sets. Since these midget M-sets lie in cut-off
patches of chaos, the patterns around the midgets are also
filled with cut-off bits and pieces of fractal stuff. Today's
picture shows one of these midgets.
I have colored the little fellow to emphasize the disconnected
nature of the features. The coloring alone determines the
pattern, as can be shown by rotating the colors a few steps in
either direction. I rate this image about a six on my personal
0-to-10 quality scale.
The parameter file renders in 2 minutes on an average Pentium.
The JPEG image file is also available from:
<alt.binaries.pictures.fractals>
and from:
<http://home.att.net/~Paul.N.Lee/FotD/FotD.html>
The fractal weather today was sunny but hugely cold -- perfect
weather for hunting fractals. The temperature of 17F (-8C) was
far far too cold for the cats to step outdoors. They passed the
day chasing the moving patches of sunlight or stretched out
before a radiator.
It was also a perfect day for organizing my pompous philosophi-
cal ponderings, which are coming together under the theme of the
three ways of obtaining knowledge. In a few days the ponderings
will be made public -- stay posted.
Oh my -- I see it's time to shut down the fractal shop, feed the
fractal cats, and call it a night. Until next time, take care,
and if your fractal becomes disconnected, it might not be all
that bad.
Jim Muth
jamth@mindspring.com
START FORMULA==============================================
MandelbrotMix4 {; Jim Muth
a=real(p1), b=imag(p1), d=real(p2), f=imag(p2),
g=1/f, h=1/d, j=1/(f-b), z=(-a*b*g*h)^j,
k=real(p3)+1, l=imag(p3)+100, c=fn1(pixel):
z=k*((a*(z^b))+(d*(z^f)))+c,
|z| < l
}
END FORMULA================================================
START PARAMETER FILE=======================================
Bits_and_Pieces { ; time=0:02:10.72, SF5 on a p200
reset=2000 type=formula formulafile=critical.frm
formulaname=MandelbrotMix4 function=ident passes=1
center-mag=-3.74526689211778600/-3.02147577593595700\
/6.83767e+007/1/57.5 params=0.333/2/-2/-2/0/0 float=y
maxiter=500 bailout=25 inside=0
logmap=27 symmetry=none periodicity=10
colors=000WGJ<3>CD9<3>6SZ5Kd4Pj<2>Tcpahrkovrvzzzz000\
mml<3>QQQKKKDDD666WRBGDc5Uz4cz37b<3>017jSy<4>UJcRH_N\
FWKESHCNDAJ<3>436213zaOwVKmPGcJCZC8764gYO<3>EB8754C9\
1<3>53032A11KAmP6jS<3>3Qc2Kf2Fi1Ab05e`qh<2>9DqmWtjTw\
VJzF9zg1z<3>Q0zM0zH0z<2>40z6Qz4Hz28z9iz6Uz3Fz7kz<7>3\
Mz2Jz2Gz<3>03zuPz<3>TCzL9zE6z73ziVz<3>FAz75zJoz<3>Da\
zCYzBVzAPz8Lz<3>3Pz2Qz1RzqSz<3>RWzKXzDYz6Zzv_z<4>9dz\
SezLfzEgz7hz2iz<2>0lz7mz<5>3sz3tz2uz<3>0yzZzn<5>IzRG\
zNDzJ<3>2z3ZzY<3>LzKHzHEzD<2>3z3fzX<2>Az8Mzc7zF<2>1z\
3czD<3>Hz5Bz35z1zzC<3>Zz6Sz5Lz4Ez27z1CznKzu<3>6zJ3z9\
Lzk<5>BzP9zM8zI<3>1z36zG5zE4zC
}
END PARAMETER FILE=========================================
START 20.0 PAR-FORMULA FILE================================
Bits_and_Pieces { ; time=0:02:10.72, SF5 on a p200
reset=2000 type=formula formulafile=critical.frm
formulaname=MandelbrotMix4 function=ident passes=1
center-mag=-3.74526689211778600/-3.02147577593595700\
/6.83767e+007/1/57.5 params=0.333/2/-2/-2/0/0 float=y
maxiter=500 bailout=25 inside=0
logmap=27 symmetry=none periodicity=10
colors=000WGJ<3>CD9<3>6SZ5Kd4Pj<2>Tcpahrkovrvzzzz000\
mml<3>QQQKKKDDD666WRBGDc5Uz4cz37b<3>017jSy<4>UJcRH_N\
FWKESHCNDAJ<3>436213zaOwVKmPGcJCZC8764gYO<3>EB8754C9\
1<3>53032A11KAmP6jS<3>3Qc2Kf2Fi1Ab05e`qh<2>9DqmWtjTw\
VJzF9zg1z<3>Q0zM0zH0z<2>40z6Qz4Hz28z9iz6Uz3Fz7kz<7>3\
Mz2Jz2Gz<3>03zuPz<3>TCzL9zE6z73ziVz<3>FAz75zJoz<3>Da\
zCYzBVzAPz8Lz<3>3Pz2Qz1RzqSz<3>RWzKXzDYz6Zzv_z<4>9dz\
SezLfzEgz7hz2iz<2>0lz7mz<5>3sz3tz2uz<3>0yzZzn<5>IzRG\
zNDzJ<3>2z3ZzY<3>LzKHzHEzD<2>3z3fzX<2>Az8Mzc7zF<2>1z\
3czD<3>Hz5Bz35z1zzC<3>Zz6Sz5Lz4Ez27z1CznKzu<3>6zJ3z9\
Lzk<5>BzP9zM8zI<3>1z36zG5zE4zC
}
frm:MandelbrotMix4 {; Jim Muth
a=real(p1), b=imag(p1), d=real(p2), f=imag(p2),
g=1/f, h=1/d, j=1/(f-b), z=(-a*b*g*h)^j,
k=real(p3)+1, l=imag(p3)+100, c=fn1(pixel):
z=k*((a*(z^b))+(d*(z^f)))+c,
|z| < l
}
END 20.0 PAR-FORMULA FILE==================================
- --------------------------------------------------------------
Thanks for using Fractint, The Fractals and Fractint Discussion List
Post Message: fractint@lists.xmission.com
Get Commands: majordomo@lists.xmission.com "help"
Administrator: twegner@swbell.net
Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
------------------------------
Date: Tue, 18 Jan 2000 08:42:31 GMT
From: "Rupert Millard" <rupertam@hotmail.com>
Subject: Re: (fractint) Colour Map Recognition, I've writen a program
Hello all,
You can download my compressed duplicate colour map program at:
http://www.geocities.com/kangarupert/page10.html
Unfortunatley it says some maps are the same when they are not, check them
carefully in Fractint before you delete them. I will improve the program
when I have time so that double-checking will not be neccesary. I have
written this with best intentoins at heart and I accept no responsibility if
something goes wrong.
From,
Rupert
______________________________________________________
Get Your Private, Free Email at http://www.hotmail.com
- --------------------------------------------------------------
Thanks for using Fractint, The Fractals and Fractint Discussion List
Post Message: fractint@lists.xmission.com
Get Commands: majordomo@lists.xmission.com "help"
Administrator: twegner@swbell.net
Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
------------------------------
Date: Wed, 19 Jan 2000 01:31:45 -0500 (EST)
From: Jim Muth <jamth@mindspring.com>
Subject: (fractint) FOTD, 19-01-00 (Schmidget Midget) (c)
FOTD -- January 19, 2000
Fractal enthusiasts and visionaries:
When I saw today's midget, I said to myself, "midget schmidget,
that's no midget. It's just an open area that will fill in when
a higher maxiter is used." Almost as a joke, I raised the
maxiter to 500,000 and recalculated the picture. But the midget
did not vanish. That is when I decided on the name: "Schmidget
Midget".
The midget is a true midget of order 1.2 or something like that.
Actually, the formula is -(Z^1.2)-1.23*(Z^1.234)+C, an obviously
whimsical arrangement. As is the case with many midgets of such
a low order, the midget in today's picture is surrounded by
strangely decorative fractal objects, all divided and sub-div-
ided to infinity. There is a smaller sub-midget near the south
shore of today's midget, which I will investigate first thing
tomorrow.
The fractal weather was cloudy and very cold again today. A few
flakes of snow drifted down during the afternoon, but did not
accumulate on the ground. The temperature of 20F (-6.5C) was
perfect for snowmaking on the area ski-resorts.
Right now, it's getting late and the cats are getting hungry.
It's time to shut down the Fractal shoppe until tomorrow, and
call it a night. Until next time, take care, and the parameter
file is slow. Don't forget to download the image from:
<alt.binaries.pictures.fractals>
or from:
<http://home.att.net/~Paul.N.Lee/FotD/FotD.html>
My philosophy foundered once again today. As always, I'll try
again tomorrow.
Jim Muth
jamth@mindspring.com
START FORMULA==============================================
MandelbrotMix4 {; Jim Muth
a=real(p1), b=imag(p1), d=real(p2), f=imag(p2),
g=1/f, h=1/d, j=1/(f-b), z=(-a*b*g*h)^j,
k=real(p3)+1, l=imag(p3)+100, c=fn1(pixel):
z=k*((a*(z^b))+(d*(z^f)))+c,
|z| < l
}
END FORMULA================================================
START PARAMETER FILE=======================================
Schmidget_Midget { ; time=0:28:29.00, SF5 on p200
reset=2000 type=formula formulafile=critical.frm
formulaname=MandelbrotMix4 function=ident
center-mag=+0.05528421807509840/+0.21748697569581040\
/259.4338/1/32.5 params=-1/1.2/-1.23/1.234/-0.4/0
float=y maxiter=500000 bailout=25 inside=0 logmap=49
symmetry=none periodicity=10
colors=000fABlB8<17>Z7CY6DY6DX6DW6DV6DV5FU5HR6JP9LND\
NLFP<2>0gu<24>1JI1IG1HF<3>1E9<3>CQ9ES9HV9JYAM`FObK<3\
>igcnhhsjmrjrrjwrjzqjw<3>pejpcgpbdp`a<2>oWaoS`<2>nL`\
<3>mU_mW_mYZm_ZmaZ<5>lhYliXljX<4>kmWknWknV<3>kpVlqWn\
qWorXqsYruY<3>pyZpzZpzZ<17>kzajzajza<2>izaizagz`<3>l\
zdmzenzf<2>qzhozi<4>nzinzjnzjmzjmzjmzj<12>hzlhzmgzm<\
2>fzmfzmgzn<6>`zg_zfZze<3>VzaRzb<3>_zXazWczVezTgzSiz\
RjzR<13>lzElzDlzD<3>lz91zY<2>0z_<5>`zF
}
END PARAMETER FILE=========================================
START 20.0 PAR-FORMULA FILE================================
Schmidget_Midget { ; time=0:28:29.00, SF5 on p200
reset=2000 type=formula formulafile=critical.frm
formulaname=MandelbrotMix4 function=ident
center-mag=+0.05528421807509840/+0.21748697569581040\
/259.4338/1/32.5 params=-1/1.2/-1.23/1.234/-0.4/0
float=y maxiter=500000 bailout=25 inside=0 logmap=49
symmetry=none periodicity=10
colors=000fABlB8<17>Z7CY6DY6DX6DW6DV6DV5FU5HR6JP9LND\
NLFP<2>0gu<24>1JI1IG1HF<3>1E9<3>CQ9ES9HV9JYAM`FObK<3\
>igcnhhsjmrjrrjwrjzqjw<3>pejpcgpbdp`a<2>oWaoS`<2>nL`\
<3>mU_mW_mYZm_ZmaZ<5>lhYliXljX<4>kmWknWknV<3>kpVlqWn\
qWorXqsYruY<3>pyZpzZpzZ<17>kzajzajza<2>izaizagz`<3>l\
zdmzenzf<2>qzhozi<4>nzinzjnzjmzjmzjmzj<12>hzlhzmgzm<\
2>fzmfzmgzn<6>`zg_zfZze<3>VzaRzb<3>_zXazWczVezTgzSiz\
RjzR<13>lzElzDlzD<3>lz91zY<2>0z_<5>`zF
}
frm:MandelbrotMix4 {; Jim Muth
a=real(p1), b=imag(p1), d=real(p2), f=imag(p2),
g=1/f, h=1/d, j=1/(f-b), z=(-a*b*g*h)^j,
k=real(p3)+1, l=imag(p3)+100, c=fn1(pixel):
z=k*((a*(z^b))+(d*(z^f)))+c,
|z| < l
}
END 20.0 PAR-FORMULA FILE==================================
- --------------------------------------------------------------
Thanks for using Fractint, The Fractals and Fractint Discussion List
Post Message: fractint@lists.xmission.com
Get Commands: majordomo@lists.xmission.com "help"
Administrator: twegner@swbell.net
Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
------------------------------
Date: Wed, 19 Jan 2000 15:00:56 GET
From: "Tony \(Anthony\) Hanmer" <a_hanmer@hotmail.com>
Subject: Re: (fractint) Colour Map Recognition, I've writen a program
Many thanks, Rupert, for developing this programme! Most gratifying. I'm
currently rather a pauper as far as free computer access goes, but in a
while I hope for this to change, and to be able to investigate what you've
produced. Rest assured that it will be very useful.
Thanks again,
Tony Hanmer
______________________________________________________
Get Your Private, Free Email at http://www.hotmail.com
- --------------------------------------------------------------
Thanks for using Fractint, The Fractals and Fractint Discussion List
Post Message: fractint@lists.xmission.com
Get Commands: majordomo@lists.xmission.com "help"
Administrator: twegner@swbell.net
Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
------------------------------
Date: Wed, 19 Jan 2000 11:41:51 -0000
From: Stephanie Dunne <sdunne@cs.may.ie>
Subject: RE: (fractint) Colour Map Recognition, I've writen a program
Hi,
Can I get a copy of that colour map recognition program please?
Stephanie Dunne
- --------------------------------------------------------------
Thanks for using Fractint, The Fractals and Fractint Discussion List
Post Message: fractint@lists.xmission.com
Get Commands: majordomo@lists.xmission.com "help"
Administrator: twegner@swbell.net
Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
------------------------------
Date: Wed, 19 Jan 2000 17:19:07 +-100
From: Severino Fernandez <severinofer@recol.es>
Subject: (fractint) change detection using multifractal analysis
I am very interested in automated change detection on images, =
particularly in small area changes
in remote sensing images, such as caused by changes in structures, new =
buildings, new roads, etc.
I have read the paper on change detection in images using multifractal =
analysis written by Cristophe Canus and Jacques Levy Vehel and their =
approach seemed extremely innovative and with a great potential.
I tried to implement it using the large deviation spectrum, estimating =
it with the simple histogram technique that Dr. Levy Vehel mentions in =
some other of his papers, and using dyadic partitions of my sample =
images, but the spectra I get are very discountinuous, due to the =
non-overlapping of the capacity measures computed at the different =
partitions. In fact, they do not resemble at all the spectra that Dr =
Canus and Dr Levy Vehel present in their paper.
I assume that I am wrongly interpreting the technique or that the method =
for spectrum estimation I am using is not the appropiate one.
If somebody has made any experience on this particular subject, I think =
I need advice or any hint on multifractal spectrum estimation algorithms =
publications, web pages or free software.
Thank you all in advance
Severino Fernandez
Instituto Nacional de Tecnica Aeroespacial
Division de Ciencias del Espacio
Departamento de Teledeteccion
Carretera de Ajalvir, Km 4
28850 Torrejon de Ardoz
Spain
Tel +34 91 677 41 30
+34 91 677 41 90
+34 91 305 16 52
Fax +34 91 677 46 46
Email severinofer@recol.es
- --------------------------------------------------------------
Thanks for using Fractint, The Fractals and Fractint Discussion List
Post Message: fractint@lists.xmission.com
Get Commands: majordomo@lists.xmission.com "help"
Administrator: twegner@swbell.net
Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
------------------------------
Date: Wed, 19 Jan 2000 18:38:02 EST
From: RParracho@aol.com
Subject: Re: (fractint) Colour Map Recognition, dos help
i was wondering if anyone remembers the dos pipe or redirect so the output of
this great little utility can go to a text file or diectly to the printer?
- --------------------------------------------------------------
Thanks for using Fractint, The Fractals and Fractint Discussion List
Post Message: fractint@lists.xmission.com
Get Commands: majordomo@lists.xmission.com "help"
Administrator: twegner@swbell.net
Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
------------------------------
Date: Thu, 20 Jan 2000 0:31 0000
From: comdotatdotcom@csi.com
Subject: RE: Re: (fractint) Colour Map Recognition, dos help
Hi,
>i was wondering if anyone remembers the dos pipe or redirect so the
>output of
>this great little utility can go to a text file or diectly to the printer?
That'll be the > sign like this:
widget > out.txt
where widget is the name of the utility, to print it try:
widget > prn:
or:
widget > lpt1:
Cheers,
Robin.
- --------------------------------------------------------------
Thanks for using Fractint, The Fractals and Fractint Discussion List
Post Message: fractint@lists.xmission.com
Get Commands: majordomo@lists.xmission.com "help"
Administrator: twegner@swbell.net
Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
------------------------------
Date: Wed, 19 Jan 2000 21:26:45 EST
From: RParracho@aol.com
Subject: Re: (fractint) Colour Map Recognition, dos help
i tried the ">" redirect and it still went to console...am i just doing this
wrong?
c:\fractint\map>multimap > dupes.txt
it didn't redirect to my printer either. can you or anyone else do it
successfully?
- --------------------------------------------------------------
Thanks for using Fractint, The Fractals and Fractint Discussion List
Post Message: fractint@lists.xmission.com
Get Commands: majordomo@lists.xmission.com "help"
Administrator: twegner@swbell.net
Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
------------------------------
Date: Wed, 19 Jan 2000 22:22:19 -0500
From: Harry Bissell <harrybissell@prodigy.net>
Subject: Re: (fractint) Colour Map Recognition, dos help
I think you have to use
copy C:\xxxxxxxx >prn or lpt1 etc...
:^) Harry
RParracho@aol.com wrote:
> i tried the ">" redirect and it still went to console...am i just doing this
> wrong?
>
> c:\fractint\map>multimap > dupes.txt
>
> it didn't redirect to my printer either. can you or anyone else do it
> successfully?
>
> --------------------------------------------------------------
> Thanks for using Fractint, The Fractals and Fractint Discussion List
> Post Message: fractint@lists.xmission.com
> Get Commands: majordomo@lists.xmission.com "help"
> Administrator: twegner@swbell.net
> Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
- --------------------------------------------------------------
Thanks for using Fractint, The Fractals and Fractint Discussion List
Post Message: fractint@lists.xmission.com
Get Commands: majordomo@lists.xmission.com "help"
Administrator: twegner@swbell.net
Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
------------------------------
Date: Thu, 20 Jan 2000 00:39:00 -0500 (EST)
From: Jim Muth <jamth@mindspring.com>
Subject: (fractint) FOTD, 20-01-00 (Cyclogenesis) (c)
FOTD -- January 20, 2000
Fractal enthusiasts and visionaries:
Here in Fractal land everyone is preparing for the big one --
the blizzard of 2000, which is due to strike Thursday and drop
as much as 6 inches on the area. Of course in these parts an
inch of snow is almost a blizzard, so 6 inches must be
considered a storm of biblical proportions.
The incredible thing is that the storm which is due to drop all
that snow does not yet exist. It is expected to develop
Thursday morning just off the coast. Since the development of a
cyclonic storm is known as Cyclogenesis, and today's fractal is
twisted into spiral bands, I named the picture "Cyclogenesis".
Yes, I realize that the central midget is turning in the wrong
direction for a cyclone in the Northern Hemisphere, but maybe
the scene is in the Southern Hemisphere.
The formula behind the cyclone is my MandelbrotMix4, with some
significant irrational numbers randomly fed into it. And yes, I
am guilty of a bit of post-processing in a graphic program.
The parameter file, if you choose to run it, is a slow one,
requiring 41 minutes on a Pentium 200mhz. The JPEG'd image file
can be downloaded in only a minute or two from:
<alt.binaries.pictures.fractals>
or from:
<http://home.att.net/~Paul.N.Lee/FotD/FotD.html>
The fractal weather today was bearable, with partly cloudy skies
and an afternoon temperature of 32F (0C), which was milder than
yesterday, but still too cold for fractal cats to venture out of
doors.
A few days ago I left Mr. Percy Smedley standing in the middle
of his four-dimensional room contemplating the job before him of
painting the six walls. Lest there be fear that I have
abandoned Mr. Smedley, let me assure you that I have not
forgotten him. I have merely been too busy to properly describe
his task as he covers, or rather fills, the six surface-volumes
of the six walls with paint. (Actually, I'm not sure how I
manage to fit all the things I do into a day with only 24 hours.)
I hope to return to Mr. Smedley tomorrow. And if the snow
develops as predicted, things will grind to a halt, giving me
plenty of time to clear the driveway and return to our intrepid
painter.
Until that glorious moment arrives, take care, and the
philosophy *will* return, but Smedley comes first.
Jim Muth
jamth@mindspring.com
START FORMULA==============================================
MandelbrotMix4 {; Jim Muth
a=real(p1), b=imag(p1), d=real(p2), f=imag(p2),
g=1/f, h=1/d, j=1/(f-b), z=(-a*b*g*h)^j,
k=real(p3)+1, l=imag(p3)+100, c=fn1(pixel):
z=k*((a*(z^b))+(d*(z^f)))+c,
|z| < l
}
END FORMULA================================================
START PARAMETER FILE=======================================
Cyclogenesis { ; time=0:41:22.64, SF5 on a p200
reset=2000 type=formula formulafile=critical.frm
formulaname=MandelbrotMix4 function=ident passes=1
center-mag=-0.23418923607130600/+0.42864036781657180\
/1.288888e+011/1/90
params=2.71828182845905/3.14159265358979/1.414213562\
38/0.70710706/0/0 float=y maxiter=6000 bailout=25
inside=0 logmap=-614 symmetry=none periodicity=10
colors=000E0eE0e<3>L0mN0oP0qR0s<2>X6xZ8z`AzZCxVEuNGm\
HHdAJX2LN0NG0P60P00T20V40X80ZA0bC0dG0eH0gJ2iN6mP8oTA\
qVCsXEu`HxbJzdLzgNziPzmTzoVzqXzuZzv`zxezsizmqogxkbzg\
ZzbTz`NzVHzPCzJ8X8kZ6z<2>Z6zZ6z`6z`6z`4z`4z`4z`4zb4z\
b4vb4Ab4Ab2Ab2ed2gd2id2id2kd2md2mb6q`8sZAuXEvVGxTHzR\
JzPNzNPzLRzJTzHXz<2>CbzGdzJezNgzPixTkuXmq`okbqgesdiu\
`mvVoxRszNvzHzzEzzAzz4zz0zz0zz0zzHzzezzdzzdzzbzzbzz`\
zz`zzZzzZzzXzzXxzVvzVuzTuzTqzNmzJkzGgzCdz8bz4Zz0Vz0T\
z0Pz0Nz0Lz0Lz0Lz0Jz0Jz0Jz0Hz0Hz0Hz0Gz2Gz4Gz6Gz6Hz6Jz\
4Lz4Nz2Pz2Rz0Tz0Vz0Xz0`z0PzAVzAzz6zz0<4>zz0zz0zz0zz0\
zz2zz6vzAuzCqzGozHkzLizPezRdzV`zZZz`VzdTzeZz`bzVezRi\
zLmzHszCvz6zz2zz0zz0zz0xz4uz8qzEmzHgzNdzR`zX<2>PziRz\
kRzkRzkRzmRzmRzmRzmTzoTzoTzoTzoTzqTzqTzqTzqVzsVzsVzu\
<3>VzvVzvVzxVzxVzz<3>VzzVzzTzzTzzRzzRzzRzzPzzPzz0zX0\
zZ0zb4ze
}
END PARAMETER FILE=========================================
START 20.0 PAR-FORMULA FILE================================
Cyclogenesis { ; time=0:41:22.64, SF5 on a p200
reset=2000 type=formula formulafile=critical.frm
formulaname=MandelbrotMix4 function=ident passes=1
center-mag=-0.23418923607130600/+0.42864036781657180\
/1.288888e+011/1/90
params=2.71828182845905/3.14159265358979/1.414213562\
38/0.70710706/0/0 float=y maxiter=6000 bailout=25
inside=0 logmap=-614 symmetry=none periodicity=10
colors=000E0eE0e<3>L0mN0oP0qR0s<2>X6xZ8z`AzZCxVEuNGm\
HHdAJX2LN0NG0P60P00T20V40X80ZA0bC0dG0eH0gJ2iN6mP8oTA\
qVCsXEu`HxbJzdLzgNziPzmTzoVzqXzuZzv`zxezsizmqogxkbzg\
ZzbTz`NzVHzPCzJ8X8kZ6z<2>Z6zZ6z`6z`6z`4z`4z`4z`4zb4z\
b4vb4Ab4Ab2Ab2ed2gd2id2id2kd2md2mb6q`8sZAuXEvVGxTHzR\
JzPNzNPzLRzJTzHXz<2>CbzGdzJezNgzPixTkuXmq`okbqgesdiu\
`mvVoxRszNvzHzzEzzAzz4zz0zz0zz0zzHzzezzdzzdzzbzzbzz`\
zz`zzZzzZzzXzzXxzVvzVuzTuzTqzNmzJkzGgzCdz8bz4Zz0Vz0T\
z0Pz0Nz0Lz0Lz0Lz0Jz0Jz0Jz0Hz0Hz0Hz0Gz2Gz4Gz6Gz6Hz6Jz\
4Lz4Nz2Pz2Rz0Tz0Vz0Xz0`z0PzAVzAzz6zz0<4>zz0zz0zz0zz0\
zz2zz6vzAuzCqzGozHkzLizPezRdzV`zZZz`VzdTzeZz`bzVezRi\
zLmzHszCvz6zz2zz0zz0zz0xz4uz8qzEmzHgzNdzR`zX<2>PziRz\
kRzkRzkRzmRzmRzmRzmTzoTzoTzoTzoTzqTzqTzqTzqVzsVzsVzu\
<3>VzvVzvVzxVzxVzz<3>VzzVzzTzzTzzRzzRzzRzzPzzPzz0zX0\
zZ0zb4ze
}
frm:MandelbrotMix4 {; Jim Muth
a=real(p1), b=imag(p1), d=real(p2), f=imag(p2),
g=1/f, h=1/d, j=1/(f-b), z=(-a*b*g*h)^j,
k=real(p3)+1, l=imag(p3)+100, c=fn1(pixel):
z=k*((a*(z^b))+(d*(z^f)))+c,
|z| < l
}
END 20.0 PAR-FORMULA FILE==================================
- --------------------------------------------------------------
Thanks for using Fractint, The Fractals and Fractint Discussion List
Post Message: fractint@lists.xmission.com
Get Commands: majordomo@lists.xmission.com "help"
Administrator: twegner@swbell.net
Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
------------------------------
Date: Thu, 20 Jan 2000 14:19:59 -0800
From: Gregory McClure <Gregory.McClure@quantum.com>
Subject: (fractint) GregsBrn.frm -- Barnsley-Style Formulas...
Here are some formulas in the same vein as the GregsMan formulas posted
earlier...
File GregsBrn.frm
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D
GregsBarnsleyC2E {; =A92000 Greg McClure -- Dual func with even =
constants
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; imag(p3) =3D type =3D 0/MOD, 1/REAL, 2/IMAG, 3/OR, 4/AND, 5/MANH, =
6/MANR
; p1 =3D 0/0, p2 =3D 1/0, p3 =3D 2/0, fn1/fn2 =3D ident/ident =3D =
standard Barnsley
z =3D p1 + Pixel:
ip =3D imag(p3)
IF(IP<0)
ip =3D abs(ip)-1
z =3D z + p1
ENDIF
IF(real(z)>=3D0)
z=3D(fn1(z)-p2)*Pixel
ELSE
z=3D(fn2(z)+p2)*Pixel
ENDIF,
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D real(z)
ELSEIF(ip<2.1)
ct =3D imag(z)
ELSEIF((ip<3.1) && (|rz|>|iz|))
ct =3D real(z)
ELSEIF((ip<3.1) && (|rz|<=3D|iz|))
ct =3D imag(z)
ELSEIF((ip<4.1) && (|rz|<|iz|))
ct =3D real(z)
ELSEIF((ip<4.1) && (|rz|>=3D|iz|))
ct =3D imag(z)
ELSEIF(ip<5.1)
ct =3D (abs(real(z)) + abs(imag(z)))
ELSEIF(ip<6.1)
ct =3D (real(z) + imag(z))
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsBarnsleyC2P {; =A92000 Greg McClure -- Dual func with pos constant
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; imag(p3) =3D type =3D 0/MOD, 1/REAL, 2/IMAG, 3/OR, 4/AND, 5/MANH, =
6/MANR
; p1 =3D 0/0, p2 =3D 1/0, p3 =3D 2/0, fn1/fn2 =3D ident/ident =3D =
standard Barnsley
z =3D p1 + Pixel:
ip =3D imag(p3)
IF(IP<0)
ip =3D abs(ip)-1
z =3D z + p1
ENDIF
IF(real(z)>=3D0)
z=3D(fn1(z)-p2)*Pixel
ELSE
z=3D(fn2(z)+1)*Pixel
ENDIF,
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D real(z)
ELSEIF(ip<2.1)
ct =3D imag(z)
ELSEIF((ip<3.1) && (|rz|>|iz|))
ct =3D real(z)
ELSEIF((ip<3.1) && (|rz|<=3D|iz|))
ct =3D imag(z)
ELSEIF((ip<4.1) && (|rz|<|iz|))
ct =3D real(z)
ELSEIF((ip<4.1) && (|rz|>=3D|iz|))
ct =3D imag(z)
ELSEIF(ip<5.1)
ct =3D (abs(real(z)) + abs(imag(z)))
ELSEIF(ip<6.1)
ct =3D (real(z) + imag(z))
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsBarnsleyC2N {; =A92000 Greg McClure -- Dual func with neg constant
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; imag(p3) =3D type =3D 0/MOD, 1/REAL, 2/IMAG, 3/OR, 4/AND, 5/MANH, =
6/MANR
; p1 =3D 0/0, p2 =3D 1/0, p3 =3D 2/0, fn1/fn2 =3D ident/ident =3D =
standard Barnsley
z =3D p1 + Pixel:
ip =3D imag(p3)
IF(IP<0)
ip =3D abs(ip)-1
z =3D z + p1
ENDIF
IF(real(z)>=3D0)
z=3D(fn1(z)-1)*Pixel
ELSE
z=3D(fn2(z)+p2)*Pixel
ENDIF,
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D real(z)
ELSEIF(ip<2.1)
ct =3D imag(z)
ELSEIF((ip<3.1) && (|rz|>|iz|))
ct =3D real(z)
ELSEIF((ip<3.1) && (|rz|<=3D|iz|))
ct =3D imag(z)
ELSEIF((ip<4.1) && (|rz|<|iz|))
ct =3D real(z)
ELSEIF((ip<4.1) && (|rz|>=3D|iz|))
ct =3D imag(z)
ELSEIF(ip<5.1)
ct =3D (abs(real(z)) + abs(imag(z)))
ELSEIF(ip<6.1)
ct =3D (real(z) + imag(z))
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsBarnsleyM2E {; =A92000 Greg McClure -- Dual func with even mult
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; imag(p3) =3D type =3D 0/MOD, 1/REAL, 2/IMAG, 3/OR, 4/AND, 5/MANH, =
6/MANR
; p1 =3D 0/0, p2 =3D 1/0, p3 =3D 2/0, fn1/fn2 =3D ident/ident =3D =
standard Barnsley
z =3D p1 + Pixel:
ip =3D imag(p3)
IF(IP<0)
ip =3D abs(ip)-1
z =3D z + p1
ENDIF
IF(real(z)>=3D0)
z=3Dp2*(fn1(z)-1)*Pixel
ELSE
z=3Dp2*(fn2(z)+1)*Pixel
ENDIF,
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D real(z)
ELSEIF(ip<2.1)
ct =3D imag(z)
ELSEIF((ip<3.1) && (|rz|>|iz|))
ct =3D real(z)
ELSEIF((ip<3.1) && (|rz|<=3D|iz|))
ct =3D imag(z)
ELSEIF((ip<4.1) && (|rz|<|iz|))
ct =3D real(z)
ELSEIF((ip<4.1) && (|rz|>=3D|iz|))
ct =3D imag(z)
ELSEIF(ip<5.1)
ct =3D (abs(real(z)) + abs(imag(z)))
ELSEIF(ip<6.1)
ct =3D (real(z) + imag(z))
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsBarnsleyM2P {; =A92000 Greg McClure -- Dual func with pos mult
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; imag(p3) =3D type =3D 0/MOD, 1/REAL, 2/IMAG, 3/OR, 4/AND, 5/MANH, =
6/MANR
; p1 =3D 0/0, p2 =3D 1/0, p3 =3D 2/0, fn1/fn2 =3D ident/ident =3D =
standard Barnsley
z =3D p1 + Pixel:
ip =3D imag(p3)
IF(IP<0)
ip =3D abs(ip)-1
z =3D z + p1
ENDIF
IF(real(z)>=3D0)
z=3Dp2*(fn1(z)-1)*Pixel
ELSE
z=3D(fn2(z)+1)*Pixel
ENDIF,
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D real(z)
ELSEIF(ip<2.1)
ct =3D imag(z)
ELSEIF((ip<3.1) && (|rz|>|iz|))
ct =3D real(z)
ELSEIF((ip<3.1) && (|rz|<=3D|iz|))
ct =3D imag(z)
ELSEIF((ip<4.1) && (|rz|<|iz|))
ct =3D real(z)
ELSEIF((ip<4.1) && (|rz|>=3D|iz|))
ct =3D imag(z)
ELSEIF(ip<5.1)
ct =3D (abs(real(z)) + abs(imag(z)))
ELSEIF(ip<6.1)
ct =3D (real(z) + imag(z))
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsBarnsleyM2N {; =A92000 Greg McClure -- Dual func with neg mult
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; imag(p3) =3D type =3D 0/MOD, 1/REAL, 2/IMAG, 3/OR, 4/AND, 5/MANH, =
6/MANR
; p1 =3D 0/0, p2 =3D 1/0, p3 =3D 2/0, fn1/fn2 =3D ident/ident =3D =
standard Barnsley
z =3D p1 + Pixel:
ip =3D imag(p3)
IF(IP<0)
ip =3D abs(ip)-1
z =3D z + p1
ENDIF
IF(real(z)>=3D0)
z=3D(fn1(z)-1)*Pixel
ELSE
z=3Dp2*(fn2(z)+1)*Pixel
ENDIF,
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D real(z)
ELSEIF(ip<2.1)
ct =3D imag(z)
ELSEIF((ip<3.1) && (|rz|>|iz|))
ct =3D real(z)
ELSEIF((ip<3.1) && (|rz|<=3D|iz|))
ct =3D imag(z)
ELSEIF((ip<4.1) && (|rz|<|iz|))
ct =3D real(z)
ELSEIF((ip<4.1) && (|rz|>=3D|iz|))
ct =3D imag(z)
ELSEIF(ip<5.1)
ct =3D (abs(real(z)) + abs(imag(z)))
ELSEIF(ip<6.1)
ct =3D (real(z) + imag(z))
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsBarnsJulC2E {; =A92000 Greg McClure -- Dual func with even =
constants
; p1 =3D point, p2 =3D multiplier, p3 =3D cutoff/type
; imag(p3) =3D type =3D 0/MOD, 1/REAL, 2/IMAG, 3/OR, 4/AND, 5/MANH, =
6/MANR
; p2 =3D 1/0, p3 =3D 2/0, fn1/fn2 =3D ident/ident =3D standard Barns =
Julia f/p1
z =3D Pixel:
IF(real(z)>=3D0)
z=3D(fn1(z)-p2)*p1
ELSE
z=3D(fn2(z)+p2)*p1
ENDIF,
ip =3D imag(p3)
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D real(z)
ELSEIF(ip<2.1)
ct =3D imag(z)
ELSEIF((ip<3.1) && (|rz|>|iz|))
ct =3D real(z)
ELSEIF((ip<3.1) && (|rz|<=3D|iz|))
ct =3D imag(z)
ELSEIF((ip<4.1) && (|rz|<|iz|))
ct =3D real(z)
ELSEIF((ip<4.1) && (|rz|>=3D|iz|))
ct =3D imag(z)
ELSEIF(ip<5.1)
ct =3D (abs(real(z)) + abs(imag(z)))
ELSEIF(ip<6.1)
ct =3D (real(z) + imag(z))
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsBarnsJulC2P {; =A92000 Greg McClure -- Dual func with pos constant
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; imag(p3) =3D type =3D 0/MOD, 1/REAL, 2/IMAG, 3/OR, 4/AND, 5/MANH, =
6/MANR
; p2 =3D 1/0, p3 =3D 2/0, fn1/fn2 =3D ident/ident =3D standard Barns =
Julia f/p1
z =3D Pixel:
IF(real(z)>=3D0)
z=3D(fn1(z)-p2)*p1
ELSE
z=3D(fn2(z)+1)*p1
ENDIF,
ip =3D imag(p3)
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D real(z)
ELSEIF(ip<2.1)
ct =3D imag(z)
ELSEIF((ip<3.1) && (|rz|>|iz|))
ct =3D real(z)
ELSEIF((ip<3.1) && (|rz|<=3D|iz|))
ct =3D imag(z)
ELSEIF((ip<4.1) && (|rz|<|iz|))
ct =3D real(z)
ELSEIF((ip<4.1) && (|rz|>=3D|iz|))
ct =3D imag(z)
ELSEIF(ip<5.1)
ct =3D (abs(real(z)) + abs(imag(z)))
ELSEIF(ip<6.1)
ct =3D (real(z) + imag(z))
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsBarnsJulC2N {; =A92000 Greg McClure -- Dual func with neg constant
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; imag(p3) =3D type =3D 0/MOD, 1/REAL, 2/IMAG, 3/OR, 4/AND, 5/MANH, =
6/MANR
; p2 =3D 1/0, p3 =3D 2/0, fn1/fn2 =3D ident/ident =3D standard Barns =
Julia f/p1
z =3D Pixel:
IF(real(z)>=3D0)
z=3D(fn1(z)-1)*p1
ELSE
z=3D(fn2(z)+p2)*p1
ENDIF,
ip =3D imag(p3)
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D real(z)
ELSEIF(ip<2.1)
ct =3D imag(z)
ELSEIF((ip<3.1) && (|rz|>|iz|))
ct =3D real(z)
ELSEIF((ip<3.1) && (|rz|<=3D|iz|))
ct =3D imag(z)
ELSEIF((ip<4.1) && (|rz|<|iz|))
ct =3D real(z)
ELSEIF((ip<4.1) && (|rz|>=3D|iz|))
ct =3D imag(z)
ELSEIF(ip<5.1)
ct =3D (abs(real(z)) + abs(imag(z)))
ELSEIF(ip<6.1)
ct =3D (real(z) + imag(z))
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsBarnsJulM2E {; =A92000 Greg McClure -- Dual func with even mult
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; imag(p3) =3D type =3D 0/MOD, 1/REAL, 2/IMAG, 3/OR, 4/AND, 5/MANH, =
6/MANR
; p2 =3D 1/0, p3 =3D 2/0, fn1/fn2 =3D ident/ident =3D standard Barns =
Julia f/p1
z =3D Pixel:
IF(real(z)>=3D0)
z=3Dp2*(fn1(z)-1)*p1
ELSE
z=3Dp2*(fn2(z)+1)*p1
ENDIF,
ip =3D imag(p3)
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D real(z)
ELSEIF(ip<2.1)
ct =3D imag(z)
ELSEIF((ip<3.1) && (|rz|>|iz|))
ct =3D real(z)
ELSEIF((ip<3.1) && (|rz|<=3D|iz|))
ct =3D imag(z)
ELSEIF((ip<4.1) && (|rz|<|iz|))
ct =3D real(z)
ELSEIF((ip<4.1) && (|rz|>=3D|iz|))
ct =3D imag(z)
ELSEIF(ip<5.1)
ct =3D (abs(real(z)) + abs(imag(z)))
ELSEIF(ip<6.1)
ct =3D (real(z) + imag(z))
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsBarnsJulM2P {; =A92000 Greg McClure -- Dual func with pos mult
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; imag(p3) =3D type =3D 0/MOD, 1/REAL, 2/IMAG, 3/OR, 4/AND, 5/MANH, =
6/MANR
; p2 =3D 1/0, p3 =3D 2/0, fn1/fn2 =3D ident/ident =3D standard Barns =
Julia f/p1
z =3D Pixel:
IF(real(z)>=3D0)
z=3Dp2*(fn1(z)-1)*p1
ELSE
z=3D(fn2(z)+1)*p1
ENDIF,
ip =3D imag(p3)
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D real(z)
ELSEIF(ip<2.1)
ct =3D imag(z)
ELSEIF((ip<3.1) && (|rz|>|iz|))
ct =3D real(z)
ELSEIF((ip<3.1) && (|rz|<=3D|iz|))
ct =3D imag(z)
ELSEIF((ip<4.1) && (|rz|<|iz|))
ct =3D real(z)
ELSEIF((ip<4.1) && (|rz|>=3D|iz|))
ct =3D imag(z)
ELSEIF(ip<5.1)
ct =3D (abs(real(z)) + abs(imag(z)))
ELSEIF(ip<6.1)
ct =3D (real(z) + imag(z))
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsBarnsJulM2N {; =A92000 Greg McClure -- Dual func with neg mult
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; imag(p3) =3D type =3D 0/MOD, 1/REAL, 2/IMAG, 3/OR, 4/AND, 5/MANH, =
6/MANR
; p2 =3D 1/0, p3 =3D 2/0, fn1/fn2 =3D ident/ident =3D standard Barns =
Julia f/p1
z =3D Pixel:
IF(real(z)>=3D0)
z=3D(fn1(z)-1)*p1
ELSE
z=3Dp2*(fn2(z)+1)*p1
ENDIF,
ip =3D imag(p3)
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D real(z)
ELSEIF(ip<2.1)
ct =3D imag(z)
ELSEIF((ip<3.1) && (|rz|>|iz|))
ct =3D real(z)
ELSEIF((ip<3.1) && (|rz|<=3D|iz|))
ct =3D imag(z)
ELSEIF((ip<4.1) && (|rz|<|iz|))
ct =3D real(z)
ELSEIF((ip<4.1) && (|rz|>=3D|iz|))
ct =3D imag(z)
ELSEIF(ip<5.1)
ct =3D (abs(real(z)) + abs(imag(z)))
ELSEIF(ip<6.1)
ct =3D (real(z) + imag(z))
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D
The Kwisatz Haderach,
=DF Gregory J. McClure
- --------------------------------------------------------------
Thanks for using Fractint, The Fractals and Fractint Discussion List
Post Message: fractint@lists.xmission.com
Get Commands: majordomo@lists.xmission.com "help"
Administrator: twegner@swbell.net
Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
------------------------------
Date: Thu, 20 Jan 2000 14:17:18 -0800
From: Gregory McClure <Gregory.McClure@quantum.com>
Subject: (fractint) GregsHyp.frm -- HyperComplex Formulas...
Here are some formulas in the same vein as the GregsMan formulas posted
earlier...
File GregsHyp.frm
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D
GregsMandelHM1 {; =A92000 Greg McClure -- HyperMandel (Type ++--)
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; z(0) =3D (xP1+xPixel,yP1+yPixel,xP2,yP2)
; z(n) =3D f[z(n-1)] + (xPixel,yPixel,xP2,yP2)
; imag(p3) =3D type =3D 0/MOD, 1/MODRI, 2/MODJK, 3/REAL, 4/IMAG,
; =3D 5/JMAG, 6/KMAG, 7/OR, 8/AND, 9/MANH, 10/MANR
z1 =3D p1 + Pixel, z2 =3D p2, z =3D 0:
f1 =3D real(z1) - imag(z2) + (0,1) * (imag(z1) + real(z2))
f2 =3D real(z1) + imag(z2) + (0,1) * (imag(z1) - real(z2))
f1 =3D fn1(f1)
f2 =3D fn2(f2)
z1 =3D (real(f1) + real(f2) + (0,1) * (imag(f1) + imag(f2))) / 2 + =
Pixel
z2 =3D (imag(f1) - imag(f2) + (0,1) * (real(f2) - real(f1))) / 2 + =
p2
z =3D sqrt(|z1| + |z2|)
ip =3D imag(p3)
rz =3D real(z1)
iz =3D imag(z1)
jz =3D real(z2)
kz =3D imag(z2)
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D z1
ELSEIF(ip<2.1)
ct =3D z2
ELSEIF(ip<3.1)
ct =3D rz
ELSEIF(ip<4.1)
ct =3D iz
ELSEIF(ip<5.1)
ct =3D jz
ELSEIF(ip<6.1)
ct =3D kz
ELSEIF((ip<7.1) && (|rz|>=3D|iz|) && (|rz|>=3D|jz|) && =
(|rz|>=3D|kz|))
ct =3D rz
ELSEIF((ip<7.1) && (|iz|>=3D|rz|) && (|iz|>=3D|jz|) && =
(|iz|>=3D|kz|))
ct =3D iz
ELSEIF((ip<7.1) && (|jz|>=3D|rz|) && (|jz|>=3D|iz|) && =
(|jz|>=3D|kz|))
ct =3D jz
ELSEIF((ip<7.1) && (|kz|>=3D|rz|) && (|iz|>=3D|kz|) && =
(|kz|>=3D|jz|))
ct =3D kz
ELSEIF((ip<8.1) && (|rz|<|iz|) && (|rz|<|jz|) && (|rz|<|kz|))
ct =3D rz
ELSEIF((ip<8.1) && (|iz|<|rz|) && (|iz|<|jz|) && (|iz|<|kz|))
ct =3D iz
ELSEIF((ip<8.1) && (|jz|<|rz|) && (|jz|<|iz|) && (|jz|<|kz|))
ct =3D jz
ELSEIF((ip<8.1) && (|kz|<|rz|) && (|kz|<|iz|) && (|kz|<|jz|))
ct =3D kz
ELSEIF(ip<9.1)
ct =3D (abs(rz) + abs(iz) + abs(jz) + abs(kz))
ELSEIF(ip<10.1)
ct =3D (rz + iz + jz + kz)
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsMandelHM2 {; =A92000 Greg McClure -- HyperMandel (Type --++)
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; z(0) =3D (xP2,yP2,xP1+xPixel,yP1+yPixel)
; z(n) =3D f[z(n-1)] + (xP2,yP2,xPixel,yPixel)
; imag(p3) =3D type =3D 0/MOD, 1/MODRI, 2/MODJK, 3/REAL, 4/IMAG,
; =3D 5/JMAG, 6/KMAG, 7/OR, 8/AND, 9/MANH, 10/MANR
z1 =3D p2, z2 =3D p1 + Pixel, z =3D 0:
f1 =3D real(z1) - imag(z2) + (0,1) * (imag(z1) + real(z2))
f2 =3D real(z1) + imag(z2) + (0,1) * (imag(z1) - real(z2))
f1 =3D fn1(f1)
f2 =3D fn2(f2)
z1 =3D (real(f1) + real(f2) + (0,1) * (imag(f1) + imag(f2))) / 2 + =
p2
z2 =3D (imag(f1) - imag(f2) + (0,1) * (real(f2) - real(f1))) / 2 + =
Pixel
z =3D sqrt(|z1| + |z2|)
ip =3D imag(p3)
rz =3D real(z1)
iz =3D imag(z1)
jz =3D real(z2)
kz =3D imag(z2)
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D z1
ELSEIF(ip<2.1)
ct =3D z2
ELSEIF(ip<3.1)
ct =3D rz
ELSEIF(ip<4.1)
ct =3D iz
ELSEIF(ip<5.1)
ct =3D jz
ELSEIF(ip<6.1)
ct =3D kz
ELSEIF((ip<7.1) && (|rz|>=3D|iz|) && (|rz|>=3D|jz|) && =
(|rz|>=3D|kz|))
ct =3D rz
ELSEIF((ip<7.1) && (|iz|>=3D|rz|) && (|iz|>=3D|jz|) && =
(|iz|>=3D|kz|))
ct =3D iz
ELSEIF((ip<7.1) && (|jz|>=3D|rz|) && (|jz|>=3D|iz|) && =
(|jz|>=3D|kz|))
ct =3D jz
ELSEIF((ip<7.1) && (|kz|>=3D|rz|) && (|iz|>=3D|kz|) && =
(|kz|>=3D|jz|))
ct =3D kz
ELSEIF((ip<8.1) && (|rz|<|iz|) && (|rz|<|jz|) && (|rz|<|kz|))
ct =3D rz
ELSEIF((ip<8.1) && (|iz|<|rz|) && (|iz|<|jz|) && (|iz|<|kz|))
ct =3D iz
ELSEIF((ip<8.1) && (|jz|<|rz|) && (|jz|<|iz|) && (|jz|<|kz|))
ct =3D jz
ELSEIF((ip<8.1) && (|kz|<|rz|) && (|kz|<|iz|) && (|kz|<|jz|))
ct =3D kz
ELSEIF(ip<9.1)
ct =3D (abs(rz) + abs(iz) + abs(jz) + abs(kz))
ELSEIF(ip<10.1)
ct =3D (rz + iz + jz + kz)
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsMandelHMA(XAXIS) {; =A92000 Greg McClure -- HyperMandel (Type =
++--)
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; z(0) =3D (xP1+xPixel,yP1+yPixel,xP2,yP2) \ FORCED SYMMETRY
; z(n) =3D f[z(n-1)] + (xPixel,yPixel,xP2,yP2) / ON X-AXIS
; imag(p3) =3D type =3D 0/MOD, 1/MODRI, 2/MODJK, 3/REAL, 4/IMAG,
; =3D 5/JMAG, 6/KMAG, 7/OR, 8/AND, 9/MANH, 10/MANR
z1 =3D p1 + Pixel, z2 =3D p2, z =3D 0:
f1 =3D real(z1) - imag(z2) + (0,1) * (imag(z1) + real(z2))
f2 =3D real(z1) + imag(z2) + (0,1) * (imag(z1) - real(z2))
f1 =3D fn1(f1)
f2 =3D fn2(f2)
z1 =3D (real(f1) + real(f2) + (0,1) * (imag(f1) + imag(f2))) / 2 + =
Pixel
z2 =3D (imag(f1) - imag(f2) + (0,1) * (real(f2) - real(f1))) / 2 + =
p2
z =3D sqrt(|z1| + |z2|)
ip =3D imag(p3)
rz =3D real(z1)
iz =3D imag(z1)
jz =3D real(z2)
kz =3D imag(z2)
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D z1
ELSEIF(ip<2.1)
ct =3D z2
ELSEIF(ip<3.1)
ct =3D rz
ELSEIF(ip<4.1)
ct =3D iz
ELSEIF(ip<5.1)
ct =3D jz
ELSEIF(ip<6.1)
ct =3D kz
ELSEIF((ip<7.1) && (|rz|>=3D|iz|) && (|rz|>=3D|jz|) && =
(|rz|>=3D|kz|))
ct =3D rz
ELSEIF((ip<7.1) && (|iz|>=3D|rz|) && (|iz|>=3D|jz|) && =
(|iz|>=3D|kz|))
ct =3D iz
ELSEIF((ip<7.1) && (|jz|>=3D|rz|) && (|jz|>=3D|iz|) && =
(|jz|>=3D|kz|))
ct =3D jz
ELSEIF((ip<7.1) && (|kz|>=3D|rz|) && (|iz|>=3D|kz|) && =
(|kz|>=3D|jz|))
ct =3D kz
ELSEIF((ip<8.1) && (|rz|<|iz|) && (|rz|<|jz|) && (|rz|<|kz|))
ct =3D rz
ELSEIF((ip<8.1) && (|iz|<|rz|) && (|iz|<|jz|) && (|iz|<|kz|))
ct =3D iz
ELSEIF((ip<8.1) && (|jz|<|rz|) && (|jz|<|iz|) && (|jz|<|kz|))
ct =3D jz
ELSEIF((ip<8.1) && (|kz|<|rz|) && (|kz|<|iz|) && (|kz|<|jz|))
ct =3D kz
ELSEIF(ip<9.1)
ct =3D (abs(rz) + abs(iz) + abs(jz) + abs(kz))
ELSEIF(ip<10.1)
ct =3D (rz + iz + jz + kz)
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsMandelHMB(XAXIS) {; =A92000 Greg McClure -- HyperMandel (Type =
- --++)
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; z(0) =3D (xP2,yP2,xP1+xPixel,yP1+yPixel) \ FORCED SYMMETRY
; z(n) =3D f[z(n-1)] + (xP2,yP2,xPixel,yPixel) / ON X-AXIS
; imag(p3) =3D type =3D 0/MOD, 1/MODRI, 2/MODJK, 3/REAL, 4/IMAG,
; =3D 5/JMAG, 6/KMAG, 7/OR, 8/AND, 9/MANH, 10/MANR
z1 =3D p2, z2 =3D p1 + Pixel, z =3D 0:
f1 =3D real(z1) - imag(z2) + (0,1) * (imag(z1) + real(z2))
f2 =3D real(z1) + imag(z2) + (0,1) * (imag(z1) - real(z2))
f1 =3D fn1(f1)
f2 =3D fn2(f2)
z1 =3D (real(f1) + real(f2) + (0,1) * (imag(f1) + imag(f2))) / 2 + =
p2
z2 =3D (imag(f1) - imag(f2) + (0,1) * (real(f2) - real(f1))) / 2 + =
Pixel
z =3D sqrt(|z1| + |z2|)
ip =3D imag(p3)
rz =3D real(z1)
iz =3D imag(z1)
jz =3D real(z2)
kz =3D imag(z2)
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D z1
ELSEIF(ip<2.1)
ct =3D z2
ELSEIF(ip<3.1)
ct =3D rz
ELSEIF(ip<4.1)
ct =3D iz
ELSEIF(ip<5.1)
ct =3D jz
ELSEIF(ip<6.1)
ct =3D kz
ELSEIF((ip<7.1) && (|rz|>=3D|iz|) && (|rz|>=3D|jz|) && =
(|rz|>=3D|kz|))
ct =3D rz
ELSEIF((ip<7.1) && (|iz|>=3D|rz|) && (|iz|>=3D|jz|) && =
(|iz|>=3D|kz|))
ct =3D iz
ELSEIF((ip<7.1) && (|jz|>=3D|rz|) && (|jz|>=3D|iz|) && =
(|jz|>=3D|kz|))
ct =3D jz
ELSEIF((ip<7.1) && (|kz|>=3D|rz|) && (|iz|>=3D|kz|) && =
(|kz|>=3D|jz|))
ct =3D kz
ELSEIF((ip<8.1) && (|rz|<|iz|) && (|rz|<|jz|) && (|rz|<|kz|))
ct =3D rz
ELSEIF((ip<8.1) && (|iz|<|rz|) && (|iz|<|jz|) && (|iz|<|kz|))
ct =3D iz
ELSEIF((ip<8.1) && (|jz|<|rz|) && (|jz|<|iz|) && (|jz|<|kz|))
ct =3D jz
ELSEIF((ip<8.1) && (|kz|<|rz|) && (|kz|<|iz|) && (|kz|<|jz|))
ct =3D kz
ELSEIF(ip<9.1)
ct =3D (abs(rz) + abs(iz) + abs(jz) + abs(kz))
ELSEIF(ip<10.1)
ct =3D (rz + iz + jz + kz)
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsJuliaHM1 {; =A92000 Greg McClure -- HyperJulia (Type ++00)
; p1 =3D point, p2 =3D multiplier, p3 =3D cutoff/type
; z(-1) =3D (xPixel,yPixel,0,0)
; z(n) =3D f[z(n-1)] + (xP1,yP1,xP2,yP2)
; imag(p3) =3D type =3D 0/MOD, 1/MODRI, 2/MODJK, 3/REAL, 4/IMAG,
; =3D 5/JMAG, 6/KMAG, 7/OR, 8/AND, 9/MANH, 10/MANR
z1 =3D Pixel, z2 =3D 0, z =3D 0,
f1 =3D real(z1) - imag(z2) + (0,1) * (imag(z1) + real(z2)),
f2 =3D real(z1) + imag(z2) + (0,1) * (imag(z1) - real(z2)),
f1 =3D fn1(f1), f2 =3D fn2(f2),
z1 =3D (real(f1) + real(f2) + (0,1) * (imag(f1) + imag(f2))) / 2 + =
p1,
z2 =3D (imag(f1) - imag(f2) + (0,1) * (real(f2) - real(f1))) / 2 + =
p2:
f1 =3D real(z1) - imag(z2) + (0,1) * (imag(z1) + real(z2))
f2 =3D real(z1) + imag(z2) + (0,1) * (imag(z1) - real(z2))
f1 =3D fn1(f1)
f2 =3D fn2(f2)
z1 =3D (real(f1) + real(f2) + (0,1) * (imag(f1) + imag(f2))) / 2 + =
p1
z2 =3D (imag(f1) - imag(f2) + (0,1) * (real(f2) - real(f1))) / 2 + =
p2
z =3D sqrt(|z1| + |z2|)
ip =3D imag(p3)
rz =3D real(z1)
iz =3D imag(z1)
jz =3D real(z2)
kz =3D imag(z2)
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D z1
ELSEIF(ip<2.1)
ct =3D z2
ELSEIF(ip<3.1)
ct =3D rz
ELSEIF(ip<4.1)
ct =3D iz
ELSEIF(ip<5.1)
ct =3D jz
ELSEIF(ip<6.1)
ct =3D kz
ELSEIF((ip<7.1) && (|rz|>=3D|iz|) && (|rz|>=3D|jz|) && =
(|rz|>=3D|kz|))
ct =3D rz
ELSEIF((ip<7.1) && (|iz|>=3D|rz|) && (|iz|>=3D|jz|) && =
(|iz|>=3D|kz|))
ct =3D iz
ELSEIF((ip<7.1) && (|jz|>=3D|rz|) && (|jz|>=3D|iz|) && =
(|jz|>=3D|kz|))
ct =3D jz
ELSEIF((ip<7.1) && (|kz|>=3D|rz|) && (|iz|>=3D|kz|) && =
(|kz|>=3D|jz|))
ct =3D kz
ELSEIF((ip<8.1) && (|rz|<|iz|) && (|rz|<|jz|) && (|rz|<|kz|))
ct =3D rz
ELSEIF((ip<8.1) && (|iz|<|rz|) && (|iz|<|jz|) && (|iz|<|kz|))
ct =3D iz
ELSEIF((ip<8.1) && (|jz|<|rz|) && (|jz|<|iz|) && (|jz|<|kz|))
ct =3D jz
ELSEIF((ip<8.1) && (|kz|<|rz|) && (|kz|<|iz|) && (|kz|<|jz|))
ct =3D kz
ELSEIF(ip<9.1)
ct =3D (abs(rz) + abs(iz) + abs(jz) + abs(kz))
ELSEIF(ip<10.1)
ct =3D (rz + iz + jz + kz)
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsJuliaHM2 {; =A92000 Greg McClure -- HyperJulia (Type 00++)
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; z(-1) =3D (0,0,xPixel,yPixel)
; z(n) =3D f[z(n-1)] + (xP2,yP2,xP1,yP1)
; imag(p3) =3D type =3D 0/MOD, 1/MODRI, 2/MODJK, 3/REAL, 4/IMAG,
; =3D 5/JMAG, 6/KMAG, 7/OR, 8/AND, 9/MANH, 10/MANR
z1 =3D 0, z2 =3D Pixel, z =3D 0,
f1 =3D real(z1) - imag(z2) + (0,1) * (imag(z1) + real(z2)),
f2 =3D real(z1) + imag(z2) + (0,1) * (imag(z1) - real(z2)),
f1 =3D fn1(f1), f2 =3D fn2(f2),
z1 =3D (real(f1) + real(f2) + (0,1) * (imag(f1) + imag(f2))) / 2 + =
p2,
z2 =3D (imag(f1) - imag(f2) + (0,1) * (real(f2) - real(f1))) / 2 + =
p1:
f1 =3D real(z1) - imag(z2) + (0,1) * (imag(z1) + real(z2))
f2 =3D real(z1) + imag(z2) + (0,1) * (imag(z1) - real(z2))
f1 =3D fn1(f1)
f2 =3D fn2(f2)
z1 =3D (real(f1) + real(f2) + (0,1) * (imag(f1) + imag(f2))) / 2 + =
p2
z2 =3D (imag(f1) - imag(f2) + (0,1) * (real(f2) - real(f1))) / 2 + =
p1
z =3D sqrt(|z1| + |z2|)
ip =3D imag(p3)
rz =3D real(z1)
iz =3D imag(z1)
jz =3D real(z2)
kz =3D imag(z2)
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D z1
ELSEIF(ip<2.1)
ct =3D z2
ELSEIF(ip<3.1)
ct =3D rz
ELSEIF(ip<4.1)
ct =3D iz
ELSEIF(ip<5.1)
ct =3D jz
ELSEIF(ip<6.1)
ct =3D kz
ELSEIF((ip<7.1) && (|rz|>=3D|iz|) && (|rz|>=3D|jz|) && =
(|rz|>=3D|kz|))
ct =3D rz
ELSEIF((ip<7.1) && (|iz|>=3D|rz|) && (|iz|>=3D|jz|) && =
(|iz|>=3D|kz|))
ct =3D iz
ELSEIF((ip<7.1) && (|jz|>=3D|rz|) && (|jz|>=3D|iz|) && =
(|jz|>=3D|kz|))
ct =3D jz
ELSEIF((ip<7.1) && (|kz|>=3D|rz|) && (|iz|>=3D|kz|) && =
(|kz|>=3D|jz|))
ct =3D kz
ELSEIF((ip<8.1) && (|rz|<|iz|) && (|rz|<|jz|) && (|rz|<|kz|))
ct =3D rz
ELSEIF((ip<8.1) && (|iz|<|rz|) && (|iz|<|jz|) && (|iz|<|kz|))
ct =3D iz
ELSEIF((ip<8.1) && (|jz|<|rz|) && (|jz|<|iz|) && (|jz|<|kz|))
ct =3D jz
ELSEIF((ip<8.1) && (|kz|<|rz|) && (|kz|<|iz|) && (|kz|<|jz|))
ct =3D kz
ELSEIF(ip<9.1)
ct =3D (abs(rz) + abs(iz) + abs(jz) + abs(kz))
ELSEIF(ip<10.1)
ct =3D (rz + iz + jz + kz)
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsJuliaHMA(ORIGIN) {; =A92000 Greg McClure -- HyperJulia (Type ++00)
; p1 =3D point, p2 =3D multiplier, p3 =3D cutoff/type
; z(-1) =3D (xPixel,yPixel,0,0) \ FORCED SYMMETRY
; z(n) =3D f[z(n-1)] + (xP1,yP2,xP2,yP2) / AT ORIGIN
; imag(p3) =3D type =3D 0/MOD, 1/MODRI, 2/MODJK, 3/REAL, 4/IMAG,
; =3D 5/JMAG, 6/KMAG, 7/OR, 8/AND, 9/MANH, 10/MANR
z1 =3D Pixel, z2 =3D 0, z =3D 0,
f1 =3D real(z1) - imag(z2) + (0,1) * (imag(z1) + real(z2)),
f2 =3D real(z1) + imag(z2) + (0,1) * (imag(z1) - real(z2)),
f1 =3D fn1(f1), f2 =3D fn2(f2),
z1 =3D (real(f1) + real(f2) + (0,1) * (imag(f1) + imag(f2))) / 2 + =
p1,
z2 =3D (imag(f1) - imag(f2) + (0,1) * (real(f2) - real(f1))) / 2 + =
p2:
f1 =3D real(z1) - imag(z2) + (0,1) * (imag(z1) + real(z2))
f2 =3D real(z1) + imag(z2) + (0,1) * (imag(z1) - real(z2))
f1 =3D fn1(f1)
f2 =3D fn2(f2)
z1 =3D (real(f1) + real(f2) + (0,1) * (imag(f1) + imag(f2))) / 2 + =
p1
z2 =3D (imag(f1) - imag(f2) + (0,1) * (real(f2) - real(f1))) / 2 + =
p2
z =3D sqrt(|z1| + |z2|)
ip =3D imag(p3)
rz =3D real(z1)
iz =3D imag(z1)
jz =3D real(z2)
kz =3D imag(z2)
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D z1
ELSEIF(ip<2.1)
ct =3D z2
ELSEIF(ip<3.1)
ct =3D rz
ELSEIF(ip<4.1)
ct =3D iz
ELSEIF(ip<5.1)
ct =3D jz
ELSEIF(ip<6.1)
ct =3D kz
ELSEIF((ip<7.1) && (|rz|>=3D|iz|) && (|rz|>=3D|jz|) && =
(|rz|>=3D|kz|))
ct =3D rz
ELSEIF((ip<7.1) && (|iz|>=3D|rz|) && (|iz|>=3D|jz|) && =
(|iz|>=3D|kz|))
ct =3D iz
ELSEIF((ip<7.1) && (|jz|>=3D|rz|) && (|jz|>=3D|iz|) && =
(|jz|>=3D|kz|))
ct =3D jz
ELSEIF((ip<7.1) && (|kz|>=3D|rz|) && (|iz|>=3D|kz|) && =
(|kz|>=3D|jz|))
ct =3D kz
ELSEIF((ip<8.1) && (|rz|<|iz|) && (|rz|<|jz|) && (|rz|<|kz|))
ct =3D rz
ELSEIF((ip<8.1) && (|iz|<|rz|) && (|iz|<|jz|) && (|iz|<|kz|))
ct =3D iz
ELSEIF((ip<8.1) && (|jz|<|rz|) && (|jz|<|iz|) && (|jz|<|kz|))
ct =3D jz
ELSEIF((ip<8.1) && (|kz|<|rz|) && (|kz|<|iz|) && (|kz|<|jz|))
ct =3D kz
ELSEIF(ip<9.1)
ct =3D (abs(rz) + abs(iz) + abs(jz) + abs(kz))
ELSEIF(ip<10.1)
ct =3D (rz + iz + jz + kz)
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
GregsJuliaHMB(ORIGIN) {; =A92000 Greg McClure -- HyperJulia (Type 00++)
; p1 =3D offset, p2 =3D multiplier, p3 =3D cutoff/type
; z(-1) =3D (0,0,xPixel,yPixel) \ FORCED SYMMETRY
; z(n) =3D f[z(n-1)] + (xP2,yP2,xP1,yP1) / AT ORIGIN
; imag(p3) =3D type =3D 0/MOD, 1/MODRI, 2/MODJK, 3/REAL, 4/IMAG,
; =3D 5/JMAG, 6/KMAG, 7/OR, 8/AND, 9/MANH, 10/MANR
z1 =3D real(Pixel), z2 =3D (0,1) * imag(Pixel), z =3D 0,
f1 =3D real(z1) - imag(z2) + (0,1) * (imag(z1) + real(z2)),
f2 =3D real(z1) + imag(z2) + (0,1) * (imag(z1) - real(z2)),
f1 =3D fn1(f1), f2 =3D fn2(f2),
z1 =3D (real(f1) + real(f2) + (0,1) * (imag(f1) + imag(f2))) / 2 + =
p2,
z2 =3D (imag(f1) - imag(f2) + (0,1) * (real(f2) - real(f1))) / 2 + =
p1:
f1 =3D real(z1) - imag(z2) + (0,1) * (imag(z1) + real(z2))
f2 =3D real(z1) + imag(z2) + (0,1) * (imag(z1) - real(z2))
f1 =3D fn1(f1)
f2 =3D fn2(f2)
z1 =3D (real(f1) + real(f2) + (0,1) * (imag(f1) + imag(f2))) / 2 + =
p2
z2 =3D (imag(f1) - imag(f2) + (0,1) * (real(f2) - real(f1))) / 2 + =
p1
z =3D sqrt(|z1| + |z2|)
ip =3D imag(p3)
rz =3D real(z1)
iz =3D imag(z1)
jz =3D real(z2)
kz =3D imag(z2)
IF(ip<0.1)
ct =3D z
ELSEIF(ip<1.1)
ct =3D z1
ELSEIF(ip<2.1)
ct =3D z2
ELSEIF(ip<3.1)
ct =3D rz
ELSEIF(ip<4.1)
ct =3D iz
ELSEIF(ip<5.1)
ct =3D jz
ELSEIF(ip<6.1)
ct =3D kz
ELSEIF((ip<7.1) && (|rz|>=3D|iz|) && (|rz|>=3D|jz|) && =
(|rz|>=3D|kz|))
ct =3D rz
ELSEIF((ip<7.1) && (|iz|>=3D|rz|) && (|iz|>=3D|jz|) && =
(|iz|>=3D|kz|))
ct =3D iz
ELSEIF((ip<7.1) && (|jz|>=3D|rz|) && (|jz|>=3D|iz|) && =
(|jz|>=3D|kz|))
ct =3D jz
ELSEIF((ip<7.1) && (|kz|>=3D|rz|) && (|iz|>=3D|kz|) && =
(|kz|>=3D|jz|))
ct =3D kz
ELSEIF((ip<8.1) && (|rz|<|iz|) && (|rz|<|jz|) && (|rz|<|kz|))
ct =3D rz
ELSEIF((ip<8.1) && (|iz|<|rz|) && (|iz|<|jz|) && (|iz|<|kz|))
ct =3D iz
ELSEIF((ip<8.1) && (|jz|<|rz|) && (|jz|<|iz|) && (|jz|<|kz|))
ct =3D jz
ELSEIF((ip<8.1) && (|kz|<|rz|) && (|kz|<|iz|) && (|kz|<|jz|))
ct =3D kz
ELSEIF(ip<9.1)
ct =3D (abs(rz) + abs(iz) + abs(jz) + abs(kz))
ELSEIF(ip<10.1)
ct =3D (rz + iz + jz + kz)
ELSE
ct =3D z
ENDIF,
|ct| <=3D |real(p3)|
}
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D
The Kwisatz Haderach,
=DF Gregory J. McClure
- --------------------------------------------------------------
Thanks for using Fractint, The Fractals and Fractint Discussion List
Post Message: fractint@lists.xmission.com
Get Commands: majordomo@lists.xmission.com "help"
Administrator: twegner@swbell.net
Unsubscribe: majordomo@lists.xmission.com "unsubscribe fractint"
------------------------------
End of fractint-digest V1 #440
******************************