search

timmyfaraday / HarmonicPowerModels.jl

An extension package of PowerModels.jl for Harmonic (Optimal) Power Network
BSD 3-Clause "New" or "Revised" License
6 stars 1 forks source link

PSCC paper: Gradient of the excitation current #1

Closed timmyfaraday closed 2 years ago

timmyfaraday commented 3 years ago

The gradient, also the Hessian but not of interest given JuMP limitations, of an interpolation is freely available which means that the autodiff setting could be avoided, possibly speeding up computation.

Issue: to get the gradient at a node x of an interpolation int, the following is required: grad = gradient(int,x...), which introduces an extra argument in the function call, making it impossible to register as is (status 29/07/21).

Possible fix: add function to the power models data file rather than the interpolation, where int is already parsed as an argument.

timmyfaraday commented 3 years ago

Both the interpolation and gradient function have been added to the hdata dictionary. However, evaluation of the gradient function does not work in JuMP/Ipopt for a coordinate in the feasible space, where the function evaluating the function itself in the terminal works for that given coordinate.

frederikgeth commented 2 years ago

was this solved?

timmyfaraday commented 2 years ago

Currently, this remains unresolved. However, similar to the excitation, extrapolation of the Splines might resolve this issue.