conv_thm_2D.m
e04f06e1e9b93D_parallel

conv_thm_2D.m
View Options

			function [ conv ] = conv_thm_2D( F, G, Pad )
			%conv_thm_2D computes the convolution between F and G with a zero pading Pad
			% kspace -> real -> product -> kspace

			[Nr, Nz] = size(F);

			f = ifft2(F,NrPad, NzPad);
			g = ifft2(G,NrPad, NzPad);

			conv_pad = fft2(f.*g); % convolution becomes product

			conv = conv_pad(1:Nr,1:Nz); % remove padding
			end