cherab.phix.inversion.Lcurve.optimize#

Lcurve.optimize(itemax=5, bounds=(-20.0, 2.0))Source#

Excute the optimization of L-curve regularization.

Warning

This method will be deprecated in the future. Please use solve() instead.

This method is used to seek the optimal regularization parameter computing curvature. The optimal regularization parameter corresponds to the index of maximum curvature. This procedure is iterated by up to the itemax times. Every time iterately calculating, the range of regularization parameters is narrowed to FWHM around the maximum curvature point. The optimal regularization parameter is cached to self._lambda_opt which can be seen lambda_opt property. And, lambdas is updated and stored to the property.

Parameters:

  • itemax (int) – iteration times, by default 5

  • bounds (tuple[float, float]) – initial bounds of log10 of regularization parameter, by default (-20.0, 2.0).

Returns:

optimized solution vector

Return type:

numpy.ndarray