On Jun 3, 2011, at 3:34 AM, Biswa Sengupta wrote:
> Since you know about the underlying-data structures, what would be the best > way to go about writing automatic differentiation rules for the solution > structure after a state solve. For example, in my script, I have a field, V. > If I were to construct the gradient of V w.r.t to let's say g_Na there are > typically two ways - Newton type approximate gradient or using automatic > differentiation (like adolc). Are the solution data structures amenable to > automatic differentiation tools? If yes, do you guys have some examples? We have not done anything along these lines, but I recall that Prof. Edwin Garcia from Purdue has done some experiments with automatic differentiation to build up FiPy equations, using Sage. He's subscribed to the list, but I don't know how often he actually checks it.
