The following code runs:
import numpy as np
from aspire.basis import FBBasis2D
L = 16
basis = FBBasis2D((L, L), dtype=np.float64)
im = np.random.randn(2*L, 2*L).astype(basis.dtype)
coef = basis.expand(im)
but should error. We should find out how widespread this behavior is (just expand or also evaluate_t or even evaluate?) and in what classes (just FBBasis2D or others as well) and then fix it.