'c': c.T, # m/s 'rho': rho.T, # density 'Q': Q, # can be inf 'nx': nx, # count 'nz': nz, # count 'freeSurf': freeSurf, # t r b l 'nPML': nPML, 'geom': geom, 'cache': cache, 'cacheDir': cacheDir, 'freqs': freqs, 'nky': nky, } sp = SeisFDFDProblem(systemConfig) jobs = sp.forwardAccumulate() def trackprogress(sp, jobs, interval=1.0): systemJobs = jobs['systemJobs'] jobkeys = systemJobs.keys() jobkeys.sort() fig = plt.figure() ax1 = fig.add_axes([0.1, 0.10, 0.15, 0.85], xlabel='Subproblem', ylabel='Source') ax1.get_xaxis().set_major_locator(ticker.MaxNLocator(integer=True)) ax2 = fig.add_axes([0.25, 0.10, 0.75, 0.85], xlabel='Receiver')
'rho': rho.T, # density 'Q': Q, # can be inf 'nx': nx, # count 'nz': nz, # count 'freeSurf': freeSurf, # t r b l 'nPML': nPML, 'geom': geom, 'cache': cache, 'cacheDir': cacheDir, 'freqs': freqs, 'nky': nky, } sp = SeisFDFDProblem(systemConfig) jobs = sp.forwardAccumulate() def trackprogress(sp, jobs, interval=1.0): systemJobs = jobs['systemJobs'] jobkeys = systemJobs.keys() jobkeys.sort() fig = plt.figure() ax1 = fig.add_axes([0.1,0.10,0.15,0.85], xlabel='Subproblem', ylabel='Source') ax1.get_xaxis().set_major_locator(ticker.MaxNLocator(integer=True)) ax2 = fig.add_axes([0.25,0.10,0.75,0.85], xlabel='Receiver') im1 = ax2.imshow(np.zeros((nsrc, nrec)), vmin=-50*nky, vmax=50*nky, cmap=cm.bwr) im2 = ax1.imshow(np.zeros((nsrc, nsp)), vmin=0, vmax=2, interpolation='nearest', aspect='auto')
'nz': nz, # count 'freeSurf': freeSurf, # t r b l 'nPML': nPML, 'geom': geom, 'cache': cache, 'cacheDir': cacheDir, 'freqs': freqs, 'nky': nky, 'parFac': parFac, 'chunksPerWorker': chunksPerWorker, } sp = SeisFDFDProblem(systemConfig) jobs, G = sp.forwardAccumulate() def colourCodeNodes(graph): def mapColours(value): if value < 0: return (0, 0, 0) elif value == 0: return (0, 0, 1) elif value == 1: return (0, 1, 0) elif value == 3: return (1, 0, 0) def assessStatus(G, node):
'nx': nx, # count 'nz': nz, # count 'freeSurf': freeSurf, # t r b l 'nPML': nPML, 'geom': geom, 'cache': cache, 'cacheDir': cacheDir, 'freqs': freqs, 'nky': nky, 'parFac': parFac, 'chunksPerWorker': chunksPerWorker, } sp = SeisFDFDProblem(systemConfig) jobs, G = sp.forwardAccumulate() def colourCodeNodes(graph): def mapColours(value): if value < 0: return (0, 0, 0) elif value == 0: return (0, 0, 1) elif value == 1: return (0, 1, 0) elif value == 3: return (1, 0, 0) def assessStatus(G, node):