desc.objectives.PlasmaVesselDistance

class desc.objectives.PlasmaVesselDistance(surface, eq=None, target=None, bounds=None, weight=1, normalize=True, normalize_target=True, surface_grid=None, plasma_grid=None, use_softmin=False, alpha=1.0, name='plasma-vessel distance')Source

Target the distance between the plasma and a surrounding surface.

Computes the minimum distance from each point on the surface grid to a point on the plasma grid. For dense grids, this will approximate the global min, but in general will only be an upper bound on the minimum separation between the plasma and the surrounding surface.

NOTE: for best results, use this objective in combination with either MeanCurvature or PrincipalCurvature, to penalize the tendency for the optimizer to only move the points on surface corresponding to the grid that the plasma-vessel distance is evaluated at, which can cause cusps or regions of very large curvature.

NOTE: When use_softmin=True, ensures that alpha*values passed in is at least >1, otherwise the softmin will return inaccurate approximations of the minimum. Will automatically multiply array values by 2 / min_val if the min of alpha*array is <1. This is to avoid inaccuracies that arise when values <1 are present in the softmin, which can cause inaccurate mins or even incorrect signs of the softmin versus the actual min.

Parameters:
  • surface (Surface) – Bounding surface to penalize distance to.

  • eq (Equilibrium) – Equilibrium that will be optimized to satisfy the Objective.

  • eq – Equilibrium that will be optimized to satisfy the Objective.

  • target ({float, ndarray}, optional) – Target value(s) of the objective. Only used if bounds is None. Must be broadcastable to Objective.dim_f.

  • bounds (tuple of {float, ndarray}, optional) – Lower and upper bounds on the objective. Overrides target. Both bounds must be broadcastable to to Objective.dim_f

  • weight ({float, ndarray}, optional) – Weighting to apply to the Objective, relative to other Objectives. Must be broadcastable to to Objective.dim_f

  • normalize (bool, optional) – Whether to compute the error in physical units or non-dimensionalize.

  • normalize_target (bool, optional) – Whether target and bounds should be normalized before comparing to computed values. If normalize is True and the target is in physical units, this should also be set to True.

  • surface_grid (Grid, optional) – Collocation grid containing the nodes to evaluate surface geometry at.

  • plasma_grid (Grid, optional) – Collocation grid containing the nodes to evaluate plasma geometry at.

  • use_softmin (bool, optional) – Use softmin or hard min.

  • alpha (float, optional) – Parameter used for softmin. The larger alpha, the closer the softmin approximates the hardmin. softmin -> hardmin as alpha -> infinity. if alpha*array < 1, the underlying softmin will automatically multiply the array by 2/min_val to ensure that alpha*array>1. Making alpha larger than this minimum value will make the softmin a more accurate approximation of the true min.

  • name (str, optional) – Name of the objective function.

Methods

build([eq, use_jit, verbose])

Build constant arrays.

compute(*args, **kwargs)

Compute plasma-surface distance.

compute_scalar(*args, **kwargs)

Compute the scalar form of the objective.

compute_scaled(*args, **kwargs)

Compute and apply weighting and normalization.

compute_scaled_error(*args, **kwargs)

Compute and apply the target/bounds, weighting, and normalization.

compute_unscaled(*args, **kwargs)

Compute the raw value of the objective.

copy([deepcopy])

Return a (deep)copy of this object.

eq(other)

Compare equivalence between DESC objects.

jit()

Apply JIT to compute methods, or re-apply after updating self.

load(load_from[, file_format])

Initialize from file.

print_value(*args, **kwargs)

Print the value of the objective.

save(file_name[, file_format, file_mode])

Save the object.

xs(eq)

Return a tuple of args required by this objective from the Equilibrium eq.

Attributes

args

Names (str) of arguments to the compute functions.

bounds

Lower and upper bounds of the objective.

built

Whether the transforms have been precomputed (or not).

constants

Constant parameters such as transforms and profiles.

derivatives

Derivatives of the function wrt the argument given by the dict keys.

dim_f

Number of objective equations.

dimensions

Dimensions of the argument given by the dict keys.

fixed

Whether the objective fixes individual parameters (or linear combo).

linear

Whether the objective is a linear function (or nonlinear).

name

Name of objective function (str).

normalization

normalizing scale factor.

scalar

Whether default "compute" method is a scalar or vector.

target

Target value(s) of the objective.

weight

Weighting to apply to the Objective, relative to other Objectives.