mandelbrot Python 3 #8 program

# The Computer Language Benchmarks Game
# https://salsa.debian.org/benchmarksgame-team/benchmarksgame/
#   
#     line-by-line from Greg Buchholz's C program



import sys

def main():

    w = h = bit_num = 0
    byte_acc = 0
    i = 0; iter = 50
    x = y = limit = 2.0
    Zr = Zi = Cr = Ci = Tr = Ti = 0.0

    w = h = int(sys.argv[1])

    sys.stdout.write(f'P4\n{w} {h}\n'); sys.stdout.flush()

    for y in range(h):

        for x in range(w):

            Zr = Zi = Tr = Ti = 0.0 
            Cr = (2.0 * x / w - 1.5); Ci = (2.0 * y / h - 1.0)  
           
            for i in range(iter):             
                if Tr+Ti <= limit*limit:
                    Zi = 2.0*Zr*Zi + Ci
                    Zr = Tr - Ti + Cr
                    Tr = Zr * Zr
                    Ti = Zi * Zi
            
            
            byte_acc = byte_acc << 1
            if Tr+Ti <= limit*limit: byte_acc = byte_acc  | 0x01            
            
            bit_num += 1         

            if bit_num == 8:          
                sys.stdout.buffer.write(bytes([byte_acc]))        
                byte_acc = 0
                bit_num = 0


            elif x == w - 1:

                byte_acc = byte_acc << (8-w%8)   
                sys.stdout.buffer.write(bytes([byte_acc]))  
                byte_acc = 0
                bit_num = 0

main() 
    

NOTES:
64-bit Ubuntu quad core
Python 3.13.0


 Fri, 17 Jan 2025 13:37:40 GMT

MAKE:
mv mandelbrot.python3-8.python3 mandelbrot.py
pyright .
0 errors, 0 warnings, 0 informations 

4.84s to complete and log all make actions

COMMAND LINE:
 /opt/src/Python-3.13.0/bin/python3 -OO mandelbrot.py 16000

(BINARY) PROGRAM OUTPUT NOT SHOWN