Computes the AB join in the sense of the matrix profile, that is, this function compares every subsequences in two input GTS. Both GTS must be bucketized and filled. The output is a GTS. Its length and metadata are determined by the first GTS. The value is the matrix profile result. The elevation is the argmin bucket index that references to the second input GTS.