Martin Roelfs

Yes that's one case, and another one I've found is when performing bivector splits: when the two bivectors in the split B = B_1 + B_2 start to become to...

I've taken the liberty of making a PR for this with a slightly different API then suggested above: I added a `.pinv()` method to really make it easy to distinguish....

@hugohadfield that sounds great! Cleaning up the code base is what GA is for ;). I have made one potentially contentious change to `_layout.py`: I've removed the checking if a...

@hugohadfield, when it comes to cleaning up the code base, should that be part of this PR or is it better to leave this one focused on adding the general...

Due to numerical errors, there are two possibilities. In theory the degenerate cases lead to a "no inverse exists" error, although in practice I've had instances where they are different...

So if I create the known_split ```python def r4_split(): alg = Layout([1, 1, 1, 1]) e1, e2, e3, e4 = alg.basis_vectors_lst return {'B': 2*e1*e2 + 2*e3*e4, 'Bs': [2*e3*e4, 2*e1*e2], 'ls':...

I've now added tests for `exp(log(R)) == R` in higher dimensions as well. This slows down the tests significantly but was well worth it because it smoked out some bugs...

I like this a lot too, just needs documentation and testing as far as I'm concerned. Other than that it's a good step in this direction, and it also sets...

One example I have is from reaction kinetics. Suppose we have an equilibrium reaction between compound A and B with rate constants `k_forward` and `k_reverse`. Then these can be expressed...

Good one @Pitje06. You are right, I don't like it ;). Integer constraints in fitting are a pain. Do you have any ideas how to solve this? Any algorithms? Naively...