This is still incredibly WIP, and it's in the context of #140. I'm not really sure on how to proceed here, so dumping this here in case @chriselrod can see a way forward.
What is the correct way of selecting the correct calls depending on the value of theta and eta? I tried macros and inline functions, but I was likely doing something wrong which is why I got error messages that were a bit too dense for me (I can paste them here if it can potentially help).
How bad (tentatively) are the extra sign multiplications going to affect things?
The final conjugation of _C/C, I'm really not sure how to incorporate that with the new unreshaped C array, so I don't really know what to ask here other than 'how do?'
This is still incredibly WIP, and it's in the context of #140. I'm not really sure on how to proceed here, so dumping this here in case @chriselrod can see a way forward.