After today's dicsussion with Mihai we realised that the Local Connectivity creator might not be giving the correct matrix, that is, is yieldining only the geodesic distance matrix.
In surfaces_scientific.py --> CortexScientific – compute_local_connectivity
self.local_connectivity.compute() # Evaluate equation based distance
This line computes the spatial profile of the local connectivity, and in principle should replace the geodesic distance matrix, with a new one.
Hower after having a look at the local connectivity creator.
It is only computing the geodesic distance via the compute+sparse matrix method.
Because of the checks that exist in the "configure" method in CortexScientific, the method "compute_local_connectivity" is never reached...