#
# PYTHON • MyCompiler
# artecomplexa
# https://[Log in to view URL]
# artecomplexa
# 
# 2026 Junho Ter, 9
# by 
# j
# 

import matplotlib.pyplot as plt

import numpy as np

from numpy import pi
from numpy import sqrt
from numpy import sin
from numpy import cos
from numpy import atan

import colorsys

MAXRESOLUTION = 1024

def arte(tipo):

    #
    # Configuração Geral
    #
    
    plt.style.use('dark_background')
    plt.figure(figsize=(7, 7))

    #
    # Parâmetros
    #

    match tipo:
        
        case 1:  
            
            a, b = 4, 5
            
        case 2:
            
            a, b = -2, 2
            
        case 3:
            
            a, b = -1, -1
            
        case 4:
            
            a, b = 5.9, 5.64
            
        case 5:
            
            a, b = 3.64, 1.71

        case 6:
            
            a, b = -5, -1

    #
    # Inicializações
    #

    H = np.zeros((MAXRESOLUTION, MAXRESOLUTION))
    L = np.zeros((MAXRESOLUTION, MAXRESOLUTION)) 
    S = np.zeros((MAXRESOLUTION, MAXRESOLUTION))
    
    xn = +1.00
    yn = -0.50
    
    max = 262144
    lim = 2

    r = 1.00

    while( r <= lim ):
        
        for t in range(0, 57600):
            
            #
            # Circunferência
            #
            
            w = sin(r*t*pi/28800) + a
            z = cos(2*r*t*pi/28800) + b
            
            #
            # Atrator de Simone
            #
        
            xnp = np.sin(xn**3 - yn**2 + w)
            ynp = np.cos(3*xn*yn + z)
            
            # 
            # Transforma o intervalo 
            # [-1.0, +1.0] para 
            # [0, MAXRESOLUTION - 1]
            #
        
            idx = int((xnp + 1.0) * 0.5 * (MAXRESOLUTION - 1))
            idy = int((ynp + 1.0) * 0.5 * (MAXRESOLUTION - 1))
            
            #
            # Filtro de segurança
            #
        
            if 0 <= idx < MAXRESOLUTION and 0 <= idy < MAXRESOLUTION:

               #
               # Cor baseada no ângulo real do ponto 
               # (Gera degradê fractal)
               # 
            
               H[idx, idy] = (np.arctan2(xnp, ynp) + np.pi) / (2 * np.pi)

               # 
               # Acumulador de brilho: 
               # Quanto mais vezes o ponto é visitado, 
               # mais ele brilha
               #
            
               L[idx, idy] = min(L[idx, idy] + 0.1, 0.8)

               #
               # Saturação Máxima
               #
            
               S[idx, idy] = 1.0

            xn = xnp
            yn = ynp
            
        #
        # Incrementa r
        #

        r += 0.5

    #
    # Conversão ultra~rápida vetorizada
    # de HLS para RGB
    # 
    
    hls_to_rgb_v = np.vectorize(colorsys.hls_to_rgb)
    R, G, B = hls_to_rgb_v(H, L, S)
    
    # 
    # Agrupa os canais e rotaciona 
    # para coincidir com os eixos 
    # cartesianos convencionais
    # 
    
    RGB = np.dstack((R, G, B))
    RGB = np.swapaxes(RGB, 0, 1)

    # 
    # Gera o gráfico
    # 
    
    tit = f"Ars ~ j ~ | a = {a} b = {b}"
    plt.imshow(RGB, extent=[-1.0, 1.0, -1.0, 1.0], origin='lower')
    plt.title(tit)
    plt.xlabel("x")
    plt.ylabel("y")
    plt.axis('off')
    plt.show()

arte(6)