setpaths;

embedding.py -s <subject> -f <output filebasename>

where input is a connectivity matrix saved as a .mat file.

[clus, randpath] = DoFindPaths(n);

where `n` is the number of permutations,
and a config file is loaded in the script.

config/config_my18_path.m

randPath = pathsPermute(clus, n);

[distDMN] = loadDistDMN_group(hemi);

[dist, zone] = distExactGeodesic(clus.peaks(find(clus.edgeNet == 13)), '32', 'L', 'zones', '1');

DoMyelin();

[distances, zones, surfi_164] = loadDistDMN_individual(sub, hemi);