Hello,
maybe I'm able to find here a solution for my current problem.
I'm trying to calculate taylor coefficients recursively. It's necessary to calculate the Jacobian, which is in my case a 2x2 matrix with interval. Then I need to multiply the Jacobian with an AScalar vector, which is not possible.
I already implemented the algorithm in Matlab with gradientinit(x) but its really slow so I'm trying to work in C++ now.
Also the solution needs to be from type AScalar because I need the value() and the derivartives for the next computed taylor coefficients.

I hope there's someone out there who wants to help me.
Best whishes