sort_()

• Sort spikes using the feature and waveform info
• Called by "jrc sort myparam.prm"

Syntax

• S0 = sort_(P)

Input

• P: Parameters struct (P)

Output

• S0: Master struct (S0)

Precondition

• Spike detection is already done using detect_(). is_detected_() checks whether the spike detection is already done. "jrc sort" command calls detect_() if [[is_detected_()] returns false.

Operations

• fet2clu_(): Cluster the spike features (trFet_spk)
  • cluster_spacetime_(): Apply spatiotemporal DPCLUS (divide and conquer density-based clustering)
    • For each site, compute Rho (density) and Delta (distance) values for each spike using DPCLUS clustering
    • calc_dc2_(): Compute distance cut-off value (dc) based on the P.dc_percent parameter.
    • cuda_rho_(): Compute local density (rho) per each spike using GPU.
    • cuda_delta_(): Compute separation distance (delta) per each spike using GPU.
• postCluster_(): Group clusters using rho and delta values (see P.rho_cut, P.delta1_cut), and apply automated merging based on the mean unit waveforms
• post_merge_(): Merge units based on the mean spike waveforms. Unit pairs with correlation above the threshold (P.maxWavCor) is automatically merged.