The isfc function currently takes a single input array and computes the pairwise ISFCs among all voxels in this array. However, an important use-case is computing ISFCs between one set of voxels and another (non-identical) set of voxels or "targets". For example, I want to compute ISFCs between the 500 voxel time series in my ROI and 360 target time series from a parcellation. This will result in an asymmetric ISFC n_voxels by n_targets ISFC matrix. I have a fix for this on deck and will make a PR soon.