def plot(): """ Example plot of "Dinkydau flake" location, classic test case for perturbation techinque and glitch correction. """ directory = "./perturb3" # A simple showcas using perturbation technique x, y = "-0.75", "0." dx = "5.e0" precision = 150 nx = 3200 x = '-1.993101465346811633595616349779370960719700854263460744227885419412572351324445321874004442403503011725492641453441484329872421586639050539267910275781311721259641283288898121495962960511188360859446141741747448336162741635241659807342073817485606900204314068415477260531866235822220430486119843399542682686903115170284286744427789259769672374264750048282753697939441835784223761880144973743249562058785490789121881765822494487680713802365108655723804325265573559505557675274602687698535315326824126504568612493712586162172913902182849175957355201038749736221172166381630971780664574945186600702295814202276821096082583371646391752258136082934974808859285874633438821365019578751567557825904520349615083013659980914914027970434909527071583051834592117828848476162531653958895351112086988431145593103631584906268842017812275327407852982198464690863881669210429524886644612378383194829443312845510612712652609161824242961337428831821452985753354390486290596759804822656081435972288149493458837417345622327424314356121019642859410599076344584439053' y = '-0.0000000000000000000001023710512431589570005203798273240286470771650981763139330824251367527539704908652681677255846096018123604758184274291544039306788236186699886780634842327075201516856958710101981983627439650646801970216473720458108184168523830340517163870169919249638047451318365933557979031258177237717753794769207546631778925762244933996868229775839453423976148160460368950010744615992543567659336441278697767936754177402161193387335896316048958999044786173769537636937690325372184302039387809032377492478883386867909568081378343026465718279702082155350634045451662043279571130329807181205168575035218579446441743895034500494162054542465662499561872308796389299543485971550119965555406797513540955180015002773089512902159298036751555750537198269153512054314778106379999904973319311775913118876664658896447125394230182911598486994660593785101070940474199338173235671779027464278763584882649334715308568951689930703319182855450455829071681801598812078546949506703052822064075723064430117087174494686425795962484572598182905448268661576794867514394048126125843131' dx = '0.000000000000000000000000000000000000000000000000000001611311677165311992036964188953416400357306950103215815939935672793028078345666262795353854791128193417180028322502563016371152356567374462535929127041725333256467561733290384009069197316210881426462964928434585521650089769009585599250595741136222346428026621232359069818161715344902647536565294138377179149230247747809792678090166161270842118179873716803137049807896171897484336423552754428273143300120233210420317384023816622150052145902931298541283952210479589483643604612244327904646956347253652556190458708451251053597744280720211750536552494355053069572018017158283906744473170534917744052363527030185810701281751448157641033852978421572439410367749824884649311769407121646490987840830692921695850309271975840408355787675238627661695949205488539553336065242972661368264433343340961667491522236015169240810915606280026794402643082165872947330057146892300382601081535552273108633947461317572323488134065825603499355505743581891071053495091972381402991134813796575186487618387985500615976611029133681690931704603523636857668558756510416666667' # Set to True if you only want to rerun the post-processing part settings.skip_calc = False # Set to True to enable multi-processing settings.enable_multiprocessing = True # xy_ratio = 1.0 # theta_deg = 0. # complex_type = np.complex128 mandelbrot = fsm.Perturbation_mandelbrot(directory) mandelbrot.zoom(precision=precision, x=x, y=y, dx=dx, nx=nx, xy_ratio=1.0, theta_deg=0., projection="cartesian", antialiasing=True) mandelbrot.calc_std_div( complex_type=np.complex128, file_prefix="dev", subset=None, max_iter=50000, M_divergence=1.e3, epsilon_stationnary=1.e-3, pc_threshold=0.1, SA_params={ "cutdeg": 64, "cutdeg_glitch": 8, "SA_err": 1.e-4, "use_Taylor_shift": True }, glitch_eps=1.e-6, interior_detect=False, #True, glitch_max_attempt=20) mandelbrot.run() mask_codes = [0, 2, 3, 4] mask = fs.Fractal_Data_array( mandelbrot, file_prefix="dev", postproc_keys=('stop_reason', lambda x: np.isin(x, mask_codes)), mode="r+raw") # cv = fs.Fractal_Data_array(mandelbrot, file_prefix="dev", # postproc_keys=('stop_reason', lambda x: x != 2), mode="r+raw") # potential_data_key = ("potential", {}) # gold = np.array([255, 210, 66]) / 255. # black = np.array([0, 0, 0]) / 255. # color_gradient = fs.Color_tools.Lch_gradient(gold, black, 200) # colormap = fs.Fractal_colormap(color_gradient) gold = np.array([255, 210, 66]) / 255. black = np.array([0, 0, 0]) / 255. purple = np.array([181, 40, 99]) / 255. colors1 = np.vstack((purple[np.newaxis, :])) colors2 = np.vstack((gold[np.newaxis, :])) colormap = fscolors.Fractal_colormap(kinds="Lch", colors1=colors1, colors2=colors2, n=200, funcs=None, extent="mirror") plotter = fs.Fractal_plotter( fractal=mandelbrot, base_data_key=("potential", {}), # potential_data_key, #glitch_sort_key, base_data_prefix="dev", base_data_function=lambda x: x, # np.sin(x*0.0001), colormap=colormap, probes_val=[ 0., 0.55 ], # 200. + 200, #* 428 - 00.,#[0., 0.5, 1.], #phi * k * 2. + k * np.array([0., 1., 1.5, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0]) / 3.5, probes_kind="qt", #"z", "qt" mask=mask) #plotter.add_calculation_layer(postproc_key=potential_data_key) # layer1_key = ("DEM_shade", {"kind": "potential", # "theta_LS": 30., # "phi_LS": 50., # "shininess": 3., # "ratio_specular": 15000.}) # plotter.add_grey_layer(postproc_key=layer1_key, intensity=0.75, # blur_ranges=[],#[[0.99, 0.999, 1.0]], # disp_layer=False, #skewness=0.2, # normalized=False, hardness=0.35, # skewness=0.0, shade_type={"Lch": 1.0, "overlay": 1., "pegtop": 4.}) # layer2_key = ("field_lines", {}) # plotter.add_grey_layer(postproc_key=layer2_key, # hardness=1.0, intensity=0.68, skewness=0.0, ## blur_ranges=[[0.50, 0.60, 1.0]], # shade_type={"Lch": 0., "overlay": 2., "pegtop": 1.}) plotter.add_grey_layer(postproc_key=("DEM_shade", { "kind": "potential", "theta_LS": 30., "phi_LS": 50., "shininess": 30., "ratio_specular": 15000. }), blur_ranges=[], hardness=0.9, intensity=0.8, shade_type={ "Lch": 1.0, "overlay": 0., "pegtop": 4. }) plotter.plot("dev", mask_color=(0., 0., 0.))
def plot(): """ Example plot of "Dinkydau flake" location, classic test case for perturbation techinque and glitch correction. """ directory = "./perturb2" # A simple showcas using perturbation technique x, y = "-1.74920463345912691e+00", "-2.8684660237361114e-04" dx = "5e-12" precision = 16 nx = 800 # Set to True if you only want to rerun the post-processing part settings.skip_calc = False # Set to True to enable multi-processing settings.enable_multiprocessing = True # xy_ratio = 1.0 # theta_deg = 0. # complex_type = np.complex128 mandelbrot = fsm.Perturbation_mandelbrot(directory) mandelbrot.zoom(precision=precision, x=x, y=y, dx=dx, nx=nx, xy_ratio=1.0, theta_deg=0., projection="cartesian", antialiasing=False) mandelbrot.calc_std_div(complex_type=np.complex128, file_prefix="dev", subset=None, max_iter=50000, M_divergence=1.e3, epsilon_stationnary=1.e-3, pc_threshold=0.1, SA_params={ "cutdeg": 8, "cutdeg_glitch": 8, "SA_err": 1.e-4, "use_Taylor_shift": True }, glitch_eps=1.e-6, interior_detect=True, glitch_max_attempt=20) mandelbrot.run() cv = fs.Fractal_Data_array(mandelbrot, file_prefix="dev", postproc_keys=('stop_reason', lambda x: x == 1), mode="r+raw") potential_data_key = ("potential", {}) citrus2 = np.array([103, 189, 0]) / 255. citrus_white = np.array([252, 251, 226]) / 255. wheat1 = np.array([244, 235, 158]) / 255. wheat2 = np.array([246, 207, 106]) / 255. wheat3 = np.array([191, 156, 96]) / 255. lavender1 = np.array([154, 121, 144]) / 255. lavender2 = np.array([140, 94, 134]) / 255. lavender3 = np.array([18, 16, 58]) / 255. def wave(x): return 0.5 + (0.4 * (x - 0.5) - 0.6 * 0.5 * np.cos(x * np.pi * 3.)) colormap = fscolors.Fractal_colormap( kinds="Lch", colors1=np.vstack( (citrus_white, wheat2, wheat1, wheat2, wheat1, wheat2, wheat3, wheat1, lavender2, wheat1, wheat2, wheat3, wheat1, lavender2, wheat1, lavender3, lavender2, lavender3, lavender1, lavender3, lavender2)), colors2=None, n=100, funcs=lambda x: wave(x), extent="mirror") plotter = fs.Fractal_plotter( fractal=mandelbrot, base_data_key= potential_data_key, # potential_data_key, #glitch_sort_key, base_data_prefix="dev", base_data_function=lambda x: x, # np.sin(x*0.0001), colormap=colormap, probes_val=np.linspace(0., 1., 22) * 0.5, # 200. + 200, #* 428 - 00.,#[0., 0.5, 1.], #phi * k * 2. + k * np.array([0., 1., 1.5, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0]) / 3.5, probes_kind="qt", #"z", "qt" mask=~cv) #plotter.add_calculation_layer(postproc_key=potential_data_key) layer1_key = ("DEM_shade", { "kind": "potential", "theta_LS": 30., "phi_LS": 50., "shininess": 3., "ratio_specular": 15000. }) plotter.add_grey_layer( postproc_key=layer1_key, intensity=0.75, blur_ranges=[], #[[0.99, 0.999, 1.0]], disp_layer=False, #skewness=0.2, normalized=False, hardness=0.35, skewness=0.0, shade_type={ "Lch": 1.0, "overlay": 1., "pegtop": 4. }) layer2_key = ("field_lines", {}) plotter.add_grey_layer(postproc_key=layer2_key, hardness=1.0, intensity=0.68, skewness=0.4, blur_ranges=[[0.50, 0.60, 1.0]], shade_type={ "Lch": 0., "overlay": 2., "pegtop": 1. }) plotter.plot("dev", mask_color=(0., 0., 1.))
def func(fractal: fsm.Perturbation_mandelbrot = fractal, file_prefix: str = "test", x: mpmath.mpf = x, y: mpmath.mpf = y, dx: mpmath.mpf = dx, xy_ratio: float = xy_ratio, dps: int = dps, max_iter: int = max_iter, nx: int = nx): # # interior_detect: bool=True): interior_detect = False # True # fractal.clean_up(file_prefix) fractal.zoom(precision=dps, x=x, y=y, dx=dx, nx=nx, xy_ratio=xy_ratio, theta_deg=0., projection="cartesian", antialiasing=False) fractal.calc_std_div( complex_type=np.complex128, file_prefix=file_prefix, subset=None, max_iter=max_iter, M_divergence=1.e3, epsilon_stationnary=1.e-4, pc_threshold=0.1, SA_params={ "cutdeg": 8, "cutdeg_glitch": 8, "SA_err": 1.e-20, # "SA_err" 1e-6, "cutdeg": 8 fails # "SA_err" 1e-20, "cutdeg": 8 fails # "SA_err" 1e-20, "cutdeg": 64 fails # "SA_err" 1e-50, "cutdeg": 64 fails # "SA_err" 1e-200, "cutdeg": 64 "use_Taylor_shift": True }, glitch_eps=1.e-6, interior_detect=interior_detect, glitch_max_attempt=4) fractal.run() gold = np.array([255, 210, 66]) / 255. black = np.array([0, 0, 0]) / 255. purple = np.array([181, 40, 99]) / 255. citrus2 = np.array([103, 189, 0]) / 255. colors1 = np.vstack((citrus2[np.newaxis, :])) colors2 = np.vstack((black[np.newaxis, :])) colormap = fscolors.Fractal_colormap(kinds="Lch", colors1=colors1, colors2=colors2, n=200, funcs=None, extent="mirror") mask_codes = [2] #, 3, 4] mask = fs.Fractal_Data_array( fractal, file_prefix=file_prefix, postproc_keys=('stop_reason', lambda x: np.isin(x, mask_codes)), mode="r+raw") plotter = fs.Fractal_plotter( fractal=fractal, base_data_key=("potential", {}), # ("field_lines", {"n_iter": 10, "swirl": 1.}), , base_data_prefix=file_prefix, base_data_function=lambda x: x, colormap=colormap, probes_val=[0.25, 0.75], probes_kind="qt", mask=mask) # plotter.add_calculation_layer(("potential", {})) plotter.add_grey_layer(postproc_key=("DEM_shade", { "kind": "potential", "theta_LS": 30., "phi_LS": 70., "shininess": 300., "ratio_specular": 15000. }), blur_ranges=[], hardness=0.9, intensity=0.8, shade_type={ "Lch": 1.0, "overlay": 0., "pegtop": 4. }) plotter.plot(file_prefix, mask_color=(0., 0., 1.))
def plot(): """ Example plot of "Dinkydau flake" location, classic test case for perturbation technique and glitch correction. """ directory = "./flake" # Dinkydau flake # http://www.fractalforums.com/announcements-and-news/pertubation-theory-glitches-improvement/msg73027/#msg73027 # Ball method 1 found period: 7884 x = "-1.99996619445037030418434688506350579675531241540724851511761922944801584242342684381376129778868913812287046406560949864353810575744772166485672496092803920095332" y = "0.00000000000000000000000000000000030013824367909383240724973039775924987346831190773335270174257280120474975614823581185647299288414075519224186504978181625478529" dx = "1.8e-157" precision = 200 # Set to True if you only want to rerun the post-processing part settings.skip_calc = False # Set to True to enable multi-processing settings.enable_multiprocessing = True nx = 600 xy_ratio = 0.5 theta_deg = 0. complex_type = np.complex128 mandelbrot = fsm.Perturbation_mandelbrot(directory) mandelbrot.zoom(precision=precision, x=x, y=y, dx=dx, nx=nx, xy_ratio=xy_ratio, theta_deg=theta_deg, projection="cartesian", antialiasing=False) mandelbrot.calc_std_div(complex_type=complex_type, file_prefix="dev", subset=None, max_iter=50000, M_divergence=1.e3, epsilon_stationnary=1.e-3, pc_threshold=0.1, SA_params={ "cutdeg": 8, "cutdeg_glitch": 8, "SA_err": 1.e-4, "use_Taylor_shift": True }, glitch_eps=1.e-6, interior_detect=False, glitch_max_attempt=20) mandelbrot.run() glitched = fs.Fractal_Data_array(mandelbrot, file_prefix="dev", postproc_keys=('stop_reason', lambda x: x == 3), mode="r+raw") potential_data_key = ("potential", {}) citrus2 = np.array([103, 189, 0]) / 255. citrus_white = np.array([252, 251, 226]) / 255. wheat1 = np.array([244, 235, 158]) / 255. wheat2 = np.array([246, 207, 106]) / 255. wheat3 = np.array([191, 156, 96]) / 255. lavender1 = np.array([154, 121, 144]) / 255. lavender2 = np.array([140, 94, 134]) / 255. lavender3 = np.array([18, 16, 58]) / 255. def wave(x): return 0.5 + (0.4 * (x - 0.5) - 0.6 * 0.5 * np.cos(x * np.pi * 3.)) colormap = fscolors.Fractal_colormap( kinds="Lch", colors1=np.vstack( (citrus_white, wheat2, wheat1, wheat2, wheat1, wheat2, wheat3, wheat1, lavender2, wheat1, wheat2, wheat3, wheat1, lavender2, wheat1, lavender3, lavender2, lavender3, lavender1, lavender3, lavender2)), colors2=None, n=100, funcs=lambda x: wave(x), extent="mirror") # colormap.extent = "mirror" #"repeat" plotter = fs.Fractal_plotter( fractal=mandelbrot, base_data_key= potential_data_key, # potential_data_key, #glitch_sort_key, base_data_prefix="dev", base_data_function=lambda x: x, # np.sin(x*0.0001), colormap=colormap, probes_val=np.linspace(0., 1., 22)** 0.2, # 200. + 200, #* 428 - 00.,#[0., 0.5, 1.], #phi * k * 2. + k * np.array([0., 1., 1.5, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0]) / 3.5, probes_kind="qt", #"z", "qt" mask=glitched) #plotter.add_calculation_layer(postproc_key=potential_data_key) layer1_key = ("DEM_shade", { "kind": "potential", "theta_LS": 30., "phi_LS": 50., "shininess": 30., "ratio_specular": 15000. }) plotter.add_grey_layer( postproc_key=layer1_key, intensity=0.75, blur_ranges=[], #[[0.99, 0.999, 1.0]], disp_layer=False, #skewness=0.2, normalized=False, hardness=0.35, skewness=0.0, shade_type={ "Lch": 1.0, "overlay": 1., "pegtop": 4. }) layer2_key = ("field_lines", {}) plotter.add_grey_layer(postproc_key=layer2_key, hardness=1.0, intensity=0.68, skewness=0.4, blur_ranges=[[0.50, 0.60, 1.0]], shade_type={ "Lch": 0., "overlay": 2., "pegtop": 1. }) plotter.plot("dev", mask_color=(0., 0., 1.))
def make_M2_img(self, x, y, dx, precision, nx, complex_type, test_name, prefix, interior_detect=False, mask_codes=[3, 4], SA_params={"cutdeg": 64, "cutdeg_glitch": 8}, antialiasing=False, colormap=None, hardness=0.75, intensity=0.95, probes_val=[0., 0.25], blur_ranges=[], grey_layer_key=None, glitch_max_attempt=1, xy_ratio=1.0): """ """ test_dir = os.path.join(self.image_dir, test_name) mandelbrot = fsmodels.Perturbation_mandelbrot(test_dir) mandelbrot.zoom( precision=precision, x=x, y=y, dx=dx, nx=nx, xy_ratio=xy_ratio, theta_deg=0., projection="cartesian", antialiasing=antialiasing) mandelbrot.prepare_calc( kind="calc_std_div", complex_type=complex_type, file_prefix=prefix, subset=None, max_iter=50000, M_divergence=1.e3, epsilon_stationnary=1.e-3, pc_threshold=0.1, SA_params=SA_params, glitch_eps=1.e-6, interior_detect=interior_detect, glitch_max_attempt=glitch_max_attempt) mandelbrot.run() mask = fs.Fractal_Data_array(mandelbrot, file_prefix=prefix, postproc_keys=('stop_reason', lambda x: np.isin(x, mask_codes)), mode="r+raw") potential = ("potential", {}) if colormap is None: colormap = self.colormap if grey_layer_key is None: grey_layer_key = ("DEM_shade", {"kind": "potential", "theta_LS": 50., "phi_LS": 40., "shininess": 40., "ratio_specular": 8.}) plotter = fs.Fractal_plotter( fractal=mandelbrot, base_data_key=potential, base_data_prefix=prefix, base_data_function=lambda x:x, colormap=colormap, probes_val=probes_val, probes_kind="qt", mask=mask) plotter.add_grey_layer( postproc_key=grey_layer_key, blur_ranges=blur_ranges, hardness=hardness, intensity=intensity, shade_type={"Lch": 1.0, "overlay": 0., "pegtop": 4.}) plotter.plot(prefix, mask_color=(0., 0., 1.)) new_file = os.path.join(self.image_dir, test_name + "_" + prefix + ".png") shutil.move(os.path.join(test_dir, prefix + ".png"), new_file) shutil.rmtree(test_dir) return new_file