we refer to h q , k q , l q as the Miller indices of q. That is, there is a subset of the integer grid Z 3 that gets mapped to Q ⋆ using the matrix M * := (a * , b * , c * ) T . In practice, real data is intrinsically noisy and thus we aim to satisfy the Laue condition up to some bounded error, i.e., we assume that there is a known maximum absolute allowed error which is a hyperparameter of our algorithm.