BSplineApprox dimension mismatch

[effigies]

Command:

fmriprep-docker sourcedata/raw . participant -u $UID \
    --output-space func T1w MNI152NLin2009cAsym -w /data/scratch/multirun \
    --skip-bids-validation --nprocs 10 --use-syn-sdc --debug fieldmaps \
    --patch sdcflows=/home/chris/Projects/nipreps/sdcflows/sdcflows \
            niworkflows=/home/chris/Projects/nipreps/niworkflows/niworkflows \
            fmriprep=/home/chris/Projects/nipreps/fmriprep/fmriprep \
            smriprep=/home/chris/Projects/nipreps/smriprep/smriprep \
    -vv --config code/nipype.cfg

Uses the 22.0.2 image with the patched packages at master.

Log:

221201-02:04:10,906 nipype.workflow INFO:                                                                                      
         [Node] Setting-up "fmriprep_0_0_wf.single_subject_EB_wf.fmap_preproc_wf.wf_auto_00000.bs_filter" in "/scratch/fmriprep_0_0_wf/single_subject_EB_wf/fmap_preproc_wf/wf_auto_00000/bs_filter".                                                         
221201-02:04:10,911 nipype.workflow INFO:    
         [Node] Executing "bs_filter" <sdcflows.interfaces.bspline.BSplineApprox>                                                                                                                                                                             
221201-02:04:11,170 nipype.interface INFO:                                                                                     
         Resampling image with resolution exceeding 'zooms_min' (1.796875x1.796875x1.8000001 → 1.95x1.95x1.95).                
221201-02:04:12,800 nipype.workflow INFO:                                                                                                                                                                                                                     
         [MultiProc] Running 2 tasks, and 0 jobs ready. Free memory (GB): 55.89/56.29, Free processors: 1/10.                                                                                                                                                 
                     Currently running:     
                       * fmriprep_0_0_wf.single_subject_EB_wf.fmap_preproc_wf.wf_auto_00000.bs_filter                                                                                                                                                         
                       * fmriprep_0_0_wf.single_subject_EB_wf.fmap_preproc_wf.wf_auto_00000.zooms_bmask                        
221201-02:04:21,708 nipype.interface INFO:                                                                                                                                                                                                                    
         Approximating B-Splines grids (5x5x6, and 17x17x18 [knots]) on a grid of 118x118x78 (1086072) voxels, of which 1086072 fall within the mask.                                                                                                         
221201-02:04:24,103 nipype.workflow INFO:                                                                                                                                                                                                                     
         [Node] Finished "zooms_bmask", elapsed time 13.197545s.                                                               
221201-02:04:24,809 nipype.workflow INFO:                                                                                                                                                                                                                     
         [Job 359] Completed (fmriprep_0_0_wf.single_subject_EB_wf.fmap_preproc_wf.wf_auto_00000.zooms_bmask).                 
221201-02:04:24,813 nipype.workflow INFO:                                                                                                                                                                                                                     
         [MultiProc] Running 1 tasks, and 1 jobs ready. Free memory (GB): 56.09/56.29, Free processors: 2/10.                  
                     Currently running:                                                                                                                                                                                                                       
                       * fmriprep_0_0_wf.single_subject_EB_wf.fmap_preproc_wf.wf_auto_00000.bs_filter                          
221201-02:04:24,897 nipype.workflow INFO:                                                                                                                                                                                                                     
         [Node] Setting-up "fmriprep_0_0_wf.single_subject_EB_wf.fmap_preproc_wf.out_merge_fmap_mask" in "/scratch/fmriprep_0_0_wf/single_subject_EB_wf/fmap_preproc_wf/out_merge_fmap_mask".                                                                 
221201-02:04:24,900 nipype.workflow INFO:                                                                                                                                                                                                                     
         [Node] Executing "out_merge_fmap_mask" <nipype.interfaces.utility.base.Merge>
221201-02:04:24,901 nipype.workflow INFO:
         [Node] Finished "out_merge_fmap_mask", elapsed time 0.000207s.                                                        
221201-02:04:26,810 nipype.workflow INFO:
         [Job 361] Completed (fmriprep_0_0_wf.single_subject_EB_wf.fmap_preproc_wf.out_merge_fmap_mask).                       
221201-02:04:26,813 nipype.workflow INFO:
         [MultiProc] Running 1 tasks, and 0 jobs ready. Free memory (GB): 56.09/56.29, Free processors: 2/10.                  
                     Currently running:
                       * fmriprep_0_0_wf.single_subject_EB_wf.fmap_preproc_wf.wf_auto_00000.bs_filter                          
221201-12:24:33,177 nipype.workflow INFO:
         [Node] Finished "bs_filter", elapsed time 37222.247152s.
221201-12:24:33,177 nipype.workflow WARNING:
         Storing result file without outputs
221201-12:24:33,178 nipype.workflow WARNING:
         [Node] Error on "fmriprep_0_0_wf.single_subject_EB_wf.fmap_preproc_wf.wf_auto_00000.bs_filter" (/scratch/fmriprep_0_0_wf/single_subject_EB_wf/fmap_preproc_wf/wf_auto_00000/bs_filter)
221201-12:24:33,670 nipype.workflow ERROR:
         Node bs_filter failed to run on host dcb28336b466.
221201-12:24:33,670 nipype.workflow ERROR:
         Saving crash info to /out/sub-EB/log/20221201-010535_1988ebcf-05b6-4a82-a374-f297f2de2c81/crash-20221201-122433-fmriprep-bs_filter-3f4730e2-0f35-4a92-a123-a3215309504f.txt
Traceback (most recent call last):
  File "/opt/conda/lib/python3.9/site-packages/nipype/pipeline/plugins/multiproc.py", line 67, in run_node
    result["result"] = node.run(updatehash=updatehash)
  File "/opt/conda/lib/python3.9/site-packages/nipype/pipeline/engine/nodes.py", line 527, in run
    result = self._run_interface(execute=True)
  File "/opt/conda/lib/python3.9/site-packages/nipype/pipeline/engine/nodes.py", line 645, in _run_interface
    return self._run_command(execute)
  File "/opt/conda/lib/python3.9/site-packages/nipype/pipeline/engine/nodes.py", line 771, in _run_command
    raise NodeExecutionError(msg)
nipype.pipeline.engine.nodes.NodeExecutionError: Exception raised while executing Node bs_filter.

Traceback:
        Traceback (most recent call last):
          File "/opt/conda/lib/python3.9/site-packages/nipype/interfaces/base/core.py", line 398, in run
            runtime = self._run_interface(runtime)
          File "/opt/conda/lib/python3.9/site-packages/sdcflows/interfaces/bspline.py", line 231, in _run_interface
            interp_data[mask] = regressors @ model.coef_
          File "/opt/conda/lib/python3.9/site-packages/scipy/sparse/_base.py", line 623, in __matmul__
            return self._mul_dispatch(other)
          File "/opt/conda/lib/python3.9/site-packages/scipy/sparse/_base.py", line 555, in _mul_dispatch
            raise ValueError('dimension mismatch')
        ValueError: dimension mismatch

In addition, note that we're still taking ~10 hours to run this.