.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "examples/batch_mode/17-perpendicular_burning_ship_DEM.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_17-perpendicular_burning_ship_DEM.py: ================================= 17 - "Perpendicular" Burning Ship ================================= This variant of the Burning Ship is interesting for the variety of its shapes. Here some decoration around a deep mini, at a depth of 1.8e-128 Reference: `fractalshades.models.Perturbation_perpendicular_burning_ship` .. GENERATED FROM PYTHON SOURCE LINES 13-182 .. image-sg:: /examples/batch_mode/images/sphx_glr_17-perpendicular_burning_ship_DEM_001.png :alt: 17 perpendicular burning ship DEM :srcset: /examples/batch_mode/images/sphx_glr_17-perpendicular_burning_ship_DEM_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 = '0.99933368699185373639022706206725770981715223678700982012404710427171219409317906966176602615220191953222863056569706047389033861855088978962654736627779241722367794435564687067' y = '0.996457170505254265452116588481786614441676204699082859164552562590920318860452552741252698252881096523313551079924623070901888949263620321040269806468649035399547439554888846827' dx = '1.89491419e-128' xy_ratio = 1.8 theta_deg = -35 dps = 140 nx = 2400 # _1b = 'Skew parameters /!\\ Re-run when modified!' has_skew = True skew_00 = 0.9393449132343535 skew_01 = 0.34255954266629574 skew_10 = 0.6121690906221282 skew_11 = 1.287817016597985 # _2 = 'Calculation parameters' max_iter = 500000 # _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.cmap_register["classic"] invert_cmap = False DEM_min = 1e-04 zmin = -9.21034049987793 zmax = 0.4 * (-9.21034049987793) + 0.6 * (-4.388575077056885) # _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=60, 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:** ( 2 minutes 59.714 seconds) .. _sphx_glr_download_examples_batch_mode_17-perpendicular_burning_ship_DEM.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: 17-perpendicular_burning_ship_DEM.py <17-perpendicular_burning_ship_DEM.py>` .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: 17-perpendicular_burning_ship_DEM.ipynb <17-perpendicular_burning_ship_DEM.ipynb>` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_