31 	if (fractal->IFS.absX) z.x = fabs(z.x);  in FormulaCode()
32 	if (fractal->IFS.absY) z.y = fabs(z.y);  in FormulaCode()
33 	if (fractal->IFS.absZ) z.z = fabs(z.z);  in FormulaCode()
37 		if (fractal->IFS.enabled[i])  in FormulaCode()
39 			z = fractal->IFS.rot[i].RotateVector(z);  in FormulaCode()
40 			double length = z.Dot(fractal->IFS.direction[i]);  in FormulaCode()
42 			if (length < fractal->IFS.distance[i])  in FormulaCode()
44 				z -= fractal->IFS.direction[i]  in FormulaCode()
45 						 * (2.0 * (length - fractal->IFS.distance[i]) * fractal->IFS.intensity[i]);  in FormulaCode()
50 	if (fractal->IFS.rotationEnabled)  in FormulaCode()
51 		z = fractal->IFS.mainRot.RotateVector(z - fractal->IFS.offset) + fractal->IFS.offset;  in FormulaCode()
53 	if (fractal->IFS.edgeEnabled)  in FormulaCode()
55 		if (fractal->IFS.edge.x > 0.0) z.x = fractal->IFS.edge.x - fabs(fractal->IFS.edge.x - z.x);  in FormulaCode()
56 		if (fractal->IFS.edge.y > 0.0) z.y = fractal->IFS.edge.y - fabs(fractal->IFS.edge.y - z.y);  in FormulaCode()
57 		if (fractal->IFS.edge.z > 0.0) z.z = fractal->IFS.edge.z - fabs(fractal->IFS.edge.z - z.z);  in FormulaCode()
60 	z *= fractal->IFS.scale;  in FormulaCode()
61 	if (fractal->IFS.mengerSpongeMode)  in FormulaCode()
63 		z.x -= fractal->IFS.offset.x * (fractal->IFS.scale - 1.0);  in FormulaCode()
64 		z.y -= fractal->IFS.offset.y * (fractal->IFS.scale - 1.0);  in FormulaCode()
65 		if (z.z > 0.5 * fractal->IFS.offset.z * (fractal->IFS.scale - 1.0))  in FormulaCode()
66 			z.z -= fractal->IFS.offset.z * (fractal->IFS.scale - 1.0);  in FormulaCode()
70 		z -= fractal->IFS.offset * (fractal->IFS.scale - 1.0);  in FormulaCode()
72 	aux.DE *= fabs(fractal->IFS.scale);  in FormulaCode()