The rank-(L,L,1,1) block term decomposition (BTD) with spatial orthonormality (BTD-O) applied to 4-way multi-subject fMRI data achieves good performance due to preserving higher spatial structure and reducing crosstalk between components. However, the high rank L value (e.g., 35) of BTD-O for fMRI data leads to high computation complexity. Moreover, multi-subject fMRI data contains high noise nature. Although an accelerated BTD-O (accBTD-O) was proposed, it showed similar performance to BTD-O. Inspired by the compression, smoothing, and spatial structure invariance features of the pooling sche...