There was a direct access to a NumPy configuration which was
removed in 1.22. This PR removes the special handling and rely
on accessing the configuration in a controlled way to avoid later
issues.

These changes remove the stride-based manual computation of the new shape, since
those are potentially sensitive to broadcasted arrays with no strides.

`Solve` has also been changed to match SciPy.

* Fix class, attribute, and function formatting
* Update the relevant glossary terms and references to them
* Rename documentation to better reflect their content
* Remove redundant documentation
* Reorder index so that it starts with the basic, necessary articles
* Grammar and wording updates

Closes #692

These changes only affect `PatternSub`, which now no longer allows constraints
on non (logic) variable terms in its patterns.  Likewise, the `values_eq_approx`
and `skip_identities_fn` options are no longer supported.

This commit introduces a `LocalOptTracker` object that performs an MRO-based
lookup of `LocalOptimizer`s that track `Op` types.

This improvement is due to the use of `x[indices] += y` instead of `np.add.at(x,
indices, y)`, but comes at the cost of not handling duplicate `indices` values.