1 /**
2 * Mandelbulber v2, a 3D fractal generator _%}}i*<. ______
3 * Copyright (C) 2020 Mandelbulber Team _>]|=||i=i<, / ____/ __ __
4 * \><||i|=>>%) / / __/ /___/ /_
5 * This file is part of Mandelbulber. )<=i=]=|=i<> / /__ /_ __/_ __/
6 * The project is licensed under GPLv3, -<>>=|><|||` \____/ /_/ /_/
7 * see also COPYING file in this folder. ~+{i%+++
8 *
9 * Modified Mandelbox (ABox) formula
10 * This formula contains aux.color and aux.actualScaleA
11 */
12
13 #include "all_fractal_definitions.h"
14
cFractalFoldBoxMod1()15 cFractalFoldBoxMod1::cFractalFoldBoxMod1() : cAbstractFractal()
16 {
17 nameInComboBox = "Fold Box - Mod 1";
18 internalName = "fold_box_mod1";
19 internalID = fractal::foldBoxMod1;
20 DEType = analyticDEType;
21 DEFunctionType = linearDEFunction;
22 cpixelAddition = cpixelEnabledByDefault;
23 defaultBailout = 100.0;
24 DEAnalyticFunction = analyticFunctionLinear;
25 coloringFunction = coloringFunctionABox;
26 }
27
FormulaCode(CVector4 & z,const sFractal * fractal,sExtendedAux & aux)28 void cFractalFoldBoxMod1::FormulaCode(CVector4 &z, const sFractal *fractal, sExtendedAux &aux)
29 {
30 double colorAdd = 0.0;
31 if (aux.i >= fractal->transformCommon.startIterations
32 && aux.i < fractal->transformCommon.stopIterations)
33 {
34 CVector4 tempA = CVector4(0.0, 0.0, 0.0, 0.0);
35 CVector4 tempB = CVector4(0.0, 0.0, 0.0, 0.0);
36
37 CVector4 oldZ = z;
38 if (fractal->transformCommon.functionEnabledx)
39 tempA.x = fabs(z.x + fractal->transformCommon.additionConstant111.x);
40 if (fractal->transformCommon.functionEnabledAx)
41 tempB.x = fabs(z.x - fractal->transformCommon.additionConstantA111.x);
42 z.x = tempA.x - tempB.x - (z.x * fractal->transformCommon.scale3D111.x);
43
44 if (fractal->transformCommon.functionEnabledy)
45 tempA.y = fabs(z.y + fractal->transformCommon.additionConstant111.y);
46 if (fractal->transformCommon.functionEnabledAy)
47 tempB.y = fabs(z.y - fractal->transformCommon.additionConstantA111.y);
48 z.y = tempA.y - tempB.y - (z.y * fractal->transformCommon.scale3D111.y);
49
50 if (fractal->transformCommon.functionEnabledz)
51 tempA.z = fabs(z.z + fractal->transformCommon.additionConstant111.z);
52 if (fractal->transformCommon.functionEnabledAz)
53 tempB.z = fabs(z.z - fractal->transformCommon.additionConstantA111.z);
54 z.z = tempA.z - tempB.z - (z.z * fractal->transformCommon.scale3D111.z);
55 if (z.x != oldZ.x) colorAdd += fractal->mandelbox.color.factor.x;
56 if (z.y != oldZ.y) colorAdd += fractal->mandelbox.color.factor.y;
57 if (z.z != oldZ.z) colorAdd += fractal->mandelbox.color.factor.z;
58 }
59
60 if (fractal->transformCommon.functionEnabledFalse
61 && aux.i >= fractal->transformCommon.startIterationsA
62 && aux.i < fractal->transformCommon.stopIterationsA) // box fold
63 {
64 if (fabs(z.x) > fractal->mandelbox.foldingLimit)
65 {
66 z.x = sign(z.x) * fractal->mandelbox.foldingValue - z.x;
67 colorAdd += fractal->mandelbox.color.factor.x;
68 }
69 if (fabs(z.y) > fractal->mandelbox.foldingLimit)
70 {
71 z.y = sign(z.y) * fractal->mandelbox.foldingValue - z.y;
72 colorAdd += fractal->mandelbox.color.factor.y;
73 }
74 double zLimit = fractal->mandelbox.foldingLimit * fractal->transformCommon.scale1;
75 double zValue = fractal->mandelbox.foldingValue * fractal->transformCommon.scale1;
76 if (fabs(z.z) > zLimit)
77 {
78 z.z = sign(z.z) * zValue - z.z;
79 colorAdd += fractal->mandelbox.color.factor.z;
80 }
81 }
82
83 if (aux.i >= fractal->transformCommon.startIterationsB
84 && aux.i < fractal->transformCommon.stopIterationsB) // spherical fold
85 {
86 double r2 = z.Dot(z);
87 // if (r2 < 1e-21 && r2 > -1e-21)
88 // r2 = (r2 > 0) ? 1e-21 : -1e-21;
89 if (r2 < fractal->mandelbox.mR2)
90 {
91 z *= fractal->mandelbox.mboxFactor1;
92 aux.DE *= fractal->mandelbox.mboxFactor1;
93 colorAdd += fractal->mandelbox.color.factorSp1;
94 }
95 else if (r2 < fractal->mandelbox.fR2)
96 {
97 double tglad_factor2 = fractal->mandelbox.fR2 / r2;
98 z *= tglad_factor2;
99 aux.DE *= tglad_factor2;
100 colorAdd += fractal->mandelbox.color.factorSp2;
101 }
102 }
103
104 // scale
105 if (fractal->transformCommon.functionEnabledXFalse
106 && aux.i >= fractal->transformCommon.startIterationsD
107 && aux.i < fractal->transformCommon.stopIterationsD)
108 {
109 double useScale = aux.actualScaleA + fractal->mandelbox.scale;
110 z *= useScale;
111 aux.DE = aux.DE * fabs(useScale) + 1.0;
112
113 // update actualScale for next iteration
114 double vary = fractal->transformCommon.scaleVary0
115 * (fabs(aux.actualScaleA) - fractal->transformCommon.scaleB1);
116 if (fractal->transformCommon.functionEnabledMFalse)
117 aux.actualScaleA = -vary;
118 else
119 aux.actualScaleA = aux.actualScaleA - vary;
120 }
121 else
122 {
123 z *= fractal->mandelbox.scale;
124 aux.DE = aux.DE * fabs(fractal->mandelbox.scale) + 1.0;
125 }
126 if (fractal->mandelbox.mainRotationEnabled && aux.i >= fractal->transformCommon.startIterationsC
127 && aux.i < fractal->transformCommon.stopIterationsC)
128 z = fractal->mandelbox.mainRot.RotateVector(z);
129
130 // color updated v2.14
131 if (fractal->foldColor.auxColorEnabled)
132 {
133 aux.color += colorAdd;
134 }
135 }
136