ConceptPythonIntermediate Finite volumes: why a conservation law is solved by bookkeeping the fluxes cp = np.pad(c, 1) # clean water beyond both endscp = np.pad(c, 1) # clean water beyond both endscp = np.pad(c_snap, 1)cp = np.pad(c, 1)point = lambda c: -uc * (np.pad(c, 1)[2:] - np.pad(c, 1)[:-2]) / (2 * h) - c * duc
ToolPythonBeginner MRI k-space reconstruction with numpy.fft: shifts, fold-over, ringing padded = np.abs(recon(np.pad(kspace_clean, N // 2))) * 4 # 512 x 512; ifft2 now divides by 4 times more samples