new julia/kaleidoscope variation
This commit is contained in:
zueuk
parent
789c0ae596
commit
489f4fa1be
416
2.10/Source/varJuliaScope.pas
Normal file
416
2.10/Source/varJuliaScope.pas
Normal file
@ -0,0 +1,416 @@
|
||||
unit varJuliaScope;
|
||||
|
||||
interface
|
||||
|
||||
uses
|
||||
BaseVariation, XFormMan;
|
||||
|
||||
const
|
||||
variation_name='juliascope';
|
||||
var_n_name='juliascope_power';
|
||||
var_c_name='juliascope_dist';
|
||||
|
||||
{$define _ASM_}
|
||||
|
||||
type
|
||||
TVariationJuliaScope = class(TBaseVariation)
|
||||
private
|
||||
N: integer;
|
||||
c: double;
|
||||
|
||||
rN: integer;
|
||||
cn: double;
|
||||
|
||||
procedure CalcPower1;
|
||||
procedure CalcPowerMinus1;
|
||||
procedure CalcPower2;
|
||||
procedure CalcPowerMinus2;
|
||||
|
||||
public
|
||||
constructor Create;
|
||||
|
||||
class function GetName: string; override;
|
||||
class function GetInstance: TBaseVariation; override;
|
||||
|
||||
class function GetNrVariables: integer; override;
|
||||
class function GetVariableNameAt(const Index: integer): string; override;
|
||||
|
||||
function SetVariable(const Name: string; var value: double): boolean; override;
|
||||
function GetVariable(const Name: string; var value: double): boolean; override;
|
||||
|
||||
procedure Prepare; override;
|
||||
procedure CalcFunction; override;
|
||||
procedure GetCalcFunction(var f: TCalcFunction); override;
|
||||
end;
|
||||
|
||||
implementation
|
||||
|
||||
uses
|
||||
math;
|
||||
|
||||
// TVariationJuliaScope
|
||||
|
||||
///////////////////////////////////////////////////////////////////////////////
|
||||
constructor TVariationJuliaScope.Create;
|
||||
begin
|
||||
N := random(5) + 2;
|
||||
c := 1.0;
|
||||
end;
|
||||
|
||||
procedure TVariationJuliaScope.Prepare;
|
||||
begin
|
||||
rN := abs(N);
|
||||
cn := c / N / 2;
|
||||
end;
|
||||
|
||||
///////////////////////////////////////////////////////////////////////////////
|
||||
procedure TVariationJuliaScope.GetCalcFunction(var f: TCalcFunction);
|
||||
begin
|
||||
if cn = 1/4 then f := CalcPower2
|
||||
else if cn = -1/4 then f := CalcPowerMinus2
|
||||
else if cn = 1/2 then f := CalcPower1
|
||||
else if cn = -1/2 then f := CalcPowerMinus1
|
||||
else f := CalcFunction;
|
||||
end;
|
||||
|
||||
///////////////////////////////////////////////////////////////////////////////
|
||||
procedure TVariationJuliaScope.CalcFunction;
|
||||
{$ifndef _ASM_}
|
||||
var
|
||||
rnd: integer;
|
||||
r: double;
|
||||
sina, cosa: extended;
|
||||
begin
|
||||
rnd := random(rN);
|
||||
if (rnd and 1) = 0 then
|
||||
sincos( (2*pi*rnd + arctan2(FTy^, FTx^)) / N, sina, cosa)
|
||||
else
|
||||
sincos( (2*pi*rnd - arctan2(FTy^, FTx^)) / N, sina, cosa);
|
||||
r := vvar * Math.Power(sqr(FTx^) + sqr(FTy^), cn);
|
||||
FPx^ := FPx^ + r * cosa;
|
||||
FPy^ := FPy^ + r * sina;
|
||||
{$else}
|
||||
asm
|
||||
mov edx, [eax + FTy]
|
||||
fld qword ptr [edx]
|
||||
fld qword ptr [eax + cn]
|
||||
mov edx, [eax + FTx]
|
||||
fld qword ptr [edx]
|
||||
fld st(2)
|
||||
fld st(1)
|
||||
fpatan
|
||||
mov ecx, eax
|
||||
mov eax, dword ptr [eax + rN]
|
||||
call System.@RandInt
|
||||
push eax
|
||||
|
||||
shr eax, 1
|
||||
jnc @even
|
||||
fchs
|
||||
@even:
|
||||
|
||||
fldpi
|
||||
fadd st, st
|
||||
fimul dword ptr [esp]
|
||||
add esp, 4
|
||||
faddp
|
||||
fidiv dword ptr [ecx + N]
|
||||
|
||||
fxch st(3)
|
||||
fmul st, st
|
||||
fxch st(1)
|
||||
fmul st, st
|
||||
faddp
|
||||
// --- x^y = 2^(y*log2(x))
|
||||
fyl2x
|
||||
fld st
|
||||
frndint
|
||||
fsub st(1), st
|
||||
fxch st(1)
|
||||
f2xm1
|
||||
fld1
|
||||
fadd
|
||||
fscale
|
||||
fstp st(1)
|
||||
// ---
|
||||
fmul qword ptr [ecx + vvar]
|
||||
fxch st(1)
|
||||
fsincos
|
||||
fmul st, st(2)
|
||||
|
||||
mov edx, [ecx + FPx]
|
||||
fadd qword ptr [edx]
|
||||
fstp qword ptr [edx]
|
||||
fmulp
|
||||
mov edx, [ecx + FPy]
|
||||
fadd qword ptr [edx]
|
||||
fstp qword ptr [edx]
|
||||
fwait
|
||||
{$endif}
|
||||
end;
|
||||
|
||||
procedure TVariationJuliaScope.CalcPower2;
|
||||
{$ifndef _ASM_}
|
||||
var
|
||||
r: double;
|
||||
sina, cosa: extended;
|
||||
begin
|
||||
if random(2) = 0 then
|
||||
sincos(arctan2(FTy^, FTx^)/2, sina, cosa)
|
||||
else
|
||||
sincos(pi - arctan2(FTy^, FTx^)/2, sina, cosa);
|
||||
|
||||
r := vvar * sqrt(sqrt(sqr(FTx^) + sqr(FTy^)));
|
||||
|
||||
FPx^ := FPx^ + r * cosa;
|
||||
FPy^ := FPy^ + r * sina;
|
||||
{$else}
|
||||
asm
|
||||
mov edx, [eax + FTy]
|
||||
fld qword ptr [edx]
|
||||
mov edx, [eax + FTx]
|
||||
fld qword ptr [edx]
|
||||
fld st(1)
|
||||
fld st(1)
|
||||
fpatan
|
||||
fld1
|
||||
fadd st, st
|
||||
fdivp st(1), st
|
||||
mov ecx, eax
|
||||
//mov eax, 2
|
||||
call System.@RandInt
|
||||
|
||||
shr eax, 1
|
||||
jnc @skip
|
||||
fldpi
|
||||
fsubrp st(1), st
|
||||
@skip:
|
||||
|
||||
{
|
||||
push eax
|
||||
|
||||
shr eax, 1
|
||||
jnc @even
|
||||
fchs
|
||||
@even:
|
||||
|
||||
fldpi
|
||||
fimul dword ptr [esp]
|
||||
add esp, 4
|
||||
faddp
|
||||
}
|
||||
fxch st(2)
|
||||
fmul st, st
|
||||
fxch st(1)
|
||||
fmul st, st
|
||||
faddp
|
||||
fsqrt
|
||||
fsqrt
|
||||
fmul qword ptr [ecx + vvar]
|
||||
fxch st(1)
|
||||
|
||||
fsincos
|
||||
|
||||
fmul st, st(2)
|
||||
mov edx, [ecx + FPx]
|
||||
fadd qword ptr [edx]
|
||||
fstp qword ptr [edx]
|
||||
fmulp
|
||||
mov edx, [ecx + FPy]
|
||||
fadd qword ptr [edx]
|
||||
fstp qword ptr [edx]
|
||||
fwait
|
||||
{$endif}
|
||||
end;
|
||||
|
||||
procedure TVariationJuliaScope.CalcPowerMinus2;
|
||||
{$ifndef _ASM_}
|
||||
var
|
||||
r: double;
|
||||
sina, cosa: extended;
|
||||
begin
|
||||
if random(2) = 0 then
|
||||
sincos(arctan2(FTy^, FTx^)/2, sina, cosa)
|
||||
else
|
||||
sincos(pi - arctan2(FTy^, FTx^)/2, sina, cosa);
|
||||
r := vvar / sqrt(sqrt(sqr(FTx^) + sqr(FTy^)));
|
||||
|
||||
FPx^ := FPx^ + r * cosa;
|
||||
FPy^ := FPy^ - r * sina;
|
||||
{$else}
|
||||
asm
|
||||
mov edx, [eax + FTy]
|
||||
fld qword ptr [edx]
|
||||
mov edx, [eax + FTx]
|
||||
fld qword ptr [edx]
|
||||
fld st(1)
|
||||
fld st(1)
|
||||
fpatan
|
||||
fld1
|
||||
fadd st, st
|
||||
fdivp st(1), st
|
||||
mov ecx, eax
|
||||
mov eax, 2
|
||||
call System.@RandInt
|
||||
|
||||
shr eax, 1
|
||||
jnc @skip
|
||||
fldpi
|
||||
fsubrp st(1), st
|
||||
@skip:
|
||||
|
||||
{ push eax
|
||||
|
||||
shr eax, 1
|
||||
jnc @even
|
||||
fchs
|
||||
@even:
|
||||
|
||||
fldpi
|
||||
fimul dword ptr [esp]
|
||||
add esp, 4
|
||||
faddp
|
||||
}
|
||||
fxch st(2)
|
||||
fmul st, st
|
||||
fxch st(1)
|
||||
fmul st, st
|
||||
faddp
|
||||
fsqrt
|
||||
fsqrt
|
||||
fdivr qword ptr [ecx + vvar]
|
||||
fxch st(1)
|
||||
|
||||
fsincos
|
||||
|
||||
fmul st, st(2)
|
||||
mov edx, [ecx + FPx]
|
||||
fadd qword ptr [edx]
|
||||
fstp qword ptr [edx]
|
||||
fmulp
|
||||
mov edx, [ecx + FPy]
|
||||
fsubr qword ptr [edx]
|
||||
fstp qword ptr [edx]
|
||||
fwait
|
||||
{$endif}
|
||||
end;
|
||||
|
||||
procedure TVariationJuliaScope.CalcPower1;
|
||||
{$ifndef _ASM_}
|
||||
begin
|
||||
FPx^ := FPx^ + vvar * FTx^;
|
||||
FPy^ := FPy^ + vvar * FTy^;
|
||||
{$else}
|
||||
asm
|
||||
mov edx, [eax + FTy]
|
||||
fld qword ptr [edx]
|
||||
mov edx, [eax + FTx]
|
||||
fld qword ptr [edx]
|
||||
fld qword ptr [eax + vvar]
|
||||
fmul st(2), st
|
||||
fmulp
|
||||
mov edx, [eax + FPx]
|
||||
fadd qword ptr [edx]
|
||||
fstp qword ptr [edx]
|
||||
mov edx, [eax + FPy]
|
||||
fadd qword ptr [edx]
|
||||
fstp qword ptr [edx]
|
||||
fwait
|
||||
{$endif}
|
||||
end;
|
||||
|
||||
procedure TVariationJuliaScope.CalcPowerMinus1;
|
||||
{$ifndef _ASM_}
|
||||
var
|
||||
r: double;
|
||||
begin
|
||||
r := vvar / (sqr(FTx^) + sqr(FTy^));
|
||||
|
||||
FPx^ := FPx^ + r * FTx^;
|
||||
FPy^ := FPy^ - r * FTy^;
|
||||
{$else}
|
||||
asm
|
||||
mov edx, [eax + FTy]
|
||||
fld qword ptr [edx]
|
||||
mov edx, [eax + FTx]
|
||||
fld qword ptr [edx]
|
||||
fld st(1)
|
||||
fmul st, st
|
||||
fld st(1)
|
||||
fmul st, st
|
||||
faddp
|
||||
fdivr qword ptr [eax + vvar]
|
||||
fmul st(2), st
|
||||
fmulp
|
||||
mov edx, [eax + FPx]
|
||||
fadd qword ptr [edx]
|
||||
fstp qword ptr [edx]
|
||||
mov edx, [eax + FPy]
|
||||
fsubr qword ptr [edx]
|
||||
fstp qword ptr [edx]
|
||||
fwait
|
||||
{$endif}
|
||||
end;
|
||||
|
||||
///////////////////////////////////////////////////////////////////////////////
|
||||
class function TVariationJuliaScope.GetInstance: TBaseVariation;
|
||||
begin
|
||||
Result := TVariationJuliaScope.Create;
|
||||
end;
|
||||
|
||||
///////////////////////////////////////////////////////////////////////////////
|
||||
class function TVariationJuliaScope.GetName: string;
|
||||
begin
|
||||
Result := variation_name;
|
||||
end;
|
||||
|
||||
///////////////////////////////////////////////////////////////////////////////
|
||||
class function TVariationJuliaScope.GetVariableNameAt(const Index: integer): string;
|
||||
begin
|
||||
case Index of
|
||||
0: Result := var_n_name;
|
||||
1: Result := var_c_name;
|
||||
else
|
||||
Result := '';
|
||||
end
|
||||
end;
|
||||
|
||||
///////////////////////////////////////////////////////////////////////////////
|
||||
function TVariationJuliaScope.SetVariable(const Name: string; var value: double): boolean;
|
||||
begin
|
||||
Result := False;
|
||||
if Name = var_n_name then begin
|
||||
N := Round(Value);
|
||||
if N = 0 then N := 1;
|
||||
Value := N;
|
||||
Result := True;
|
||||
end
|
||||
else if Name = var_c_name then begin
|
||||
c := value;
|
||||
end;
|
||||
end;
|
||||
|
||||
///////////////////////////////////////////////////////////////////////////////
|
||||
class function TVariationJuliaScope.GetNrVariables: integer;
|
||||
begin
|
||||
Result := 2;
|
||||
end;
|
||||
|
||||
///////////////////////////////////////////////////////////////////////////////
|
||||
function TVariationJuliaScope.GetVariable(const Name: string; var value: double): boolean;
|
||||
begin
|
||||
Result := False;
|
||||
if Name = var_n_name then begin
|
||||
Value := N;
|
||||
Result := true;
|
||||
end
|
||||
else if Name = var_c_name then begin
|
||||
Value := c;
|
||||
Result := true;
|
||||
end;
|
||||
end;
|
||||
|
||||
///////////////////////////////////////////////////////////////////////////////
|
||||
initialization
|
||||
RegisterVariation(TVariationJuliaScope);
|
||||
end.
|
||||
Loading…
Reference in New Issue
Block a user