def METHOD_NAME(
    datasets,
    ds001419_data,
    surface_files,
    tmp_path_factory,
):
    """Test surface-warping workflow with mesh surfaces available, but not in standard space.
    The transforms should be applied and all of the standard-space outputs should be generated.
    """
    tmpdir = tmp_path_factory.mktemp("test_warp_surfaces_to_template_wf")
    subject_id = "01"
    wf = anatomical.init_warp_surfaces_to_template_wf(
        fmri_dir=datasets["ds001419"],
        subject_id=subject_id,
        output_dir=tmpdir,
        omp_nthreads=1,
        mem_gb=0.1,
    )
    wf.inputs.inputnode.lh_pial_surf = surface_files["native_lh_pial"]
    wf.inputs.inputnode.rh_pial_surf = surface_files["native_rh_pial"]
    wf.inputs.inputnode.lh_wm_surf = surface_files["native_lh_wm"]
    wf.inputs.inputnode.rh_wm_surf = surface_files["native_rh_wm"]
    # transforms (only used if warp_to_standard is True)
    wf.inputs.inputnode.anat_to_template_xfm = ds001419_data["anat_to_template_xfm"]
    wf.inputs.inputnode.template_to_anat_xfm = ds001419_data["template_to_anat_xfm"]
    wf.base_dir = tmpdir
    wf.run()
    # All of the possible fsLR surfaces should be available.
    out_anat_dir = os.path.join(tmpdir, "xcp_d", "sub-01", "anat")
    for key, filename in surface_files.items():
        if "fsLR" in key:
            out_fname = os.path.basename(filename)
            out_file = os.path.join(out_anat_dir, out_fname)
            assert os.path.isfile(out_file), "\n".join(sorted(os.listdir(out_anat_dir)))