Skip to content

Instantly share code, notes, and snippets.

@vadimkantorov
Last active September 22, 2021 07:51
Show Gist options
  • Select an option

  • Save vadimkantorov/d9b56f9b85f1f4ce59ffecf893a1581a to your computer and use it in GitHub Desktop.

Select an option

Save vadimkantorov/d9b56f9b85f1f4ce59ffecf893a1581a to your computer and use it in GitHub Desktop.
Compact Bilinear Pooling in PyTorch using the new FFT support
import torch
import torch.nn as nn
class CompactBilinearPooling(nn.Module):
def forward(self, bottom1, bottom2):
assert bottom1.size(1) == self.input_dim1 and bottom2.size(1) == self.input_dim2
batch_size, _, height, width = bottom1.size()
bottom1_flat = bottom1.permute(0, 2, 3, 1).contiguous().view(-1, self.input_dim1)
bottom2_flat = bottom2.permute(0, 2, 3, 1).contiguous().view(-1, self.input_dim2)
sketch_1 = bottom1_flat.mm(self.sparse_sketch_matrix1)
sketch_2 = bottom2_flat.mm(self.sparse_sketch_matrix2)
fft1_real, fft1_imag = torch.rfft(sketch_1, 1).permute(2, 0, 1)
fft2_real, fft2_imag = torch.rfft(sketch_2, 1).permute(2, 0, 1)
fft_product_real = fft1_real * fft2_real - fft1_imag * fft2_imag
fft_product_imag = fft1_real * fft2_imag - fft1_imag * fft2_real
cbp_flat = torch.irfft(torch.stack([fft_product_real, fft_product_imag], dim = -1)
cbp = cbp_flat.view(batch_size, height, width, self.output_dim) * self.output_dim
if self.sum_pool:
return cbp.sum(dim = 1).sum(dim = 1)
return cbp.permute(0, 3, 1, 2)
@pangjh3
Copy link

pangjh3 commented Apr 20, 2018

Thanks for your code, how to install the new fft support?

@vadimkantorov
Copy link
Author

Just install PyTorch from master branch or even 0.4 version probably has FFT

@ayumiymk
Copy link

Thanks for your code first. I have a question that in the other implements, like Torch version and Tensorflow version, there is a zero_padding before feeding the tensor into the fft. But in this code, I don't see the zero_padding.

Thanks very much!

@hj0921
Copy link

hj0921 commented Mar 19, 2021

hello,

torch.stack([torch.arange(in_features), rand_h]) where in_features is not defined. How to fix it?

thanks!

@vadimkantorov
Copy link
Author

Thanks for noting this. Fixed! It should have been in_channels

@vadimkantorov
Copy link
Author

Some ways to improve the code: make use of the new PyTorch fft module, complex support. Figure out dense x sparse matmul (currently I'm materializing the sparse sketch)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment