быстрый и эффективный алгоритм аппроксимации методом наименьших квадратов в C?

Я пытаюсь реализовать линейный метод наименьших квадратов, помещающийся на 2 массивах данные: время против амплитуды. Единственная техника, которую я знаю до сих пор, - это проверить все возможные точки m и b в (y = m * x + b), а затем выяснить, какая комбинация лучше всего подходит для моих данных, чтобы иметь наименьшую ошибку. Однако я думаю, что повторение такого количества комбинаций иногда бесполезно, потому что оно проверяет все. Есть ли способы ускорить процесс, которые я не использую? не знаю о чем? Спасибо.