Use permanent temporaries in descriptor for minimum cache misses

[mattsignorelli]

No description provided.

[mattsignorelli]

https://github.com/jump-dev/MutableArithmetics.jl

This (hopefully) will make it very easy to handle temporaries and increase speed

[mattsignorelli]

https://github.com/jump-dev/MutableArithmetics.jl

This (hopefully) will make it very easy to handle temporaries and increase speed

After a while of trying to learn the very not well documented API for this interface, I found this only works for simple arithmetic (+,-,*,/) and no other functions, therefore I think it we will have to write our own temporary management. One idea I have is to have a macro @usetemps which appends some symbol to each of the operators ( similar to @. ) and converts the product of the expression to heap allocated TPS or ComplexTPS. Those special operators with appended symbols return temporary types which have special behaviors / avoid heap allocation if possible. If a temporary is included in the evaluation, it can be reused or mutated.

[mattsignorelli]

This actually doesn't work, because you'd have to basically redefine PEMDAS

[mattsignorelli]

The solution for this will be to define a macro that replaces all +, -, *, /, ^ symbols with other overloaded infix symbols of same precendence (see here) . These overloaded symbols will be transparent to all types except for TPSs, which will use special temporary functions. All other known overloaded functions (sin, cos, etc) will be converted in the macro to something like _sin_t which will also be transparent to all types except TPSs. The macro will then wrap the expression in a function that copies the resulting TPS temporary type into a TPS, or be transparent if not. The macro will have to be appended to all expressions

[mattsignorelli]

Implemented with @FastGTPSA