autoarray.inversion.regularization.AdaptSplit#

class AdaptSplit[source]#

Bases: Adapt

Regularization which uses the derivatives at a cross of four points around each pixel centre and values adapted to the data being fitted to smooth an inversion’s solution.

An adaptive regularization scheme which splits every source pixel into a cross of four regularization points and interpolates to these points in order to smooth an inversion’s solution.

The size of this cross is determined via the size of the source-pixel, for example if the source pixel is a Delaunay pixel the area of the pixel is computed and the distance of each point of the cross is given by the area times 0.5.

For the weighted regularization scheme, each pixel is given an ‘effective regularization weight’, which is applied when each set of pixel neighbors are regularized with one another. The motivation of this is that different regions of a pixelization’s mesh require different levels of regularization (e.g., high smoothing where the no signal is present and less smoothing where it is, see (Nightingale, Dye and Massey 2018)).

Unlike Constant regularization, neighboring pixels must now be regularized with one another in both directions (e.g. if pixel 0 regularizes pixel 1, pixel 1 must also regularize pixel 0). For example:

B = [-1, 1] [0->1]: [-1, -1] 1 now also regularizes 0

For Constant regularization this would NOT produce a positive-definite matrix. However, for the weighted scheme, it does!

The regularize weight_list change the B matrix as shown below - we simply multiply each pixel’s effective regularization weight by each row of B it has a -1 in, so:

regularization_weights = [1, 2, 3, 4]

B = [-1, 1, 0 ,0] # [0->1]: [0, -2, 2 ,0] # [1->2] [0, 0, -3 ,3] # [2->3] [4, 0, 0 ,-4] # [3->0]

If our -1’s werent down the diagonal this would look like:

B = [4, 0, 0 ,-4] # [3->0]: [0, -2, 2 ,0] # [1->2] [-1, 1, 0 ,0] # [0->1] [0, 0, -3 ,3] # [2->3] This is valid!

A full description of regularization and this matrix can be found in the parent AbstractRegularization class.

JAX & gradient support (2026-07 gradient sweep): as for ConstantSplit — differentiable and FD-certified on the Delaunay mesh family (e.g. the KNN meshes), structurally incompatible with the rectangular meshes. Note the defaults inner_coefficient == outer_coefficient == 1.0 make the weighting uniform — numerically identical to ConstantSplit(coefficient=1.0).

Parameters:

coefficients – The regularization coefficients which controls the degree of smoothing of the inversion reconstruction in high and low signal regions of the reconstruction.
signal_scale (float) – A factor which controls how rapidly the smoothness of regularization varies from high signal regions to low signal regions.

Methods

`log_det_regularization_matrix_term_from`	Returns `log det H` of this scheme's regularization matrix computed from a factorization the scheme itself knows about, or `None` when no such shortcut exists (the default).
`regularization_matrix_from`	Returns the regularization matrix with shape [pixels, pixels].
`regularization_term_from`	Returns this scheme's contribution to the regularization term `s^T H s` computed from a factorization the scheme itself knows about, or `None` when no such shortcut exists (the default).
`regularization_weights_from`	Returns the regularization weights of this regularization scheme.

Attributes

is_split_regularization

Whether this scheme is a "split" regularization variant, which regularizes using a split-cross calculation of the mesh's mappings rather than the mappings themselves.

is_split_regularization = True#

Whether this scheme is a “split” regularization variant, which regularizes using a split-cross calculation of the mesh’s mappings rather than the mappings themselves.

Split schemes require the mesh’s interpolator to provide _mappings_sizes_weights_split, which only the adaptive meshes (e.g. Delaunay, DelaunayNN, KNNBarycentric) do. Pixelization uses this flag together with AbstractMesh.supports_split_regularization to reject unsupported combinations at construction.

regularization_matrix_from(linear_obj, xp=<module 'numpy' from '/home/docs/checkouts/readthedocs.org/user_builds/pyautogalaxy/envs/latest/lib/python3.12/site-packages/numpy/__init__.py'>)[source]#

Returns the regularization matrix with shape [pixels, pixels].

Parameters:: linear_obj (LinearObj) – The linear object (e.g. a Mapper) which uses this matrix to perform regularization.
Return type:: The regularization matrix.