Пример #1
0
def _matmul_linear_operator(linop_a, linop_b):
    """Generic matmul of two `LinearOperator`s."""
    is_square = registrations_util.is_square(linop_a, linop_b)
    is_non_singular = None
    is_self_adjoint = None
    is_positive_definite = None

    if is_square:
        is_non_singular = registrations_util.combined_non_singular_hint(
            linop_a, linop_b)
    elif is_square is False:  # pylint:disable=g-bool-id-comparison
        is_non_singular = False
        is_self_adjoint = False
        is_positive_definite = False

    return linear_operator_composition.LinearOperatorComposition(
        operators=[linop_a, linop_b],
        is_non_singular=is_non_singular,
        is_self_adjoint=is_self_adjoint,
        is_positive_definite=is_positive_definite,
        is_square=is_square,
    )
Пример #2
0
def _matmul_linear_operator(linop_a, linop_b):
    """Generic matmul of two `LinearOperator`s."""
    is_square = _is_square(linop_a, linop_b)
    is_non_singular = None
    is_self_adjoint = None
    is_positive_definite = None

    if is_square:
        is_non_singular = _combined_non_singular_hint(linop_a, linop_b)
        is_self_adjoint = _combined_self_adjoint_hint(linop_a, linop_b)
    elif is_square is False:
        is_non_singular = False
        is_self_adjoint = False
        is_positive_definite = False

    return linear_operator_composition.LinearOperatorComposition(
        operators=[linop_a, linop_b],
        is_non_singular=is_non_singular,
        is_self_adjoint=is_self_adjoint,
        is_positive_definite=is_positive_definite,
        is_square=is_square,
    )