Newton-Raphson method for state estimation

[TonyXiang8787]

Newton-Raphson method should solve the issue in #237

[mgovers]

Update 2024/02/08:

Since #438 , Newton-Raphson for state estimation is now in beta as an experimental feature.

Because we do not officially support it yet, nor guarantee its correctness, it is not yet possible to call this feature via the regular API's PowerGridModel.calculate_state_estimation).

Instead, to use the experimental feature, call the protected method PowerGridModel._calculate_state_estimation with as arguments:

• calculation_method=CalculationMethod.newton_raphson
• experimental_features=power_grid_model.enum._ExperimentalFeatures.enabled
• any other arguments which are the same arguments as for the iterative linear case.

E.g. someone who would usually call

from power_grid_model import PowerGridModel
from power_grid_model.enum import CalculationMethod

# create input data
model = PowerGridModel(input_data)
model.calculate_state_estimation(calculation_method=CalculationMethod.iterative_linear)

may run the Newton-Raphson state estimation by running

from power_grid_model import PowerGridModel
from power_grid_model.enum import CalculationMethod, _ExperimentalFeatures

# create input data
model = PowerGridModel(input_data)
model.calculate_state_estimation(calculation_method=CalculationMethod.newton_raphson, experimental_features=_ExperimentalFeatures.enabled)