Skip to content

Add numerical stabilization for difference of exponentials #1303

@rlouf

Description

@rlouf

From aesara-devs/aeppl#131

def logdiffexp(a, b):
    """log(exp(a) - exp(b))"""
    return a + at.log1mexp(b - a)

Metadata

Metadata

Assignees

No one assigned

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions