【RPi PICO】计算Mandelbrot

dcexpert

【RPi PICO】计算Mandelbrot [复制链接]

 

将micropython源码中的examples/mandel.py例子修改了一下，用于测试计算性能

 

from time import ticks_ms, ticks_diff

def run(func, param = None):
    t1 = ticks_ms()
    if param == None:
        func()
    else:
        func(param)
    t2 = ticks_ms()
    print('calc time:', ticks_diff(t2, t1), 'ms')

def mandelbrot():
    # returns True if c, complex, is in the Mandelbrot set
    # @micropython.native
    def in_set(c):
        z = 0
        for i in range(40):
            z = z * z + c
            if abs(z) > 4:
                return False
        return True

    for v in range(31):
        for u in range(81):
            if in_set((u / 30 - 2) + (v / 15 - 1) * 1j):
                print(' ', end='')
            else:
                print('#', end='')
        print()

run(mandelbrot)

 

 

运行结果

 

############################################################ ####################
#################################################################################
########################################################  #######################
#####################################################       #####################
######################################################      #####################
################################################### ###    ### # ################
############################################   ##                ################
############################################                           ##########
###########################################                          ############
##########################################                             ##########
########################################                                 ########
########################## ### ########                                 #########
##########################          ###                                 #########
########################             #                                   ########
#######################                                                ##########
                                                                    #############
#######################                                                ##########
########################             #                                   ########
##########################          ###                                 #########
########################## ### ########                                 #########
########################################                                 ########
##########################################                             ##########
###########################################                          ############
############################################                           ##########
############################################   ##                ################
################################################### ###    ### # ################
######################################################      #####################
#####################################################       #####################
########################################################  #######################
#################################################################################
############################################################ ####################

 

dcexpert

circuitpython版本程序

 

from time import monotonic

def run(func, param = None):
    t1 = monotonic()
    if param == None:
        func()
    else:
        func(param)
    t2 = monotonic()
    print('calc time:', round((t2 - t1)*1000), 'ms')

def mandelbrot():
    # returns True if c, complex, is in the Mandelbrot set
    # @micropython.native
    def in_set(c):
        z = 0
        for i in range(40):
            z = z * z + c
            if abs(z) > 4:
                return False
        return True

    for v in range(31):
        for u in range(81):
            if in_set((u / 30 - 2) + (v / 15 - 1) * 1j):
                print(' ', end='')
            else:
                print('>', end='')
        print()

run(mandelbrot)

 

w494143467

感谢分享！！！