fixef: Implement Kacmarz method

[s3alfisc]

Context

Currently, pyfixest can compute fixed effect estimates in the fixef() method by using a sparse solver.

A potentially faster alternative is to implement a Kacmarz solver. See section 3 in the lfe paper.

To Do:

• Read up on how the Kacmarz solver actually works.
• Implement it.