()
| 33 | from MB_fortran import mb |
| 34 | |
| 35 | def main(): |
| 36 | imax = 255 |
| 37 | for N in [500,1000,2000,5000]: |
| 38 | T_s = time.time() |
| 39 | nR, nC = N, N |
| 40 | Re = np.linspace(-0.7440, -0.7433, nC) |
| 41 | Im = np.linspace( 0.1315, 0.1322, nR) |
| 42 | img = np.zeros((nR*nC,), dtype=np.int16) |
| 43 | img = mb(Re, Im, imax).reshape(nR,nC) |
| 44 | print(N, time.time() - T_s) |
| 45 | if __name__ == '__main__': main() |