Apophysis-AV/Variations/varInversion3D.pas
2022-03-08 20:25:51 +03:00

169 lines
4.0 KiB
ObjectPascal

{ Apophysis AV "Phoenix Edition" Copyright (C) 2021 Alice V. Koryagina }
unit varInversion3D;
interface
uses
BaseVariation, XFormMan;
const
sr = 'inversion3D_radius';
sx0 = 'inversion3D_x0';
sy0 = 'inversion3D_y0';
sz0 = 'inversion3D_z0';
type
TVariationInversion3D = class(TBaseVariation)
private
r, r2, x0, y0, z0: double;
public
constructor Create;
class function GetName: string; override;
class function GetInstance: TBaseVariation; override;
function GetNrVariables: integer; override;
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;
function ResetVariable(const Name: string): boolean; override;
procedure Prepare; override;
procedure CalcFunction; override;
end;
implementation
uses
Math;
{ TVariationInversion3D }
//////////////////////////////////////////
procedure TVariationInversion3D.Prepare;
begin
r2 := sqr(r);
end;
//////////////////////////////////////////
procedure TVariationInversion3D.CalcFunction;
var x, y, z, dist: double;
begin
x := FTx^ - x0;
y := FTy^ - y0;
z := FTz^ - z0;
dist := sqr(x) + sqr(y) + sqr(z) + 1E-30;
FPx^ := FPx^ + vvar * (r2 * x / dist + x0);
FPy^ := FPy^ + vvar * (r2 * y / dist + y0);
FPz^ := FPz^ + vvar * (r2 * z / dist + z0);
end;
///////////////////////////////////////////////////////////////////////////////
constructor TVariationInversion3D.Create;
begin
r := 1;
x0 := 0;
y0 := 0;
z0 := 0;
end;
///////////////////////////////////////////////////////////////////////////////
class function TVariationInversion3D.GetInstance: TBaseVariation;
begin
Result := TVariationInversion3D.Create;
end;
///////////////////////////////////////////////////////////////////////////////
class function TVariationInversion3D.GetName: string;
begin
Result := 'inversion3D';
end;
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
function TVariationInversion3D.GetVariableNameAt(const Index: integer): string;
begin
case Index of
0: Result := sr;
1: Result := sx0;
2: Result := sy0;
3: Result := sz0;
else
Result := '';
end
end;
///////////////////////////////////////////////////////////////////////////////
function TVariationInversion3D.SetVariable(const Name: string; var value: double): boolean;
begin
Result := False;
if Name = sr then begin
if (Value < 1E-6) then Value := 1E-6;
r := value;
Result := True;
end else if Name = sx0 then begin
x0 := Value;
Result := True;
end else if Name = sy0 then begin
y0 := Value;
Result := True;
end else if Name = sz0 then begin
z0 := Value;
Result := True;
end;
end;
function TVariationInversion3D.ResetVariable(const Name: string): boolean;
begin
Result := False;
if Name = sr then begin
r := 1;
Result := True;
end else if Name = sx0 then begin
x0 := 0;
Result := True;
end else if Name = sy0 then begin
y0 := 0;
Result := True;
end else if Name = sz0 then begin
z0 := 0;
Result := True;
end;
end;
/////////////////////////////////////////////////////////////////////
function TVariationInversion3D.GetNrVariables: integer;
begin
Result := 4;
end;
///////////////////////////////////////////////////////////////////////////////
function TVariationInversion3D.GetVariable(const Name: string; var value: double): boolean;
begin
Result := False;
if Name = sr then begin
Value := r;
Result := True;
end else if Name = sx0 then begin
Value := x0;
Result := True;
end else if Name = sy0 then begin
Value := y0;
Result := True;
end else if Name = sz0 then begin
Value := z0;
Result := True;
end;
end;
///////////////////////////////////////////////////////////////////////////////
initialization
RegisterVariation(TVariationClassLoader.Create(TVariationInversion3D), true, false);
end.