.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "examples/batch_mode/18-perpendicular_burning_ship_glynn.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note Click :ref:`here ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_examples_batch_mode_18-perpendicular_burning_ship_glynn.py: ====================================================== 18 - "Perpendicular" Burning Ship: hidden Glynn spiral ====================================================== An example of hidden feature in this fractal: a Glynn spiral in a very skewed area at a depth of 4.e-40. Reference: `fractalshades.models.Perturbation_perpendicular_burning_ship` .. GENERATED FROM PYTHON SOURCE LINES 13-193 .. image-sg:: /examples/batch_mode/images/sphx_glr_18-perpendicular_burning_ship_glynn_001.png :alt: 18 perpendicular burning ship glynn :srcset: /examples/batch_mode/images/sphx_glr_18-perpendicular_burning_ship_glynn_001.png :class: sphx-glr-single-img .. code-block:: default import os import numpy as np import fractalshades as fs import fractalshades.models as fsm import fractalshades.colors as fscolors from fractalshades.postproc import ( Postproc_batch, Continuous_iter_pp, DEM_normal_pp, DEM_pp, Raw_pp, ) from fractalshades.colors.layers import ( Color_layer, Bool_layer, Normal_map_layer, Virtual_layer, Blinn_lighting, ) def plot(plot_dir): fs.settings.enable_multithreading = True fs.settings.inspect_calc = True # A simple showcase using perturbation technique calc_name = 'test' # _1 = 'Zoom parameters' x = '-1.6221172452279831275586824847368230989301274844265' y = '-0.0043849065564689427951877101597546609652950526531633' dx = '4.646303299697506e-40' xy_ratio = 1.8 theta_deg = -2 dps = 55 nx = 2400 # _1b = 'Skew parameters /!\\ Re-run when modified!' has_skew = True skew_00 = 1.011753723519244 skew_01 = -1.157539989768796 skew_10 = -0.5299787188179303 skew_11 = 1.5947275737676074 # _2 = 'Calculation parameters' max_iter = 20000 # _3 = 'Bilinear series parameters' eps = 1e-06 # _4 = 'Plotting parameters: base field' base_layer = 'distance_estimation' interior_color = (0.6627451181411743, 0.4313725531101227, 0.0) colormap = fscolors.Fractal_colormap( colors=[ [1. , 1. , 0. ], [0.05098039, 0.03921569, 0.3137255 ], [0.10588235, 0.78039217, 0.78039217], [0.33333334, 1. , 1. ] ], kinds=['Lch', 'Lch', 'Lch', 'Lch'], grad_npts=[ 3, 3, 32, 32], grad_funcs=['x', 'x**2', 'x', 'x'], extent='mirror' ) invert_cmap = True DEM_min = 1e-08 zmin = 6.540582180023193 zmax = 18.42068099975586 * 1.1 - 6.540582180023193 * 0.1 # _5 = 'Plotting parameters: shading' shade_kind = 'glossy' gloss_intensity = 100.0 light_angle_deg = 35.0 light_color = (1.0, 1.0, 1.0) gloss_light_color = (1.0, 1.0, 1.0) # Run the calculation fractal = fsm.Perturbation_burning_ship( plot_dir, flavor="Perpendicular burning ship" ) fractal.zoom(precision=dps, x=x, y=y, dx=dx, nx=nx, xy_ratio=xy_ratio, theta_deg=theta_deg, projection="cartesian", has_skew=has_skew, skew_00=skew_00, skew_01=skew_01, skew_10=skew_10, skew_11=skew_11 ) fractal.calc_std_div( calc_name=calc_name, subset=None, max_iter=max_iter, M_divergence=1.e3, BLA_eps=eps, ) pp = Postproc_batch(fractal, calc_name) if base_layer == "continuous_iter": pp.add_postproc(base_layer, Continuous_iter_pp()) elif base_layer == "distance_estimation": pp.add_postproc("continuous_iter", Continuous_iter_pp()) pp.add_postproc(base_layer, DEM_pp()) pp.add_postproc("interior", Raw_pp("stop_reason", func=lambda x: x != 1)) if shade_kind != "None": pp.add_postproc("DEM_map", DEM_normal_pp(kind="potential")) plotter = fs.Fractal_plotter(pp) plotter.add_layer(Bool_layer("interior", output=False)) if shade_kind != "None": plotter.add_layer(Normal_map_layer( "DEM_map", max_slope=40, output=False )) if base_layer != 'continuous_iter': plotter.add_layer( Virtual_layer("continuous_iter", func=None, output=False) ) sign = {False: 1., True: -1.}[invert_cmap] if base_layer == 'distance_estimation': cmap_func = lambda x: sign * np.where( np.isinf(x), np.log(DEM_min), np.log(np.clip(x, DEM_min, None)) ) else: cmap_func = lambda x: sign * np.log(x) plotter.add_layer(Color_layer( base_layer, func=cmap_func, colormap=colormap, probes_z=[zmin, zmax], output=True)) plotter[base_layer].set_mask( plotter["interior"], mask_color=interior_color ) if shade_kind != "None": light = Blinn_lighting(0.6, np.array([1., 1., 1.])) light.add_light_source( k_diffuse=0.8, k_specular=.0, shininess=350., polar_angle=light_angle_deg, azimuth_angle=10., color=np.array(light_color)) if shade_kind == "glossy": light.add_light_source( k_diffuse=0.2, k_specular=gloss_intensity, shininess=400., polar_angle=light_angle_deg, azimuth_angle=10., color=np.array(gloss_light_color)) plotter[base_layer].shade(plotter["DEM_map"], light) plotter.plot() if __name__ == "__main__": # Some magic to get the directory for plotting: with a name that matches # the file or a temporary dir if we are building the documentation try: realpath = os.path.realpath(__file__) plot_dir = os.path.splitext(realpath)[0] plot(plot_dir) except NameError: import tempfile with tempfile.TemporaryDirectory() as plot_dir: fs.utils.exec_no_output(plot, plot_dir) .. rst-class:: sphx-glr-timing **Total running time of the script:** ( 0 minutes 51.564 seconds) .. _sphx_glr_download_examples_batch_mode_18-perpendicular_burning_ship_glynn.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: 18-perpendicular_burning_ship_glynn.py <18-perpendicular_burning_ship_glynn.py>` .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: 18-perpendicular_burning_ship_glynn.ipynb <18-perpendicular_burning_ship_glynn.ipynb>` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_